SENSigma:=proc(n)局部ifs,i,a,r,p;ifs:=ifactors(n)[2];a:=1;对于i从1到nops(ifs)dor:=op(2,op(i,ifs));p:=op(1,op(i,ifs));a:=a*(p*(1-p^r)/(1-p)-(-1)^r);od;返回(a);结束时间:125141英镑:=proc(nmax)局部a;a:=[2];而nops(a)<nmax do a:=[op(a),SENSigma(op(-1,a))];od;返回(a);结束时间:A125141号(40) ; #R.J.马塔尔2007年5月18日