(PARI){对于素数(p=3,10^3,对于(i=precprime(p-1)+1,p,(2*i)%3==0&isprime(2*i-1)&next(2);i%2==0&(i/2)%3!=2&isprim(i/2
nextA124123(p)={while(p=nextprime(p+1),对于(i=precprime(p-1)+1,p,(2*i)%3==0&isprime(2*i-1)&next(2);i%2==0&(i/2)%3!=2&isprim(i/2
t=2;向量(200,i,t=nextA124123(t))\\前200项的60%位于1+3Z:
t=[0,0];向量(#%,i,t[%i]%3]++);t产生[120,80]
t=10^11;向量(200,i,t=nextA124123(t))\\这些项中有50%位于1+3Z:
t=[0,0];向量(#%,i,t[%i]%3]++);t产生[100100]
t=10^30;向量(200,i,t=nextA124123(t+1));t-10^30得出31773=10^30以外第200项的距离
t=10^30;向量(200,i,t=下一素数(t+1));(t-1e30)/%产生0.52…,约为底漆的局部密度。(结束)
|