#来自在线整数序列百科全书的问候!http://oeis.org/搜索:id:a248015显示第1-1页,共1页%一A248015%S A248015 8,18,28,30,34,44,46,48,50,58,60,64,68,70,76,78,86,88,96,98100104,%电话:A248015 108114118128136144148158164166168178186188190194196,%美国A248015 198200%nA248015个正数N使得N^2+1是复合的,并且没有正整数c和z,使得N=c*z^2+z+c。%C A134407的A248015子集。%C A248015如果f(x)=x^2+1且g(C,y)=C*y^2+y+C,则g对x的代数代换给出了一个因式分解:f(g(C,y))=(y^2+1)*(C^2*y^2+C^2+2*C*y+1)。因为f(g(c,y))的两个因子都是大于1的整数,所以f(g(c,y))是一个复合数。%C A248015这些数字必然是来自A134407的偶数项,因为对于奇数n=2c+1,有z=1的“禁止”分解_M、 F.Hasler_2014年10月4日%H A248015埃里克·韦斯坦的数学世界,兰道的问题%最大A2200:48025牛顿:%p A248015 mb:=过程(n::整数)::整数;%p A248015如果isprime(n^2+1)=false,则返回n else返回-1 fi;%p A248015结束程序:%p A248015 A134407:=矢量(最大值):%p A248015适用于a从1到最大值A134407[a]:=mb(a):结束do:%p A248015 A134407s:=转换(A134407,'set')减{-1}:%p A248015 A134407l:=转换(A134407s,“列表”):%p A248015,适用于c从1到200 do%p A248015,用于1至20 do的z%p A248015 A134407s:=A134407s减去{c*z^2+z+c}:%p A248015端do:%p A248015端do:%p A248015 A134407s;%o A248015(PARI)is(n)={!比特(n,0)&&!isprime(n^2+1)&&!for(z=2,sqrtint(n),(n-z)%(z^2+1)| | return)}\\\\\\\\\\\\哈斯勒,2014年10月4日%Y A248015,参见A134407。%K A248015无%O A248015 1,1%2014年9月29日,A A248015马特C.安德森#根据OEIS最终用户许可协议提供内容:http://oeis.org/LICENSE