A111821
Number of partitions of 4*5^n into powers of 5, also equals column 1 of triangle A111820, which shifts columns left and up under matrix 5th power.
Original entry on oeis.org
1, 5, 55, 2055, 291430, 165397680, 390075741430, 3927972221522680, 172358768282285194555, 33479766506261422878944555, 29150234311482124092454001991430
Offset: 0
-
a(n,q=5)=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 || j==1,B[i,j]=1,B[i,j]=(A^q)[i-1,j-1]);));A=B); return(A[n+2,2]))
A111825
Triangle P, read by rows, that satisfies [P^6](n,k) = P(n+1,k+1) for n>=k>=0, also [P^(6*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(0,k)=1 and P(k,k)=1 for all k>=0.
Original entry on oeis.org
1, 1, 1, 1, 6, 1, 1, 96, 36, 1, 1, 6306, 3816, 216, 1, 1, 1883076, 1625436, 139536, 1296, 1, 1, 2700393702, 3121837776, 360839016, 5036256, 7776, 1, 1, 19324893252552, 28794284803908, 4200503990976, 78293629296, 181382976, 46656, 1
Offset: 0
Let q=6; 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/(1-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 + 42/3!*x^3 + 7296/4!*x^4 +... (A111829).
Thus the g.f. of column 0 of matrix power P^m is:
1 + m*L(x) + m^2/2!*L(x)*L(6*x) + m^3/3!*L(x)*L(6*x)*L(6^2*x) +
m^4/4!*L(x)*L(6*x)*L(6^2*x)*L(6^3*x) + ...
Triangle P begins:
1;
1,1;
1,6,1;
1,96,36,1;
1,6306,3816,216,1;
1,1883076,1625436,139536,1296,1;
1,2700393702,3121837776,360839016,5036256,7776,1; ...
where P^6 shifts columns left and up one place:
1;
6,1;
96,36,1;
6306,3816,216,1; ...
A111829
Column 0 of the matrix logarithm (A111828) of triangle A111825, which shifts columns left and up under matrix 6th power; these terms are the result of multiplying the element in row n by n!.
Original entry on oeis.org
0, 1, -4, 42, 7296, -7931976, -45557382240, 3064554175021200, 801993619807364206080, -2618439032548254776387771520, -30580166025709706974876961026475520, 4440597519115996836838709580481861376121600
Offset: 0
A(x) = x - 4/2!*x^2 + 42/3!*x^3 + 7296/4!*x^4 - 7931976/5!*x^5 +...
where e.g.f. A(x) satisfies:
x/(1-x) = A(x) + A(x)*A(6*x)/2! + A(x)*A(6*x)*A(6^2*x)/3! +
A(x)*A(6*x)*A(6^2*x)*A(6^3*x)/4! + ...
Let G(x) be the g.f. of A111826 (column 1 of A111825), then
G(x) = 1 + 6*A(x) + 6^2*A(x)*A(6*x)/2! +
6^3*A(x)*A(6*x)*A(6^2*x)/3! +
6^4*A(x)*A(6*x)*A(6^2*x)*A(6^3*x)/4! + ...
-
{a(n,q=6)=local(A=x/(1-x+x*O(x^n)));for(i=1,n, A=x/(1-x)/(1+sum(j=1,n,prod(k=1,j,subst(A,x,q^k*x))/(j+1)!))); return(n!*polcoeff(A,n))}
A111831
Number of partitions of 6*7^n into powers of 7, also equals column 1 of triangle A111830, which shifts columns left and up under matrix 7th power.
Original entry on oeis.org
1, 7, 154, 16275, 9106461, 28543862991, 521136519414483, 56980036448207052005, 38084892600214854893482284, 158081960770204032330986210466109, 4125860571927530263431055188002578191656
Offset: 0
-
a(n,q=7)=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 || j==1,B[i,j]=1,B[i,j]=(A^q)[i-1,j-1]);));A=B); return(A[n+2,2]))
A111836
Number of partitions of 7*8^n into powers of 8, also equals column 1 of triangle A111835, which shifts columns left and up under matrix 8th power.
Original entry on oeis.org
1, 8, 232, 36968, 35593832, 219379963496, 9003699178010216, 2530260913162860295784, 4970141819535151534947497576, 69322146154435681317709098939119208
Offset: 0
-
a(n,q=8)=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 || j==1,B[i,j]=1,B[i,j]=(A^q)[i-1,j-1]);));A=B); return(A[n+2,2]))
A111817
Number of partitions of 3*4^n into powers of 4, also equals column 1 of triangle A078536, which shifts columns left and up under matrix 4th power.
Original entry on oeis.org
1, 4, 28, 524, 29804, 5423660, 3276048300, 6744720496300, 48290009081437356, 1221415413140406958252, 110523986015743458745929900, 36150734459755630877180158951596
Offset: 0
-
a(n,q=4)=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 || j==1,B[i,j]=1,B[i,j]=(A^q)[i-1,j-1]);));A=B); return(A[n+2,2]))
Showing 1-6 of 6 results.
Comments