(PARI)eva(n)=subst(Pol(n),x,10)
插入(n,len,pos)=我的(d=数字(n),v=[],w=[]);对于(y=1,pos,v=concat(v,d[y]));v=concat(v,向量(len));对于(z=pos+1,#d,v=concat(v,d[z]));eva(v)
a306920(n)=表示素数(p=10,表示(k=1,#数字(p)-1,my(zins=插入(p,n,k));if(ispseudoprime(zins),返回(p))
删除零(n)=我的(v=[],w=[]、d=数字(n));对于(x=1,#d,如果(d[x]==0,对于(y=1,x-1,w=concat(w,d[y]));对于(z=x+1,#d,w=concat(w,d[z]));如果(#w>0,v=concat(v,[eva(w)]);w=[]);向量排序(v,8)
是(n)=my(x=a306920(n),y=a306820(n+1),rz=remove_zeros(x));如果(#setintersect([0],vecsort(digits(x)))==0,则返回(0));对于(k=1,#rz,如果(y==rz[k],返回(0));1
|