A045917 -id:A045917 - cp's OEIS frontend

A061357 Number of 0

0, 0, 0, 1, 1, 1, 1, 2, 2, 2, 2, 3, 2, 2, 3, 2, 3, 4, 1, 3, 4, 3, 3, 5, 4, 3, 5, 3, 3, 6, 2, 5, 6, 2, 5, 6, 4, 5, 7, 4, 4, 8, 4, 4, 9, 4, 4, 7, 3, 6, 8, 5, 5, 8, 6, 7, 10, 6, 5, 12, 3, 5, 10, 3, 7, 9, 5, 5, 8, 7, 7, 11, 5, 5, 12, 4, 8, 11, 4, 8, 10, 5, 5, 13, 9, 6, 11, 7, 6, 14, 6, 8, 13, 5, 8, 11, 6, 9

Offset: 1

Amarnath Murthy, Apr 28 2001

Number of prime pairs (p,q) with p < n < q and q-n = n-p.
The same as the number of ways n can be expressed as the mean of two distinct primes.
Conjecture: for n>=4 a(n)>0. - Benoit Cloitre, Apr 29 2003
Conjectures from Rick L. Shepherd, Jun 24 2003: (Start)
1) For each integer N>=1 there exists a positive integer m(N) such that for n>=m(N) a(n)>a(N). (After the first m(N)-1 terms, a(N) does not reappear). In particular, for N=1 (or 2 or 3), m(N)=4 and a(N)=0, giving Benoit Cloitre's conjecture. (cont.)
(cont.) Conjectures based upon observing a(1),...,a(10000):
m(4)=m(5)=m(6)=m(7)=m(19)=20 for a(4)=a(5)=a(6)=a(7)=a(19)=1,
m(8)=...(7 others)...=m(34)=35 for a(8)=...(7 others)...=a(34)=2,
m(12)=...(10 others)...=m(64)=65 for a(12)=...(10 others)...=a(64)=3,
m(18)=...(10 others)...=m(79)=80 for a(18)=...(10 others)...=a(79)=4,
m(24)=...(14 others)...=m(94)=95 for a(24)=...(14 others)...=a(94)=5,
m(30)=...(17 others)...=m(199)=200 for a(30)=...(17 others)...=a(199)=6, etc.
2) Each nonnegative integer appears at least once in the current sequence.
3) Stronger than 2): A001477 (nonnegative integers) is a subsequence of the current sequence. (Supporting evidence: I've observed that 0,1,2,...,175 is a subsequence of a(1),...,a(10000)).
(End)
a(n) is also the number of k such that 2*k+1=p and 2*(n-k-1)+1=q are both odd primes with p < q with p*q = n^2 - m^2. [Pierre CAMI, Sep 01 2008]
Also: Number of ways n^2 can be written as b^2+pq where 0a(n) = sum (A010051(2*n - p): p prime < n). [Reinhard Zumkeller, Oct 19 2011]
a(n) is also the number of partitions of 2*n into two distinct primes. See the first formula by T. D. Noe, and the Alois P. Heinz, Nov 14 2012, crossreference. - Wolfdieter Lang, May 13 2016
All 0Jamie Morken, Jun 02 2017
a(n) is the number of appearances of n in A143836. - Ya-Ping Lu, Mar 05 2023

			a(10)= 2: there are two such pairs (3,17) and (7,13), as 10 = (3+17)/2 = (7+13)/2.

Pierre CAMI, Table of n, a(n) for n = 1..60000

Cf. A071681 (subsequence for prime n only).
Cf. A092953.
Bisection of A117929 (even part). - Alois P. Heinz, Nov 14 2012

a061357 n = sum $
   zipWith (\u v -> a010051 u * a010051 v) [n+1..] $ reverse [1..n-1]
-- Reinhard Zumkeller, Nov 10 2012, Oct 19 2011

Table[Count[Range[n - 1], k_ /; And[PrimeQ[n - k], PrimeQ[n + k]]], {n, 98}] (* Michael De Vlieger, May 14 2016 *)

a(n)=my(s);forprime(p=2,n-1,s+=isprime(2*n-p));s \\ Charles R Greathouse IV, Mar 08 2013

from sympy import primerange, isprime
def A061357(n): return sum(1 for p in primerange(n) if isprime((n<<1)-p)) # Chai Wah Wu, Sep 03 2024

a(n) = A045917(n) - A010051(n). - T. D. Noe, May 08 2007
a(n) = sum(A010051(n-k)*A010051(n+k): 1 <= k < n). - Reinhard Zumkeller, Nov 10 2012
a(n) = sum_{i=2..n-1} A010051(i)*A010051(2n-i). [Wesley Ivan Hurt, Aug 18 2013]

More terms from Larry Reeves (larryr(AT)acm.org), May 15 2001

A035026 Number of times that i and 2n-i are both prime, for i = 1, ..., 2n-1.

Original entry on oeis.org

0, 1, 1, 2, 3, 2, 3, 4, 4, 4, 5, 6, 5, 4, 6, 4, 7, 8, 3, 6, 8, 6, 7, 10, 8, 6, 10, 6, 7, 12, 5, 10, 12, 4, 10, 12, 9, 10, 14, 8, 9, 16, 9, 8, 18, 8, 9, 14, 6, 12, 16, 10, 11, 16, 12, 14, 20, 12, 11, 24, 7, 10, 20, 6, 14, 18, 11, 10, 16, 14, 15, 22, 11, 10, 24, 8, 16, 22, 9, 16, 20, 10

Offset: 1

Gordon R. Bower (siegmund(AT)mosquitonet.com)

a(n) is the convolution of terms 1 to 2n of the characteristic function of the primes, A010051, with itself. Related to Goldbach's conjecture that every even number can be expressed as the sum of two primes. - T. D. Noe, Aug 01 2002
The following sequences all appear to have the same parity (with an extra zero term at the start of A010051): A010051, A061007, A035026, A069754, A071574. - Jeremy Gardiner, Aug 09 2002
Total number of printer jobs in all possible schedules for n time slots in the first-come-first-served (FCFS) policy.
a(n) = Sum_{p prime < 2*n} A010051(2*n - p). - Reinhard Zumkeller, Oct 19 2011
For n > 1: length of n-th row of triangle A171637. - Reinhard Zumkeller, Mar 03 2014
a(n) = A001221(A238711(n)) = A238778(n) / n. - Reinhard Zumkeller, Mar 06 2014
From Robert G. Wilson v, Dec 15 2016: (Start)
First occurrence of k: 1, 2, 4, 5, 8, 11, 12, 17, 18, 37, 24, 53, 30, 89, 39, 71, 42, 101, 45, 179, 57, 137, 72, 193, 60, 233, ..., .
Conjectured last occurrence of k: 1, 3, 6, 19, 34, 31, 64, 61, 76, 79, 94, 83, 166, 199, 136, 181, 184, 229, 244, 271, 316, 277, 346, 313, 301, 293, ..., .
Conjectured number occurrences of k: 1, 2, 2, 3, 6, 3, 8, 4, 7, 5, 11, 5, 11, 8, 10, 3, 17, 7, 16, 3, 13, 8, 21, 4, 12, 3, 22, 7, 20, 8, 15, ..., .
Records: 0, 1, 2, 3, 4, 5, 6, 7, 8, 10, 12, 14, 16, 18, 20, 24, 26, 28, 38, 42, 48, 54, 60, 64, 82, 88, 102, 104, 114, 116, 136, 146, 152, 166, 182, ..., .
(End)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Rémy Sigrist, Colored scatterplot of the first 500000 terms
Index entries for sequences related to Goldbach conjecture

Cf. A010051. Essentially the same as A002372.

Cf. A073610.

a035026 n = sum $ map (a010051 . (2 * n -)) $
   takeWhile (< 2 * n) a000040_list
-- Reinhard Zumkeller, Oct 19 2011

A035026 := proc(n)
    local a,i ;
    a := 0 ;
    for i from 1 to 2*n-1 do
        if isprime(i) and isprime(2*n-i) then
            a := a+1 ;
        end if;
    end do:
    a ;
end proc: # R. J. Mathar, Jul 01 2013

For[lst={}; n=1, n<=100, n++, For[cnt=0; i=1, i<=2n-1, i++ If[PrimeQ[i]&&PrimeQ[2n-i], cnt++ ]]; AppendTo[lst, cnt]]; lst
f[n_] := Block[{c = Boole@ PrimeQ[ n/2], p = 2}, While[ 2p < n, If[ PrimeQ[n - p], c += 2]; p = NextPrime@ p]; c];; Array[ f[ 2#] &, 90] (* Robert G. Wilson v, Dec 15 2016 *)

For n > 1, a(n) = 2*A045917(n) - A010051(n).
a(n) = A010051(n) + 2*A061357(n). - Wesley Ivan Hurt, Aug 21 2013
a(n) = A073610(2*n). - Ridouane Oudra, Sep 06 2023

Corrected by T. D. Noe, May 05 2002

A047160 For n >= 2, a(n) = smallest number m >= 0 such that n-m and n+m are both primes, or -1 if no such m exists.

Original entry on oeis.org

0, 0, 1, 0, 1, 0, 3, 2, 3, 0, 1, 0, 3, 2, 3, 0, 1, 0, 3, 2, 9, 0, 5, 6, 3, 4, 9, 0, 1, 0, 9, 4, 3, 6, 5, 0, 9, 2, 3, 0, 1, 0, 3, 2, 15, 0, 5, 12, 3, 8, 9, 0, 7, 12, 3, 4, 15, 0, 1, 0, 9, 4, 3, 6, 5, 0, 15, 2, 3, 0, 1, 0, 15, 4, 3, 6, 5, 0, 9, 2, 15, 0, 5, 12, 3, 14, 9, 0, 7, 12, 9, 4, 15, 6, 7, 0, 9, 2, 3

Offset: 2

Lior Manor

I have confirmed there are no -1 entries through integers to 4.29*10^9 using PARI. - Bill McEachen, Jul 07 2008
From Daniel Forgues, Jul 02 2009: (Start)
Goldbach's Conjecture: for all n >= 2, there are primes (distinct or not) p and q s.t. p+q = 2n. The primes p and q must be equidistant (distance m >= 0) from n: p = n-m and q = n+m, hence p+q = (n-m)+(n+m) = 2n.
Equivalent to Goldbach's Conjecture: for all n >= 2, there are primes p and q equidistant (distance >= 0) from n, where p and q are n when n is prime.
If this conjecture is true, then a(n) will never be set to -1.
Twin Primes Conjecture: there is an infinity of twin primes.
If this conjecture is true, then a(n) will be 1 infinitely often (for which each twin primes pair is (n-1, n+1)).
Since there is an infinity of primes, a(n) = 0 infinitely often (for which n is prime).
(End)
If n is composite, then n and a(n) are coprime, because otherwise n + a(n) would be composite. - Jason Kimberley, Sep 03 2011
From Jianglin Luo, Sep 22 2023: (Start)
a(n) < primepi(n)+sigma(n,0);
a(n) < primepi(primepi(n)+n);
a(n) < primepi(n), for n>344;
a(n) = o(primepi(n)), as n->+oo. (End)
If -1 < a(n) < n-3, then a(n) is divisible by 3 if and only if n is not divisible by 3, and odd if and only if n is even. - Robert Israel, Oct 05 2023

			16-3=13 and 16+3=19 are primes, so a(16)=3.

T. D. Noe, Table of n, a(n) for n = 2..10000
Jason Kimberley, Symmetrical plot of A047160

Cf. A001031, A002092, A002372, A002373, A002374, A002375, A014092, A025583, A035026, A047949, A071406, A082467, A102084, A103147, A112823, A155764, A155765, A177461, A078611, A010051, A045917, A325142.

a047160 n = if null ms then -1 else head ms
            where ms = [m | m <- [0 .. n - 1],
                            a010051' (n - m) == 1, a010051' (n + m) == 1]
-- Reinhard Zumkeller, Aug 10 2014

A047160:=func;[A047160(n):n in[2..100]]; // Jason Kimberley, Sep 02 2011

Table[k = 0; While[k < n && (! PrimeQ[n - k] || ! PrimeQ[n + k]), k++]; If[k == n, -1, k], {n, 2, 100}]
smm[n_]:=Module[{m=0},While[AnyTrue[n+{m,-m},CompositeQ],m++];m]; Array[smm,100,2] (* Harvey P. Dale, Nov 16 2024 *)

a(n)=forprime(p=n,2*n, if(isprime(2*n-p), return(p-n))); -1 \\ Charles R Greathouse IV, Jun 23 2017

10 N=2// 20 M=0// 30 if and{prmdiv(N-M)=N-M,prmdiv(N+M)=N+M} then print M;:goto 50// 40 inc M:goto 30// 50 inc N: if N>130 then stop// 60 goto 20

a(n) = n - A112823(n).

a(n) = A082467(n) * A005171(n), for n > 3. - Jason Kimberley, Jun 25 2012

More terms from Patrick De Geest, May 15 1999
Deleted a comment. - T. D. Noe, Jan 22 2009
Comment corrected and definition edited by Daniel Forgues, Jul 08 2009

A117929 Number of partitions of n into 2 distinct primes.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 2, 0, 2, 1, 2, 1, 2, 0, 3, 1, 2, 0, 2, 0, 3, 1, 2, 1, 3, 0, 4, 0, 1, 1, 3, 0, 4, 1, 3, 1, 3, 0, 5, 1, 4, 0, 3, 0, 5, 1, 3, 0, 3, 0, 6, 1, 2, 1, 5, 0, 6, 0, 2, 1, 5, 0, 6, 1, 4, 1, 5, 0, 7, 0, 4, 1, 4, 0, 8, 1, 4, 0, 4, 0, 9, 1, 4, 0, 4, 0, 7, 0, 3, 1, 6, 0, 8, 1, 5, 1

Offset: 1

Emeric Deutsch, Apr 03 2006

Number of distinct rectangles with prime length and width such that L + W = n, W < L. For example, a(16) = 2; the two rectangles are 3 X 13 and 5 X 11. - Wesley Ivan Hurt, Oct 29 2017

			a(24) = 3 because we have [19,5], [17,7] and [13,11].

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

Cf. A010051, A045917, A061358, A073610, A166081 (positions of 0), A077914 (positions of 2), A080862 (positions of 6).

Column k=2 of A219180. - Alois P. Heinz, Nov 13 2012

g:=sum(sum(x^(ithprime(i)+ithprime(j)),i=1..j-1),j=1..35): gser:=series(g,x=0,130): seq(coeff(gser,x,n),n=1..125);
# alternative
A117929 := proc(n)
    local a,i,p ;
    a := 0 ;
    p := 2 ;
    for i from 1 do
        if 2*p >= n then
            return a;
        end if;
        if isprime(n-p) then
            a := a+1 ;
        end if;
        p := nextprime(p) ;
    end do:
end proc:
seq(A117929(n),n=1..80) ; # R. J. Mathar, Oct 01 2021

l = {}; For[n = 1, n <= 1000, n++, c = 0; For[k = 1, Prime[k] < n/2, k++, If[PrimeQ[n - Prime[k]], c = c + 1] ]; AppendTo[l, c] ] l (* Jake Foster, Oct 27 2008 *)
Table[Count[IntegerPartitions[n,{2}],?(AllTrue[#,PrimeQ]&&#[[1]]!= #[[2]] &)],{n,120}] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale, Jul 26 2020 *)

a(n)=my(s);forprime(p=2,(n-1)\2,s+=isprime(n-p));s \\ Charles R Greathouse IV, Feb 26 2014

from sympy import sieve
from collections import Counter
from itertools import combinations
def aupton(max):
    sieve.extend(max)
    a = Counter(c[0]+c[1] for c in combinations(sieve._list, 2))
    return [a[n] for n in range(1, max+1)]
print(aupton(105)) # Michael S. Branicky, Feb 16 2024

G.f.: Sum_{j>0} Sum_{i=1..j-1} x^(p(i)+p(j)), where p(k) is the k-th prime.
G.f.: A(x)^2/2 - A(x^2)/2 where A(x) = Sum_{p in primes} x^p. - Geoffrey Critzer, Nov 21 2012
a(n) = [x^n*y^2] Product_{i>=1} (1+x^prime(i)*y). - Alois P. Heinz, Nov 22 2012
a(n) = Sum_{i=2..floor((n-1)/2)} A010051(i) * A010051(n-i). - Wesley Ivan Hurt, Oct 29 2017

A067187 Numbers that can be expressed as the sum of two primes in exactly one way.

Original entry on oeis.org

4, 5, 6, 7, 8, 9, 12, 13, 15, 19, 21, 25, 31, 33, 39, 43, 45, 49, 55, 61, 63, 69, 73, 75, 81, 85, 91, 99, 103, 105, 109, 111, 115, 129, 133, 139, 141, 151, 153, 159, 165, 169, 175, 181, 183, 193, 195, 199, 201, 213, 225, 229, 231, 235, 241, 243, 253, 259, 265, 271

Offset: 1

Amarnath Murthy, Jan 10 2002

nonn

All primes + 2 are terms of this sequence. Is 12 the last even term? - Frank Ellermann, Jan 17 2002
A048974, A052147, A067187 and A088685 are very similar after dropping terms less than 13. - Eric W. Weisstein, Oct 10 2003
Values of n such that A061358(n)=1. - Emeric Deutsch, Apr 03 2006

			4 is a term as 4 = 2+2, 15 is a term as 15 = 13+2.

Robert Price, Table of n, a(n) for n = 1..5002
Index entries for sequences related to Goldbach conjecture

Cf. A023036, A045917, A061358.
Subsequence of A014091.
Numbers that can be expressed as the sum of two primes in k ways for k=0..10: A014092 (k=0), this sequence (k=1), A067188 (k=2), A067189 (k=3), A067190 (k=4), A067191 (k=5), A066722 (k=6), A352229 (k=7), A352230 (k=8), A352231 (k=9), A352233 (k=10).

g:=sum(sum(x^(ithprime(i)+ithprime(j)),i=1..j),j=1..80): gser:=series(g,x=0,280): a:=proc(n) if coeff(gser,x^n)=1 then n else fi end: seq(a(n),n=1..272); # Emeric Deutsch, Apr 03 2006

cQ[n_]:=Module[{c=0},Do[If[PrimeQ[n-i]&&PrimeQ[i],c++],{i,2,n/2}]; c==1]; Select[Range[4,271],cQ[#]&] (* Jayanta Basu, May 22 2013 *)
y = Select[Flatten@Table[Prime[i] + Prime[j], {i, 60}, {j, 1, i}], # < Prime[60] &]; Select[Union[y], Count[y, #] == 1 &] (* Robert Price, Apr 21 2025 *)

Edited by Frank Ellermann, Jan 17 2002

A000954 Conjecturally largest even integer which is an unordered sum of two primes in exactly n ways.

Original entry on oeis.org

2, 12, 68, 128, 152, 188, 332, 398, 368, 488, 632, 692, 626, 992, 878, 908, 1112, 998, 1412, 1202, 1448, 1718, 1532, 1604, 1682, 2048, 2252, 2078, 2672, 2642, 2456, 2936, 2504, 2588, 2978, 3092, 3032, 3218, 3272, 3296, 3632, 3548, 3754, 4022, 4058, 4412

Offset: 0

Bill Gosper

The Goldbach conjecture is that every even number is the sum of two primes.

			2 is largest even integer which is the sum of two primes in 0 ways, 12 is largest even integer which is the unordered sum of two primes in 1 way (5+7), etc.

T. D. Noe, Table of n, a(n) for n = 0..5000

Cf. A045917, A023036, A000974, A001172, A002375.

f[n_] := Block[{c = 0, k = 3}, While[k <= n/2, If[PrimeQ[k] && PrimeQ[n - k], c++ ]; k++ ]; c]; a = Table[0, {50}]; a[[1]] = 2; a[[2]] = 4; Do[m = n; b = f[n]; If[b < 100, a[[b + 1]] = n], {n, 6, 20000, 2}] (* Robert G. Wilson v, Dec 20 2003 *)

A067188 Numbers that can be expressed as the (unordered) sum of two primes in exactly two ways.

Original entry on oeis.org

10, 14, 16, 18, 20, 28, 32, 38, 68

Offset: 1

Amarnath Murthy, Jan 10 2002

Corresponds to numbers 2m such that A045917(m)=2. Subsequence of A014091. - Lekraj Beedassy, Apr 22 2004

			18 is a term as 18 = 13+5 = 11+7 are the only two ways to express 18 as a sum of two primes.

Index entries for sequences related to Goldbach conjecture

Cf. A023036.

Numbers that can be expressed as the sum of two primes in k ways for k=0..10: A014092 (k=0), A067187 (k=1), this sequence (k=2), A067189 (k=3), A067190 (k=4), A067191 (k=5), A066722 (k=6), A352229 (k=7), A352230 (k=8), A352231 (k=9), A352233 (k=10).

y = Select[Flatten@Table[Prime[i] + Prime[j], {i, 100}, {j, 1, i}], # < Prime[100] &]; Select[Union[y], Count[y, #] == 2 &] (* Robert Price, Apr 22 2025 *)

Corrected by Peter Bertok (peter(AT)bertok.com), who finds (Jan 13 2002) that there are no other terms below 10000 and conjectures there are no further terms in this sequence and A067189, A067190, etc.

R. K. Guy (Jan 14 2002) remarks: "I believe that these conjectures follow from a more general one by Hardy & Littlewood (probably in Some problems of 'partitio numerorum' III, on the expression of a number as a sum of primes, Acta Math. 44(1922) 1-70)."

Previous Showing 11-20 of 109 results. Next

cp's OEIS Frontend

A068307 From Goldbach problem: number of decompositions of n into a sum of three primes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

Extensions

A001031 Goldbach conjecture: a(n) = number of decompositions of 2n into sum of two primes (counting 1 as a prime).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

Extensions

A023036 Smallest positive even integer that is an unordered sum of two primes in exactly n ways.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A061357 Number of 0

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

Extensions

A035026 Number of times that i and 2n-i are both prime, for i = 1, ..., 2n-1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Formula

Extensions

A047160 For n >= 2, a(n) = smallest number m >= 0 such that n-m and n+m are both primes, or -1 if no such m exists.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Mathematica