login
A242648
Dirichlet g.f.: Product_{n>=2} 1/(1-1/n^s)^sigma(n).
1
1, 3, 4, 13, 6, 24, 8, 46, 23, 36, 12, 116, 14, 48, 48, 161, 18, 156, 20, 174, 64, 72, 24, 484, 52, 84, 112, 232, 30, 360, 32, 526, 96, 108, 96, 841, 38, 120, 112, 726, 42, 480, 44, 348, 312, 144, 48, 1864, 93, 357, 144, 406, 54, 888, 144, 968, 160, 180
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