A216427 Numbers of the form a^2*b^3, where a >= 2 and b >= 2.

32, 72, 108, 128, 200, 243, 256, 288, 392, 432, 500, 512, 576, 648, 675, 800, 864, 968, 972, 1024, 1125, 1152, 1323, 1352, 1372, 1568, 1600, 1728, 1800, 1944, 2000, 2048, 2187, 2304, 2312, 2592, 2700, 2888, 2916, 3087, 3125, 3136, 3200, 3267, 3456, 3528, 3872, 3888, 4000, 4096, 4232, 4500, 4563, 4608

Offset: 1

Author

V. Raman, Sep 07 2012

Keywords

• nonn

Comments

Powerful numbers (A001694) that are not squares of cubefree numbers (A004709), cubes of squarefree numbers (A062838), or 6th powers of primes (A030516). - Amiram Eldar, Feb 07 2023

Links

• Charles R Greathouse IV, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A143610, A216426.

Cf. A001694, A002117, A004709, A013661, A013664, A062838, A085966.

Programs

Mathematica

With[{max = 5000}, Union[Table[i^2*j^3, {j, 2, max^(1/3)}, {i, 2, Sqrt[max/j^3]}] // Flatten]] (* Amiram Eldar, Feb 07 2023 *)

PARI

list(lim)=my(v=List()); for(b=2, sqrtnint(lim\4, 3), for(a=2, sqrtint(lim\b^3), listput(v, a^2*b^3))); Set(v) \\ Charles R Greathouse IV, Jan 03 2014

Python

from math import isqrt
from sympy import mobius, integer_nthroot, primepi
def A216427(n):
    def squarefreepi(n): return int(sum(mobius(k)*(n//k**2) for k in range(1, isqrt(n)+1)))
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x):
        j, b = isqrt(x), integer_nthroot(x,6)[0]
        l, c = 0, n+x-1+primepi(b)+sum(mobius(k)*(j//k**3) for k in range(1, b+1))
        while j>1:
            k2 = integer_nthroot(x//j**2,3)[0]+1
            w = squarefreepi(k2-1)
            c -= j*(w-l)
            l, j = w, isqrt(x//k2**3)
        return c+l
    return bisection(f,n,n) # Chai Wah Wu, Sep 13 2024

Formula

Sum_{n>=1} 1/a(n) = 1 + ((zeta(2)-1)*(zeta(3)-1)-1)/zeta(6) - P(6) = 0.12806919584708298724..., where P(s) is the prime zeta function. - Amiram Eldar, Feb 07 2023