%I#8 2017年6月13日22:33:49

%S 0,1,0，-1,3,0，-3，-3,9,0150，-9，-9,27,01236450，-27，-27,81,0，-2555748，

%电话：37081350，-81，-81243,0，-64342116，-7667244111244050，-243，-243729，

%电话：05885700899760，-193026348，-230017323337212150，-729，-7292187,0

%N三角形A078122的矩阵对数，在矩阵立方体下左右移动列；这些项是第n行和第k列中的每个元素乘以（n-k）！的结果！。

%C忽略对角线上方的零时，k列等于3^k乘以0列（A111816）。

%对于n>=k>=0，F T（n，k）=3^k*T（n-k，0）=A111816（n-k）。

%e带有阶乘分母的A078122矩阵日志开始于：

%e 0；

%e 1/1！，0;

%e-1/2！，3/1!, 0;

%e-3/3-3/2!, 9/1!, 0;

%e 150/4-9/3!, -9/2！，27/1!, 0;

%e 1236/5！，450/4!, -27/3!, -27/2!, 81/1！，0;

%e-2555748/6！，3708/5!, 1350/4!, -81/3!, -81/2!, 243/1!, 0; ...

%o（PARI）T（n，k，q=3）=局部（A=Mat（1），B）；如果（n<k | | k<0,0，对于（m=1，n+1，B=矩阵（m，m））；对于（i=1，m，对于（j=1，i，如果（j==i|j==1，B[i，j]=1，B[i，j]=（A^q）[i-1，j-1]）；）；A=B）；B=总和（i=1，#A，-（A^0-A）^i/i）；返回（n-k）*B（n+1，k+1））

%Y参见A078122、A111816（第0列）、A111840（变体）；对数矩阵：A110504（q=-1）、A111813。

%K压裂、签名、表格

%0、5

%A Gottfried Helms_和Paul D.Hanna，2005年8月22日