OFFSET
1,2
COMMENTS
Terms p^p, p prime, form the subsequence A051674. - Bernard Schott, Sep 21 2022
Terms p^q * q^p with distinct primes p and q form the subsequence A082949. - Bernard Schott, Feb 01 2023
EXAMPLE
576 = 2^6 * 3^2, lcm(2,3) = 6 = lcm(6,2), hence 576 is a term.
MATHEMATICA
Select[Range[10^6], Equal @@ LCM @@ FactorInteger[#] &] (* Amiram Eldar, Aug 07 2022 *)
PROG
(PARI) isok(k) = my(f=factor(k)); lcm(f[, 1]) == lcm(f[, 2]); \\ Michel Marcus, Aug 07 2022
(Python)
from math import lcm
from sympy import factorint
def ok(n): f = factorint(n); return lcm(*f.keys()) == lcm(*f.values())
print([k for k in range(10**6) if ok(k)]) # Michael S. Branicky, Aug 07 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Jean-Marc Rebert, Aug 07 2022
STATUS
approved