%I#17 2023年11月5日09:03:31
%S 1128218616384781262798088235422097152478078310000128，
%电话1948717035815424627481054133761707833436286435456410338674，
%电话：611940224893871738128001638418002628122494357760340482544645843742726103593751
%N a（N）=Sum_{d|N，N/d==1模4}d^7-和{d|N，N/d=3模4{d^7。
%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^7*x^k/（1+x^（2*k））_伊利亚·古特科夫斯基，2018年11月26日
%F与a（p^e）相乘=四舍五入（p^（7e+7）/（p^7+p%4-2）），其中p%4是p模4的余数。（在A321833中跟随R.Israel）-M.F.Hasler_，2018年11月26日
%F和{k=1..n}a（k）~c*n^8/8，其中c=A258815.-_Amiram Eldar，2023年11月4日
%t s[n_，r_]：=除数和[n，#^7&，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^（7*e+7）-s[p]^（e+1））/（p^7-s[p]）；a[1]＝1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2023年11月4日*）
%o（PARI）适用（A321831（n）=退单（适用（f->f[1]^（7*f[2]+7）//（f[1]_7+f[1]%4-2），列（系数（n））），[1..30]）\\M.f.哈斯勒，2018年11月26日
%Y参见A101455、A258815。
%Y参考A321543-A321565、A321807-A321836了解类似序列。
%Y Glaisher的E'_i（i=0..12）：A002654、A050469、A050470、A050441、A050488、A321829、A321830，该序列，A321832、A3218.33、A321844、A321935、A321866。
%K nonn，简单，多
%O 1,2号机组
%A _N.J.A.Sloane，2018年11月24日
