[go: up one dir, main page]

login
A291208
Number of noncube divisors of n.
2
0, 1, 1, 2, 1, 3, 1, 2, 2, 3, 1, 5, 1, 3, 3, 3, 1, 5, 1, 5, 3, 3, 1, 6, 2, 3, 2, 5, 1, 7, 1, 4, 3, 3, 3, 8, 1, 3, 3, 6, 1, 7, 1, 5, 5, 3, 1, 8, 2, 5, 3, 5, 1, 6, 3, 6, 3, 3, 1, 11, 1, 3, 5, 4, 3, 7, 1, 5, 3, 7, 1, 10, 1, 3, 5, 5, 3, 7, 1, 8, 3, 3, 1, 11, 3, 3, 3, 6, 1, 11, 3, 5, 3, 3, 3, 10, 1, 5, 5, 8, 1, 7, 1, 6, 7
OFFSET
1,4
FORMULA
G.f.: Sum_{k>=1} x^A007412(k)/(1 - x^A007412(k)).
G.f.: Sum_{k>=1} (x^k - x^(k^3))/((1 - x^k)*(1 - x^(k^3))).
a(n) = A000005(n) - A061704(n).
EXAMPLE
a(8) = 2 because 8 has 4 divisors {1, 2, 4, 8} among which 2 are noncubes {2, 4}.
MATHEMATICA
nmax = 105; Rest[CoefficientList[Series[Sum[(x^k - x^k^3)/((1 - x^k) (1 - x^k^3)), {k, 1, nmax}], {x, 0, nmax}], x]]
PROG
(PARI) a(n) = sumdiv(n, d, !ispower(d, 3)); \\ Michel Marcus, Aug 21 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Aug 21 2017
STATUS
approved