OFFSET
1,1
COMMENTS
Numbers in this sequence are neither cubefree nor cubeful (in A362148), hence are in A389558, a proper subset of A286708.
A389558(48) = 5400 = 2^3 * 3^3 * 5, divisible by 2 prime cubes; a(48) = 5408.
Numbers whose prime factorization exponents include at least one 2, and only one exponent that is larger than 2. - Amiram Eldar, Dec 21 2025
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
FORMULA
Sum_{n>=1} 1/a(n) = (15/Pi^2) * Sum_{p prime} 1/((p-1)*(p^2+1)) - Sum_{p prime} 1/((p-1)*p^2) = 0.081675144208728966934... . - Amiram Eldar, Dec 21 2025
EXAMPLE
Table of n, a(n) for select n:
n a(n)
-----------------------------------
1 72 = 2^3 * 3^2
2 108 = 2^2 * 3^3
3 144 = 2^4 * 3^2
4 200 = 2^3 * 5^2
5 288 = 2^5 * 3^2
6 324 = 2^2 * 3^4
11 675 = 3^3 * 5^2
23 1800 = 2^3 * 3^2 * 5^2
29 2700 = 2^2 * 3^3 * 5^2
37 3528 = 2^3 * 3^2 * 7^2
38 3600 = 2^4 * 3^2 * 5^2
241 88200 = 2^3 * 3^2 * 5^2 * 7^2
MATHEMATICA
With[{nn = 5400}, Select[Rest@ Union@ Flatten@ Table[a^2*b^3, {b, Surd[nn, 3]}, {a, Sqrt[nn/b^3]}], And[Length[#] > 1, Count[#, _?(# > 2 &)] == 1] &@ FactorInteger[#][[;; , -1]] &] ]
PROG
(PARI) isok(k) = {my(e = vecsort(factor(k)[, 2])); #e > 1 && e[1] > 1 && e[#e-1] == 2 && e[#e] > 2; } \\ Amiram Eldar, Dec 21 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael De Vlieger, Dec 10 2025
STATUS
approved
