A111843
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)!.
Original entry on oeis.org
0, 1, 0, 3, 3, 0, 27, 9, 9, 0, 486, 81, 27, 27, 0, 7776, 1458, 243, 81, 81, 0, -2423196, 23328, 4374, 729, 243, 243, 0, -97338996, -7269588, 69984, 13122, 2187, 729, 729, 0, 5883879500784, -292016988, -21808764, 209952, 39366, 6561, 2187, 2187, 0
Offset: 0
Matrix log of A111840, with factorial denominators, begins:
0;
1/1!, 0;
3/2!, 3/1!, 0;
27/3!, 9/2!, 9/1!, 0;
486/4!, 81/3!, 27/2!, 27/1!, 0;
7776/5!, 1458/4!, 243/3!, 81/2!, 81/1!, 0;
-2423196/6!, 23328/5!, 4374/4!, 729/3!, 243/2!, 243/1!, 0;
A111844
Column 0 of the matrix logarithm (A111843) of triangle A111840, which shifts columns left and up under matrix cube; these terms are the result of multiplying the element in row n by n!.
Original entry on oeis.org
0, 1, 3, 27, 486, 7776, -2423196, -97338996, 5883879500784, 548540050402080, -1737375315124971951360, -405928706169160555680960, 60788545124934395018363657569920, 36207408592259278909089966337224960, -237458310218887960183820317532070376189904640
Offset: 0
E.g.f. A(x) = x + 3/2!*x^2 + 27/3!*x^3 + 486/4!*x^4 + 7776/5!*x^5
- 2423196/6!*x^6 - 97338996/7!*x^7 +...
where A(x) satisfies:
x = A(x) - A(x)*A(3*x)/2! + A(x)*A(3*x)*A(3^2*x)/3!
- A(x)*A(3*x)*A(3^2*x)*A(3^3*x)/4! + ...
also:
Let G(x) be the g.f. of A111841 (column 0 of A111840), then
G(x) = 1 + x + 3*x^2 + 18*x^3 + 216*x^4 + 5589*x^5 + 336555*x^6 +...
= 1 + A(x) + A(x)*A(3*x)/2! + A(x)*A(3*x)*A(3^2*x)/3!
+ A(x)*A(3*x)*A(3^2*x)*A(3^3*x)/4! +...
-
{a(n,q=3)=local(A=Mat(1),B);if(n<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!*B[n+1,1]))}
A111841
Number of partitions of 3^n-1 into powers of 3, also equals column 0 of triangle A111840, which shifts columns left and up under matrix cube.
Original entry on oeis.org
1, 1, 3, 18, 216, 5589, 336555, 49768101, 18707873562, 18299531019402, 47379925800261099, 328983441917303863134, 6190598463101580564238419, 318441251661562459898972204796, 45106336219710244780433937129788943
Offset: 0
-
{a(n,q=3)=local(A=Mat(1),B);if(n<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);return(A[n+1,1]))}
A111842
Row sums of triangle A111840, which shifts columns left and up under matrix cube.
Original entry on oeis.org
1, 2, 7, 46, 595, 16444, 1048303, 162728110, 63746277967, 64594795730680, 172419318632651104, 1229463017642626881490, 23684690483668583872503679, 1244115601652916934000237966330, 179585081405174505374545193721101377
Offset: 0
-
{a(n,q=3)=local(A=Mat(1),B);if(n<0,0, for(m=1,n+2,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); return(sum(k=0,n,A[n+1,k+1])))}
A111815
Matrix log of triangle A078122, 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)!.
Original entry on oeis.org
0, 1, 0, -1, 3, 0, -3, -3, 9, 0, 150, -9, -9, 27, 0, 1236, 450, -27, -27, 81, 0, -2555748, 3708, 1350, -81, -81, 243, 0, -64342116, -7667244, 11124, 4050, -243, -243, 729, 0, 5885700899760, -193026348, -23001732, 33372, 12150, -729, -729, 2187, 0
Offset: 0
Matrix log of A078122, with factorial denominators, begins:
0;
1/1!, 0;
-1/2!, 3/1!, 0;
-3/3!, -3/2!, 9/1!, 0;
150/4!, -9/3!, -9/2!, 27/1!, 0;
1236/5!, 450/4!, -27/3!, -27/2!, 81/1!, 0;
-2555748/6!, 3708/5!, 1350/4!, -81/3!, -81/2!, 243/1!, 0; ...
A111845
Triangle P, read by rows, that satisfies [P^4](n,k) = P(n+1,k+1) for n>=k>=0, also [P^(4*m)](n,k) = [P^m](n+1,k+1) for all m, where [P^m](n,k) denotes the element at row n, column k, of the matrix power m of P, with P(k,k)=1 and P(k+1,1)=P(k+1,0) for k>=0.
Original entry on oeis.org
1, 1, 1, 4, 4, 1, 40, 40, 16, 1, 1040, 1040, 544, 64, 1, 78240, 78240, 48960, 8320, 256, 1, 18504256, 18504256, 13110400, 2878720, 131584, 1024, 1, 14463224448, 14463224448, 11192599808, 2982187520, 180270080, 2099200, 4096, 1
Offset: 0
Let q=4; the g.f. of column k of matrix power P^m is:
1 + (m*q^k)*L(x) + (m*q^k)^2/2!*L(x)*L(q*x) +
(m*q^k)^3/3!*L(x)*L(q*x)*L(q^2*x) +
(m*q^k)^4/4!*L(x)*L(q*x)*L(q^2*x)*L(q^3*x) + ...
where L(x) satisfies:
x = L(x) - L(x)*L(q*x)/2! + L(x)*L(q*x)*L(q^2*x)/3! +- ...
and L(x) = x + 4/2!*x^2 + 56/3!*x^3 + 1728/4!*x^4 + ... (A111849).
Thus the g.f. of column 0 of matrix power P^m is:
1 + m*L(x) + m^2/2!*L(x)*L(4*x) + m^3/3!*L(x)*L(4*x)*L(4^2*x) +
m^4/4!*L(x)*L(4*x)*L(4^2*x)*L(4^3*x) + ...
Triangle P begins:
1;
1,1;
4,4,1;
40,40,16,1;
1040,1040,544,64,1;
78240,78240,48960,8320,256,1;
18504256,18504256,13110400,2878720,131584,1024,1; ...
where P^4 shifts columns left and up one place:
1;
4,1;
40,16,1;
1040,544,64,1;
78240,48960,8320,256,1; ...
-
{P(n,k,q=4) = my(A=Mat(1),B);if(nPaul D. Hanna, Jul 11 2025):
for(n=0,10, for(k=0,n, print1(P(n,k),", ")); print(""))
Showing 1-6 of 6 results.
Comments