OFFSET
1,4
COMMENTS
PROG
(PARI)
A175992(n, k)=if(k==n, 0, if(n%k==0, 1, 0));
N=40; /* that many terms */
M=matrix(N, N, r, c, A175992(r, c) );
E=sum(j=0, N, M^j/j!); /* matrix exponential */
c1=E[, 1]; /* first column */
vector(N, j, numerator(c1[j])) /* show terms */ /* Joerg Arndt, Apr 16 2011 */
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Mats Granvik, Apr 04 2011
EXTENSIONS
Sequence terms added, title changed, and comment lines changed by Mats Granvik, Apr 14 2011
STATUS
approved
