%I#41 2021年11月10日01:16:18

%S 1,2,2,3,2,9,2,4,3,4,4,2,17,2,4,6,5,2,19,2,10,4,4，4,2,37,3,4，12,2,45,2，

%电话：6,4,4,57,2,4,4,18,2,29,2,6,16,4,15,3,6,6,2,35,4,23,4,2301,2，

%U 4,6,7,4,28,2,6,4,19,2255,2,4,10,6,20,2,61

%调和平均数为整数的N的除数的非空子集的数目。

%H Chai Wah Wu，n的表格，n=1..3000的a（n）</a>

%H Eric Weisstein的数学世界，<a href=“http://mathworld.wolfram.com/HarmonicMean.html“>谐波平均值</a>

%H<a href=“/index/Di#divisors”>与数字除数相关的序列的索引项</a>

%e a（6）=9个子集：{1}、{2}、}、6}、2,6}，{3,6}和{1,2,3,6neneneep。

%t a[n_]：=计数[子集[除数[n]]，_？（长度[#]>0&&IntegerQ[HarmonicMean[#]]&）]；阵列[a，100]（*_Amiram Eldar_，2021年11月9日*）

%o（PARI）h（s，d）=#s/sum（k=1，#s，1/d[s[k]]）；

%o a（n）=my（d=除数（n），nb=0）；对于子集（#d，s，如果（#s&&（分母（h（s，d））==1，nb++））；编号；\\_米歇尔·马库斯（Michel Marcus），2020年12月15日

%o（Python）

%o从itertools导入组合

%o来自sympy导入除数

%o定义A339665（n）：

%o ds=元组（除数（n，生成器=True））

%o返回范围（1，len（ds）+1）内i的总和（如果n*i%总和（d）==0，则组合（ds，i）中d的总和（1）#_Chai Wah Wu_，2021年11月9日

%Y参见A027750、A100587、A339453、A339663、A339、664、A339和666。

%K nonn公司

%O 1,2号机组

%A _Ilya Gutkovskiy_，2020年12月11日