OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..18030
EXAMPLE
n a(n)
------------------------------------
1 144 = 12^2 = 2^4 * 3^2
2 324 = 18^2 = 2^2 * 3^4
3 576 = 24^2 = 2^6 * 3^2
4 1728 = 12^3 = 2^6 * 3^3
5 2304 = 48^2 = 2^8 * 3^2
6 2916 = 54^2 = 2^2 * 3^6
7 3600 = 60^2 = 2^4 * 3^2 * 5^2
8 5832 = 18^3 = 2^3 * 3^6
9 8100 = 90^2 = 2^2 * 3^4 * 5^2
10 9216 = 96^2 = 2^10 * 3^2
11 13824 = 24^3 = 2^9 * 3^3
12 14400 = 120^2 = 2^6 * 3^2 * 5^2
MATHEMATICA
nn = 2^20; mm = Sqrt[nn]; i = 1; k = 2; fQ[x_] := And[#[[1, 1]] == 2, Length[#] > 1, Union@ Differences@ Map[PrimePi, #[[;; , 1]] ] == {1}, 1 == Min[#] < Max[#] &[#[[;; , -1]] ] ] &[FactorInteger[x] ]; 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
