%I#20 2017年1月6日16:01:45
%S 1,5,2510138513974921169495740919190963477720814776775873,
%电话:21921941705487932259952857210325772292237893264134169,
%电话:229546639737235077632127512682165713919106105223590049714169901312790521817681693157
%N矩阵行列式的绝对值,其中元素a(i,j)=如果i+j>N,则2*(i+j-N)-1,否则2*(N+1-i-j)。
%C主对角线是1,左上角元素增加较大的偶数,右下角元素增加较多的奇数。
%H Harvey P.Dale,n表,n=1..1000的a(n)</a>
%H<a href=“/index/Rec#order_03”>带常系数的线性重复出现的索引条目,签名(5,-3,-9)。
%F a(n)=1/12*[(4n-1)3^n-3(-1)^n]。
%F a(n)=5*a(n-1)-3*a(n-2)-9*a(n-3)。通用格式:x*(3*x^2+1)/((x+1)*(3*1)^2)_科林·巴克,2014年9月28日
%e a(5)=测定值{8 6 4 2 1/6 4 2 1 3/4 2 1 3 5/2 1 3 5 7/1 3 5 7 9}=385。
%tf[i_,j_,n_]:=块[{a=2*(i+j)-2*n-1},如果[i+j<=n,a=Abs[a-1]];a] ;表[Abs[Det[表[f[i,j,n],{i,1,n},{j,1,n}]],{n,1,24}]
%t线性递归[{5,-3,-9},{1,5,25},30](*哈维·P·戴尔,2017年1月6日*)
%o(PARI)Vec(x*(3*x^2+1)/((x+1)*(3*1)^2)+o(x^100))\\科林·巴克,2014年9月28日
%o(PARI)a(n)=abs(matdet(矩阵(n,n,i,j,如果(i+j>n,2*(i+j-n)-1,2*(n+1-i-j))));\\_米歇尔·马库斯,2014年9月29日
%K nonn,简单
%O 1,2号机组
%A _罗伯特·G·威尔逊,2003年5月7日
%E来自Colin Barker的更多条款,2014年9月28日
%E定义由Michel Marcus修正,2014年10月1日
|