OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
n a(n)
-------------------------------------
1 144 = 12^2 = 2^4 * 3^2
2 576 = 24^2 = 2^6 * 3^2
3 1728 = 12^3 = 2^6 * 3^3
4 2304 = 48^2 = 2^8 * 3^2
5 3600 = 60^2 = 2^4 * 3^2 * 5^2
6 9216 = 96^2 = 2^10 * 3^2
7 13824 = 24^3 = 2^9 * 3^3
8 14400 = 120^2 = 2^6 * 3^2 * 5^2
9 20736 = 12^4 = 2^8 * 3^4
10 32400 = 180^2 = 2^4 * 3^4 * 5^2
11 36864 = 192^2 = 2^12 * 3^2
12 57600 = 240^2 = 2^8 * 3^2 * 5^2
MATHEMATICA
nn = 2^31; mm = Sqrt[nn]; i = 1; k = 2; fQ[x_] := And[Length[#] > 1, 1 == Min[#] < Max[#], Times @@ MapIndexed[Prime[First[#2]]^#1 &, ReverseSort[#]] == x] &[FactorInteger[x][[;; , -1]] ]; MapIndexed[Set[S[First[#2]], #1] &, Select[Range@ Sqrt[nn], fQ]]; Union@ Reap[While[j = 2; While[S[i]^j < nn, Sow[S[i]^j]; j++]; j > 2, k++; i++] ][[-1, 1]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael De Vlieger, Oct 10 2025
STATUS
approved
