%I#34 2024年6月2日14:40:14

%S 1,1，-26,1126，-26，-342,1703126，-1330，-262198，-342，-3276，14914，

%电话703，-68581268892，-1330，-12166，-26157512198，-18980，-34224390，

%U-3276、-29790、1345804914、-4309270350654、-6858、-57148126689228892、-79506、-1330

%N a（N）=Sum_{d|N，d==1模4}d^3-和{d|N，d==3模4{d^3。

%C乘法，因为它是[10-3^305^30-7^3…]的逆Möbius变换，它是乘法的_Christian G.Bower，2005年5月18日

%H Seiichi Manyama，n的表，n的a（n）=1..10000</a>

%H J.W.L.Glaisher，<a href=“https://books.google.com/books？id=bLs9AQAAMAAAJ&amp；pg=RA1-PA1“>关于数字表示为二、四、六、八、十和十二个方块之和，Quart.J.Math.38（1907），1-62（见第4页和第8页）。

%H<a href=“/index/Ge#Glaisher”>为Glaisher</a>提到的序列索引条目。

%F a（n）=A050451（n）-A050454（n）。

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

%F与a（2^e）=1相乘，对于奇素数p，如果p==1（mod 4），则为（p^3）^（e+1）-1）/（p^3-1）；如果p==3_阿米拉姆·埃尔达尔，2023年9月27日

%F a（n）=和{d|n}d^3*sin（d*Pi/2）_Ridouane Oudra，2024年6月2日

%p A050459：=proc（n）局部a；a:=0；对于numtheory[除数]（n）中的d，如果d模4=1，那么a:=a+d^3；elif d mod 4=3，则a：=a-d ^3；结束条件：；结束do；a；结束进程：

%p序列（A050459（n），n=1..100）；#_R.J.Mathar_，2011年1月7日

%t s[n_，r_]：=除数和[n，#^3&，Mod[#，4]==r&]；a[n]：=s[n，1]-s[n，3]；阵列[a，30]（*_Amiram Eldar_，2018年12月6日*）

%t f[p_，e_]：=如果[模式[p，4]==1，（p^3）^（e+1）-1）/（p^3-1）；f[2，e_]：=1；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，60]（*_Amiram Eldar_，2023年9月27日*）

%A322143的Y列k=3。

%Y Glaisher的E_i（i=0..12）：A002654、A050457、A002173、A050449、A050466、A321821、A32182、A321H23、A3218024、A321855、A321866、A322827、A321888。

%Y参见A050451、A050454。

%K符号，简单，多，更改

%氧1,3

%A _N.J.A.Sloane，1999年12月23日