%I#38 2019年3月2日02:32:03
%S 1,10,25,34,58,64,85,91121130169196214226289324370400526,
%电话:52962567670673077178411024108912251255135114141444,
%电话:152116811849190619362116220923042500280929863136336434813600
%不是平方和素数之和的数字。
%C序列是无限的:如果2n-1是复合的,那么n^2在序列中。(证明:如果n^2=x^2+p带有p素数,那么p=(n-x)(n+x),那么n-x=1和n+x=p。因此2n-1=p是素数,而不是复合数。)_Dean Hickerson,2002年11月27日
%C 21679是这个序列中最后一个已知的非方形。参见A020495_T.D.Noe_,2006年8月5日
%C A002471(a(n))=0;A014089.-的补充_Reinhard Zumkeller_,2008年9月7日
%这个序列中没有素数,因为它们至少可以表示为p+0^2_阿隆索·德尔·阿特(Alonso del Arte),2012年5月26日
%C术语数量<10^k,k=0..8:1,8,27,75,223,719,2361,7759,…,.-_Robert G.Wilson v_,2012年5月26日
%C到目前为止,只有21项不是正方形,它们是A020495的项。那些是正方形的,它们的平方根是A104275的成员_Robert G.Wilson v_,2012年5月26日
%H Robert G.Wilson v,<a href=“/A014090/b014090.txt”>n表,n=1..100000的a(n)
%e From _Alonso del Arte,2012年5月26日:(开始)
%e 10在序列中,因为10-p_i中没有一个是正方形(8,7,5,3),10-b^2中没有一个是素数(10,9,6,1);i从1变为pi(10),或者b从0变为floor(sqrt(10))。
%e11不在序列中,因为它可以表示为3^2+2或0^2+11。(结束)
%t={};Do[k=0;While[k^2<n&&!PrimeQ[n-k^2],k++];当[k^2<n&&!PrimeQ[n-k^2],k++]时;如果[k^2>=n,附加到[t,n]],{n,25000}];t(*t.D.Noe_,2006年8月5日*)
%t最大值=5000;补码[Range[max],Flatten[Table[Prime[p]+b^2,{p,PrimePi[max]},{b,0,Ceiling[Sqrt[max]]](*_Alonso del Arte_,2012年5月26日*)
%t fQ[n_]:=块[{j=Sqrt[n],k},如果[IntegerQ[j]&&!PrimeQ[2j-1],真,k=地板[j];而[k>-1&&!素数Q[n-k^2],k--];如果[k==-1,真,假]]];选择[范围[3600],fQ](*RobertG.Wilson v_1012年5月26日*)
%Y参考A020495,A104275。
%Y参考A064233(不允许0^2)。
%不,简单,好
%O 1,2号机组
%A.N.J.A.Sloane,R.K.盖伊_
|