login
A125802
Column 4 of table A125800; also equals row sums of matrix power A078122^4.
5
1, 5, 35, 485, 15200, 1144664, 215155493, 103674882878, 130648799730635, 437302448840089232, 3936208033244539574405, 96244898501021613327012635, 6446494058446469307795159512465, 1191218783863555524342034469450207222
OFFSET
0,2
COMMENTS
Triangle A078122 shifts left one column under matrix cube and is related to partitions into powers of 3.
PROG
(PARI) a(n)=local(p=4, q=3, A=Mat(1), B); for(m=1, n+1, 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(sum(c=0, n, (A^p)[n+1, c+1]))
CROSSREFS
Cf. A125800, A078122; other columns: A078125, A078124, A125801, A125803.
Sequence in context: A247596 A325408 A097492 * A034217 A011556 A361055
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Dec 10 2006
STATUS
approved