%I#14 2021年7月15日15:05:50
%S 1,3,4618991634762533359064123698602473269931755,
%电话137981705682772010691509081194746281113779885682143602440,
%电话:151890124741061619463551026519936538762415715172280581764069557390720735010689620
%N-log(总和{N>=0}(-x)^N/N^3)=和{n>=1}a(n)*x^n/n^三。
%F等于三角形T(n,k)=(-1)^(n-k)*C(n,k)^3矩阵对数的第0列。
%F a(n)=-(-1)^n+(1/n)*和{k=1..n-1}(-1)*(n-k-1)*二项式(n,k)^3*k*a(k).-_伊利亚·古特科夫斯基,2021年7月15日
%e L(x)=-log(1-x+x^2/2!^3-x^3/3!^3+x^4/4!^3-x^5/5!^3+-…)
%e其中
%e L(x)=x+3*x^2/2^3+46*x^3/3^3+1899*x^4/4^3+163476*x^5/5^3+。。。
%e备用生成方法。
%e A181543(n,k)=C(n,k)^3的签名版本开始于:
%e 1;
%e 1,1;
%e 1、8、1;
%e 1、27、27、1;
%e 1、64、216、64、1;
%e 1、125、1000、1000、125、1。。。
%e三角形A181543的矩阵对数开始于:
%e 0;
%e 1,0;
%e-3、8、0;
%e 46,-81,27,0;
%e-1899、2944、-648、64、0;
%e 163476、-237375、46000、-3000、125、0。。。
%e,其中该序列(有符号)位于列0中。
%o(PARI){a(n)=n!^3*polcoeff(-log(总和(m=0,n,(-x)^m/m!^3)+x*o(x^n)),n)}
%o(PARI)/*作为有符号三角形A181543的矩阵对数的第0列*/
%o{a(n)=局部(L,M=矩阵(n+1,n+1,r,c,如果(r>=c,(-1)^(r-c)*二项式(r-1,c-1)^3));
%o L=总和(n=1,#M,(M^0-M)^n/n);如果(n<0.0,L[n+1,1])}
%Y参考A002190(变型),A181543。
%K nonn公司
%O 1,2号机组
%A _保罗·D·汉纳,2011年7月26日
|