A121431 - OEIS

%I #3 Mar 30 2012 18:36:58

%S 1,1,1,2,3,4,9,15,22,30,69,118,178,250,335,769,1317,1995,2820,3810,

%T 4984,11346,19311,29126,41061,55410,72492,92652,208914,352636,528097,

%U 740035,993678,1294776,1649634,2065146,4613976,7722840,11476963,15971180

%N Number of subpartitions of partition P=[0,0,1,1,1,2,2,2,2,3,3,3,3,3,4,...] (A052146).

%C See A115728 for the definition of subpartitions of a partition.

%F G.f.: 1/(1-x) = Sum_{n>=0} a(n)*x^n*(1-x)^A052146(n).

%e The g.f. may be illustrated by:

%e 1/(1-x) = (1 + 1*x)*(1-x)^0 + (x^2 + 2*x^3 + 3*x^4)*(1-x)^1 +

%e (4*x^5 + 9*x^6 + 15*x^7 + 22*x^8)*(1-x)^2 +

%e (30*x^9 + 69*x^10 + 118*x^11 + 178*x^12 + 250*x^13)*(1-x)^3 +

%e (335*x^14 + 769*x^15 + 1317*x^16 + 1995*x^17 + 2820*x^18 + 3810*x^19)*(1-x)^4 +...

%e When the sequence is put in the form of a triangle:

%e 1, 1,

%e 1, 2, 3,

%e 4, 9, 15, 22,

%e 30, 69, 118, 178, 250,

%e 335, 769, 1317, 1995, 2820, 3810,

%e 4984, 11346, 19311, 29126, 41061, 55410, 72492,

%e 92652, 208914, 352636, 528097, 740035, 993678, 1294776, ...

%e then the columns of this triangle form column 1 (with offset)

%e of successive matrix powers of triangle H=A121412.

%e This sequence is embedded in table A121426 as follows.

%e Column 1 of successive powers of matrix H begin:

%e H^1: [1,1,4,30,335,4984,92652,2065146,53636520,...];

%e H^2: [1,2,9,69,769,11346,208914,4613976,118840164,...];

%e H^3: 1, [3,15,118,1317,19311,352636,7722840,197354133,...];

%e H^4: 1,4, [22,178,1995,29126,528097,11476963,291124693,...];

%e H^5: 1,5,30, [250,2820,41061,740035,15971180,402319275,...];

%e H^6: 1,6,39,335, [3810,55410,993678,21310710,533345745,...];

%e H^7: 1,7,49,434,4984, [72492,1294776,27611970,686872893,...];

%e H^8: 1,8,60,548,6362,92652, [1649634,35003430,865852191,...];

%e H^9: 1,9,72,678,7965,116262,2065146, [43626510,1073540871,...];

%e the terms enclosed in brackets form this sequence.

%o (PARI) {a(n)=local(A); if(n==0,1,A=x+x*O(x^n); for(k=0, n, A+=polcoeff(A, k)*x^k*(1-(1-x)^( (sqrtint(8*k+9)+1)\2 - 1 ) )); polcoeff(A, n))}

%Y Cf. A121412 (triangle H), A121416 (H^2), A121420 (H^3); A121426, A121427; column 1 of H^n: A121414, A121418, A121422; variants: A121430, A121432, A121433.

%K nonn

%O 0,4

%A _Paul D. Hanna_, Jul 30 2006