cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A362148 Numbers that are neither cubefree nor cubefull.

This page as a plain text file.
%I A362148 #61 Dec 06 2024 13:17:55
%S A362148 24,40,48,54,56,72,80,88,96,104,108,112,120,135,136,144,152,160,162,
%T A362148 168,176,184,189,192,200,208,224,232,240,248,250,264,270,272,280,288,
%U A362148 296,297,304,312,320,324,328,336,344,351,352,360,368,375,376,378,384,392,400
%N A362148 Numbers that are neither cubefree nor cubefull.
%C A362148 In fact, every cubefull number > 1 is noncubefree, but the converse is false.
%C A362148 This sequence = A046099 \ A036966 and lists these counterexamples.
%C A362148 Numbers k such that for some primes p and q, k is divisible by p^3*q but not by q^3. - _Robert Israel_, Apr 28 2023
%C A362148 The asymptotic density of this sequence is 1 - 1/zeta(3) = 0.168092... - _Charles R Greathouse IV_, Apr 28 2023
%C A362148 From _Amiram Eldar_, Sep 17 2023: (Start)
%C A362148 Numbers k such that A360539(k) > 1 and A360540(k) > 1.
%C A362148 Equivalently, numbers that have in their prime factorization at least one exponent that is smaller than 3 and at least one exponent that is larger than 2. (End)
%H A362148 Michael De Vlieger, <a href="/A362148/b362148.txt">Table of n, a(n) for n = 1..10000</a>
%F A362148 Equals A362147 \ A004709.
%F A362148 Sum_{n>=1} 1/a(n) = 1 + zeta(s) - zeta(s)/zeta(3*s) - Product_{p prime}(1 + 1/(p^(2*s)*(p^s-1))), s > 1. - _Amiram Eldar_, Sep 17 2023
%e A362148 24 = 2^3 * 3 is noncubefree as it is divisible by the cube 2^3, but it is not cubefull because 3 divides 24 but 3^3 does not divide 24, hence 24 is a term.
%e A362148 648 = 2^4 * 3^3 is noncubefree as it is divisible by the cube 3^3, but it is also cubefull because primes 2 and 3 divide 648, and 2^3 and 3^3 divide also 648, so 648 is not a term.
%p A362148 filter:= proc(n) local F;
%p A362148 F:= ifactors(n)[2][..,2];
%p A362148   min(F) < 3 and max(F) >= 3
%p A362148 end proc:
%p A362148 select(filter, [$1..400]); # _Robert Israel_, Apr 28 2023
%t A362148 Select[Range[500], Min[(e = FactorInteger[#][[;; , 2]])] < 3 && Max[e] > 2 &] (* _Amiram Eldar_, Apr 09 2023 *)
%o A362148 (PARI) isok(k) = (k>1) && (vecmax(factor(k)[, 2])>2) && (vecmin(factor(k)[, 2])<=2); \\ _Michel Marcus_, Apr 19 2023
%Y A362148 Intersection of A046099 (not cubefree) and A362147 (not cubefull)
%Y A362148 Cf. A004709 (cubefree), A036966 (cubefull), A360539, A360540.
%K A362148 nonn,easy
%O A362148 1,1
%A A362148 _Bernard Schott_, Apr 09 2023