%I#12 2023年3月5日13:26:01
%S 1,1,-1,12,-16,12160,-3312381,-16048000,-101376001603680,-10496,1,
%电话266716800000,-47670336000092708406000,-1022881200307505,-1,
%电话:18631342033920000,-3499358557504000078981336366912000,-12427237734400750409713900,-9316560,1
%逆希尔伯特矩阵特征多项式的三角函数。
%C按行读取的三角形:对于0<=k<=n,T(n,k)是λ^k在det中的系数(H^(-1)-lambda I),其中H是n x n Hilbert矩阵。
%C行总和为:1、0、-3、-772、-2496415、-11830727696、-8588275706265059、-9728358620103886374316、-17352032030130344466515679407359、-4898477775570499454780372858733881836257416、-219554569246037985920541114453120558720536422853379
%H Robert Israel,n表,n=0..902时的a(n)
%H Eric Weisstein的数学世界,<a href=“http://mathworld.wolfram.com/HilbertMatrix.html“>Hilbert矩阵</a>
%F t(n,m)=系数表[Characteristic Polynomial[Inverse[HilbertMatrix[n]],x],x'
%e{1},
%e{1,-1},
%e{12,-16,1},
%电子邮箱{2160,-3312,381,-1},
%电子邮箱{6048000,-10137600,1603680,-10496,1},
%电子邮箱{266716800000,-476703360000,92708406000,-1022881200,307505,-1},
%电子邮箱:{18631342033920000,-34993585575040000,78981336366912000,-1242627237734400,750409713900,-9316560,1}
%pf:=proc(n)使用线性代数;
%p局部λ,p,j;
%p p:=特征多项式(Hilbert矩阵(n),λ)/行列式(HilbertMatrix(n));
%p seq(系数(p,λ,n-j),j=0..n);
%p端程序:
%p序列(f(n),n=0..10);#_罗伯特·伊斯雷尔(Robert Israel),2016年10月5日
%t<<线性代数`矩阵运算`;a=连接[{{1}},表[CoefficientList[CharacteristicPolynomial[Inverse[HilbertMatrix[n]],x],x]{n,1,10}]];压扁[a]
%Y参考A005249。
%K tabl,签名
%0、4
%A _Roger L.Bagula,2007年12月14日
%E编辑:罗贝尔·伊斯雷尔,2016年10月5日
|