%I#10 2021年9月15日05:11:30

%S 1,1,3,10,25,70182476122031227883197944934012223730114，

%电话：73792317995974369204105638002544137761048713145988775，

%电话：3479817138269219219593637784629903905109117574322562950459601658313611407922150373287503982757660419301781452975346

%N乘积{k>=1}（1+x^k）^（2^k-1）的展开。

%C A081362和A102866的卷积。

%C A000225的称重转换。

%H Vaclav Kotesovec，n=0..3000时的n表，a（n）</a>

%H N.J.A.Sloane，转换</a>

%F G.F.：exp（总和{k>=1}（-1）^（k+1）*x^k/（k*（1-x^k）*（1-2*x^k。

%F a（n）~c*exp（2*sqrt（n）-1/2）*2^（n-1）/（A079555*sqort（Pi）*n^（3/4）），其中c=exp（总和{k>=2}（-1）^（k-1）/

%p a：=级数（mul（（1+x^k）^（2^k-1），k=1..100），x=0,33）：seq（系数（a，x，n），n=0..32）；#_保罗·拉瓦（Paolo P.Lava），2019年4月2日

%t nmax=32；系数列表[系列[乘积[（1+x^k）^（2^k-1），{k，1，nmax}]，{x，0，nmax{]，x]

%t nmax=32；系数列表[级数[Exp[Sum[（-1）^（k+1）x^k/（k（1-x^k）（1-2 x^k

%ta[n_]：=a[n]=如果[n==0，1，Sum[Sum[（-1）^（k/d+1）d（2^d-1），{d，除数[k]}]a[n-k]，{k，1，n}]/n]；表[a[n]，{n，0，32}]

%Y参见A000225、A007070、A081362、A098407、A102866、A319918。

%K nonn公司

%0、3

%A _Ilya Gutkovskiy_，2018年10月1日