Problem 58 - Project Euler
gmpのisprimeはベクトル渡せるから便利.
# library(gmp) plus <- 4 r.btm <- 9 primes <- 3 all.diag <- 5 while(primes/all.diag > 0.1){ candy <- r.btm + plus*(1:4) primes <- primes + length(candy[isprime(candy) > 0]) all.diag <- all.diag + 4 plus <- plus + 2 r.btm <- candy[4] } plus - 1