%I#41 2022年10月4日12:36:02

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

%电话：3,40,56,30,72,32,63,48,54,48,91,38,60,4,90,42,96,44,84,78,72,48124，

%U 57,93,72,7,54120,72120,80,90,60168,62,96104127,6144,68单位

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

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

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

%F与a（13^e）=1和a（p^e）=（p^（e+1）-1）/（p-1）相乘_Amiram Eldar，2020年9月17日

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

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

%t f[p_，e_]：=如果[p==13，1，（p^（e+1）-1）/（p-1）]；a[1]=1；a[n_]：=倍@@f@@FactorInteger[n]；阵列[a，100]（*_Amiram Eldar_，2020年9月17日*）

%o（PARI）a（n）=汇总（n，d，（（d%13）>0）*d）；\\_安德鲁·霍罗伊，2018年7月20日

%Y参考A000203，A133099。

%Y不能被k整除的n的除数之和：A046913（k=3），A046897（k=4），A116073（k=5），A284326（k=6），A113957（k=7），A284 341（k=8），A116607（k=9），A28.4344（k=10），这个序列（k=13），A227131（k=25）。

%K nonn，多个

%O 1,2号机组

%A _Seichi Manyama_，2017年3月29日

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