%I#32 2022年10月4日12:36:21

%S 1,3,4,7,6,12,8,7,13,18,12,28,14,24,7,18,39,20,42,32,36,24,28,31，

%T 42,40,56,30,72,32,7,48,54,48,91,38,60,56,42,42,96,44,84,78,72,48,28，

%U 57,93,72,98,54120,72,56,80,90,60168,62,96104,7,84144,68

%N不能被8整除的N的除数之和。

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

%计算公式：和{k>=1}k*x^k/（1-x^k）-8*k*x*（8*k）/（1-x^（8*k））_伊利亚·古特科夫斯基，2017年3月25日

%F与a（2^e）=7相乘，如果e>=3，则a（p^e）=（p^（e+1）-1）/（p-1）相乘。-_Amiram Eldar，2020年9月17日

%F和{k=1..n}a（k）~（7*Pi^2/96）*n^2.-_Amiram Eldar，2022年10月4日

%t表[Sum[Boole[Mod[d，8]>0]d，{d，Divisors[n]}]，{n，100}]（*_Indranil Ghosh_，2017年3月25日*）

%t表[Total[DeleteCases[Divisions[n]，_？（可分割[#，8]&）]]，{n，120}]（*_哈维·P·戴尔，2018年3月18日*）

%t f[p_，e_]：=如果[p==2&&e>=3,7，（p^（e+1）-1）/（p-1）]；a[1]=1；a[n_]：=次数@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2020年9月17日*）

%o（PARI）针对（n=1100，print1（sumdiv（n，d，（d%8）>0）*d），“，”））\\_Indranil Ghosh_，2017年3月25日

%o（Python）

%o来自sympy导入除数

%o打印（[sum（[i for i in divisors（n）in divisor（n）if i%8]）for n in range（1101）]）#_Indranil Ghosh_，2017年3月25日

%Y Cf.不可被k整除的n的除数之和：A046913（k=3），A046897（k=4），A116073（k=5），A284326（k=6），A113957（k=7），这个序列（k=8），A116607（k=9），A284 344（k=10）。

%K nonn，多个

%O 1,2号机组

%A _Seiichi Manyama，2017年3月25日

%E关键词：由Andrew Howroyd_2018年7月20日添加的mult