%I#24 2023年10月13日06:52:31

%S 0,16,8125662513122401409666421001614641209922856138432，

%电话：513316553683521106288130321160256196963234272279841335872，

%电话：391250456992538083614912707281821312923521104857612006431336352150365117006081874161

%N的奇数正除数的除数补码的四次方的和。

%H Robert Israel，n的表，a（n）表示n=1..10000</a>

%F a（n）=n^4*和{d|n，d<n，d奇数}1/d^4。

%F G.F.：和{k>=2}k^4*x^k/（1-x^（2*k））_伊利亚·古特科夫斯基，2023年5月14日

%F From _Amiram Eldar_，2023年10月13日：（开始）

%F a（n）=A051001（n）*A006519（n）^4-A000035（n）。

%F和{k=1..n}a（k）=c*n^5/5，其中c=31*zeta（5）/32=1.00452376。（结束）

%e a（10）=10^4*和{d|10，d<10，d奇数}1/d^4=10^4*（1/1^4+1/5^4）=10016。

%p f：=proc（n）局部m，d；

%pm:=n/2^padic:-ordp（n，2）；

%p加法（（n/d）^4，d=select（`<`，数字理论：-除数（m），n））

%p结束进程：映射（f，[$1..40]）；#_罗伯特·伊斯雷尔，2023年4月3日

%t A352050[n_]：=除数和[n，1/#^4&，#<n&&OddQ[#]&]n^4；阵列[A352050,50]（*_Paolo Xausa_，2023年8月9日*）

%t a[n_]：=除数Sigma[-4，n/2^整数指数[n，2]]*n^4-模[n，2]；数组[a，100]（*_Amiram Eldar_，2023年10月13日*）

%o（PARI）a（n）=n^4*sigma（n>>估值（n，2），-4）-n%2；\\_Amiram Eldar，2023年10月13日

%Y k=0..10:A091954（k=0），A352047（k=1），A35.248（k=2），A35249（k=3），这个序列（k=4），A35.2（k=5），A352（k=6），A35 253（k=7），A3 52054（k=8），A3.5 255（k=9），A35.2056（k=10），n的奇真除数的除数补的k次幂之和。

%Y参考A000035、A006519、A013663、A051001。

%K nonn，简单

%O 1,2号机组

%A _韦斯利·伊万·赫特，2022年3月1日