| 数据
|
10, 15, 17, 57, 68, 520, 525, 547, 789, 894, 3149, 3700, 3845, 4272, 8304, 15792, 63555, 119926, 121267
|
| MAPLE公司
|
使用(数字理论):P:=proc(q,h)局部a,b,k,n,t,v;v: =阵列(1..h);
对于从2到q的n,做a:=n*(n+1)/2-sigma(n);b: =ilog10(a)+1;如果b>1,则
对于从1到b的k,dov[b-k+1]:=(a mod 10);a: =trunc(a/10);od;t: =b+1;v[t]:=加(v[k],k=1..b);
而v[t]<n dot:=t+1;v[t]:=加(v[k],k=t-b..t-1);od;
如果v[t]=n,则打印(n);fi;fi;od;结束:P(10^61000);
|