/*求σ(n)的和,其中σ(n)是n的除数之和。(非常快的程序)由P L Patodia于2008年1月9日编写2008年6月25日修改(有小错误,但这并不影响A072692的值)*/A024916(z)={局部(s,u,d,n,a,p);s=z*z;u=平方(z);p=0;对于(d=1,u,n=z\d-z\(d+1);如果(n<=1,p=d;break(),a=z%d;s-=(2*a+(n-1)*d)*n/2););如果(p==0,u=z\(u+1),u=z\p);对于(d=2,u,s-=z%d);申报表;}/*A072692-1<=j<=10^n的σ(j)之和,其中σ(j)是j的除数之和。*/A072692(n)={对于(i=1,n,打印(i,“:”,A024916(10^i))}