OFFSET
1,6
FORMULA
a(1) = 1, a(p) = 1, a(p*q) = p*q, a(p*q*...*z) = (p*q*...*z)^(2^(k-1)-1), for p, q = primes, k = natural numbers, p*q*...*z = product of k (k > 2) distinct primes p, q, ..., z. [corrected by Jason Yuen, Dec 17 2025]
EXAMPLE
For n = 16, set of such divisors is {1, 8}; a(16) = 1*8 = 8.
PROG
(SageMath) A178646 = lambda n: prod(d for d in divisors(n) if not n.is_power_of(d)) # D. S. McNeil, Dec 28 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Dec 25 2010
EXTENSIONS
More terms from Michel Marcus, Dec 17 2025
STATUS
approved
