OFFSET
1,3
LINKS
FORMULA
Multiplicative with a(p^e) = p^e if (p^(1+e) - 1)/(p-1) = 1 + p + p^2 + ... + p^e is composite, and 1 otherwise.
a(n) = n / A324892(n).
EXAMPLE
For n = 150 = 2 * 3 * 5^2, sigma(2) = 3 is prime, sigma(3) = 4 is not prime, and sigma(25) = 31 is prime, thus a(150) = 3.
PROG
(PARI) A324894(n) = { my(f=factor(n)); prod(i=1, #f~, (f[i, 1]^!isprime(sigma(f[i, 1]^f[i, 2])))^f[i, 2]); };
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Antti Karttunen, Mar 29 2019
STATUS
approved