range <- (1:20)*2 win <- sapply(range, function(N) 1-pbinom(N/2, N, 12/25)) max(win) range[win==max(win)]