使用(数字理论):P:=proc(q,h,w)局部a,b,k,n,t,v;v: =阵列(1..h);
对于从1到q的n,做a:=w*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,6);
|