Problem 71 - Project Euler
13秒。
の中に答があると決めてかかったけど、それほど自明でもない気がしてきた。
library(gmp) divisor <- 1:10^6 numrtor <- floor(divisor*3/7) frac <- numrtor/divisor frac.o <- order(frac, decreasing=TRUE) for(i in frac.o[-1]){ if(gcd.bigz(divisor[i], numrtor[i])==1) break }