%I#11 2016年1月12日00:11:56
%S 2,5,13,37,41,53,73109137157173181193197233349373401421,
%电话:4575095415575695795776136617097337577698218779411009,
%电话:10331069111711291193120112371301137314531493931549159716211657166916971721
%形式为x^2+y^2的N个素数,其中x^2+x*y+y^2也是带有正x和y的素数。
%C序列关注x和y的正值。否则,例如17=4^2+(-1)^2和4^2+4*(-1)+(-1)*2=13也是素数。
%C相应的广义古巴素数是3,7,19,43,61,67,97,139,181,199,211,223,271,277,337,421,439,499,541,601,619,631,751,787,811,823,829,919。。。
%e5是一个项,因为2^2+1^2=5是素数,2^2+2*1^2=7是素数。
%e13是一个项,因为3^2+2^2=13是素数,3^2+3*2+2^2=19是素数。
%e37是一个项,因为6^2+1^2=37是质数,6^2+6*1+1^2=43是质数。
%t极限=50;取[Select[Union@Flatten@Table[If[PrimeQ[Abs[x^2+xy+y^2],x^2+y^2,Nothing],{x,lim},{y,lim{],PrimeQ,53](*_Michael De Vlieger_,2016年1月7日*)
%o(PARI)列表(lim)=我的(v=列表(),t);lim=1;对于(x=1,平方(lim),对于(y=1,min(平方(lim-x^2),x),如果(i素数(t=x^2+y^2)&&i素数;向量排序(Vec(v),8)
%Y参考A002313、A007645。
%K nonn公司
%O 1,1号机组
%A Altug Alkan,2016年1月7日
|