A242648 - OEIS

%I #16 Nov 24 2023 12:04:23

%S 1,3,4,13,6,24,8,46,23,36,12,116,14,48,48,161,18,156,20,174,64,72,24,

%T 484,52,84,112,232,30,360,32,526,96,108,96,841,38,120,112,726,42,480,

%U 44,348,312,144,48,1864,93,357,144,406,54,888,144,968,160,180

%N Dirichlet g.f.: Product_{n>=2} 1/(1-1/n^s)^sigma(n).

%o (PARI) \\ Based on _Michael Somos_'s code for A007896

%o n=101;

%o v = vector(n, k, k==1);

%o    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));

%o v

%Y Cf. A000203, A050367, A007896, A242649.

%K nonn

%O 1,2

%A _N. J. A. Sloane_, May 26 2014