|
（PARI）admissibleMod（M=1000，t=[3,9]，debug=0）={local（tt=1）；而（M>tt*=10，t=concat（[t，t+向量（#t，i，t[1]）*tt，t+矢量（#t、i，t[2]）*tt]））；if（debug，print（“t=”t））；t=Set（t）；tt=[]；
对于（i=1，M，如果（setsearch（t，i^2%M），tt=concat（tt，i））；concat（tt，M+tt[1]）}A058433号（Nmax=1e10，N=1，addMod=100000，debug=1）={local（a=[]，addNext=admissibleMod（addMod=10^round（log（addMod）/log（10）），[3，9]），d=1，加=向量（addMod，i，if（i-1>addNext[d]，d++）；addNext[d]-i+1），pow10=向量（d=#Str（（（（Nmax=round（Nmax）））^2），i，10^（i-1）））），nextOK=[0，2、1、0、5、4、3、2、1和0]、N下一个、N2、数字位数、位置）；
nextOK=向量（#nextOK，i，if（nextOK[i]，nextK[i]*pow10））；N=圆形（N）；而（Nmax>=N，numDigits=#Str（N2=N^2）；如果（place=nextOK[1+d=N2\pow10[numDigits]]，N=max（平方（place[numBigits]+d*pow10[数字位数]），N+1）；下一步）；Nnext=最小值（Nmax，sqrtint（round（（d+1）*pow10[numDigits]））；
if（调试，打印（“从“N”检查到“Nnext”：<=”，1+max（0，Nnext-N）*（#addNext-1）\addMod，“候选”）；N+=添加[1+N%addMod]；位置=1；而（Nnext>=N，dr=divrem（N2=N^2，pow10[位置=numDigits]）；while（place--&！d=nextOK[1+（dr=divrem（dr[2]，pow10[place]））[1]，）；
if（！place，break）；N=平方（N2-dr[2]+d[位置]）+1；N+=添加[1+N%addMod]；）；if（！place，if（调试，打印（N，“^2=”，N^2））；a=凹面（a，N））；N=Nnext*3\2+1）；一个}-M.F.哈斯勒2007年5月14日
| 