OFFSET
1,4
COMMENTS
LINKS
PROG
(PARI) A352896(n) = if(n<=2, n-1, my(m=0); while(n>2, n = A352892(n); m = max(m, bigomega(n))); (m)); \\ Needs also code from A352892.
(PARI) A352896(n) = if(n<=2, n-1, my(m=0); while(n>2, n = A341515(n); m = max(m, bigomega(n))); (m)); \\ Slower, but equivalent.
(PARI)
\\ Faster:
A156552(n) = { my(f = factor(n), p, 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 };
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Antti Karttunen, Apr 08 2022
STATUS
approved