OFFSET
1,2
COMMENTS
Numbers whose largest cubefree divisor (A007948) and cubefull part (A360540) have the same number of divisors (A000005).
If k and m are coprime terms, then k*m is also a term.
The characteristic function of this sequence depends only on prime signature.
1 is the only cubefree (A004709) term.
Includes the 4th powers of squarefree numbers (1 and A113849).
Numbers of the for m*p^(3*2^k+1), where m is squarefree, p is prime, gcd(m, p) = 1 and omega(m) = k, are all terms. In particular, this sequence includes numbers of the form p^7*q, where p != q are primes (A179664), and numbers of the form p^13*q*r where p, q, and r are distinct primes.
The corresponding numbers of cubefree (or 3-full) divisors are 1, 3, 3, 6, 3, 6, 6, 9, 6, 6, 6, 3, 6, 6, ... .
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
q[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, Times @@ (Min[#, 3] & /@ (e + 1)) == Times @@ (Max[#, 1] & /@ (e - 1))]; q[1] = True; Select[Range[10^4], q]
PROG
(PARI) is(k) = {my(e = factor(k)[, 2]); prod(i = 1, #e, min(e[i] + 1, 3)) == prod(i = 1, #e, max(e[i] - 1, 1)); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Feb 25 2023
STATUS
approved