OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
Table of n, a(n) for select n:
n a(n)
-----------------------------------
1 400 = 20^2 = 2^4 * 5^2
2 784 = 28^2 = 2^4 * 7^2
3 1600 = 40^2 = 2^6 * 5^2
4 1936 = 44^2 = 2^4 * 11^2
5 2025 = 45^2 = 3^4 * 5^2
6 2500 = 50^2 = 2^2 * 5^4
7 2704 = 52^2 = 2^4 * 13^2
8 3136 = 56^2 = 2^6 * 7^2
9 3969 = 63^2 = 3^4 * 7^2
11 5625 = 75^2 = 3^2 * 5^4
14 7056 = 84^2 = 2^4 * 3^2 * 7^2
47 40000 = 200^2 = 2^6 * 5^4
MATHEMATICA
a053669[x_] := Block[{qx}, qx = 2; While[Divisible[x, qx], qx = NextPrime[qx]]; qx]; nn = 2^15; mm = Sqrt[nn]; i = 1; k = 2; fQ[x_] := And[Length[#] > 1, CountDistinct[#[[;; , -1]] ] > 1, a053669[x] < #[[-1, 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, Nov 03 2025
STATUS
approved
