OFFSET
1,1
COMMENTS
Powers k^m, m > 1, where k is in A386224.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
n a(n)
-------------------------------------
1 324 = 18^2 = 2^2 * 3^4
2 2916 = 54^2 = 2^2 * 3^6
3 5832 = 18^3 = 2^3 * 3^6
4 8100 = 90^2 = 2^2 * 3^4 * 5^2
5 22500 = 150^2 = 2^2 * 3^2 * 5^4
6 26244 = 162^2 = 2^2 * 3^8
7 72900 = 270^2 = 2^2 * 3^6 * 5^2
8 90000 = 300^2 = 2^4 * 3^2 * 5^4
9 104976 = 18^4 = 2^4 * 3^8
10 157464 = 54^3 = 2^3 * 3^9
11 202500 = 450^2 = 2^2 * 3^4 * 5^4
12 236196 = 486^2 = 2^2 * 3^10
MATHEMATICA
nn = 6000000; mm = Sqrt[nn]; i = 1; k = 2; fQ[x_] := And[#[[1, 1]] == 2, Length[#] > 1, Union@ Differences@ Map[PrimePi, #[[;; , 1]] ] == {1}, ReverseSort[#[[;; , -1]] ] != #[[;; , -1]], 1 == Min[#] < Max[#] &[#[[;; , -1]] ] ] &[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, Oct 17 2025
STATUS
approved
