(PARI)\\数值函数的连续分式有理逼近
cfrac(m,f)=x=f;对于(n=0,m,i=楼层(x));x=1/(x-i);打印1(i,“,”)
cfracnumprime(m,f)={cf=向量(100000);x=f;对于(n=0,m,i=楼层(x);x=1/(x-i);cf[n+1]=i;);对于(m1=0,m,r=cf[m1+1];对于步进(n=m1,1,-1,r=1/r;r+=cf[n];);numer=分子(r);分母(r):如果是素数(分子),则打印1(分子,“,”);)}
(PARI)
默认值(realprecision,10^5);
cf=控制(Pi);
n=0;
{对于(k=1,#cf,\\生成b文件
pq=contfracpnqn(向量(k,j,cf[j]));
p=pq[1,1];q=pq[2,1];
如果(ispseudoprime(p),n+=1;打印(n,“,p))\\A086785号
\\如果(ispseudoprime(q),n+=1;打印(n,“,q))\\A086788号
); }
|