A125802 Column 4 of table A125800; also equals row sums of matrix power A078122^4.

1, 5, 35, 485, 15200, 1144664, 215155493, 103674882878, 130648799730635, 437302448840089232, 3936208033244539574405, 96244898501021613327012635, 6446494058446469307795159512465, 1191218783863555524342034469450207222

Offset: 0

Author

Paul D. Hanna, Dec 10 2006

Keywords

• nonn

Comments

Triangle A078122 shifts left one column under matrix cube and is related to partitions into powers of 3.

Crossrefs

Cf. A125800, A078122; other columns: A078125, A078124, A125801, A125803.

Programs

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]))