OFFSET
1,1
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
Table of n, a(n), and q(a(n)) for n = 1..12:
n a(n) q(a(n))
--------------------------------
1 200 = 2^3 * 5^2 3
2 392 = 2^3 * 7^2 3
3 500 = 2^2 * 5^3 3
4 675 = 3^3 * 5^2 2
5 800 = 2^5 * 5^2 3
6 968 = 2^3 * 11^2 3
7 1125 = 3^2 * 5^3 2
8 1323 = 3^3 * 7^2 2
9 1352 = 2^3 * 13^2 3
10 1372 = 2^2 * 7^3 3
11 1568 = 2^5 * 7^2 3
12 2000 = 2^4 * 5^3 3
The number 12 is not a term since it is not powerful (i.e., not in A001694).
The number 36, though powerful, is not a term since it is a perfect square.
s(1) = 72 is not in the sequence since q(72) > gpf(72), i.e., 5 > 3.
s(2) = 108 is not in the sequence since q(108) > gpf(108), i.e., 5 > 3.
a(1) = s(3) = 200 because q(200) < gpf(200), i.e., 3 < 5.
a(2) = s(4) = 392 because q(392) < gpf(392), i.e., 3 < 7, etc.
MATHEMATICA
nn = 6; mm = Times @@ Prime@ Range[nn]; Complement[Select[Union@ Flatten@ Table[a^2*b^3, {b, Surd[mm, 3]}, {a, Sqrt[mm/b^3]}], And[Length[#2] > 1, GCD @@ #2 == 1] & @@ {#, FactorInteger[#][[;; , -1]]} &], Union@ Flatten[f[nn][[3 ;; -1, 2 ;; -1]] ] ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael De Vlieger, Jul 21 2025
STATUS
approved
