%我#15 2021年4月22日04:44:59

%S 1,1,2,2,5,2,25,2,34,19129,21046,27424501597,244254,227517，

%电话332129967，2187275757120039018560854850，215404511，23524578，

%电话：2260209262157518863353855285,26334659820448951255304727,2185493291001,21285159458520352972038875,26006303471178

%N将N的组成（有序分区）转换为N的无平方因子的数目。

%H Alois P.Heinz，n的表格，n=0..2000的a（n）</a>

%H Eric Weisstein的数学世界，<a href=“http://mathworld.wolfram.com/Squarefree.html“>无方块</a>

%H<a href=“/index/Com#comp”>与合成相关的序列的索引条目</a>

%F a（n）=[x^n]1/（1-求和{d|n，|mu（d）|=1}x^d），其中mu（d）是Moebius函数（A008683）。

%如果n是素数，则F a（n）=2。

%e a（4）=5，因为4有3个除数{1，2，4}，其中2个是无平方的{1，2}，所以我们有[2，2]，[2，1，1]，[1，2，2]和[1，1,1,1]。

%p（数字理论）：

%p a:=proc（n）选项记忆；局部b，l；

%p l，b:=选择（issqrfree，除数（n）），

%p过程（m）选项记忆`如果`（m=0，1，

%p加（`if`（j>m，0，b（m-j）），j=l）

%p端；b（n）

%p结束：

%p序列（a（n），n=0..50）；#_Alois P.Heinz_，2017年3月30日

%t表[d=除数[n]；系数[级数[1/（1-总和[MoebiusMu[d[[k]]^2 x^d[[k]]，{k，长度[d]}]），{x，0，n}]，x，n]，{n，0，48}]

%o（Python）

%o来自sympy导入除数

%o来自sympy.theory.factor导入核心

%o从sympy.core.cache导入缓存

%o@缓存

%o定义a（n）：

%o l=[x代表除数（n）中的x，如果核（x）==x]

%o@缓存

%o def b（m）：如果m==0，则返回1，否则求和（如果j<=m，则b（m-j）用于l中的j）

%o返回b（n）

%o打印（[a（n）代表范围（51）中的n）]#_Indranil Ghosh，2017年8月1日，在Maple代码之后

%Y参见A005117、A008683、A100346、A225244、A225145、A280194。

%K nonn公司

%0、3

%A _Ilya Gutkovskiy_，2017年3月27日