OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
n a(n)
-----------------------------------
1 100 = 10^2 = 2^2 * 5^2
2 196 = 14^2 = 2^2 * 7^2
3 225 = 15^2 = 3^2 * 5^2
4 324 = 18^2 = 2^2 * 3^4
5 400 = 20^2 = 2^4 * 5^2
6 441 = 21^2 = 3^2 * 7^2
7 484 = 22^2 = 2^2 * 11^2
8 676 = 26^2 = 2^2 * 13^2
9 784 = 28^2 = 2^4 * 7^2
10 1000 = 10^3 = 2^3 * 5^3
11 1089 = 33^2 = 3^2 * 11^2
17 1764 = 42^2 = 2^2 * 3^2 * 7^2
MATHEMATICA
nn = 10000; mm = Sqrt[nn]; i = 1; k = 2; fQ[x_] := And[CompositeQ[x], GCD @@ # == 1, Times @@ MapIndexed[Prime[First[#2]]^#1 &, ReverseSort[#]] != x] &@ FactorInteger[x][[;; , -1]]; 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 09 2025
STATUS
approved
