A112022 -id:A112022 - cp's OEIS frontend

A000586 Number of partitions of n into distinct primes.

1, 0, 1, 1, 0, 2, 0, 2, 1, 1, 2, 1, 2, 2, 2, 2, 3, 2, 4, 3, 4, 4, 4, 5, 5, 5, 6, 5, 6, 7, 6, 9, 7, 9, 9, 9, 11, 11, 11, 13, 12, 14, 15, 15, 17, 16, 18, 19, 20, 21, 23, 22, 25, 26, 27, 30, 29, 32, 32, 35, 37, 39, 40, 42, 44, 45, 50, 50, 53, 55, 57, 61, 64, 67, 70, 71, 76, 78, 83, 87, 89, 93, 96

Offset: 0

N. J. A. Sloane

			n=16 has a(16) = 3 partitions into distinct prime parts: 16 = 2+3+11 = 3+13 = 5+11.

H. Gupta, Certain averages connected with partitions. Res. Bull. Panjab Univ. no. 124 1957 427-430.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence in two entries, N0004 and N0039).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from T. D. Noe)
Edray Herber Goins and Talitha M. Washington, On the generalized climbing stairs problem, Ars Combin. 117 (2014), 183-190. MR3243840 (Reviewed), arXiv:0909.5459 [math.CO], 2009.
H. Gupta, Partitions into distinct primes, Proc. Nat. Acad. Sci. India, 21 (1955), 185-187. [broken link]
BongJu Kim, Partition number identities which are true for all set of parts, arXiv:1803.08095 [math.CO], 2018.
Vaclav Kotesovec, Plot of log(a(n)) / log(Qas(n)) for n = 2..10^8, for Qas see the formula (25) from the article by Murthy, Brack and Bhaduri, p. 7.
M. V. N. Murthy, M. Brack, R. K. Bhaduri, On the asymptotic distinct prime partitions of integers, arXiv:1904.02776 [math.NT], Mar 22 2019.
K. F. Roth and G. Szekeres, Some asymptotic formulae in the theory of partitions, Quart. J. Math., Oxford Ser. (2) 5 (1954), 241-259.

Cf. A000041, A070215, A000607 (parts may repeat), A112022, A000009, A046675, A319264, A319267.

a000586 = p a000040_list where
   p _  0 = 1
   p (k:ks) m = if m < k then 0 else p ks (m - k) + p ks m
-- Reinhard Zumkeller, Aug 05 2012

b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
      b(n, i-1)+`if`(ithprime(i)>n, 0, b(n-ithprime(i), i-1))))
    end:
a:= n-> b(n, numtheory[pi](n)):
seq(a(n), n=0..100);  # Alois P. Heinz, Nov 15 2012

CoefficientList[Series[Product[(1+x^Prime[k]), {k, 24}], {x, 0, Prime[24]}], x]
b[n_, i_] := b[n, i] = If[n == 0, 1, If[i < 1, 0, b[n, i-1] + If[Prime[i] > n, 0, b[n - Prime[i], i-1]]]]; a[n_] := b[n, PrimePi[n]]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Apr 09 2014, after Alois P. Heinz *)
nmax = 100; pmax = PrimePi[nmax]; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 0; poly[[3]] = 1; Do[p = Prime[k]; Do[poly[[j]] += poly[[j - p]], {j, nmax + 1, p + 1, -1}];, {k, 2, pmax}]; poly (* Vaclav Kotesovec, Sep 20 2018 *)

a(n,k=n)=if(n<1, !n, my(s);forprime(p=2,k,s+=a(n-p,p-1));s) \\ Charles R Greathouse IV, Nov 20 2012

from sympy import isprime, primerange
from functools import cache
@cache
def a(n, k=None):
    if k == None: k = n
    if n < 1: return int(n == 0)
    return sum(a(n-p, p-1) for p in primerange(1, k+1))
print([a(n) for n in range(83)]) # Michael S. Branicky, Sep 03 2021 after Charles R Greathouse IV

G.f.: Product_{k>=1} (1+x^prime(k)).
a(n) = A184171(n) + A184172(n). - R. J. Mathar, Jan 10 2011
a(n) = Sum_{k=0..A024936(n)} A219180(n,k). - Alois P. Heinz, Nov 13 2012
log(a(n)) ~ Pi*sqrt(2*n/(3*log(n))) [Roth and Szekeres, 1954]. - Vaclav Kotesovec, Sep 13 2018

Entry revised by N. J. A. Sloane, Jun 10 2012

A109611 Chen primes: primes p such that p + 2 is either a prime or a semiprime.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 47, 53, 59, 67, 71, 83, 89, 101, 107, 109, 113, 127, 131, 137, 139, 149, 157, 167, 179, 181, 191, 197, 199, 211, 227, 233, 239, 251, 257, 263, 269, 281, 293, 307, 311, 317, 337, 347, 353, 359, 379, 389, 401, 409

Offset: 1

Paul Muljadi, Jul 31 2005

nonn

43 is the first prime which is not a member (see A102540).
Contains A001359 = lesser of twin primes.
A063637 is a subsequence. - Reinhard Zumkeller, Mar 22 2010
In 1966 Chen proved that this sequence is infinite; his proof did not appear until 1973 due to the Cultural Revolution. - Charles R Greathouse IV, Jul 12 2016
Primes p such that p + 2 is a term of A037143. - Flávio V. Fernandes, May 08 2021
Named after the Chinese mathematician Chen Jingrun (1933-1996). - Amiram Eldar, Jun 10 2021

			a(4) = 7 because 7 + 2 = 9 and 9 is a semiprime.
a(5) = 11 because 11 + 2 = 13, a prime.

R. J. Mathar, Table of n, a(n) for n = 1..34076
Jing Run Chen, On the representation of a larger even integer as the sum of a prime and the product of at most two primes, Sci. Sinica 16 (1973), pp. 157-176.
Ben Green and Terence Tao, Restriction theory of the Selberg sieve, with applications, arXiv:math/0405581 [math.NT], 2004-2005, pp. 5, 14, 18-19, 21.
Ben Green and Terence Tao, Restriction theory of the Selberg sieve, with applications, J. Théor. Nombres Bordeaux, Vol. 18, No. 1 (2006), pp. 147-182.
Eric Weisstein's World of Mathematics, Chen's Theorem.
Eric Weisstein's World of Mathematics, Chen Prime.
Wikipedia, Chen prime.
Binbin Zhou, The Chen primes contain arbitrarily long arithmetic progressions, Acta Arithmetica, Vol. 138 (2009), pp. 301-315.

Union of A001359 and A063637.
Cf. A001358, A112021, A112022, A139689, A269256.
Cf. A037143.

A109611 := proc(n)
    option remember;
    if n =1 then
        2;
    else
        a := nextprime(procname(n-1)) ;
        while true do
            if isprime(a+2) or numtheory[bigomega](a+2) = 2 then
                return a;
            end if;
            a := nextprime(a) ;
        end do:
    end if;
end proc: # R. J. Mathar, Apr 26 2013

semiPrimeQ[x_] := TrueQ[Plus @@ Last /@ FactorInteger[ x ] == 2]; Select[Prime[Range[100]], PrimeQ[ # + 2] || semiPrimeQ[ # + 2] &] (* Alonso del Arte, Aug 08 2005 *)
SequencePosition[PrimeOmega[Range[500]], {1, , 1|2}][[All, 1]] (* _Jean-François Alcover, Feb 10 2018 *)

isA001358(n)= if( bigomega(n)==2, return(1), return(0) );
isA109611(n)={ if( ! isprime(n), return(0), if( isprime(n+2), return(1), return( isA001358(n+2)) ); ); }
{ n=1; for(i=1,90000, p=prime(i); if( isA109611(p), print(n," ",p); n++; ); ); } \\ R. J. Mathar, Aug 20 2006

list(lim)=my(v=List([2]),semi=List(),L=lim+2,p=3); forprime(q=3,L\3, forprime(r=3,min(L\q,q), listput(semi,q*r))); semi=Set(semi); forprime(q=7,lim, if(setsearch(semi,q+2), listput(v,q))); forprime(q=5,L, if(q-p==2, listput(v,p)); p=q); Set(v) \\ Charles R Greathouse IV, Aug 25 2017

from sympy import isprime, primeomega
def ok(n): return isprime(n) and (primeomega(n+2) < 3)
print(list(filter(ok, range(1, 410)))) # Michael S. Branicky, May 08 2021

a(n)+2 = A139690(n).

Sum_{n>=1} 1/a(n) converges (Zhou, 2009). - Amiram Eldar, Jun 10 2021

Corrected by Alonso del Arte, Aug 08 2005

A112020 Number of partitions of n into distinct semiprimes.

Original entry on oeis.org

1, 0, 0, 0, 1, 0, 1, 0, 0, 1, 2, 0, 0, 1, 2, 2, 1, 0, 1, 3, 2, 2, 1, 2, 3, 5, 2, 2, 3, 5, 4, 5, 3, 4, 6, 9, 6, 5, 6, 10, 10, 9, 7, 9, 12, 14, 12, 11, 14, 18, 17, 16, 16, 19, 21, 24, 21, 23, 26, 29, 30, 32, 31, 33, 39, 40, 39, 41, 45, 49, 54, 53, 54, 59, 68, 66, 68, 70, 78, 82, 88, 86, 93, 101

Offset: 0

Reinhard Zumkeller, Aug 26 2005

nonn

			For n=4 one partition: {2*2}.
For n=6 one partition: {2*3}.
For n=10 two partitions: {2*2+2*3,2*5}.

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

Cf. A101048, A112022, A001358, A000586.

h:= proc(n) option remember; `if`(n=0, 0,
     `if`(numtheory[bigomega](n)=2, n, h(n-1)))
    end:
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
      b(n-i, h(min(n-i, i-1)))+b(n, h(i-1))))
    end:
a:= n-> b(n, h(n)):
seq(a(n), n=0..100);  # Alois P. Heinz, Mar 19 2024

nmax = 100;
CoefficientList[Series[Product[1+x^(Prime[j] Prime[k]), {j, 1, nmax}, {k, j, nmax}], {x, 0, nmax}], x] (* Jean-François Alcover, Nov 10 2021 *)

A112021 Number of partitions of n into Chen primes.

Original entry on oeis.org

0, 1, 1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 9, 10, 12, 14, 17, 19, 23, 26, 30, 35, 40, 46, 52, 60, 67, 77, 87, 98, 111, 124, 140, 157, 175, 197, 219, 244, 272, 302, 336, 372, 412, 456, 503, 556, 613, 675, 742, 816, 896, 983, 1078, 1180, 1291, 1411, 1542, 1683, 1836, 2001, 2178

Offset: 1

Reinhard Zumkeller, Aug 26 2005

nonn

a(n) = A000607(n) for n <= 42.

Andrew Howroyd, Table of n, a(n) for n = 1..1000

Cf. A112022, A101048, A109611.

fQ[n_] := PrimeQ@n && (PrimeQ[n + 2] || 2 == Plus @@ Last /@ FactorInteger[n + 2]); f[n_] := Block[{c = k = 0, l = PartitionsP@n, p = Union /@ IntegerPartitions@n}, While[k++; k < l, If[Union[fQ /@ p[[k]]] == {True}, c++ ]]; c]; lst = {}; Do[ AppendTo[lst, f[n]], {n, 61}]; lst (* or *)
Rest@ CoefficientList[ Series[1/Times @@ (1 - x^Select[ Range@100, fQ@# &]), {x, 0, 61}], x] (* Robert G. Wilson v, Jun 16 2006 *)

ok(n)={isprime(n) && bigomega(n+2)<3}
{my(n=80); Vec(prod(k=1, n, if(ok(k), 1/(1-x^k) + O(x*x^n), 1))-1,-n)} \\ Andrew Howroyd, Dec 28 2017

G.f.: Product_{k>=1} 1/(1 - x^A109611(k)). - Andrew Howroyd, Dec 28 2017

cp's OEIS Frontend

A000586 Number of partitions of n into distinct primes.

Views

Author

Keywords

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

Extensions

A109611 Chen primes: primes p such that p + 2 is either a prime or a semiprime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A112020 Number of partitions of n into distinct semiprimes.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A112021 Number of partitions of n into Chen primes.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula