A300547 a(n) = [x^n] Product_{d|n} (1 - x^d).
1, -1, -1, -1, -1, -1, -2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 0, -1, -1, -1, 0, -1, -1, -1, -2, -1, -2, -1, -1, -1, -1, -1, -2, -1, -1, -1, 0, -1, -1, -1, -1, -1, -1, -1, -3, -1, -1, -1, -1, -1, -2, -1, -1, -1, -1, -1, -5, -1, -1, -1, -1, -1, -3, -1, -1, -1, -1, -1, 0, -1, -1, -1, -1, -1, -2, -1, 1, -1, -1, -1, -2, -1
Offset: 0
Links
Programs
-
Mathematica
Table[SeriesCoefficient[Product[(1 - Boole[Mod[n, k] == 0] x^k), {k, 1, n}], {x, 0, n}], {n, 0, 85}]
-
PARI
A300547(n) = { if(!n,return(1)); my(p=1); fordiv(n,d, p *= (1 - 'x^d)); polcoeff(p,n); }; \\ Antti Karttunen, Sep 25 2018