（PARI）search（u）={/*慢，u必须是一个小整数*/
my（log2=log（2），q，t，t0，L1=List（））；
对于素数（y=3，素数（10^u），
t=对数（y+1）\log2；
而（t>t0，
q=4*t^2+y+1；
如果（q%（t*（2^（2*t+1）））＝＝0，
列表输入（L1，[t，y]）；
t0=t；
打破
,
t吨--
)));
第一层
}
（PARI）a（n）=我的（k=n<<（2*n+1），p=-4*n^2-1）；while（！i素数（p+=k），）；第页\\查尔斯·格里特豪斯四世2014年9月18日
