OFFSET
1,2
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)) + 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