OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
Table of n, a(n) for n = 1..12:
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 5184 = 72^2 = 2^6 * 3^4
7 9216 = 96^2 = 2^10 * 3^2
8 13824 = 24^3 = 2^9 * 3^3
9 14400 = 120^2 = 2^6 * 3^2 * 5^2
10 20736 = 12^4 = 2^8 * 3^4
11 32400 = 180^2 = 2^4 * 3^4 * 5^2
12 36864 = 192^2 = 2^12 * 3^2
MATHEMATICA
nn = 2^21; mm = Sqrt[nn]; i = 1; k = 2; fQ[x_] := And[Length[#] > 1, CountDistinct[#[[;; , -1]] ] > 1, Times @@ MapIndexed[Prime[First[#2] ]^#1 &, ReverseSort[#[[;; , -1]] ] ] == x] &[FactorInteger[x] ]; MapIndexed[Set[S[First[#2]], #1] &, Select[Range[mm], 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, Nov 02 2025
STATUS
approved
