OFFSET
1,1
COMMENTS
Powers k^m, with k in A378002 and m > 1.
Numbers j such that omega(j) > 1, whose prime power factor exponents have a greatest common factor g > 1, but m^(1/g) has prime power factor exponents exceed 1 and are coprime.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
EXAMPLE
Table of n, a(n) for n <= 12:
n a(n)
-----------------------------------------
1 5184 = 72^2 = 2^6 * 3^4
2 82944 = 288^2 = 2^10 * 3^4
3 186624 = 432^2 = 2^8 * 3^6
4 373248 = 72^3 = 2^9 * 3^6
5 746496 = 864^2 = 2^10 * 3^6
6 1327104 = 1152^2 = 2^14 * 3^4
7 3240000 = 1800^2 = 2^6 * 3^4 * 5^4
8 6718464 = 2592^2 = 2^10 * 3^8
9 11943936 = 3456^2 = 2^14 * 3^6
10 21233664 = 4608^2 = 2^18 * 3^4
11 23887872 = 288^3 = 2^15 * 3^6
12 26873856 = 72^4 = 2^12 * 3^8
MATHEMATICA
nn = 2^31; mm = Sqrt[nn]; i = 1; k = 2; fQ[x_] := And[Length[#] > 1, GCD @@ # == 1, Times @@ MapIndexed[Prime[First[#2]]^#1 &, ReverseSort[#]] == x] &[FactorInteger[x][[;; , -1]] ]; MapIndexed[Set[S[First[#2]], #1] &, Select[Union@ Flatten@ Table[a^2*b^3, {b, Surd[mm, 3]}, {a, Sqrt[mm/b^3]}], 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 01 2025
STATUS
approved
