%I#14 2022年7月9日11:09:27

%S 1,3,4,3,6,12,8,3,4,18,12,16,14,24,24,3,18,12,20,22,32,36,24,16,6,42，

%电话：4.28,30,72,32,3,48,54,48,16,38,60,56,30,42,96,44,40,33,72,48,16，8,18，

%U 72,46,54,12,72,36,80,90,60108,62,96,41,84144,68,58,96144,72,16,74

%N不大于N的无平方核的N的除数之和。

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

%F a（n）>=A048250（n）。

%e 56=7*2^3的无平方核是14=7*2，56的除数<=14是{1,2,4,7,8,14}，因此a（56）=1+2+4+7+8+14=36。

%t a[n_]：=除数和[n，#&，#<=倍@@FactorInteger[n][[；；，1]]&]；阵列[a，100]（*_Amiram Eldar_，2022年7月9日*）

%o（PARI）

%o A007947（n）=因子回收（因子（n）[，1]）；\\此函数来自_Andrew Lelechenko_

%o A073181（n）=sumdiv（n，d，d*（d<=A007947（n）））；\\_Antti Karttunen，2017年9月10日，以Michel Marcus的A073183代码命名。

%Y参见A000203、A007947、A048250、A073180、A073183。

%K nonn，看

%O 1,2号机组

%A _Reinhard Zumkeller，2002年7月19日