%I#18 2024年4月18日06:11:00
%S 14720439165356283471471040741208614363886033463081172360113,
%电话:136815694115741008891924496102198925349920218602436774546339,
%电话:97705416101223085596312311606487125408424461451372377726433294893864872419847638558043501958869165081131993156449248367
%N可以用两种方式写成连续素数平方和的数字。
%H Michael S.Branicky,n表,n=1..991的a(n)</a>
%H Cathal O'Sullivan、Jonathan P.Sorenson和Aryn Stahl,<a href=“https://arxiv.org/abs/2204.10930“>求素数连续幂和的算法</a>,arXiv:2204.10930[math.NT],2022-2023。见4.2副本第8-9页。
%H Michael S.Branicky,Python程序</a>
%o(Python)#查看适用于生成b文件的版本的链接
%o来自sympy import primerage,integer_nthroot
%o定义缺陷(极限):
%o adict=dict()
%o根限制=integer_ntroot(限制,2)[0]
%o对于素数范围(2,根极限+1)中的x:
%o s=x**2
%o根[s]=1
%o表示素数范围内的y(x+1,根极限+1):
%o s+=y**2
%o如果s<=极限:
%o如果s不在adic中:
%o根[s]=1
%o其他:
%o修正[s]+=1
%o其他:
%o中断
%o如果adict[s]==2,则返回adict中s的排序结果
%o打印(aupto(6*10**10))#_Michael S.Branicky_,2022年4月26日
%Y参考A001248,A340771。
%K nonn公司
%O 1,1号机组
%马库斯,2022年4月26日
|