登录
n×n素元矩阵特征多项式系数的三角序列:M=A.B.A^(-1);(A(3)为单数):示例;A(4)={{2,3,5,7,11},{3,5,17,11,13},}5,7,7,11,13,17},[7,11,11,13,19},[11,13,17,19,23}}B(4)={3,5,7,11,13},[25,7],11,13,17},[2],{7,11,13,13,17,17,18},f11,13,17,19},#11,17,13,16,17,12}。
0

%I#3 2012年3月30日17:34:26

%S 1,3,-1,-4,-10,1,12,44,23,-1,0576,-864,-288316,71,-1,-7104,-608,

%电话58001168,-670,-108,14524844096,-19712,-21712161768151,-1,

%U 450432135424,-346176,-66560690087344,-3204,12240512,-5071104,-14220802168128188672,-212928,-109245960265,-1

%由N×N素元矩阵的特征多项式系数得到的N三角序列:M=A.B.A^(-1);(A(3)为单数):示例;A(4)={{2,3,5,7,11},{3,5,17,11,13},}5,7,7,11,13,17},[7,11,11,13,19},[11,13,17,19,23}}B(4)={3,5,7,11,13},[25,7],11,13,17},[2],{7,11,13,13,17,17,18},f11,13,17,19},#11,17,13,16,17,12}。

%C行和是(对于奇异矩阵级别,我将其置零):{1,2,-13,78,0,-190,-1521,49954,246065,-2113500,9798535,74936314,…}。

%C据我所知,这种方法的唯一其他用途是推导A086515。

%对于n<=m,=d:A(d)=素数[n+m+1]B(d)=素数[n+m+2]p(x,d)=特征多项式[A.B.A^(-1)]out_n,m=系数(p(x、d))

%e{1},

%e{3,-1},

%e{-4,-10,1},

%e{12,44,23,-1},

%e{0},

%电子邮箱{576,-864,-288,316,71,-1},

%电子邮箱{-7104,-608,5800,1168,-670,-108,1},

%电子邮箱{45248,44096,-19712,-21712,116,1768151,-1},

%电子邮箱{450432135424,-346176,-66560,69008,7344,-3204,1},

%电子邮箱{2240512,-5071104,-1422080,2168128,188672,-212928,-10924,5960,265,-1},

%电子邮箱{5292544,-27469312,33984128,4845952,-6953248,-363232,455688,16064,-9714,-336,1},

%电子邮箱{-88076288、158813184、142065920、-134528512、-2576128、18750592、1438864、-925536、-41344、15148、415、-1}

%t a=表[素数[n],{n,1500}];aa[d_]:=表[a[[n+m+1]],{n,0,d},{m,0,d}];bb[d_]:=表[a[[n+m+2]],{n,0,d},{m,0,d}];M[d_]:=aa[d].bb[d]。矩阵幂[aa[d],-1]g=Join[{1},Table[If[d==3,0,CharacteristicPolynomial[M[d]、x]],{d,0,10}]];a0=连接[{{1}},表[系数表[If[d==3,0,特征多项式[M[d],x],{d,0,10}]];压扁[{{1}、{3、-1}、}-4、-10、1},{12、44、23、-1}、{0},}576、-864、-288、316、71、-1},{-7104、-608、5800、1168、-670、-108、1}、{45248、44096、-19712、-21712、116、1768、151、-1}.、{450432、135424、-346176、-66560、69008、7344、-3204、-204,1},{2240512,-5071104,-1422080,2168128,188672,-212928,-10924,5960,265,-1}, {5292544, -27469312, 33984128, 4845952, -6953248, -363232, 455688, 16064, -9714, -336, 1}, {-88076288, 158813184, 142065920, -134528512, -22576128, 18750592, 1438864, -925536, -41344, 15148, 415, -1}}]

%Y参考A086515。

%K tabf,未编辑,签名

%O 1,2号机组

%A _Roger L.Bagula,2008年4月14日