OFFSET
1,1
LINKS
PROG
(PARI)
A156552(n) = {my(f = factor(n), p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res}; \\ From A156552 by David A. Corneth
for(n=1, oo, if(A324824(n), print1(n, ", ")));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Mar 16 2019
STATUS
approved