(PARI)
默认值(质数限制,(2^31)+(2^30));
is_persistently_squarefere(n,base)={while(n>1,如果(!issquarefree(n),return(0));n=base);(1);};
A005940号(n) ={my(p=2,t=1);n--;直到(!n\=2,if((n%2),(t*=p),p=nextprime(p+1));t};\\根据的代码修改M.F.哈斯勒
isA293430(n)=is_persistently_squarefere(n,2);
n=0;k=1;while(n<=2^26,n=n+1;if(isA293430(n)&&isA29343(A163511号(n) ),写入(“b293437.txt”,k,“”,n);k=k+1));