Pythagorean prime
オンライン整数列大辞典の
A002144（http://oeis.org/A002144/list）
と比較し、答え合わせしてみる。
require 'prime'
def A002144(n)
Prime.each(10 ** 6).select{|i| i % 4 == 1}[0..n - 1]
end
ary = A002144(54)
# OEIS A002144のデータ
ary0 =
[5,13,17,29,37,41,53,61,73,89,97,101,109,113,137,
149,157,173,181,193,197,229,233,241,257,269,277,
281,293,313,317,337,349,353,373,389,397,401,409,
421,433,449,457,461,509,521,541,557,569,577,593,
601,613,617]
# 一致の確認
p ary == ary0
0 件のコメント:
コメントを投稿
注: コメントを投稿できるのは、このブログのメンバーだけです。