登录
A125141号
a(1)=2;对于n>1,a(n)=SENSigma(a(n-1)),其中SENSigga(m)=(-1)^((Sum_i r_i)+Omega(m(rj)是pj除以m的最大幂。
2, 3, 4, 5, 6, 12, 20, 30, 72, 165, 288, 693, 1056, 3024, 9280, 22500, 42845, 60480, 240000, 794580, 1814400, 7040040, 26352000, 98654400, 321552000, 1260230400, 5311834416, 17570520000, 75087810000, 325180275840, 1526817600000
抵消
1,1
评论
“Max(r_j)”的含义如下:如果d|m,d=p^e*q^f,m=p^x*q^y*r^z,则Max(e)=x,Max(f)=y。
MAPLE公司
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日
数学
SENSigma[n_]:=模块[{Ifs,i,a,r,p},Ifs=系数整数[n];a=1;对于[i=1,i<=长度[Ifs],i++,r=Ifs[[i,2]];p=Ifs[[i,1]];a=a(p(1-p^r)/(1-p)-(-1)^r)];返回[a]];
A125141号[nmax_]:=模块[{a},a={2};而[Length[a]<nmax,a=附加[a,SENSigma[a[-1]]];返回[a]];
关键词
非n
作者
Yasutoshi Kohmoto公司2007年1月12日
扩展
编辑人N.J.A.斯隆根据…的建议安德鲁·普列韦(Andrew S.Plewe)2007年5月14日
更多术语来自R.J.马塔尔2007年5月18日
状态
经核准的