simulateP <- 0.5
answer <- 0
for (i in 0:30){
   for (j in i:min(44,25+i)){
      for (k in j:min(58,25+j)){
         pri <- dbinom(i,50,simulateP)
         prj <- dbinom(j-i,25,simulateP)
         prk <- dbinom(k-j,25,simulateP)
         answer <- answer + pri*prj*prk
      }
   }
}