%I#10 2021年4月13日04:43:24
%S 1,1,3,12,39138469160354271837261869207909653722328039,
%电话77622662582614285749969284171598940027872310428088510234808334,
%电话:3369254724911073171784363561071175119186048756139023506274341276156548741685086306917136012008938158
%N乘积{k>=1}(1+x^k)^(3^(k-1))的展开。
%F a(n)~exp(2*sqrt(n/3)-1/6-c/3)*3^(n-1/4)/(2*m2(Pi)*n^(3/4)),其中c=Sum_{j>=2}(-1)^j/(j*(3^)-1))_瓦茨拉夫·科特索维奇,2021年4月13日
%ph:=proc(n,i)选项记忆`如果`(n=0,1,`如果`(i<1,0,
%p加(h(n-i*j,i-1)*二项式(3^(i-1),j),j=0..n/i))
%p端:
%价格:=n->小时(n$2):
%p序列(a(n),n=0..28);#_阿洛伊斯·海因茨,2021年4月12日
%t nmax=28;系数列表[系列[积[(1+x^k)^(3^(k-1))),{k,1,nmax}],{x,0,nmax{],x]
%ta[n_]:=a[n]=如果[n==0,1,(1/n)Sum[Sum[(-1)^(k/d+1)d3^(d-1),{d,除数[k]}]a[n-k],{k,1,n}]];表[a[n],{n,0,28}]
%o(PARI)seq(n)={Vec(prod(k=1,n,(1+x^k+o(x*x^n))^(3^(k-1)))}\\ Andrew Howroyd_,2021年4月12日
%Y请参阅A098407、A104460、A256142、A343361、A343365、A343、A3362、A34336、A343和A343364。
%K nonn公司
%0、3
%A _Ilya Gutkovskiy,2021年4月12日
|