OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
n a(n)
----------------------------------
1 5184 = 72^2 = 2^6 * 3^4
2 11664 = 108^2 = 2^4 * 3^6
3 82944 = 288^2 = 2^10 * 3^4
4 186624 = 432^2 = 2^8 * 3^6
5 373248 = 72^3 = 2^9 * 3^6
6 419904 = 648^2 = 2^6 * 3^8
7 746496 = 864^2 = 2^10 * 3^6
8 944784 = 972^2 = 2^4 * 3^10
9 1259712 = 108^3 = 2^6 * 3^9
10 1327104 = 1152^2 = 2^14 * 3^4
11 3240000 = 1800^2 = 2^6 * 3^4 * 5^4
12 3779136 = 1944^2 = 2^6 * 3^10
MATHEMATICA
nn = 2^30; mm = Sqrt[nn]; i = 1; k = 2;
MapIndexed[Set[S[First[#2]], #1] &,
Select[
Union@ Flatten@ Table[a^2*b^3, {b, Surd[mm, 3]}, {a, Sqrt[mm/b^3]}],
And[#[[1, 1]] == 2, Length[#] > 1,
Union@ Differences@ Map[PrimePi, #[[;; , 1]] ] == {1},
Apply[GCD, #[[;; , -1]] ] == 1] &[FactorInteger[#]] &] ];
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, Sep 30 2025
STATUS
approved
