OFFSET
1,1
EXAMPLE
63 is here because for 126 = 2*63, lcm(dd(126)) = lcm(1, 1, 3, 1, 2, 5, 4, 3, 21, 21, 63) = 1260, so gcd(126, lcm(dd(126))) = gcd(126, 1260) = 126.
PROG
(PARI) f(n) = {my(d = divisors(n), dd = vector(#d-1, k, d[k+1] - d[k])); gcd(n, lcm(dd)); }
isok(n) = (f(2*n) == 2*n); \\ Michel Marcus, Mar 29 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Apr 25 2001
STATUS
approved