0,3

P = A136220 is a triangular matrix where column k of P^3 equals column 0 of P^(3k+3) such that column 0 of P^3 equals column 0 of P shift one place left. Tables A136217 and A136218 are defined by recurrences seemingly unrelated to the matrix product recurrence of A136220 and yet they all generate this same sequence in column 0.

Paul D. Hanna, Table of n, a(n) for n = 0..90

Equals column 0 of triangle P=A136220, which begins:

1;

1, 1;

3, 2, 1;

15, 10, 3, 1;

108, 75, 21, 4, 1;

1036, 753, 208, 36, 5, 1;

12569, 9534, 2637, 442, 55, 6, 1;

185704, 146353, 40731, 6742, 805, 78, 7, 1; ...

where column k of P^3 = column 0 of P^(3k+3) such that

column 0 of P^3 = column 0 of P shift one place left.

Surprisingly, column 0 of P is also found in square A136218:

(1),(1),1,(1),1,(1),1,(1),1,1,(1),1,1,(1),1,1,(1),1,1,1,(1),...;

(1),(2),3,(4),5,(6),7,(8),9,10,(11),12,13,(14),15,16,(17),...;

(3),(8),15,(24),34,(46),59,(74),90,108,(127),147,169,(192),...;

(15),(49),108,(198),306,(453),622,(838),1080,1377,(1704),...;

(108),(414),1036,(2116),3493,(5555),8040,(11477),15483,...;

(1036),(4529),12569,(28052),48800,(82328),124335,(186261),...;

(12569),(61369),185704,(446560),811111,(1438447),2250731,...;

...

and has a recurrence similar to that of square array A136212

which generates the triple factorials.

(PARI) /* Generate using matrix product recurrences of triangle A136220: */ {a(n)=local(P=Mat(1), U, PShR); if(n>0, for(i=0, n, PShR=matrix(#P, #P, r, c, if(r>=c, if(r==c, 1, if(c==1, 0, P[r-1, c-1])))); U=P*PShR^2; U=matrix(#P+1, #P+1, r, c, if(r>=c, if(r<#P+1, U[r, c], if(c==1, (P^3)[ #P, 1], (P^(3*c-1))[r-c+1, 1])))); P=matrix(#U, #U, r, c, if(r>=c, if(r<#R, P[r, c], (U^c)[r-c+1, 1]))))); P[n+1, 1]}

(PARI) /* Generated as column 0 in triangle A136218 (faster): */ {a(n)=local(A=[1], B); if(n>0, for(i=1, n, m=1; B=[0]; for(j=1, #A, if(j+m-1==(m*(m+7))\6, m+=1; B=concat(B, 0)); B=concat(B, A[j])); A=Vec(Polrev(Vec(Pol(B)/(1-x+O(x^#B))))))); A[1]}

Cf. A136220 (P), A136228 (U), A136231 (W=P^3).

Cf. other columns of P: A136222, A136223, A136224.

Cf. related tables: A136217, A136218.

Cf. variants: A091352, A135881.

Adjacent sequences: A136218 A136219 A136220 this_sequence A136222 A136223 A136224

Sequence in context: A105618 A120732 A090351 this_sequence A110328 A054201 A090355

nonn

Paul D. Hanna (pauldhanna(AT)juno.com), Dec 25 2007