OFFSET
0,12
COMMENTS
Column k=1 appears to be A032190(n), n=2,3,...
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..1274
EXAMPLE
Triangle begins:
1;
1, 1;
1, 1, 1;
1, 1, 1, 1;
1, 2, 1, 1, 1;
1, 2, 2, 1, 1, 1;
1, 4, 4, 2, 1, 1, 1;
1, 4, 6, 4, 2, 1, 1, 1;
1, 7, 11, 8, 4, 2, 1, 1, 1;
1, 9, 19, 14, 8, 4, 2, 1, 1, 1;
1, 14, 33, 27, 16, 8, 4, 2, 1, 1, 1;
...
PROG
(PARI) \\ here R(n) is A048887 transposed
R(n)={Mat(vector(n, k, Col((1-x)/(1-2*x+x^(k+1)) - 1 + O(x*x^n))))}
S(M)={matrix(#M-1, #M-1, n, k, if(k<n-1, (k+1)*M[n-k-1, k+1] + sum(j=1, k, j*(M[n-j, k+1]-M[n-j, k])), if(k<n, n)))}
T(n)={my(M=S(R(n+1))); matid(n) + matrix(n, n, n, k, sumdiv(n, d, if(k<d, eulerphi(n/d)*M[d, k]))/n) }
my(M=T(10)); for(n=0, #M, for(k=0, n, print1(if(k==0, 1, M[n, k]), ", ")); print) \\ Andrew Howroyd, Oct 15 2017
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Christopher Lenard (c.lenard(AT)bendigo.latrobe.edu.au), Aug 22 2001
EXTENSIONS
a(0,0)=1 from Andrew Howroyd, Oct 15 2017
STATUS
approved