A007304 -id:A007304 - cp's OEIS frontend

A014311 Numbers with exactly 3 ones in binary expansion.

7, 11, 13, 14, 19, 21, 22, 25, 26, 28, 35, 37, 38, 41, 42, 44, 49, 50, 52, 56, 67, 69, 70, 73, 74, 76, 81, 82, 84, 88, 97, 98, 100, 104, 112, 131, 133, 134, 137, 138, 140, 145, 146, 148, 152, 161, 162, 164, 168, 176, 193, 194, 196, 200, 208, 224, 259, 261, 262, 265, 266, 268, 273, 274, 276, 280, 289, 290, 292, 296, 304

Offset: 1

Al Black (gblack(AT)nol.net)

Equivalently, sums of three distinct powers of 2.
Appears to give all n such that 64 is the highest power of 2 dividing A005148(n). - Benoit Cloitre, Jun 22 2002
From Gus Wiseman, Oct 05 2020: (Start)
These are numbers k such that the k-th composition in standard order has length 3. The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions. The sequence together with the corresponding standard compositions begins:
(1,1,1)     44: (2,1,3)     97: (1,5,1)
(2,1,1)     49: (1,4,1)     98: (1,4,2)
(1,2,1)     50: (1,3,2)    100: (1,3,3)
(1,1,2)     52: (1,2,3)    104: (1,2,4)
(3,1,1)     56: (1,1,4)    112: (1,1,5)
(2,2,1)     67: (5,1,1)    131: (6,1,1)
(2,1,2)     69: (4,2,1)    133: (5,2,1)
(1,3,1)     70: (4,1,2)    134: (5,1,2)
(1,2,2)     73: (3,3,1)    137: (4,3,1)
(1,1,3)     74: (3,2,2)    138: (4,2,2)
(4,1,1)     76: (3,1,3)    140: (4,1,3)
(3,2,1)     81: (2,4,1)    145: (3,4,1)
(3,1,2)     82: (2,3,2)    146: (3,3,2)
(2,3,1)     84: (2,2,3)    148: (3,2,3)
(2,2,2)     88: (2,1,4)    152: (3,1,4)
(End)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Robert Baillie, Summing the curious series of Kempner and Irwin, arXiv:0806.4410 [math.CA], 2008-2015. See p. 18 for Mathematica code irwinSums.m.
Stephen Morley, HAKMEM Item 175 (Gosper).
Tilman Piesk, First 56 elements in a tetrahedral array.

Cf. A038465 (base 3), A038471 (base 4), A038475 (base 5).
Cf. A081091 (primes), A212190 (squares), A212192 (triangular numbers), A173589 (Fibbinary).
Cf. A057168.
Cf. A000079, A018900, A014311, A014312, A014313, A023688, A023689, A023690, A023691 (Hammingweight = 1, 2, ..., 9).
Cf. A056558, A194848, A194847.
A000217(n-2) counts compositions into three parts.
A001399(n-3) = A069905(n) = A211540(n+2) counts the unordered case.
A001399(n-6) = A069905(n-3) = A211540(n-1) counts the unordered strict case.
A001399(n-6)*6 = A069905(n-3)*6 = A211540(n-1)*6 counts the strict case.
A014612 is an unordered version, with strict case A007304.
A337453 is the strict case.
A337461 counts the coprime case.
A033992 lists numbers divisible by exactly three different primes.
A323024 lists numbers with exactly three different prime multiplicities.
Cf. A220377, A307534, A337459, A337460, A337561, A337603, A337604.

unsigned hakmem175(unsigned x) {
    unsigned s, o, r;
    s = x & -x;  r = x + s;
    o = r ^ x;  o = (o >> 2) / s;
    return r | o;
}
unsigned A014311(int n) {
    if (n == 1) return 7;
    return hakmem175(A014311(n - 1));
}  // Peter Luschny, Jan 01 2014

a014311 n = a014311_list !! (n-1)
a014311_list = [2^x + 2^y + 2^z |
                x <- [2..], y <- [1..x-1], z <- [0..y-1]]
-- Reinhard Zumkeller, May 03 2012

Select[Range[200], (Count[IntegerDigits[#, 2], 1] == 3)&]
nn = 8; Flatten[Table[2^i + 2^j + 2^k, {i, 2, nn}, {j, 1, i - 1}, {k, 0, j - 1}]] (* T. D. Noe, Nov 05 2013 *)

for(n=0,10^3,if(hammingweight(n)==3,print1(n,", "))); \\ Joerg Arndt, Mar 04 2014

print1(t=7);for(i=2,50,print1(","t=A057168(t))) \\ M. F. Hasler, Aug 27 2014

A014311_list = [2**a+2**b+2**c for a in range(2,6) for b in range(1,a) for c in range(b)] # Chai Wah Wu, Jan 24 2021

from itertools import islice
def A014311_gen(): # generator of terms
    yield (n:=7)
    while True: yield (n:=n^((a:=-n&n+1)|(a>>1)) if n&1 else ((n&~(b:=n+(a:=n&-n)))>>a.bit_length())^b)
A014311_list = list(islice(A014311_gen(),20)) # Chai Wah Wu, Mar 10 2025

from math import isqrt, comb
from sympy import integer_nthroot
def A014311(n): return (1<<(r:=n-1-comb((m:=integer_nthroot(6*n,3)[0])+(t:=(n>comb(m+2,3)))+1,3))-comb((k:=isqrt(b:=r+1<<1))+(b>k*(k+1)),2))+(1<<(a:=isqrt(s:=n-comb(m-(t^1)+2,3)<<1))+((s<<2)>(a<<2)*(a+1)+1))+(1<Chai Wah Wu, Mar 10 2025

A000120(a(n)) = 3. - Reinhard Zumkeller, May 03 2012
Start with A084468. If n is in sequence, then 2n is too. - Ralf Stephan, Aug 16 2013
a(n+1) = A057168(a(n)). - M. F. Hasler, Aug 27 2014
a(n) = 2^A056558(n-1) + 2^A194848(n-1) + 2^A194847(n-1). - Ridouane Oudra, Sep 06 2020
Sum_{n>=1} 1/a(n) = A367110 = 1.428591545852638123996854844400537952781688750906133068397189529775365950039... (calculated using Baillie's irwinSums.m, see Links). - Amiram Eldar, Feb 14 2022

Extension and program by Olivier Gérard

A030059 Numbers that are the product of an odd number of distinct primes.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 30, 31, 37, 41, 42, 43, 47, 53, 59, 61, 66, 67, 70, 71, 73, 78, 79, 83, 89, 97, 101, 102, 103, 105, 107, 109, 110, 113, 114, 127, 130, 131, 137, 138, 139, 149, 151, 154, 157, 163, 165, 167, 170, 173, 174, 179, 181, 182, 186, 190, 191, 193

Offset: 1

N. J. A. Sloane

From Enrique Pérez Herrero, Jul 06 2012: (Start)
This sequence and A030229 partition the squarefree numbers: A005117.
Also solutions to the equation mu(n) = -1.
Sum_{n>=1} 1/a(n)^s = (zeta(s)^2 - zeta(2*s))/(2*zeta(s)*zeta(2*s)). (End) [See A088245 and the Hardy reference. - Wolfdieter Lang, Oct 18 2016]
The lexicographically least sequence of integers > 1 such that for each entry, the number of proper divisors occurring in the sequence is equal to 0 modulo 3. - Masahiko Shin, Feb 12 2018
The asymptotic density of this sequence is 3/Pi^2 (A104141). - Amiram Eldar, May 22 2020
Solutions to the equation Sum_{d|n} mu(d)*sigma(d) = -n, where sigma(n) is the sum of divisors function (A000203). - Robert D. Rosales, May 20 2024

B. C. Berndt & R. A. Rankin, Ramanujan: Letters and Commentary, see p. 23; AMS Providence RI 1995.
G. H. Hardy, Ramanujan, AMS Chelsea Publishing, 2002, pp. 64 - 65, (misprint on p. 65, line starting with Hence: it should be ... -1/Zeta(s) not ... -Zeta(s)).
S. Ramanujan, Collected Papers, Ed. G. H. Hardy et al., Cambridge 1927; Chelsea, NY, 1962, p. xxiv, 21.

T. D. Noe, Table of n, a(n) for n = 1..1000
Debmalya Basak, Nicolas Robles, and Alexandru Zaharescu, Exponential sums over Möbius convolutions with applications to partitions, arXiv:2312.17435 [math.NT], 2023. Mentions this sequence.
S. Ramanujan, Irregular numbers, J. Indian Math. Soc. 5 (1913) 105-106.
Eric Weisstein's World of Mathematics, Prime Factor
Eric Weisstein's World of Mathematics, Moebius Function
Eric Weisstein's World of Mathematics, Prime Sums
H. S. Wilf, A Greeting; and a view of Riemann's Hypothesis, Amer. Math. Monthly, 94:1 (1987), 3-6.

Cf. A000040, A001221, A005408, A007304, A008683, A030231, A051270, A123321, A030229, A005117, A088245, A104141.

a := n -> `if`(numtheory[mobius](n)=-1,n,NULL); seq(a(i),i=1..193); # Peter Luschny, May 04 2009
# alternative
A030059 := proc(n)
    option remember;
    local a;
    if n = 1 then
        2;
    else
        for a from procname(n-1)+1 do
            if numtheory[mobius](a) = -1 then
                return a;
            end if;
        end do:
    end if;
end proc: # R. J. Mathar, Sep 22 2020

Select[Range[300], MoebiusMu[#] == -1 &] (* Enrique Pérez Herrero, Jul 06 2012 *)

is(n)=my(f=factor(n)[,2]); #f%2 && vecmax(f)==1 \\ Charles R Greathouse IV, Oct 16 2015

is(n)=moebius(n)==-1 \\ Charles R Greathouse IV, Jan 31 2017

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A030059(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b+1,isqrt(x//c)+1),a+1)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b+1,integer_nthroot(x//c,m)[0]+1),a+1) for d in g(x,a2,b2,c*b2,m-1)))
    def f(x): return int(n+x-primepi(x)-sum(sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,i)) for i in range(3,x.bit_length(),2)))
    kmin, kmax = 0,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 # Chai Wah Wu, Aug 29 2024

omega(a(n)) = A001221(a(n)) gives A005408. {primes A000040} UNION {sphenic numbers A007304} UNION {numbers that are divisible by exactly 5 different primes A051270} UNION {products of 7 distinct primes (squarefree 7-almost primes) A123321} UNION {products of 9 distinct primes; also n has exactly 9 distinct prime factors and n is squarefree A115343} UNION.... - Jonathan Vos Post, Oct 19 2007

a(n) < n*Pi^2/3 infinitely often; a(n) > n*Pi^2/3 infinitely often. - Charles R Greathouse IV, Sep 07 2017

More terms from David W. Wilson

A085986 Squares of the squarefree semiprimes (p^2*q^2).

Original entry on oeis.org

36, 100, 196, 225, 441, 484, 676, 1089, 1156, 1225, 1444, 1521, 2116, 2601, 3025, 3249, 3364, 3844, 4225, 4761, 5476, 5929, 6724, 7225, 7396, 7569, 8281, 8649, 8836, 9025, 11236, 12321, 13225, 13924, 14161, 14884, 15129, 16641, 17689, 17956, 19881

Offset: 1

Alford Arnold, Jul 06 2003

This sequence is a member of a family of sequences directly related to A025487. First terms and known sequences are listed below: 1, A000007; 2, A000040; 4, A001248; 6, A006881; 8, A030078; 12, A054753; 16, A030514; 24, A065036; 30, A007304; 32, A050997; 36, this sequence; 48, ?; 60, ?; 64, ?; ....
Subsequence of A077448. The numbers in A077448 but not in here are 1, the squares of A046386, the squares of A067885, etc. - R. J. Mathar, Sep 12 2008
a(4)-a(3)=29 and a(3)+a(4)=421 are both prime.  There are no other cases where the sum and difference of two members of this sequence are both prime. - Robert Israel and J. M. Bergot, Oct 25 2019

			A006881 begins 6 10 14 15 ... so this sequence begins 36 100 196 225 ...

Subsequence of A036785 and of A077448.
Subsequence of A062503.
Cf. A025487.
Cf. A085548, A085964.

[k^2:k in [1..150]| IsSquarefree(k) and #PrimeDivisors(k) eq 2]; // Marius A. Burtea, Oct 24 2019

f[n_]:=Sort[Last/@FactorInteger[n]]=={2,2}; Select[Range[20000], f] (* Vladimir Joseph Stephan Orlovsky, Aug 14 2009 *)
Select[Range[200],PrimeOmega[#]==2&&SquareFreeQ[#]&]^2 (* Harvey P. Dale, Mar 07 2013 *)

list(lim)=my(v=List(), x=sqrtint(lim\=1), t); forprime(p=2, x\2, t=p; forprime(q=2, min(x\t,p-1), listput(v, (t*q)^2))); Set(v) \\ Charles R Greathouse IV, Sep 22 2015

is(n)=factor(n)[,2]==[2,2]~ \\ Charles R Greathouse IV, Oct 19 2015

from math import isqrt
from sympy import primepi, primerange
def A085986(n):
    def f(x): return int(n+x+(t:=primepi(s:=isqrt(x)))+(t*(t-1)>>1)-sum(primepi(x//k) for k in primerange(1, s+1)))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return m**2 # Chai Wah Wu, Aug 18 2024

a(n) = A006881(n)^2.

Sum_{n>=1} 1/a(n) = (P(2)^2 - P(4))/2 = (A085548^2 - A085964)/2 = 0.063767..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020

A046387 Products of exactly 5 distinct primes.

Original entry on oeis.org

2310, 2730, 3570, 3990, 4290, 4830, 5610, 6006, 6090, 6270, 6510, 6630, 7410, 7590, 7770, 7854, 8610, 8778, 8970, 9030, 9282, 9570, 9690, 9870, 10010, 10230, 10374, 10626, 11130, 11310, 11730, 12090, 12210, 12390, 12558, 12810, 13090, 13110

Offset: 1

Patrick De Geest, Jun 15 1998

Subsequence of A051270. 4620 = 2^2*3*5*7*11 is in A051270 but not in here, for example. - R. J. Mathar, Nov 10 2014

			a(1) = 2310 = 2 * 3 * 5 * 7 * 11 = A002110(5) = 5#.
a(2) = 2730 = 2 * 3 * 5 * 7 * 13.
a(3) = 3570 = 2 * 3 * 5 * 7 * 17.
a(10) = 6006 = 2 * 3 * 7 * 11 * 13.

T. D. Noe, Table of n, a(n) for n = 1..10000

Products of exactly k distinct primes, for k = 1 to 6: A000040, A006881. A007304, A046386, A046387, A067885.
Cf. A000040, A000961, A001221, A005117, A000977, A002110, A006881, A007304, A007774, A033992, A033993, A046386.
Cf. A014614, A046403, A051270.

A046387 := proc(n)
    option remember;
    local a;
    if n = 1 then
        2*3*5*7*11 ;
    else
        for a from procname(n-1)+1 do
            if A001221(a)= 5 and issqrfree(a) then
                return a;
            end if;
        end do:
    end if;
end proc: # R. J. Mathar, Oct 13 2019

f5Q[n_]:=Last/@FactorInteger[n]=={1, 1, 1, 1, 1}; lst={};Do[If[f5Q[n], AppendTo[lst, n]], {n, 8!}];lst (* Vladimir Joseph Stephan Orlovsky, Aug 26 2008 *)

is(n)=factor(n)[,2]==[1,1,1,1,1]~ \\ Charles R Greathouse IV, Sep 17 2015

is(n)= omega(n)==5 && bigomega(n)==5 \\ Hugo Pfoertner, Dec 18 2018

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A046387(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b+1,isqrt(x//c)+1),a+1)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b+1,integer_nthroot(x//c,m)[0]+1),a+1) for d in g(x,a2,b2,c*b2,m-1)))
    def f(x): return int(n+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,5)))
    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
    return bisection(f) # Chai Wah Wu, Aug 30 2024

Entry revised by N. J. A. Sloane, Apr 10 2006

A220377 Number of partitions of n into three distinct and mutually relatively prime parts.

Original entry on oeis.org

1, 0, 2, 1, 3, 1, 6, 1, 7, 3, 7, 3, 14, 3, 15, 6, 14, 6, 25, 6, 22, 10, 25, 9, 42, 8, 34, 15, 37, 15, 53, 13, 48, 22, 53, 17, 78, 17, 65, 30, 63, 24, 99, 24, 88, 35, 84, 30, 126, 34, 103, 45, 103, 38, 166, 35, 124, 57, 128, 51, 184, 44, 150, 67, 172, 52, 218

Offset: 6

Carl Najafi, Dec 13 2012

nonn

The Heinz numbers of these partitions are the intersection of A005117 (strict), A014612 (triples), and A302696 (coprime). - Gus Wiseman, Oct 14 2020

			For n=10 we have three such partitions: 1+2+7, 1+4+5 and 2+3+5.
From _Gus Wiseman_, Oct 14 2020: (Start)
The a(6) = 1 through a(20) = 15 triples (empty column indicated by dot, A..H = 10..17):
321  .  431  531  532  731  543  751  743  753  754  971  765  B53  875
        521       541       651       752  951  853  B51  873  B71  974
                  721       732       761  B31  871  D31  954  D51  A73
                            741       851       952       972       A91
                            831       941       B32       981       B54
                            921       A31       B41       A71       B72
                                      B21       D21       B43       B81
                                                          B52       C71
                                                          B61       D43
                                                          C51       D52
                                                          D32       D61
                                                          D41       E51
                                                          E31       F41
                                                          F21       G31
                                                                    H21
(End)

Fausto A. C. Cariboni, Table of n, a(n) for n = 6..10000 (terms 6..1000 from Seiichi Manyama)

Cf. A015617, A300815.
A023022 is the 2-part version.
A101271 is the relative prime instead of pairwise coprime version.
A220377*6 is the ordered version.
A305713 counts these partitions of any length, with Heinz numbers A302797.
A307719 is the non-strict version.
A337461 is the non-strict ordered version.
A337563 is the case with no 1's.
A337605 is the pairwise non-coprime instead of pairwise coprime version.
A001399(n-6) counts strict 3-part partitions, with Heinz numbers A007304.
A008284 counts partitions by sum and length, with strict case A008289.
A318717 counts pairwise non-coprime strict partitions.
A326675 ranks pairwise coprime sets.
A327516 counts pairwise coprime partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.
Cf. A000217, A007360, A023023, A051424, A078374, A087087, A302696, A333227, A337485, A337561.

Table[Length@Select[ IntegerPartitions[ n, {3}], #[[1]] != #[[2]] != #[[3]] && GCD[#[[1]], #[[2]]] == 1 && GCD[#[[1]], #[[3]]] == 1 && GCD[#[[2]], #[[3]]] == 1 &], {n, 6, 100}]
Table[Count[IntegerPartitions[n,{3}],?(CoprimeQ@@#&&Length[ Union[#]] == 3&)],{n,6,100}] (* _Harvey P. Dale, May 22 2020 *)

a(n)=my(P=partitions(n));sum(i=1,#P,#P[i]==3&&P[i][1]Charles R Greathouse IV, Dec 14 2012

a(n > 2) = A307719(n) - 1. - Gus Wiseman, Oct 15 2020

A067885 Products of exactly 6 distinct primes.

Original entry on oeis.org

30030, 39270, 43890, 46410, 51870, 53130, 62790, 66990, 67830, 71610, 72930, 79170, 81510, 82110, 84630, 85470, 91770, 94710, 98670, 99330, 101010, 102102, 103530, 106590, 108570, 110670, 111930, 114114, 115710, 117390, 122430, 123690, 124410, 125970, 128310

Offset: 1

Benoit Cloitre, Mar 02 2002

nonn

Subsequence of A074969. - R. J. Mathar, Nov 24 2009

Products of exactly k distinct primes, for k = 1 to 6: A000040, A006881. A007304, A046386, A046387, A067885.

Select[Range[125000],PrimeNu[#]==PrimeOmega[#]==6&] (* Harvey P. Dale, May 14 2014 *)

is(n)=factor(n)[,2]==[1,1,1,1,1,1]~ \\ Charles R Greathouse IV, Sep 14 2015

is(n)=omega(n)==6 && bigomega(n)==6 \\ Hugo Pfoertner, Dec 18 2018

list(lim)=lim\=1; my(v=List(), L1,L2,L3,L4,P4,P5); forprime(p=13,lim\2310, L1=lim\p; forprime(q=11,min(L1\210,p-2), L2=L1\q; forprime(r=7, min(L2\30,q-2), L3=L2\r; forprime(s=5,min(L3\6,r-2), L4=L3\s; P4=p*q*r*s; forprime(t=3, min(L4\2,s-2), P5=P4*t; forprime(u=2, min(L4\t,t-1), listput(v,P5*u))))))); Set(v) \\ Charles R Greathouse IV, Aug 27 2021

from math import prod, isqrt
from sympy import primerange, integer_nthroot, primepi
def A067885(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b+1,isqrt(x//c)+1),a+1)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b+1,integer_nthroot(x//c,m)[0]+1),a+1) for d in g(x,a2,b2,c*b2,m-1)))
    def f(x): return int(n+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,6)))
    kmin, kmax = 0,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 # Chai Wah Wu, Aug 29 2024

{k: A001221(k) = A001222(k) = 6}. - R. J. Mathar, Jul 18 2023

A085987 Product of exactly four primes, three of which are distinct (p^2qr).

Original entry on oeis.org

60, 84, 90, 126, 132, 140, 150, 156, 198, 204, 220, 228, 234, 260, 276, 294, 306, 308, 315, 340, 342, 348, 350, 364, 372, 380, 414, 444, 460, 476, 490, 492, 495, 516, 522, 525, 532, 550, 558, 564, 572, 580, 585, 620, 636, 644, 650, 666, 693, 708, 726

Offset: 1

Alford Arnold, Jul 08 2003

nonn

A014613 is completely determined by A030514, A065036, A085986, A085987 and A046386 since p(4) = 5. (cf. A000041). More generally, the first term of sequences which completely determine the k-almost primes can be found in A036035 (a resorted version of A025487).

A050326(a(n)) = 4. - Reinhard Zumkeller, May 03 2013

			a(1) = 60 since 60 = 2*2*3*5 and has three distinct prime factors.

Cf. A001248, A006881, A030078, A054753, A007304, A050997, A046387, A036035, A086974.

Subsequence of A014613, A307341, A178212.

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,1,2}; Select[Range[2000], f] (* Vladimir Joseph Stephan Orlovsky, May 03 2011 *)
pefp[{a_,b_,c_}]:={a^2 b c,a b^2 c,a b c^2}; Module[{upto=800},Select[ Flatten[ pefp/@Subsets[Prime[Range[PrimePi[upto/6]]],{3}]]//Union,#<= upto&]] (* Harvey P. Dale, Oct 02 2018 *)

list(lim)=my(v=List(),t,x,y,z);forprime(p=2,lim^(1/4),t=lim\p^2;forprime(q=p+1,sqrtint(t),forprime(r=q+1,t\q,x=p^2*q*r;y=p*q^2*r;listput(v,x);if(y<=lim,listput(v,y);z=p*q*r^2;if(z<=lim,listput(v,z))))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 15 2011

is(n)=vecsort(factor(n)[,2]~)==[1,1,2] \\ Charles R Greathouse IV, Oct 19 2015

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A085987(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x+sum((t:=primepi(s:=isqrt(y:=x//r**2)))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)) for r in primerange(isqrt(x)+1))+sum(primepi(x//p**3) for p in primerange(integer_nthroot(x,3)[0]+1))-primepi(integer_nthroot(x,4)[0])
    return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025

More terms from Reinhard Zumkeller, Jul 25 2003

A307719 Number of partitions of n into 3 mutually coprime parts.

Original entry on oeis.org

0, 0, 0, 1, 1, 1, 2, 1, 3, 2, 4, 2, 7, 2, 8, 4, 8, 4, 15, 4, 16, 7, 15, 7, 26, 7, 23, 11, 26, 10, 43, 9, 35, 16, 38, 16, 54, 14, 49, 23, 54, 18, 79, 18, 66, 31, 64, 25, 100, 25, 89, 36, 85, 31, 127, 35, 104, 46, 104, 39, 167, 36, 125, 58, 129, 52, 185, 45

Offset: 0

Wesley Ivan Hurt, Apr 24 2019

The Heinz numbers of these partitions are the intersection of A014612 (triples) and A302696 (pairwise coprime). - Gus Wiseman, Oct 16 2020

			There are 2 partitions of 9 into 3 mutually coprime parts: 7+1+1 = 5+3+1, so a(9) = 2.
There are 4 partitions of 10 into 3 mutually coprime parts: 8+1+1 = 7+2+1 = 5+4+1 = 5+3+2, so a(10) = 4.
There are 2 partitions of 11 into 3 mutually coprime parts: 9+1+1 = 7+3+1, so a(11) = 2.
There are 7 partitions of 12 into 3 mutually coprime parts: 10+1+1 = 9+2+1 = 8+3+1 = 7+4+1 = 6+5+1 = 7+3+2 = 5+4+3, so a(12) = 7.

Fausto A. C. Cariboni, Table of n, a(n) for n = 0..10000 (terms 0..2000 from Robert Israel)
Index entries for sequences related to partitions

A023022 is the version for pairs.
A220377 is the strict case, with ordered version A220377*6.
A327516 counts these partitions of any length, with strict version A305713 and Heinz numbers A302696.
A337461 is the ordered version.
A337563 is the case with no 1's.
A337599 is the pairwise non-coprime instead of pairwise coprime version.
A337601 only requires the distinct parts to be pairwise coprime.
A001399(n-3) = A069905(n) = A211540(n+2) counts 3-part partitions.
A002865 counts partitions with no 1's, with strict case A025147.
A007359 and A337485 count pairwise coprime partitions with no 1's.
A200976 and A328673 count pairwise non-coprime partitions.
Cf. A007304, A014612, A078374, A082024, A284825, A304709, A337462, A337605.

N:= 200: # to get a(0)..a(N)
A:= Array(0..N):
for a from 1 to N/3 do
  for b from a to (N-a)/2 do
    if igcd(a,b) > 1 then next fi;
    ab:= a*b;
    for c from b to N-a-b do
       if igcd(ab,c)=1 then A[a+b+c]:= A[a+b+c]+1 fi
od od od:
convert(A,list); # Robert Israel, May 09 2019

Table[Sum[Sum[Floor[1/(GCD[i, j] GCD[j, n - i - j] GCD[i, n - i - j])], {i, j, Floor[(n - j)/2]}], {j, Floor[n/3]}], {n, 0, 100}]
Table[Length[Select[IntegerPartitions[n,{3}],CoprimeQ@@#&]],{n,0,100}] (* Gus Wiseman, Oct 15 2020 *)

a(n) = Sum_{j=1..floor(n/3)} Sum_{i=j..floor((n-j)/2)} [gcd(i,j) * gcd(j,n-i-j) * gcd(i,n-i-j) = 1], where [] is the Iverson bracket.

a(n > 2) = A220377(n) + 1. - Gus Wiseman, Oct 15 2020

A046389 Products of exactly three distinct odd primes.

Original entry on oeis.org

105, 165, 195, 231, 255, 273, 285, 345, 357, 385, 399, 429, 435, 455, 465, 483, 555, 561, 595, 609, 615, 627, 645, 651, 663, 665, 705, 715, 741, 759, 777, 795, 805, 861, 885, 897, 903, 915, 935, 957, 969, 987, 1001, 1005, 1015, 1023, 1045, 1065, 1085

Offset: 1

Patrick De Geest, Jun 15 1998

nonn

The old name was "Odd numbers with exactly 3 distinct prime factors", which is slightly ambiguous, since it might be interpreted to include 315 = 3^2*5*7. Cf. A278569. - N. J. A. Sloane, Nov 27 2016

T. D. Noe, Table of n, a(n) for n=1..1000

Cf. A046316, A046405, A007304, A278569.

f[n_] := OddQ[n] && Last/@FactorInteger[n]=={1,1,1}; lst={}; Do[If[f[n], AppendTo[lst,n]], {n, 2000}]; lst (* Vladimir Joseph Stephan Orlovsky, Nov 23 2009 *)

list(lim)=my(v=List(), t); forprime(p=3, lim^(1/3), forprime(q=p+1, sqrt(lim\p), t=p*q; forprime(r=q+1, lim\t, listput(v, t*r)))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 26 2011

from math import isqrt
from sympy import primepi, integer_nthroot, primerange
def A046389(n):
    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): return int(n+x-sum(primepi(x//(k*m))-b for a,k in enumerate(primerange(3,integer_nthroot(x,3)[0]+1),2) for b,m in enumerate(primerange(k+1,isqrt(x//k)+1),a+1)))
    return bisection(f,n,n) # Chai Wah Wu, Sep 10 2024

Name clarified by N. J. A. Sloane, Nov 27 2016

A030626 Numbers with exactly 8 divisors.

Original entry on oeis.org

24, 30, 40, 42, 54, 56, 66, 70, 78, 88, 102, 104, 105, 110, 114, 128, 130, 135, 136, 138, 152, 154, 165, 170, 174, 182, 184, 186, 189, 190, 195, 222, 230, 231, 232, 238, 246, 248, 250, 255, 258, 266, 273, 282, 285, 286, 290, 296, 297, 310, 318, 322, 328, 344, 345, 351, 354, 357, 366, 370, 374, 375, 376, 385, 399, 402

Offset: 1

Jeff Burch

nonn

Since A119479(8)=7, there are never more than 7 consecutive terms. Runs of 7 consecutive terms start at 171897, 180969, 647385, ... (subsequence of A049053). - Ivan Neretin, Feb 08 2016

Jon E. Schoenfield, Table of n, a(n) for n = 1..10000 (first 1000 terms from R. J. Mathar)
Jérôme Germoni, Nombres à huit diviseurs, Images des Mathématiques, CNRS, 2017 (in French).
Eric Weisstein's World of Mathematics, Divisor Product.

Essentially the same as A111398.

Cf. A000005, A007304, A049053, A065036, A092759, A119479.

[n: n in [1..400] | DivisorSigma(0, n) eq 8]; // Vincenzo Librandi, Oct 05 2017

select(numtheory:-tau=8, [$1..1000]); # Robert Israel, Dec 17 2014

Select[Range[400], DivisorSigma[0, #]== 8 &] (* Vincenzo Librandi, Oct 05 2017 *)

Vec(select(x->x==8,vector(500, i, numdiv(i)),1)) \\ Michel Marcus, Dec 17 2014

from sympy import divisor_count
isok = lambda n: divisor_count(n) == 8
print([n for n in range(1, 400) if isok(n)]) # Darío Clavijo, Oct 17 2023

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A030626(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return int(n+x-sum(primepi(x//(k*m))-b for a,k in enumerate(primerange(integer_nthroot(x,3)[0]+1),1) for b,m in enumerate(primerange(k+1,isqrt(x//k)+1),a+1))-sum(primepi(x//p**3) for p in primerange(integer_nthroot(x,3)[0]+1))+primepi(integer_nthroot(x,4)[0])-primepi(integer_nthroot(x,7)[0]))
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

A000005(a(n))=8. - Juri-Stepan Gerasimov, Oct 10 2009

Equals A065036 (p*q^3) U A007304 (p*q*r) U A092759 (p^7). - Amarnath Murthy, Apr 21 2001

cp's OEIS Frontend

A014311 Numbers with exactly 3 ones in binary expansion.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

C

Haskell

Mathematica

PARI

PARI

Python

Python

Python

Formula

Extensions

A030059 Numbers that are the product of an odd number of distinct primes.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A085986 Squares of the squarefree semiprimes (p^2*q^2).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Python

Formula

A046387 Products of exactly 5 distinct primes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Extensions

A220377 Number of partitions of n into three distinct and mutually relatively prime parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A067885 Products of exactly 6 distinct primes.

Views

Author

A085987 Product of exactly four primes, three of which are distinct (p^2qr).