OFFSET
1,2
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10001
PROG
(PARI) \\ Based on Michael Somos's code for A007896
n=101;
v = vector(n, k, k==1);
for(k=2, n, m = #digits(n, k) - 1; A = (1 - x)^ -(sigma(k, 0)) + x * O(x^m); w = vector(n); for(i=0, m, w[k^i] = polcoeff(A, i)); v = dirmul(v, w));
v
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, May 26 2014
STATUS
approved