OFFSET
1,2
COMMENTS
PROG
(PARI) a188999(n) = my(f = factor(n)); for (i=1, #f~, p = f[i, 1]; e = f[i, 2]; f[i, 1] = if (e % 2, (p^(e+1)-1)/(p-1), (p^(e+1)-1)/(p-1) -p^(e/2)); f[i, 2] = 1; ); factorback(f) \\ after Michel Marcus in A189000
is_a189000(n) = ! frac(a188999(n)/n) \\ after Michel Marcus in A189000
for(n=1, oo, if(is_a189000(n), print1(a188999(n)/n, ", "))) \\ Felix Fröhlich, Sep 03 2018
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Michel Marcus, Sep 03 2018, following a suggestion from Felix Fröhlich
EXTENSIONS
a(33)-a(42) from Giovanni Resta, Sep 03 2018
STATUS
approved