面白い問題おしえて～な 28問目

面白い問題おしえて～ ..

110:１３２人目の素数さん
18/11/06 08:57:44.56 9NNsjRpE.net
>>102
ようやくコードの意味が理解できたのでRに移植。
実数計算なので誤差がでます。
f = function(N){
p=numeric()
p[1]=p[2]=p[3]=p[4]=p[5]=p[6]=1
for(n in 7:(N+1)){
p[n]=p[n-1]-p[n-6]/(6^6)
}
1-p[N+1]
}
> f(1000)
[1] 0.02110296

次ページ

続きを表示

1を表示