%I#32 2023年11月5日15:55:42
%S 1,8,26,641262083425127031008133016642198273632764096，
%电话：4914562468588064889210641215751175841898021888，
%电话：2439026208297903276834580393124309244992506545486457148
%N a（N）=和{d|N，N/d=1模4}d^3-和{d| N，N/d=3模4{d^3。
%C相乘，因为它是A000578=n^3和A101455=[10-10-1…]的Dirichlet卷积，两者都是相乘的_Christian G.Bower，2005年5月17日
%H Seiichi Manyama，n=1..10000的n，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.：和{n>=1}n^3*x^n/（1+x^（2*n））_Vladeta Jovovic_，2002年10月16日
%F From _Amiram Eldar_，2023年11月4日：（开始）
%F与a（p^e）的乘积=（p^（3*e+3）-A101455（p）^（e+1））/（p^3-A101455（p））。
%F和{k=1..n}a（k）~c*n^4/4，其中c=A175572。（结束）
%t最大值=40；s=总和[n^3*x^（n-1）/（1+x^，2*n）），{n，1，max}]+O[x]^max；系数列表[s，x]（*_Jean-François Alcover_，2015年12月2日，在_Vladeta Jovovic_*之后）
%t s[n_]：=如果[OddQ[n]，（-1）^（（n-1）/2），0]；（*A101455*）
%tf[p，e]：=（p^（3*e+3）-s[p]^（e+1））/（p^3-s[p]）；a[1]＝1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2023年11月4日*）
%o（PARI）a（n）=sumdiv（n，d，d^3*（（（n/d）%4）==1））-sumdiv_米歇尔·马库斯，2015年2月16日
%Y参考A000578、A101455、A175572。
%Y Glaisher的E'_i（i=0..12）：A002654、A050469、A050470，该序列，A050468、A321829、A321830、A3218.31、A32183、A321933、A321844、A321735、A321866。
%K非n，简单，多
%O 1,2号机组
%A _N.J.A.Sloane，1999年12月23日
%2010年7月15日，R.J.Mathar_将E偏移从0更改为1
