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 400 = 20^2 2^4 * 5^2 3
2 784 = 28^2 2^4 * 7^2 3
3 1600 = 40^2 2^6 * 5^2 3
4 1936 = 44^2 2^4 * 11^2 3
5 2025 = 45^2 3^4 * 5^2 2
6 2500 = 50^2 2^2 * 5^4 3
7 2704 = 52^2 2^4 * 13^2 3
8 3136 = 56^2 2^6 * 7^2 3
9 3969 = 63^2 3^4 * 7^2 2
14 7056 = 84^2 = 2^4 * 3^2 * 7^2 5
16 8000 = 20^3 = 2^6 * 5^3 3
26 17424 = 132^2 = 2^4 * 3^2 * 11^2 5
MATHEMATICA
nn = 2^16; mm = Sqrt[nn]; i = 1; k = 2; a053669[x_] := Module[{q = 2}, While[Divisible[x, q], q = NextPrime[q]]; q]; fQ[x_] := And[Length[#] > 1, a053669[x] < #[[-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 13 2025
STATUS
approved
