底辺私立医大を卒業した医者って頭悪いよね？ Part22

n=12345
f=\(x) 3.135<(n-x)/x & (n-x)/x<3.145
de=(1:n)[sapply(1:n,f)]
nu=n-de
nu/de
ans=NULL
for(i in 1:length(nu)) ans=c(ans,numbers::GCD(nu[i],de[i])==1)
NU=nu[ans]
DE=de[ans]
for(j in 1:length(NU)) cat(paste0(NU[j],'/',DE[j],' ',round(NU[j]/DE[j],4),'\n'))
NU/DE