A014612 -id:A014612 - cp's OEIS frontend

A069278 17-almost primes (generalization of semiprimes).

131072, 196608, 294912, 327680, 442368, 458752, 491520, 663552, 688128, 720896, 737280, 819200, 851968, 995328, 1032192, 1081344, 1105920, 1114112, 1146880, 1228800, 1245184, 1277952, 1492992, 1507328, 1548288, 1605632, 1622016

Offset: 1

Rick L. Shepherd, Mar 13 2002

nonn

Product of 17 not necessarily distinct primes.
Divisible by exactly 17 prime powers (not including 1).
For n = 1..2628 a(n)=2*A069277(n). - Zak Seidov, Jun 25 2017

D. W. Wilson, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Almost Prime.

Sequences listing r-almost primes, that is, the n such that A001222(n) = r: A000040 (r = 1), A001358 (r = 2), A014612 (r = 3), A014613 (r = 4), A014614 (r = 5), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), A069272 (r = 11), A069273 (r = 12), A069274 (r = 13), A069275 (r = 14), A069276 (r = 15), A069277 (r = 16), this sequence (r = 17), A069279 (r = 18), A069280 (r = 19), A069281 (r = 20). - Jason Kimberley, Oct 02 2011

Select[Range[2*10^6],PrimeOmega[#]==17&] (* Harvey P. Dale, Sep 28 2016 *)

k=17; start=2^k; finish=2000000; v=[]
for(n=start,finish, if(bigomega(n)==k,v=concat(v,n))); v

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A069278(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) 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,17)))
    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 31 2024

Product p_i^e_i with Sum e_i = 17.

A069280 19-almost primes (generalization of semiprimes).

Original entry on oeis.org

524288, 786432, 1179648, 1310720, 1769472, 1835008, 1966080, 2654208, 2752512, 2883584, 2949120, 3276800, 3407872, 3981312, 4128768, 4325376, 4423680, 4456448, 4587520, 4915200, 4980736, 5111808, 5971968, 6029312, 6193152

Offset: 1

Rick L. Shepherd, Mar 13 2002

nonn

Product of 19 not necessarily distinct primes.

Divisible by exactly 19 prime powers (not including 1).

D. W. Wilson, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Almost Prime.

Sequences listing r-almost primes, that is, the n such that A001222(n) = r: A000040 (r = 1), A001358 (r = 2), A014612 (r = 3), A014613 (r = 4), A014614 (r = 5), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), A069272 (r = 11), A069273 (r = 12), A069274 (r = 13), A069275 (r = 14), A069276 (r = 15), A069277 (r = 16), A069278 (r = 17), A069279 (r = 18), this sequence (r = 19), A069281 (r = 20). - Jason Kimberley, Oct 02 2011

k=19; start=2^k; finish=8000000; v=[]; for(n=start,finish, if(bigomega(n)==k,v=concat(v,n))); v

from math import prod, isqrt
from sympy import primerange, integer_nthroot, primepi
def A069280(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
    def f(x): return int(n-1+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,19)))
    kmin, kmax = 1,2
    while f(kmax) >= kmax:
        kmax <<= 1
    while True:
        kmid = kmax+kmin>>1
        if f(kmid) < kmid:
            kmax = kmid
        else:
            kmin = kmid
        if kmax-kmin <= 1:
            break
    return kmax # Chai Wah Wu, Aug 23 2024

Product p_i^e_i with Sum e_i = 19.

A062069 a(n) = sigma(d(n)), where d(k) is the number of divisors function (A000005) and sigma(k) is the sum of divisors function (A000203).

Original entry on oeis.org

1, 3, 3, 4, 3, 7, 3, 7, 4, 7, 3, 12, 3, 7, 7, 6, 3, 12, 3, 12, 7, 7, 3, 15, 4, 7, 7, 12, 3, 15, 3, 12, 7, 7, 7, 13, 3, 7, 7, 15, 3, 15, 3, 12, 12, 7, 3, 18, 4, 12, 7, 12, 3, 15, 7, 15, 7, 7, 3, 28, 3, 7, 12, 8, 7, 15, 3, 12, 7, 15, 3, 28, 3, 7, 12, 12, 7, 15, 3, 18, 6, 7, 3, 28, 7, 7, 7, 15, 3

Offset: 1

Amarnath Murthy, Jun 13 2001

nonn

a(1) = 1, a(p) = 3 for p = primes (A000040), a(pq) = 7 for pq = product of two distinct primes (A006881), a(pq...z) = 2^(k+1)-1 = A000225(k+1) for pq...z = product of k (k > 2) distinct primes p,q,...,z (A120944), a(p^k) = sigma(k+1) = A000203(k+1) for p^k = prime powers (A000961(n) for n > 1). Sequence {1,3,4,12} is finite sequence of numbers n such that sigma(tau(n)) = n. [Jaroslav Krizek, Jul 16 2009]
For semiprime n, a(n) is either 4 or 7. Also a(n) = d(n) + omega(n) + mu(n), the sum of three core sequences A000005, A001221 and A008683. When n is semiprime, a(n) is completely defined by the Mobius function as: a(n) = 4 + 3*mu(n). a(n) also has the fractal-like identities a(d(n)) = d(n) and a(n) = sigma(a(d(n))). - Wesley Ivan Hurt, Sep 02 2013
If n is a triprime (A014612), d(n) is 4, 6, or 8 and a(n) = sigma(d(n)) is 7, 12, or 15 respectively. Then a(n) = -d(n)^2/4 + 5*d(n) - 9. - Wesley Ivan Hurt, Sep 08 2013

			sigma(d(12)) = sigma(6) = 12.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harry J. Smith)

Cf. A000005, A000203, A001221, A008683, A062068.

A062069:= (n-> numtheory[sigma](numtheory[tau](n))):
seq (A062069(n), n=1..40); # Jani Melik, Jan 25 2011

Table[DivisorSigma[1, DivisorSigma[0, n]], {n, 1, 80}] (* Carl Najafi, Aug 16 2011 *)

v=[]; for(n=1,150,v=concat(v, sigma(numdiv(n)))); v

{ for (n=1, 1000, write("b062069.txt", n, " ", sigma(numdiv(n))) ) } \\ Harry J. Smith, Jul 31 2009

a(n) = A000203(A000005(n)). - Wesley Ivan Hurt, Sep 09 2013

More terms from Jason Earls, Jun 19 2001

A101605 a(n) = 1 if n is a product of exactly 3 (not necessarily distinct) primes, otherwise 0.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0

Offset: 1

Jonathan Vos Post, Dec 09 2004

			a(28) = 1 because 28 = 2 * 2 * 7 is the product of exactly 3 primes, counted with multiplicity.

Antti Karttunen, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Almost Prime.
Index entries for characteristic functions
Index entries for sequences computed from exponents in factorization of n

Cf. A010051, A064911, (char funct. of) A014612, A101637, A123074.

A101605 := proc(n)
    if numtheory[bigomega](n) = 3 then
        1;
    else
        0;
    end if;
end proc: # R. J. Mathar, Mar 13 2015

Table[Boole[PrimeOmega[n] == 3], {n, 100}] (* Jean-François Alcover, Mar 23 2020 *)

is(n)=bigomega(n)==3 \\ Charles R Greathouse IV, Apr 25 2016

a(n) = 1 if n has exactly three prime factors (not necessarily distinct), else a(n) = 0. a(n) = 1 if n is an element of A014612, else a(n) = 0.

a(n) = floor(Omega(n)/3) * floor(3/Omega(n)). - Wesley Ivan Hurt, Jan 10 2013

Description clarified by Antti Karttunen, Jul 23 2017

A211540 Number of ordered triples (w,x,y) with all terms in {1..n} and 2w = 3x + 4y.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 2, 3, 4, 5, 7, 8, 10, 12, 14, 16, 19, 21, 24, 27, 30, 33, 37, 40, 44, 48, 52, 56, 61, 65, 70, 75, 80, 85, 91, 96, 102, 108, 114, 120, 127, 133, 140, 147, 154, 161, 169, 176, 184, 192, 200, 208, 217, 225, 234, 243, 252, 261, 271, 280, 290

Offset: 0

Clark Kimberling, Apr 15 2012

For a guide to related sequences, see A211422.
Also the number of partitions of n+1 into three parts, where each part > 1. - Peter Woodward, May 25 2015
a(n) is also equal to the number of partitions of n+4 into three distinct parts, where each part > 1. - Giovanni Resta, May 26 2015
Number of different distributions of n+1 identical balls in 3 boxes as x,y,z where 0 < x < y < z. - Ece Uslu and Esin Becenen, Dec 31 2015
After the first three terms, partial sums of A008615. - Robert Israel, Dec 31 2015
For n >= 2, also the number of partitions of n - 2 into 3 parts. The Heinz numbers of these partitions are given by A014612. - Gus Wiseman, Oct 11 2020

			a(5) = a(6) = 1 with only one ordered triple (5,2,1). - _Michael Somos_, Feb 02 2015
a(11) = 5 Number of different distributions of 11 identical balls in 3 boxes as x,y and z where 0 < x < y < z. - _Ece Uslu_, Esin Becenen, Dec 31 2015
a(1) = a(2) = a(3) = a(4) = a(5) = 0, since with fewer than 6 identical balls there is no such distribution with 3 boxes that holds for 0 < x < y < z. - _Ece Uslu_, Esin Becenen, Dec 31 2015
G.f.: x^5 + x^6 + 2*x^7 + 3*x^8 + 4*x^9 + 5*x^10 + 7*x^11 + 8*x^12 + ...
From _Gus Wiseman_, Oct 11 2020: (Start)
The a(5) = 1 through a(15) = 14 partitions of n + 1 into three parts > 1 [Woodward] are the following (A = 10, B = 11, C = 12). The ordered version is A000217(n - 4) and the Heinz numbers are A046316.
  222  322  332  333  433  443  444  544  554  555  655
            422  432  442  533  543  553  644  654  664
                 522  532  542  552  643  653  663  754
                      622  632  633  652  662  744  763
                           722  642  733  743  753  772
                                732  742  752  762  844
                                822  832  833  843  853
                                     922  842  852  862
                                          932  933  943
                                          A22  942  952
                                               A32  A33
                                               B22  A42
                                                    B32
                                                    C22
The a(5) = 1 through a(15) = 14 partitions of n + 4 into three distinct parts > 1 [Resta] are the following (A = 10, B = 11, C = 12, D = 13, E = 14). The ordered version is A211540*6 and the Heinz numbers are A046389.
  432  532  542  543  643  653  654  754  764  765  865
            632  642  652  743  753  763  854  864  874
                 732  742  752  762  853  863  873  964
                      832  842  843  862  872  954  973
                           932  852  943  953  963  982
                                942  952  962  972  A54
                                A32  A42  A43  A53  A63
                                     B32  A52  A62  A72
                                          B42  B43  B53
                                          C32  B52  B62
                                               C42  C43
                                               D32  C52
                                                    D42
                                                    E32
The a(5) = 1 through a(15) = 14 partitions of n + 1 into three distinct parts [Uslu and Becenen] are the following (A = 10, B = 11, C = 12, D = 13). The ordered version is A211540(n)*6 and the Heinz numbers are A007304.
  321  421  431  432  532  542  543  643  653  654  754
            521  531  541  632  642  652  743  753  763
                 621  631  641  651  742  752  762  853
                      721  731  732  751  761  843  862
                           821  741  832  842  852  871
                                831  841  851  861  943
                                921  931  932  942  952
                                     A21  941  951  961
                                          A31  A32  A42
                                          B21  A41  A51
                                               B31  B32
                                               C21  B41
                                                    C31
                                                    D21
(End)

Robert Israel, Table of n, a(n) for n = 0..10000
Clark Kimberling, A Combinatorial Classification of Triangle Centers on the Line at Infinity, J. Int. Seq., Vol. 22 (2019), Article 19.5.4.
Ece Uslu and Esin Becenen, Identical Object Distributions.
Index entries for linear recurrences with constant coefficients, signature (1,1,0,-1,-1,1).

Cf. A001399, A069905, A211422.
All of the following pertain to 3-part strict partitions.
- A000009 counts these partitions of any length, with non-strict version A000041.
- A007304 gives the Heinz numbers, with non-strict version A014612.
- A101271 counts the relatively prime case, with non-strict version A023023.
- A220377 counts the pairwise coprime case, with non-strict version A307719.
- A337605 counts the pairwise non-coprime case, with non-strict version A337599.
Cf. A000217, A001840, A156040, A284825, A337453, A337483, A337484, A337563.

I:=[0,0,0,0,0,1]; [n le 6 select I[n] else Self(n-1)+Self(n-2)-Self(n-4)-Self(n-5)+Self(n-6): n in [1..70]]; // Vincenzo Librandi, Dec 31 2015

f:= gfun:-rectoproc({a(n) = a(n-1)+a(n-2)-a(n-4)-a(n-5)+a(n-6),seq(a(i)=0,i=0..4),a(5)=1},a(n),remember):
seq(f(i),i=0..100); # Robert Israel, Dec 31 2015

t[n_] := t[n] = Flatten[Table[-2 w + 3 x + 4 y, {w, n}, {x, n}, {y, n}]]
c[n_] := Count[t[n], 0]
t = Table[c[n], {n, 0, 80}]  (* A211540 *)
FindLinearRecurrence[t]
LinearRecurrence[{1, 1, 0, -1, -1, 1}, {0, 0, 0, 0, 0, 1}, 70] (* Vincenzo Librandi, Dec 31 2015 *)
Table[Length[Select[IntegerPartitions[n+1,{3}],UnsameQ@@#&]],{n,0,30}] (* Gus Wiseman, Oct 05 2020 *)

{a(n) = round( (n-2)^2 / 12 )}; / * Michael Somos, Feb 02 2015 */

concat(vector(5), Vec(x^5/(1-x-x^2+x^4+x^5-x^6) + O(x^100))) \\ Altug Alkan, Jan 10 2016

a(n) = a(n-1) + a(n-2) - a(n-4) - a(n-5) + a(n-6).
a(n) = A069905(n-2) = A001399(n-5) for n >= 5. - Alois P. Heinz, Nov 03 2012
a(n) = 3*k^2-6*k+3 (for n = 6*k-3), 3*k^2-5*k+2 (for n = 6*k-2), 3*k^2-4*k+1 (for n = 6*k-1), 3*k^2-3*k+1 (for n = 6*k), 3*k^2-2*k (for n = 6*k+1), 3*k^2-k (for n = 6*k+2). - Ece Uslu, Esin Becenen, Dec 31 2015
a(n) = A004526(n-2) + a(n-2) for n > 2. - Ece Uslu, Esin Becenen, Dec 31 2015
G.f.: x^5/(1 - x - x^2 + x^4 + x^5 - x^6). - Robert Israel, Dec 31 2015
a(n) = Sum_{k=1..floor(n/3)} floor((n-k)/2)-k. - Wesley Ivan Hurt, Apr 27 2019
From Gus Wiseman, Oct 11 2020: (Start)
a(n+2) = A069905(n) = A001399(n-3) counts 3-part partitions.
a(n-1) = A069905(n-3) = A001399(n-6) counts 3-part strict partitions.
a(n-1) = A069905(n-3) = A001399(n-6) counts 3-part partitions with no 1's.
a(n-4) = A069905(n-6) = A001399(n-9) counts 3-part strict partitions with no 1's.
A000217(n-2) counts 3-part compositions.
a(n-1)*6 = A069905(n-3)*6 = A001399(n-6)*6 counts 3-part strict compositions.
A000217(n-5) counts 3-part compositions with no 1's.
a(n-4)*6 = A069905(n-6)*6 = A001399(n-9)*6 counts 3-part strict compositions with no 1's.
(End)

A284825 Number of partitions of n into 3 parts without common divisors such that every pair of them has common divisors.

Original entry on oeis.org

1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 2, 0, 1, 0, 0, 0, 3, 0, 5, 0, 0, 0, 1, 0, 5, 0, 1, 0, 6, 0, 6, 0, 0, 0, 4, 0, 6, 0, 0, 0, 9, 0, 2, 1, 2, 0, 9, 0, 8, 1, 1, 0, 5, 0, 14, 0, 1, 0, 15, 0, 14, 0, 0, 1, 14, 0, 14, 0, 2, 0, 15, 0, 6, 1, 2, 1, 11, 0, 18, 1, 1, 0, 10, 0, 23

Offset: 31

Alois P. Heinz, Apr 03 2017

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

			a(31) = 1: [6,10,15] = [2*3,2*5,3*5].
a(41) = 2: [6,14,21], [6,15,20].
From _Gus Wiseman_, Oct 14 2020: (Start)
Selected terms and the corresponding triples:
  a(31)=1: a(41)=2: a(59)=3:  a(77)=4:  a(61)=5:  a(71)=6:
-------------------------------------------------------------
  15,10,6  20,15,6  24,20,15  39,26,12  33,22,6   39,26,6
           21,14,6  24,21,14  42,20,15  40,15,6   45,20,6
                    35,14,10  45,20,12  45,10,6   50,15,6
                              50,15,12  28,21,12  35,21,15
                                        36,15,10  36,20,15
                                                  36,21,14
(End)

Alois P. Heinz, Table of n, a(n) for n = 31..10000

Cf. A082024, A230034, A230035.
A023023 does not require pairwise non-coprimality, with strict case A101271.
A202425 and A328672 count these partitions of any length, ranked by A328868.
A284825*6 is the ordered version.
A307719 is the pairwise coprime instead of non-coprime version.
A337599 does not require relatively primality, with strict case A337605.
A200976 and A328673 count pairwise non-coprime partitions.
A289509 gives Heinz numbers of relatively prime partitions.
A327516 counts pairwise coprime partitions, ranked by A333227.
A337694 gives Heinz numbers of pairwise non-coprime partitions.
Cf. A000217, A000741, A001399, A007304, A014612, A220377, A318716, A328679, A337604, A337666.

a:= proc(n) option remember; add(add(`if`(igcd(i, j)>1
      and igcd(i, j, n-i-j)=1 and igcd(i, n-i-j)>1 and
      igcd(j, n-i-j)>1, 1, 0), j=i..(n-i)/2), i=2..n/3)
    end:
seq(a(n), n=31..137);

a[n_] := a[n] = Sum[Sum[If[GCD[i, j] > 1 && GCD[i, j, n - i - j] == 1 && GCD[i, n - i - j] > 1 && GCD[j, n - i - j] > 1, 1, 0], {j, i, (n - i)/2} ], {i, 2, n/3}];
Table[a[n], {n, 31, 137}] (* Jean-François Alcover, Jun 13 2018, from Maple *)
stabQ[u_,Q_]:=And@@Not/@Q@@@Tuples[u,2];
Table[Length[Select[IntegerPartitions[n,{3}],GCD@@#==1&&stabQ[#,CoprimeQ]&]],{n,31,100}] (* Gus Wiseman, Oct 14 2020 *)

a(n) > 0 iff n in { A230035 }.

a(n) = 0 iff n in { A230034 }.

A046363 Composite numbers whose sum of prime factors (with multiplicity) is prime.

Original entry on oeis.org

6, 10, 12, 22, 28, 34, 40, 45, 48, 52, 54, 56, 58, 63, 75, 76, 80, 82, 88, 90, 96, 99, 104, 108, 117, 118, 136, 142, 147, 148, 153, 165, 172, 175, 176, 184, 198, 202, 207, 210, 214, 224, 245, 248, 250, 252, 268, 273, 274, 279, 294, 296, 298, 300, 316, 320, 325

Offset: 1

Patrick De Geest, Jun 15 1998

nonn

If prime numbers were included the sequence would be 2, 3, 5, 6, 7, 10, 11, 12, 13, 17, 19, 22, 23, 28, 29, ... which is A100118. - Hieronymus Fischer, Oct 20 2007
Conjecture: a(n) can be approximated with the formula c*n^k, where c is approximately 0.46 and k is approximately 1.05. - Elijah Beregovsky, May 01 2019
The ternary Goldbach Conjecture implies that this sequence contains infinitely many terms of A014612 (triprimes). - Elijah Beregovsky, Dec 17 2019
A proof that this sequence is infinite: There are infinitely many odd primes, let p2 > p1 > 2 be two odd primes, p2-p1=2*k  then (2^k)*p1 is a term because 2*k+p1=p2 is prime. For example: 5+6=11, 6=2*3, 2^3*5=40 is a term. - Metin Sariyar, Dec 17 2019
Regarding the 2019 conjecture, with k the same, the correct value of "c" is greater than 5, based on data to n = 10^7. - Bill McEachen, Feb 17 2024

			214 = 2 * 107 -> Sum of factors is 109 -> 109 is prime.

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A046364, A046365, A100118, A000040, A002808, A066038, A001414.

f:=func; [k:k in [2..350]| not IsPrime(k) and IsPrime(f(k))]; // Marius A. Burtea, Dec 17 2019

ifac := proc (n) local L, x: L := ifactors(n)[2]: map(proc (x) options operator, arrow: seq(x[1], j = 1 .. x[2]) end proc, L) end proc: a := proc (n) if isprime(n) = false and isprime(add(t, t = ifac(n))) = true then n else end if end proc: seq(a(n), n = 1 .. 350); # with help from W. Edwin Clark - Emeric Deutsch, Jan 21 2009

PrimeFactorsAdded[n_] := Plus @@ Flatten[Table[ #[[1]]*#[[2]], {1}] & /@ FactorInteger[n]]; GenerateA046363[n_] := Select[Range[n], PrimeQ[PrimeFactorsAdded[ # ]] && PrimeQ[ # ] == False &]; (* GenerateA046363[100] would give all elements of this sequence below 100. - Ryan Witko (witko(AT)nyu.edu), Mar 08 2004 *)
Select[Range[325], !PrimeQ[#] && PrimeQ[Total[Times@@@FactorInteger[#]]]&] (* Jayanta Basu, May 29 2013 *)

is(n)=if(isprime(n),return(0)); my(f=factor(n)); isprime(sum(i=1,#f~,f[i,1]*f[i,2])) \\ Charles R Greathouse IV, Sep 21 2013

A100118 INTERSECT A002808. - R. J. Mathar, Sep 09 2015

Edited by R. J. Mathar, Nov 02 2009

A046316 Numbers of the form pqr where p,q,r are (not necessarily distinct) odd primes.

Original entry on oeis.org

27, 45, 63, 75, 99, 105, 117, 125, 147, 153, 165, 171, 175, 195, 207, 231, 245, 255, 261, 273, 275, 279, 285, 325, 333, 343, 345, 357, 363, 369, 385, 387, 399, 423, 425, 429, 435, 455, 465, 475, 477, 483, 507, 531, 539, 549, 555, 561, 575, 595, 603, 605

Offset: 1

Patrick De Geest, Jun 15 1998

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

A369979 sorted into ascending order.
Subsequence of A014612 and of A046340.
Cf. A255646 (final digits), A369054, A369058 (characteristic function), A369252 [= A003415(a(n))].
Subsequences: A046389, A046373, A046405, A075814, A338469, A338556, A338557, A369246.

a046316 n = a046316_list !! (n-1)
a046316_list = filter ((== 3) . a001222) [1, 3 ..]
-- Reinhard Zumkeller, May 05 2015

list(lim)=my(v=List(),pq); forprime(p=3,lim\9, forprime(q=3,min(lim\3\p,p), pq=p*q; forprime(r=3,lim\pq, listput(v, pq*r)))); Set(v) \\ Charles R Greathouse IV, Aug 23 2017

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

Definition clarified by N. J. A. Sloane, Dec 19 2017

A023023 Number of partitions of n into 3 unordered relatively prime parts.

Original entry on oeis.org

1, 1, 2, 2, 4, 4, 6, 6, 10, 8, 14, 12, 16, 16, 24, 18, 30, 24, 32, 30, 44, 32, 50, 42, 54, 48, 70, 48, 80, 64, 80, 72, 96, 72, 114, 90, 112, 96, 140, 96, 154, 120, 144, 132, 184, 128, 196, 150, 192, 168, 234, 162, 240, 192, 240, 210, 290, 192, 310, 240, 288, 256, 336, 240, 374

Offset: 3

David W. Wilson

nonn

			From _Gus Wiseman_, Oct 08 2020: (Start)
The a(3) = 1 through a(13) = 14 triples (A = 10, B = 11):
  111   211   221   321   322   332   432   433   443   543   544
              311   411   331   431   441   532   533   552   553
                          421   521   522   541   542   651   643
                          511   611   531   631   551   732   652
                                      621   721   632   741   661
                                      711   811   641   831   733
                                                  722   921   742
                                                  731   A11   751
                                                  821         832
                                                  911         841
                                                              922
                                                              931
                                                              A21
                                                              B11
(End)

Fausto A. C. Cariboni, Table of n, a(n) for n = 3..10000
Mohamed El Bachraoui, Relatively Prime Partitions with Two and Three Parts, Fibonacci Quart. 46/47 (2008/2009), no. 4, 341-345.

Cf. A023024-A023030, A000742-A000743, A023032-A023035.
A000741 is the ordered version.
A000837 counts these partitions of any length.
A001399(n-3) does not require relative primality.
A023022 is the 2-part version.
A101271 is the strict case.
A284825 counts the case that is also pairwise non-coprime.
A289509 intersected with A014612 gives the Heinz numbers.
A307719 is the pairwise coprime instead of relatively prime version.
A337599 is the pairwise non-coprime instead of relative prime version.
A008284 counts partitions by sum and length.
A078374 counts relatively prime strict partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.
Cf. A000010, A000217, A007434, A055684, A078374, A200976, A220377, A302698, A327516, A337563, A337600, A337605.

Table[Length[Select[IntegerPartitions[n,{3}],GCD@@#==1&]],{n,3,50}] (* Gus Wiseman, Oct 08 2020 *)

G.f. for the number of partitions of n into m unordered relatively prime parts is Sum(moebius(k)*x^(m*k)/Product(1-x^(i*k), i=1..m), k=1..infinity). - Vladeta Jovovic, Dec 21 2004
a(n) = (n^2/12)*Product_{prime p|n} (1 - 1/p^2) = A007434(n)/12 for n > 3 (proved by Mohamed El Bachraoui). [Jonathan Sondow, May 27 2009]
a(n) = Sum_{k=1..floor(n/3)} Sum_{i=k..floor((n-k)/2)} floor(1/gcd(i,k,n-i-k)). - Wesley Ivan Hurt, Jan 02 2021

A078843 Where 3^n occurs in n-almost primes, starting at a(0)=1.

Original entry on oeis.org

1, 2, 3, 5, 8, 14, 23, 39, 64, 103, 169, 269, 427, 676, 1065, 1669, 2628, 4104, 6414, 10023, 15608, 24281, 37733, 58503, 90616, 140187, 216625, 334527, 516126, 795632, 1225641, 1886570, 2901796, 4460359, 6851532, 10518476, 16138642, 24748319

Offset: 0

Benoit Cloitre and Paul D. Hanna, Dec 10 2002

nonn

			a(3) = 5 since 3^3 is the 5th 3-almost-prime: 8,12,18,20,27,....., A014612.

Max Alekseyev, Table of n, a(n) for n = 0..50
Eric Weisstein's World of Mathematics, Almost Prime.

Cf. A078840, A078841, A078842, A078844, A078845, A078846.

AlmostPrimePi[k_Integer /; k > 1, n_] := Module[{a, i}, a[0] = 1; Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
Table[ AlmostPrimePi[n, 3^n], {n, 2, 37}] (* Robert G. Wilson v, Feb 09 2006 *)

a(n)=sum(i=1,3^n,if(bigomega(i)-n,0,1))

{ appi(k,n,m=2) = local(r=0);
if(k==0,return(1));
if(k==1,return(primepi(n)));
forprime(p=m, floor(sqrtn(n,k)+1e-20),
r+=appi(k-1,n\p,p)-(k==2)*(primepi(p)-1));
r }
{ appi3(k,n) = appi(k,n) - if(k>=1,appi(k-1,n\3)) }
a=1; for(n=1,50, k=ceil(n*log(5/3)/log(5/2)); a+=appi3(n-k,3^n\2^k); print1(a,", "))
\\ Max Alekseyev, Jan 06 2008

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A078843(n):
    def almostprimepi(n,k):
        def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
        return int(sum(primepi(n//prod(c[1] for c in a))-a[-1][0] for a in g(n,0,1,1,k)) if k>1 else primepi(n))
    return almostprimepi(3**n,n) if n else 1 # Chai Wah Wu, Sep 01 2024

a(n) = a(n-1) + appi3(n-k, floor(3^n/2^k)), where k = ceiling(n*c) with c = log(5/3)/log(5/2) = 0.55749295065024006729857073190835923443... and appi3(k,n) is the number of k-almost primes not divisible by 3 and not exceeding n. - Max Alekseyev, Jan 06 2008

a(14)-a(37) from Robert G. Wilson v, Feb 09 2006

cp's OEIS Frontend

A069278 17-almost primes (generalization of semiprimes).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A069280 19-almost primes (generalization of semiprimes).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Python

Formula

A062069 a(n) = sigma(d(n)), where d(k) is the number of divisors function (A000005) and sigma(k) is the sum of divisors function (A000203).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

Extensions

A101605 a(n) = 1 if n is a product of exactly 3 (not necessarily distinct) primes, otherwise 0.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A211540 Number of ordered triples (w,x,y) with all terms in {1..n} and 2w = 3x + 4y.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Formula

A284825 Number of partitions of n into 3 parts without common divisors such that every pair of them has common divisors.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A046363 Composite numbers whose sum of prime factors (with multiplicity) is prime.

A046316 Numbers of the form pqr where p,q,r are (not necessarily distinct) odd primes.