%I#11 2018年8月27日10:04:22
%S 1,1,7,19,711735831443435106482892971159184740445626,
%电话:111012226383286369490148701943503162780465028185556696,
%电话:4199614995078558021214717784727971847104122306982287688652979799776871862107974178843232302695301
%N产品{k>=1}1/(1-k*x^k)^sigma(k)的展开,其中sigma=A000203。
%H Vaclav Kotesovec,n的表格,n=0..5920的a(n)</a>
%F a(n)~c*n^3*3^(n/3),其中
%F c=280631952508395331838833549352368682635.581151020…如果mod(n,3)=0
%F c=280631952508395331838833549352368682635.059082354…如果mod(n,3)=1
%F c=2806319525083953312838833549352333682635.088610121…如果mod(n,3)=2
%F在封闭形式中,c=(Product_{k>=4}((1-k/3^(k/3))^(-sigma(k)))/(1/3)^3*(-3+(-3)^(2/3))))*(1-2*(-1/3)^(2/3))^3)
%t nmax=40;系数列表[系列[乘积[1/(1-k*x^k)^DivisorSigma[1,k],{k,1,nmax}],{x,0,nmax{],x]
%t最大值nmax=40;s=1-x;Do[s*=和[二项式[DivisorSigma[1,k],j]*(-1)^j*k^j*x^(j*k),{j,0,nmax/k}];s=展开[s];s=取[s,Min[nmax+1,指数[s,x]+1,长度[s]];,{k,2,nmax}];系数列表[系列[1/s,{x,0,nmax}],x]
%Y参见A000203、A006906、A266941、A318415、A006171、A061256、A318484。
%K非n
%0、3
%A _卡拉夫·科特索维奇,2018年8月27日
|