OFFSET
1,2
MATHEMATICA
a295666[n_]:=Product[GCD[i, n/i]^Boole[PrimeQ[GCD[i, n/i]]], {i, Divisors[n]}]; Select[Range[14500], a295666[#]==# &] (* Stefano Spezia, Feb 20 2024 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Dec 18 2017
STATUS
approved