A327731 - OEIS

A327731

Expansion of Product_{i>=1, j>=1} (1 + x^(i*(2*j - 1))).

1, 1, 1, 3, 3, 5, 8, 10, 14, 21, 28, 36, 51, 65, 86, 117, 148, 190, 251, 316, 402, 519, 647, 814, 1032, 1282, 1593, 1994, 2457, 3029, 3754, 4591, 5617, 6895, 8381, 10193, 12411, 14999, 18125, 21919, 26359, 31672, 38074, 45556, 54468, 65134, 77576, 92322

(list; graph; refs; listen; history; text; internal format)

OFFSET

0,4

COMMENTS

Weigh transform of A001227.

LINKS

Table of n, a(n) for n=0..47.

FORMULA

G.f.: Product_{k>=1} (1 + x^k)^A001227(k).

MATHEMATICA

nmax = 47; CoefficientList[Series[Product[(1 + x^k)^DivisorSum[k, Mod[#, 2] &], {k, 1, nmax}], {x, 0, nmax}], x]

a[n_] := a[n] = If[n == 0, 1, Sum[Sum[(-1)^(k/d + 1) d DivisorSum[d, Mod[#, 2] &], {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 47}]