Q: =proc(n)局部a,e,p,f;a: =1;对于ifactors(n)[2]中的f,做e:=op(2,f);p: =op(1,f);
如果类型为(e,奇数),则a:=a*(p^(e+1)-1)/(p-1);否则a:=a*((p^(e+1)-1)/(p-1)-p^(e/2));fi;od:a;结束时间:
P: =proc(h)局部a,b,k,n;
对于从1到h的n,做a:=除数(n);b: =0;
对于从1到nops(a)的k,如果gcd(a[k],n/a[k])=1,那么b:=b+a[k';fi;od;
如果Q(n)-b=n,则打印(n);fi;od;结束:P(10^6);
|