OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
Let q = A053669(a(n)).
Table of n, a(n) for select n:
n a(n) q
----------------------------------------
1 40000 = 200^2 = 2^6 * 5^4 3
2 153664 = 392^2 = 2^6 * 7^4 3
3 250000 = 500^2 = 2^4 * 5^6 3
4 455625 = 675^2 = 3^6 * 5^4 2
5 640000 = 800^2 = 2^10 * 5^4 3
6 937024 = 968^2 = 2^6 * 11^4 3
7 1265625 = 1125^2 = 3^4 * 5^6 2
8 1750329 = 1323^2 = 3^6 * 7^4 2
9 1827904 = 1352^2 = 2^6 * 13^4 3
10 1882384 = 1372^2 = 2^4 * 7^6 3
11 2458624 = 1568^2 = 2^10 * 7^4 3
14 8000000 = 200^3 = 2^9 * 5^6 3
MATHEMATICA
nn = 2^26; mm = Sqrt[nn]; i = 1; k = 2;
a053669[x_] := Module[{qx}, qx = 2; While[Divisible[x, qx], qx = NextPrime[qx]]; qx];
MapIndexed[Set[S[First[#2]], #1] &,
Select[Union@ Flatten@ Table[a^2*b^3, {b, Surd[mm, 3]}, {a, Sqrt[mm/b^3]}],
And[GCD @@ #2[[;; , -1]] == 1,
a053669[#1] < #2[[-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, Oct 02 2025
STATUS
approved
