A111843 - OEIS

%I #7 Jun 13 2017 22:38:28

%S 0,1,0,3,3,0,27,9,9,0,486,81,27,27,0,7776,1458,243,81,81,0,-2423196,

%T 23328,4374,729,243,243,0,-97338996,-7269588,69984,13122,2187,729,729,

%U 0,5883879500784,-292016988,-21808764,209952,39366,6561,2187,2187,0

%N Matrix log of triangle A111840, which shifts columns left and up under matrix cube; these terms are the result of multiplying each element in row n and column k by (n-k)!.

%C Column k equals 3^k multiplied by column 0 (A111844) when ignoring zeros above the diagonal.

%F T(n, k) = 3^k*T(n-k, 0) = 3^k*A111844(n-k) for n>=k>=0.

%e Matrix log of A111840, with factorial denominators, begins:

%e 0;

%e 1/1!, 0;

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

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

%e 486/4!, 81/3!, 27/2!, 27/1!, 0;

%e 7776/5!, 1458/4!, 243/3!, 81/2!, 81/1!, 0;

%e -2423196/6!, 23328/5!, 4374/4!, 729/3!, 243/2!, 243/1!, 0;

%o (PARI) T(n,k,q=3)=local(A=Mat(1),B);if(n<k || k<0,0, for(m=1,n+1,B=matrix(m,m);for(i=1,m, for(j=1,i, if(j==i,B[i,j]=1,if(j==1,B[i,j]=(A^q)[i-1,1], B[i,j]=(A^q)[i-1,j-1]));));A=B); B=sum(i=1,#A,-(A^0-A)^i/i);return((n-k)!*B[n+1,k+1]))

%Y Cf. A111840 (triangle), A111844 (column 0), A111815 (variant), A111941 (q=-1), A111810 (q=2), A111848 (q=4).

%K frac,sign,tabl

%O 0,4

%A _Paul D. Hanna_, Aug 23 2005