其中（numtheory）：P:=proc（q）局部a，b，c，d，j，k，n，ok；对于从3到q的n，做a:=[]；b:=n；确定：=0；
d： =n*加法（op（2，p）/op（1，p），p=ifactors（n）[2]）；a:=[]；b:=n；确定：=0；
对于k从1到ilog10（n）+1，如果（b mod 10）>1，则确定：=1；fi；a:=[（b mod 10），op（a）]；b:=trunc（b/10）；od；b:=-1；c:=0；
如果ok=1，则当c<d时，执行b:=b+1；
如果b>0，则c:=加（a[k]^b，k=1..nops（a））；否则，对于k从1到nops（a），如果a[k]=0，那么c:=0；断裂；
否则c:=c+1；fi；od；fi；od；如果c=d，则lprint（n，b）；fi；fi；od；结束：P（10^9）；
