%I M3797 N1549#48 2023年3月9日05:04:21

%S 1,5,10,17,16,32,22,41,37,50,34,80,40,68,70,89,52113,58122,94104，

%电话70176,91122118164,88212,94185142158142269112176166266，

%电话：124284130248232212142368169275214290160356214356238266178500184284310377250428202374286

%N a（N）=和{d|N，d<=3}d^2+3*和{d| N，d>3}d。

%D P.A.MacMahon，除数平方和与给定数的分区数之间的联系，Messenger Math。，54 (1924), 113-116. 《论文集》，麻省理工出版社，1978年，第一卷，第1364-1367页。

%D N.J.A.Sloane，《整数序列手册》，学术出版社，1973年（包括该序列）。

%D N.J.A.Sloane和Simon Plouffe，《整数序列百科全书》，学术出版社，1995年（包括该序列）。

%F猜想：g.F.的逆Moebius变换：（x+2x^2+2x^3-2x^4）/（1-x）^2.-_Sean A.Irvine_，2014年5月16日

%F a（n）=3*σ（n）-k，其中，如果n是偶数，k=4；如果n是奇数，k=2_Sean A.Irvine_，2014年5月16日

%p（数字理论）：

%p A:=proc（s，n）局部d，s1，s2；

%ps1:=0；s2:=0；

%除数（n）do中d的p

%p如果d<=s，则s1:=s1+d^2，否则s2:=s2+d；fi；日期：

%ps1+s*s2；结束；

%p f：=s->[序列（A（s，n），n=1..80）]；f（3）；

%t a[n_]：=除数和[n，如果[#<=3，#^2，3#]&]；

%t表[a[n]，{n，1100}]（*Jean-François Alcover_，2023年3月9日*）

%o（PARI）a（n）=总和（n，d，如果（d<=3，d^2）+3*，如果（d>3，d））；\\_米歇尔·马库斯，2023年3月9日

%Y参见A000203、A002659、A002791。

%Y A242639中数组的一行。

%K nonn公司

%O 1,2号机组

%A _N.J.A.斯隆_

%E条目由N.J.A.Sloane修订，2014年5月21日