login
A296204
Numbers k such that Product_{d|k, gcd(d,k/d) is prime} gcd(d,k/d) = k; the fixed points of A295666.
1
1, 36, 100, 144, 196, 225, 324, 400, 441, 484, 676, 784, 1089, 1156, 1225, 1296, 1444, 1521, 1936, 2025, 2116, 2500, 2601, 2704, 3025, 3249, 3364, 3844, 3969, 4225, 4624, 4761, 5476, 5625, 5776, 5929, 6724, 7225, 7396, 7569, 8281, 8464, 8649, 8836, 9025, 9604, 9801, 10000, 11236, 12321, 13225, 13456, 13689, 13924, 14161
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
Cf. A295666, A296205 (the square roots).
Sequence in context: A072413 A131605 A376119 * A063734 A069057 A342402
KEYWORD
nonn
AUTHOR
Antti Karttunen, Dec 18 2017
STATUS
approved