提出
经核准的
编辑
打印 (A类 ) #彼得·卢什尼2014年6月10日
F.查波顿：为py3调整sage代码
okQ[k_]：=模块[{ff=FactorInteger[k^2+1]}，长度[ff]==2&ff[[All，2]]=={1，1}&&IntegerQ[Sqrt[ff[[2，1]]-ff[[1，1]]]]；
选择[范围[2，40000，2]，okQ](*让-弗朗索瓦·奥尔科弗2019年6月25日*）
（鼠尾草）
A=[]
对于范围（22000,2）中的k：
K=K^2+1
f=素数除数（K）
如果len（f）==2：
如果mul（f）==K：
如果是平方（abs（f[0]-f[1]））：
A.附录（k）
打印A#彼得·卢什尼2014年6月10日
