OFFSET
1,4
COMMENTS
a(p) = 1, a(p*q) = 6, a(p^2*q) = 19, a(p^4)= 10 etc. where p and q are primes. Question: To find an expression for a(N) where N = p^a*q^b*r^c...p,q,r are primes.
The positions of records are: 1, 2, 4, 6, 12, 24, 30, 36, 48, 60, 120, 180, 210, 240, ... - Antti Karttunen, May 19 2017
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..359
FORMULA
a(n) = A207329(n) - 1. - Michel Marcus, Sep 17 2013
EXAMPLE
a(6) = 6 as the divisors of 6 are 1, 2, 3 and 6 and the distinct products of these divisors are 2, 3, 6, 12, 18 and 36.
PROG
(PARI) a(n) = {d = divisors(n); s = Set(); for (i = 1, 2^#d, b = binary(i); if (sum(j = 1, #b, b[j]) > 1, s = Set(concat(s, prod(j = 1, #b, if (b[j] == 1, d[j], 1)))); ); ); #s; } \\ Michel Marcus, Sep 17 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Mar 27 2002
EXTENSIONS
Corrected and extended by David Wasserman, Apr 23 2003
STATUS
approved