(PARI)默认值(primelimit,10^7)istri(n)=t=楼层(sqrt(2*n));如果(2*n==t*(t+1),1,0)表示素数(p=2,5*10^6,if(istri(p^2-1),print1(p“”))
(PARI)istriang(n)=发行方(8*n+1);
对于素数(p=2,10^10,if(istriang(p^2-1),print1(p,“,”));
(PARI)/*效率更高:*/
N=1166;f=(1+x-4*x^2-2*x^3)/((x^2+2*x-1)*(x^2-2*1))+O(x^N);
对于(n=0,n-1,my(c=polceoff(f,n));如果(i首字母(c),打印1(c,“,”));