|
黄体脂酮素
|
(PARI)
\\数值常数f.m=步长的连续分式有理逼近。
cfracnumprime(m,f)={默认值(realprecision,3000);cf=矢量(m+10);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];);分子=分子(r);分母=分母(r),如果是udoprime(分子),打印1(数字,“,”);)}
(PARI)
默认值(realprecision,10^5);
cf=控制(exp(1));
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))\\A086791号
\\如果(ispseudoprime(q),n+=1;打印(n,“,q))\\A094008号
); }
|