%I#28 2023年11月26日03:11:02

%S 1,1,1,1,1,1,8,1,1,1,1,14,8,11,1,1,2,20,1,81,1,1,,1,26,14,1,8,1,1,1,32,1，

%T 1,1,8,1,38,20,14,1,1,8,44,1,1,1,1,57,26,1,14,1,56,8,20,1,1,62，

%U 32,8,1,14,1,68,1,1,8,1,74,38,26,20,8,14,80,1,1

%N a（N）=和{d|N，d==1（mod 6）}d。

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

%F G.F.：和{k>=0}（6*k+1）*x^（6xk+1）/（1-x^_伊利亚·古特科夫斯基，2017年3月21日

%F G.F.：和{n>=1}x^n*（1+5*x^（6*n））/（1-x^_彼得·巴拉（Peter Bala），2021年12月19日

%F和{k=1..n}a（k）=c*n^2+O（n*log（n）），其中c=Pi^2/72=0.13777…（A086729）_阿米拉姆·埃尔达尔，2023年11月26日

%t表[Sum[If[Mod[d，6]==1，d，0]，{d，Divisors[n]}]，{n，80}]（*_Indranil Ghosh_，2017年3月21日*）

%o（PARI）用于（n=1,82，print1（sumdiv（n，d，if（Mod（d，6）==1,d，0）），“，”））\\_Indranil Ghosh，2017年3月21日

%o（Python）

%o来自sympy导入除数

%o def a（n）：返回和（如果d%6==1，则[d代表除数（n）中的d）#_Indranil Ghosh_，2017年3月21日

%Y参考A086729，A109701。

%Y参考Sum_{d|n，d==1（mod k）}d:A000593（k=2），A078181（k=3），P050449（k=4），A284097（k=5），该序列（k=6），A28.4099（k=7），A284 100（k=8）。

%K nonn，简单

%O 1,7型

%A _Seiichi Manyama，2017年3月20日