|
%I#19 2023年11月5日09:03:42
%第1102459048104857697665626604651522824752481073741824页，
%电话3486725353100000010242593742460061916315648137858491850，
%电话：289254653952576640684048109951162777620159939004503570406761472613106625780010240001048576
%N a（N）=Sum_{d|N，N/d==1模4}d^10-和{d|N，N/d=3模4{d^10。
%H Seiichi Manyama，n的表格，n=1..10000的a（n）</a>
%H J.W.L.Glaisher，<a href=“https://books.google.com/books？id=bLs9AQAAMAAAJ&pg=RA1-PA1“>关于数字表示为二、四、六、八、十和十二个方块之和，Quart.J.Math.38（1907），1-62（见第4页和第8页）。
%H<a href=“/index/Ge#Glaisher”>为Glaisher</a>提到的序列索引条目。
%F G.F.：和{k>=1}k^10*x^k/（1+x^（2*k））_伊利亚·古特科夫斯基，2018年11月26日
%F From _Amiram Eldar_，2023年11月4日：（开始）
%F与a（p^e）的乘积=（p^（10*e+10）-A101455（p）^（e+1））/（p^10-A101455（p））。
%F Sum_{k=1..n}a（k）~c*n^11/11，其中c=beta（11）=50521*Pi^11/14863564800=0.99999 4374973…，beta是Dirichlet beta函数。（结束）
%t s[n_，r_]：=除数和[n，#^10&，Mod[n/#，4]==r&]；a[n]：=s[n，1]-s[n，3]；阵列[a，30]（*_Amiram Eldar_，2018年11月26日*）
%t s[n_]：=如果[OddQ[n]，（-1）^（（n-1）/2），0]；（*A101455*）
%tf[p，e]：=（p^（10*e+10）-s[p]^（e+1））/（p^10-s[p]）；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2023年11月4日*）
%o（PARI）适用（a（n）=sumdiv（n，d，if（bittest（n，0），（2-n，d%4）*d^10）），[1..30]）\\ M.F.Hasler_，2018年11月26日
%Y参考A101455。
%Y参考A321807-A321836了解类似序列。
%Y Glaisher的E'_i（i=0..12）：A002654、A050469、A050470、A050441、A050488、A321829、A321830、A3218.31、A32183、A321935、A321863，该序列。
%K nonn，简单，多
%O 1,2号机组
%A _N.J.A.Sloane，2018年11月24日
| 