某氏より、これは有名な秘書問題だよと教えていただいた。 https://ja.wikipedia.org/wiki/%E7%A7%98%E6%9B%B8%E5%95%8F%E9%A1%8Cこれで簡単に解析的に解けると。 n=1000 f <- function(r){ sum(1 / (r:n - 1)) * (r - 1) / n } f2 <- Vectorize(f) ans <- f…