使用（数字理论）：T:=进程（w）局部x，y，z；x： =w；y： =0；对于z从1到ilog10（x）+1做y:=10*y+（xmod10）；x： =trunc（x/10）；od；年；结束时间：
P： =proc（q）局部a，j，k，n，t；打印（1）；对于从1到q的n，做k：=0；j： =n；
而jmod2<>1做k:=k+1；j： =j/2；od；
a： =西格玛（2*n+1）+sigma（2*n-1）+sigrama（n/2^k）*2^（k+1）-6*n-2；
如果T（a）=a且σ（n）=T（σ（n）），则打印（n）；fi；od；结束：P（10^6）；
