fish <- function(n) choose(50,3)*choose(n-50,47)/choose(n,50)
n <- 50:2000
m <- max(fish(n))
m
n[fish(n)==m]