2017年12月9日土曜日

171209

Ruby


Number of Gaussian primes a+b*i such that sqrt(a^2 + b^2) <= 10^n

A091098 とA091099 のデータを用いて計算してみた。

ary1 = [0] +
[1,11,80,609,4783,39175,332180,2880504,25423491,
 227523275,2059020280,18803924340,173032709183,
 1602470783672,14922284735484,139619168787795,
 1311778575685086,12369977142579584,
 117028833597800689,1110409801150582707]
ary3 = [0] +
[2,13,87,619,4808,39322,332398,2880950,25424042,
 227529235,2059034532,18803987677,173032827655,
 1602470967129,14922285687184,139619172246129,
 1311778581969146,12369977145161275,
 117028833678543917,1110409801410336132]
p (1..10).map{|i| 8 * ary1[2 * i] + 4 * ary3[i] + 4}

出力結果
[100, 4928, 313752, 23046512, 1820205436, 150431552012, 12819767598972, 1116953361826164, 98959817242332844, 8883278410114778600]

0 件のコメント:

コメントを投稿

注: コメントを投稿できるのは、このブログのメンバーだけです。