登录
A328667型
a(n)=Sum_{d除以n}(-1)^(n+1+d+n/d)*d^2。
2
1, 5, 10, 13, 26, 50, 50, 45, 91, 130, 122, 130, 170, 250, 260, 173, 290, 455, 362, 338, 500, 610, 530, 450, 651, 850, 820, 650, 842, 1300, 962, 685, 1220, 1450, 1300, 1183, 1370, 1810, 1700, 1170, 1682, 2500, 1850, 1586, 2366, 2650, 2210, 1730, 2451, 3255
抵消
1,2
链接
公式
与a(2^e)相乘=(2^(2*e+1)+7)/3=A321358型(e) 如果e>0,则a(p^e)=(p^(2*e+2)-1)/(p^2-1)。
G.f.:总和=1}k^2*x^k/(1+(-x)^k)=总和=1}x^k*(1-(-x)^k)/(1+(-x)^k)^3。
a(n)=-(-1)^n*A321558型(n) ●●●●。a(2*n-1)=A001157号(2*n-1)=A099978号(n) ●●●●。a(4*n+2)=A001157号(4*n+2)。
求和{k=1..n}a(k)~c*n^3,其中c=7*zeta(3)/24=0.350599-阿米拉姆·埃尔达尔2022年11月1日
例子
G.f.=x+5*x ^2+10*x ^3+13*x ^4+26*x ^5+50*x ^6+50*x ^7+45*x ^8+。。。
数学
a[n_]:=如果[n<1,0,DivisorSum[n,(-1)^(n+1+#+n/#)#^2&];
黄体脂酮素
(PARI){a(n)=sumdiv(n,d,(-1)^(n+1+n\d+d)*d^2)};
关键词
非n,多重
作者
迈克尔·索莫斯2019年10月24日
状态
经核准的