A061358 -id:A061358 - cp's OEIS frontend

A141340 Positive integers n such that A061358(n) = #{primes p | n/2 <= p < n-1}.

1, 2, 3, 4, 5, 6, 7, 8, 10, 12, 14, 16, 18, 24, 30, 36, 42, 48, 60, 90, 210

Offset: 1

Rick L. Shepherd, Jun 25 2008

According to Brouwers et al., Deshouillers et al. showed that the maximum term of this sequence is 210. A141341 is a subsequence.

Positive integers k such that, for each prime p with k/2 <= p <= k - 2, k - p is prime. - Charles R Greathouse IV, May 28 2017

			For each prime 210/2 <= p <= 210, 210 - p is prime, and so 210 is in this sequence: 210 - 107 = 103, 210 - 109 = 101, 210 - 113 = 97, 210 - 127 = 83, 210 - 131 = 79, 210 - 137 = 73, 210 - 139 = 71, 210 - 149 = 61, 210 - 151 = 59, 210 - 157 = 53, 210 - 163 = 47, 210 - 167 = 43, 210 - 173 = 37, 210 - 179 = 31, 210 - 181 = 29, 210 - 191 = 19, 210 - 193 = 17, 210 - 197 = 13, 210 - 199 = 11. - _Charles R Greathouse IV_, May 28 2017

J-M. Deshouillers, A. Granville, W. Narkiewicz and C. Pomerance, An upper bound in Goldbach's problem, Math. Comp. 61 (1993), 209-213.
Brady Haran and Carl Pomerance, 210 is VERY Goldbachy, Numberphile video (2017)
David van Golstein Brouwers, John Bamberg and Grant Cairns, Totally Goldbach numbers and related conjectures, The Australian Mathematical Society, Gazette, Volume 31 Number 4, September 2004.

Cf. A061358, A141341.

Block[{r = {}}, Do[ If[ AllTrue[i - #, PrimeQ] &@ NextPrime[i/2, Range[ PrimePi[i - 2] - PrimePi[i/2]]], AppendTo[r, i]], {i, 210}]; r] (* Mikk Heidemaa, May 29 2024 *)

is(n)=forprime(p=n/2,n-2, if(!isprime(n-p),return(0))); 1 \\ Charles R Greathouse IV, May 28 2017; corrected by Michel Marcus, May 30 2024

A162197 a(n) = A161912(n) - A061358(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 6, 3, 3, 2, 5, 5, 7, 5, 8, 3, 6, 4, 11, 8, 11, 7, 10, 4, 8, 11, 15, 14, 13, 16, 18, 11, 19, 17, 20, 16, 19, 17, 20, 16, 14, 19, 33, 22, 25, 20, 26, 19, 30, 23, 29, 25, 26, 24, 30, 14, 27, 35, 35, 26

Offset: 0

Omar E. Pol, Jul 04 2009

Cf. A006218, A040976, A061358, A161911, A161912, A161913.

a(n) = A040976(n+1) - A006218(n) - A061358(n).

A002375 From Goldbach conjecture: number of decompositions of 2n into an unordered sum of two odd primes.

Original entry on oeis.org

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

Offset: 1

N. J. A. Sloane

A weaker form of this conjecture, the ternary form, was proved by Helfgott (see link below). - T. D. Noe, May 14 2013
The Goldbach conjecture is that for n >= 3, this sequence is always positive.
This has been checked up to at least 10^18 (see A002372).
With the exception of the n=2 term, identical to A045917.
The conjecture has been verified up to 3 * 10^17 (see MathWorld link). - Dmitry Kamenetsky, Oct 17 2008
Languasco and Zaccagnini proved that, where Lambda is the von Mangoldt function, and R(n) = Sum_{i + j = n} Lambda(i)*Lambda(j) is the counting function for the Goldbach numbers, and for N >= 2 and assume the Riemann hypothesis (RH) holds, then Sum_{n = 1..N} R(n) = (N^2)/2 - 2*Sum_{rho} ((N^(rho+1))/(rho*(rho+1))) + O(N * log^3 N).
If 2n is the sum of two distinct primes, then neither prime divides 2n. - Christopher Heiling, Feb 28 2017

			2 and 4 are not the sum of 2 odd primes, so a(1) = a(2) = 0; 6 = 3 + 3 (one way, so a(3) = 1); 8 = 3 + 5 (so a(4) = 1); 10 = 3 + 7 = 5 + 5 (so a(5) = 2); etc.

Calvin C. Clawson, "Mathematical Mysteries, the beauty and magic of numbers," Perseus Books, Cambridge, MA, 1996, Chapter 12, Pages 236-257.
Apostolos K. Doxiadis, Uncle Petros and Goldbach's Conjecture, Bloomsbury Pub. PLC USA, 2000.
D. A. Grave, Traktat z Algebrichnogo Analizu (Monograph on Algebraic Analysis). Vol. 2, p. 19. Vidavnitstvo Akademiia Nauk, Kiev, 1938.
H. Halberstam and H. E. Richert, 1974, "Sieve methods", Academic press, London, New York, San Francisco.
D. H. Lehmer, Guide to Tables in the Theory of Numbers. Bulletin No. 105, National Research Council, Washington, DC, 1941, p. 80.
N. V. Maslova, On the coincidence of Grünberg-Kegel graphs of a finite simple group and its proper subgroup, Proceedings of the Steklov Institute of Mathematics April 2015, Volume 288, Supplement 1, pp 129-141; Original Russian Text: Trudy Instituta Matematiki i Mekhaniki UrO RAN, 2014, Vol. 20, No. 1.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

H. J. Smith, Table of n, a(n) for n = 1..20000 (first 10000 terms from T. D. Noe)
J.-M. Deshouillers, H. J. J. te Riele and Y. Saouter, New Experimental Results Concerning the Goldbach Conjecture, Modelling, Analysis and Simulation [MAS], R 9804, p.1-12, Technical report, 1998.
James A. Farrugia, Brun's 1920 Theorem on Goldbach's Conjecture, Master's Thesis, Utah State University, All Graduate Theses and Dissertations (2018). 7153.
H. A. Helfgott, Minor arcs for Goldbach's problem, arXiv:1205.5252 [math.NT], 2012-2013.
H. A. Helfgott, Major arcs for Goldbach's theorem, arXiv:1305.2897 [math.NT], 2013-2014.
H. A. Helfgott, The ternary Goldbach conjecture is true, arxiv:1312.7748 [math.NT], 2013.
H. A. Helfgott, The ternary Goldbach problem, arXiv:1404.2224 [math.NT], 2014.
M. Herkommer, Goldbach Conjecture Research
A. V. Kumchev and D. I. Tolev, An invitation to additive number theory, arXiv:math/0412220 [math.NT], 2004.
Alessandro Languasco and Alessandro Zaccagnini, The number of Goldbach representations of an integer, Proc. Amer. Math. Soc. 140 (2012), 795-804 (preliminary version, arXiv:1011.3198 [math.NT], 2010).
Jörg Richstein, Verifying the Goldbach conjecture up to 4 * 10^14, Math. Comput., 70 (2001), 1745-1749.
Vladimir Shevelev, Binary additive problems: recursions for numbers of representations, arXiv:0901.3102 [math.NT], 2009-2013.
Matti K. Sinisalo, Checking the Goldbach conjecture up to 4*10^11, Math. Comp. 61 (1993), pp. 931-934.
Eric Weisstein's World of Mathematics, Goldbach Partition
Wikipedia, Goldbach's conjecture
G. Xiao, WIMS server, Goldbach
Index entries for sequences related to Goldbach conjecture

See also A061358. Cf. A002372 (ordered sums), A002373, A002374, A045917.
A023036 is (essentially) the first appearance of n and A000954 is the last (assumed) appearance of n.
Cf. A065091, A010051, A001031 (a weaker form of the conjecture).

a002375 n = sum $ map (a010051 . (2 * n -)) $ takeWhile (<= n) a065091_list
-- Reinhard Zumkeller, Sep 02 2013

A002375 := func; [A002375(n):n in[1..98]];

A002375 := proc(n) local s, p; s := 0; p := 3; while p<2*n do s := s+x^p; p := nextprime(p) od; (coeff(s^2, x, 2*n)+coeff(s,x,n))/2 end; [seq(A002375(n), n=1..100)];
a:=proc(n) local c,k; c:=0: for k from 1 to floor((n-1)/2) do if isprime(2*k+1)=true and isprime(2*n-2*k-1)=true then c:=c+1 else c:=c fi od end: A:=[0,0,seq(a(n),n=3..98)]; # Emeric Deutsch, Aug 27 2007
g:=sum(sum(x^(ithprime(i)+ithprime(j)),i=2..j),j=2..50): seq(coeff(g,x,2*n), n =1..98); # Emeric Deutsch, Aug 27 2007

f[n_] := Length[ Select[2n - Prime[ Range[2, PrimePi[n]]], PrimeQ]]; Table[ f[n], {n, 100}] (* Paul Abbott, Jan 11 2005 *)
nn = 10^2; ps = Boole[PrimeQ[Range[1,2*nn,2]]]; Table[Sum[ps[[i]] ps[[n-i+1]], {i, Ceiling[n/2]}], {n, nn}] (* T. D. Noe, Apr 13 2011 *)
Table[Count[IntegerPartitions[2n,{2}],?(AllTrue[#,PrimeQ]&&FreeQ[#,2]&)],{n,100}] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale, Mar 01 2018 *)
j[n_] := If[PrimeQ[2 n - 1], 2 n - 1, 0]; A085090 = Array[j, 98];
r[n_] := Table[A085090[[k]] + A085090[[n - k + 1]], {k, 1, n}];
countzeros[l_List] := Sum[KroneckerDelta[0, k], {k, l}];
Table[((x = n - 2 countzeros[A085090[[1 ;; n]]] + countzeros[r[n]]) +
KroneckerDelta[OddQ[x], True])/2, {n, 1, 98}] (* Fred Daniel Kline, Aug 30 2018 *)

A002375 := proc(n) local s,p; begin s := 0; p := 3; repeat if isprime(2*n-p) then s := s+1 end_if; p := nextprime(p+2); until p>n end_repeat; s end_proc:

A002375(n)=sum(i=2,primepi(n),isprime(2*n-prime(i))) /* ...i=1... gives A045917 */

apply( {A002375(n,s=0,N=2*n)=forprime(p=n, N-3, isprime(N-p)&&s++);s}, [1..100]) \\ M. F. Hasler, Jan 03 2023

from sympy import primerange, isprime
def A002375(n): return sum(1 for p in primerange(3,n+1) if isprime((n<<1)-p)) # Chai Wah Wu, Feb 20 2025

def A002375(n):
    P = primes(3, n+1)
    M = (2*n - p for p in P)
    F = [k for k in M if is_prime(k)]
    return len(F)
[A002375(n) for n in (1..98)] # Peter Luschny, May 19 2013

From Halberstam and Richert: a(n) < (8+0(1))*c(n)*n/log(n)^2 where c(n) = Product_{p > 2} (1-1/(p-1)^2)*Product_{p|n, p > 2} (p-1)/(p-2). It is conjectured that the factor 8 can be replaced by 2. Is a(n) > n/log(n)^2 for n large enough? - Benoit Cloitre, May 20 2002
a(n) = ceiling(A002372(n)/2). - Emeric Deutsch, Jul 14 2004
G.f.: Sum_{j>=2} Sum_{i=2..j} x^(p(i) + p(j)), where p(k) is the k-th prime. - Emeric Deutsch, Aug 27 2007
Not very efficient: a(n) = (Sum_{i=1..n} (pi(i) - pi(i-1)) * (pi(2n-i) - pi(2n-i-1))) - floor(2/n)*floor(n/2). - Wesley Ivan Hurt, Jan 06 2013
For n >= 2, a(n) = Sum_{3 <= p <= n, p is prime} A(2*n - p) - binomial(A(n), 2) - a(n-1) - a(n-2) - ... - a(1), where A(n) = A033270(n) (see Example 1 in link of V. Shevelev). - Vladimir Shevelev, Jul 08 2013

Beginning corrected by Paul Zimmermann, Mar 15 1996
More terms from James Sellers
Edited by Charles R Greathouse IV, Apr 20 2010

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

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 1, 1, 2, 1, 2, 2, 2, 1, 3, 2, 4, 2, 3, 2, 5, 2, 5, 3, 5, 3, 7, 3, 7, 2, 6, 3, 9, 2, 8, 4, 9, 4, 10, 2, 11, 3, 10, 4, 12, 3, 13, 4, 12, 5, 15, 4, 16, 3, 14, 5, 17, 3, 16, 4, 16, 6, 19, 3, 21, 5, 20, 6, 20, 2, 22, 5, 21, 6, 22, 5, 28, 5, 24, 7

Offset: 1

Naohiro Nomoto, Feb 24 2002

For even n > 2, a(n) = A061358(n-2). - Reinhard Zumkeller, Aug 08 2009
Vinogradov proved that every sufficiently large odd number is the sum of three primes. - T. D. Noe, Mar 27 2013
The two Helfgott papers show that every odd number greater than 5 is the sum of three primes (this is the Odd Goldbach Conjecture). - T. D. Noe, May 14 2013, N. J. A. Sloane, May 18 2013

			a(6) = 1 because 6 = 2+2+2,
a(9) = 2 because 9 = 2+2+5 = 3+3+3,
a(15) = 3 because 15 = 2+2+11 = 3+5+7 = 5+5+5,
a(17) = 4 because 17 = 2+2+13 = 3+3+11 = 3+7+7 = 5+5+7.
- _Zak Seidov_, Jun 29 2017

Robert G. Wilson v, Table of n, a(n) for n = 1..36000 (first 10000 terms from T. D. Noe)
H. A. Helfgott, Minor arcs for Goldbach's problem, arXiv:1205.5252 [math.NT], 2012.
H. A. Helfgott, Major arcs for Goldbach's theorem, arXiv:1305.2897 [math.NT], 2013.
H. A. Helfgott, The ternary Goldbach conjecture is true, arxiv:1312.7748 [math.NT], 2013.
H. A. Helfgott, The ternary Goldbach problem, arXiv:1404.2224 [math.NT], 2014.
Yannick Saouter, Checking the odd Goldbach conjecture up to 10^20, Math. Comp. 67 (222) (1998) 863-866.
Eric Weisstein's World of Mathematics, Vinogradov's Theorem
Wikipedia, Goldbach's conjecture.

Bisections: A045917, A054860. Cf. A002375, A007963, A061358, A059998.
First occurrence: A139321. Records: A139322.
Column k=3 of A117278.

f[n_] := Block[{c = 0, lmt = PrimePi@ Floor[n/2], p, q}, Do[p = Prime@ i; q = Prime@ j; r = n - p - q; If[ PrimeQ@ r && r >= p, c++ ], {i, lmt}, {j, i}]; c]; Array[f, 91] (* Robert G. Wilson v, Apr 13 2008 *)
Table[Count[IntegerPartitions[n,{3}],?(AllTrue[#,PrimeQ]&)],{n,50}] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale, Sep 10 2019 *)

a(n)=my(s); forprime(p=(n+2)\3,n-4, forprime(q=(n-p+1)\2,min(n-p-2,p), if(isprime(n-p-q), s++))); s \\ Charles R Greathouse IV, Jun 29 2017

from sympy import isprime, primerange, floor
def a(n):
    s=0
    for p in primerange(((n + 2)//3), n - 3):
        for q in primerange(((n - p + 1)//2), min(n - p - 2, p) + 1):
            if isprime(n - p - q): s+=1
    return s
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jul 01 2017, after PARI code by Charles R Greathouse IV

a(n) = Sum_{k=1..floor(n/3)} Sum_{i=k..floor((n-k)/2)} A010051(i) * A010051(k) * A010051(n-i-k). - Wesley Ivan Hurt, Mar 26 2019

a(n) = [x^n y^3] Product_{k>=1} 1/(1 - y*x^prime(k)). - Ilya Gutkovskiy, Apr 18 2019

More terms from Vladeta Jovovic, Mar 10 2002

A014092 Numbers that are not the sum of 2 primes.

Original entry on oeis.org

1, 2, 3, 11, 17, 23, 27, 29, 35, 37, 41, 47, 51, 53, 57, 59, 65, 67, 71, 77, 79, 83, 87, 89, 93, 95, 97, 101, 107, 113, 117, 119, 121, 123, 125, 127, 131, 135, 137, 143, 145, 147, 149, 155, 157, 161, 163, 167, 171, 173, 177, 179, 185, 187, 189, 191, 197, 203, 205, 207, 209

Offset: 1

N. J. A. Sloane

Suggested by the Goldbach conjecture that every even number larger than 2 is the sum of 2 primes.
Since (if we believe the Goldbach conjecture) all the entries > 2 in this sequence are odd, they are equal to 2 + an odd composite number (or 1).
Otherwise said, the sequence consists of 2 and odd numbers k such that k-2 is not prime. In particular there is no element from A006512, greater of a twin prime pair. - M. F. Hasler, Sep 18 2012
Values of k such that A061358(k) = 0. - Emeric Deutsch, Apr 03 2006
Values of k such that A073610(k) = 0. - Graeme McRae, Jul 18 2006

G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, Section 2.8 (for Goldbach conjecture).

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Index entries for sequences related to Goldbach conjecture

Cf. A002372, A002373, A002374, A048974, A061358.
Cf. A010051, A000040, A051035 (composites).
Equivalent sequence for prime powers: A071331.
Numbers that can be expressed as the sum of two primes in k ways for k=0..10: this sequence (k=0), A067187 (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).

a014092 n = a014092_list !! (n-1)
a014092_list = filter (\x ->
   all ((== 0) . a010051) $ map (x -) $ takeWhile (< x) a000040_list) [1..]
-- Reinhard Zumkeller, Sep 28 2011

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

s1falsifiziertQ[s_]:= Module[{ip=IntegerPartitions[s, {2}], widerlegt=False},Do[If[PrimeQ[ip[[i,1]] ] ~And~ PrimeQ[ip[[i,2]] ], widerlegt = True; Break[]],{i,1,Length[ip]}];widerlegt]; Select[Range[250],s1falsifiziertQ[ # ]==False&] (* Michael Taktikos, Dec 30 2007 *)
Join[{1,2},Select[Range[3,300,2],!PrimeQ[#-2]&]] (* Zak Seidov, Nov 27 2010 *)
Select[Range[250],Count[IntegerPartitions[#,{2}],?(AllTrue[#,PrimeQ]&)]==0&] (* _Harvey P. Dale, Jun 08 2022 *)

isA014092(n)=local(p,i) ; i=1 ; p=prime(i); while(pA014092(a), print(n," ",a); n++)) \\ R. J. Mathar, Aug 20 2006

from sympy import prime, isprime
def ok(n):
    i=1
    x=prime(i)
    while xIndranil Ghosh, Apr 29 2017

Odd composite numbers + 2 (essentially A014076(n) + 2 ).

Equals {2} union A005408 \ A052147, i.e., essentially the complement of A052147 (or rather A048974) within the odd numbers A005408. - M. F. Hasler, Sep 18 2012

A047845 a(n) = (m-1)/2, where m is the n-th odd nonprime (A014076(n)).

Original entry on oeis.org

0, 4, 7, 10, 12, 13, 16, 17, 19, 22, 24, 25, 27, 28, 31, 32, 34, 37, 38, 40, 42, 43, 45, 46, 47, 49, 52, 55, 57, 58, 59, 60, 61, 62, 64, 66, 67, 70, 71, 72, 73, 76, 77, 79, 80, 82, 84, 85, 87, 88, 91, 92, 93, 94, 97, 100, 101, 102, 103, 104, 106, 107, 108, 109, 110, 112, 115

Offset: 1

Enoch Haga

Also (starting with 2nd term) numbers of the form 2xy+x+y for x and y positive integers. This is also the numbers of sticks needed to construct a two-dimensional rectangular lattice of unit squares. See A090767 for the three-dimensional generalization. - John H. Mason, Feb 02 2004
Note that if k is not in this sequence, then 2*k+1 is prime. - Jose Brox (tautocrona(AT)terra.es), Dec 29 2005
Values of k for which A073610(2k+3)=0; values of k for which A061358(2k+3)=0. - Graeme McRae, Jul 18 2006
This sequence also arises in the following way: take the product of initial odd numbers, i.e., the product (2n+1)!/(n!*2^n) and factor it into prime numbers. The result will be of the form 3^f(3)*5^f(5)*7^f(7)*11^f(11)... . Then f(3)/f(5) = 2, f(3)/f(7) = 3, f(3)/f(11) = 5, ... and this sequence forms (for sufficiently large n, of course) the sequence of natural numbers without 4,7,10,12,..., i.e., these numbers are what is lacking in the present sequence. - Andrzej Staruszkiewicz (uszkiewicz(AT)poczta.onet.pl), Nov 10 2007
Also "flag short numbers", i.e., number of dots that can be arranged in successive rows of K, K+1, K, K+1, K, ..., K+1, K (assuming there is a total of L > 1 rows of size K > 0). Adapting Skip Garibaldi's terms, sequence A053726 would be "flag long numbers" because those patterns begin and end with the long lines. If you convert dots to sticks, you get the lattice that John H. Mason mentioned. - Juhani Heino, Oct 11 2014
Numbers k such that (2*k)!/(2*k + 1) is an integer. - Peter Bala, Jan 24 2017
Except for a(1)=0: numbers of the form k == j (mod 2j+1), j >= 1, k > 2j+1. - Bob Selcoe, Nov 07 2017

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Slate, Article about USA flag patterns -- this is where Skip Garibaldi gave definitions.

Complement of A005097.

a047845 = (`div` 2) . a014076  -- Reinhard Zumkeller, Jan 02 2013

[(n-1)/2 : n in [1..350] | (n mod 2) eq 1 and not IsPrime(n)]; // G. C. Greubel, Oct 16 2023

for n from 0 to 120 do
    if irem(factorial(2*n), 2*n+1) = 0 then print(n); end if;
end do:
# Peter Bala, Jan 24 2017

(Select[Range[1, 231, 2], PrimeOmega[#] != 1 &] - 1)/2 (* Jayanta Basu, Aug 11 2013 *)

print1(0,", ");
forcomposite(n=1,250,if(1==n%2,print1((n-1)/2,", "))); \\ Joerg Arndt, Oct 16 2023

from sympy import primepi
def A047845(n):
    if n == 1: return 0
    m, k = n-1, primepi(n) + n - 1 + (n>>1)
    while m != k:
        m, k = k, primepi(k) + n - 1 + (k>>1)
    return m-1>>1 # Chai Wah Wu, Jul 31 2024

[(n-1)/2 for n in (1..350) if n%2==1 and not is_prime(n)] # G. C. Greubel, Oct 16 2023

A193773(a(n)) > 1 for n > 1. - Reinhard Zumkeller, Jan 02 2013

Name edited by Jon E. Schoenfield, Oct 16 2023

A073610 Number of primes of the form n-p where p is a prime.

Original entry on oeis.org

0, 0, 0, 1, 2, 1, 2, 2, 2, 3, 0, 2, 2, 3, 2, 4, 0, 4, 2, 4, 2, 5, 0, 6, 2, 5, 0, 4, 0, 6, 2, 4, 2, 7, 0, 8, 0, 3, 2, 6, 0, 8, 2, 6, 2, 7, 0, 10, 2, 8, 0, 6, 0, 10, 2, 6, 0, 7, 0, 12, 2, 5, 2, 10, 0, 12, 0, 4, 2, 10, 0, 12, 2, 9, 2, 10, 0, 14, 0, 8, 2, 9, 0, 16, 2, 9, 0, 8, 0, 18, 2, 8, 0, 9, 0, 14, 0, 6

Offset: 1

Amarnath Murthy, Aug 05 2002

nonn

a(p) = 2 if p-2 is a prime else a(p) = 0. If n = 2p, p is a prime then a(n) is odd else a(n) is even. As p is counted only once and if q and n-q both are prime then the count is increased by 2. ( Analogous to the fact that perfect squares have odd number of divisors).
a(2k+1) = 2 if (2k-1) is prime, else a(2k+1)=0 (for any k). This sequence can be used to re-describe a couple of conjectures: the Goldbach conjecture == a(2n) > 0 for all n>=2; twin primes conjecture == for any n, there is a prime p>n s.t. a(p)>0.
Number of ordered ways of writing n as the sum of two primes.

			a(16) = 4 as there are 4 primes 3,5,11 and 13 such that 16-3,16-5,16-11and 16-13 are primes.

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

Cf. A061358, A065577, A107318, A098238

for i from 1 to 500 do a[i] := 0:j := 1:while(ithprime(j)

nn=20;a[x]:=Sum[x^i,{i,Table[Prime[n],{n,1,nn}]}];Drop[CoefficientList[a[x]^2,x],1]  (* Geoffrey Critzer, Nov 22 2012 *)

Vec(sum(i=1,100,x^prime(i),O(x^prime(101)))^2) \\ Charles R Greathouse IV, Jan 21 2015

G.f.: (Sum_{k>0} x^prime(k))^2. - Vladeta Jovovic, Mar 12 2005

Self-convolution of characteristic function of primes (A010051). - Graeme McRae, Jul 18 2006

Corrected and extended by Vladeta Jovovic and Sascha Kurz, Aug 06 2002

A048138 a(n) = number of m such that sum of proper divisors of m (A001065(m)) is n.

Original entry on oeis.org

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

Offset: 2

Naohiro Nomoto

The offset is 2 since there are infinitely many numbers (all the primes) for which A001065 = 1.
The graph of this sequence, shifted by 1, looks similar to that of A061358, which counts Goldbach partitions of n. - T. D. Noe, Dec 05 2008
For n > 2, a(n) <= A000009(n) as all divisor lists must have distinct values. - Roderick MacPhee, Sep 13 2016
The smallest k > 0 such that there are exactly n numbers whose sum of proper divisors is k is A125601(n). - Bernard Schott, Mar 23 2023

			a(6) = 2 since 6 is the sum of the proper divisors of 6 and 25.

T. D. Noe, Table of n, a(n) for n = 2..10000
Carl Pomerance, The first function and its iterates, pp. 125-138 in Connections in Discrete Mathematics, ed. S. Butler et al., Cambridge, 2018.

Cf. A001065, A064440, A125601, A238895, A238896 (records).

Cf. A005114, A057709, A057710, A160133.

with(numtheory): for n from 2 to 150 do count := 0: for m from 1 to n^2 do if sigma(m) - m = n then count := count+1 fi: od: printf(`%d,`,count): od:

list(n)=my(v=vector(n-1),k); for(m=4,n^2, k=sigma(m)-m; if(k>1 & k<=n, v[k-1]++)); v \\ Charles R Greathouse IV, Apr 21 2011

From Bernard Schott, Mar 23 2023: (Start)
a(n) = 0 iff n is in A005114 (untouchable numbers).
a(n) = 1 iff n is in A057709 ("hermit" numbers).
a(n) = 2 iff n is in A057710.
a(n) > 1 iff n is in A160133. (End)

More terms from James Sellers, Feb 19 2001

A071330 Number of decompositions of n into sum of two prime powers.

Original entry on oeis.org

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

Offset: 1

Reinhard Zumkeller, May 19 2002

a(2*n) > 0 (Goldbach's conjecture).

a(A071331(n)) = 0; A095840(n) = a(A000961(n)).

			10 = 1 + 3^2 = 2 + 2^3 = 3 + 7 = 5 + 5, therefore a(10) = 4;
11 = 2 + 3^2 = 3 + 2^3 = 4 + 7, therefore a(11) = 3;
12 = 1 + 11 = 3 + 3^2 = 2^2 + 2^3 = 5 + 7, therefore a(12) = 4;
a(149)=0, as for all x<149: if x is a prime power then 149-x is not.

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

Cf. A000961, A002375, A010055, A061358, A071331, A109829.

a071330 n = sum $
   map (a010055 . (n -)) $ takeWhile (<= n `div` 2) a000961_list
-- Reinhard Zumkeller, Jan 11 2013

primePowerQ[n_] := Length[ FactorInteger[n]] == 1; a[n_] := (r = 0; Do[ If[ primePowerQ[k] && primePowerQ[n-k], r++], {k, 1, Floor[n/2]}]; r); Table[a[n], {n, 1, 95}](* Jean-François Alcover, Nov 17 2011, after Michael B. Porter *)

ispp(n) = (omega(n)==1 || n==1)
A071330(n) = {local(r);r=0;for(i=1,floor(n/2),if(ispp(i) && ispp(n-i),r++));r} \\ Michael B. Porter, Dec 04 2009

a(n)=my(s); forprime(p=2,n\2,if(isprimepower(n-p), s++)); for(e=2,log(n)\log(2), forprime(p=2, sqrtnint(n\2,e), if(isprimepower(n-p^e), s++))); s+(!!isprimepower(n-1))+(n==2) \\ Charles R Greathouse IV, Nov 21 2014

A259201 Number of partitions of n into ten primes.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 4, 4, 5, 7, 8, 9, 11, 11, 14, 16, 18, 20, 25, 24, 31, 33, 38, 39, 48, 47, 59, 59, 69, 69, 87, 80, 102, 98, 118, 114, 143, 131, 168, 154, 191, 179, 227, 200, 261, 236, 297, 268, 344, 300, 396, 345, 442, 390, 509, 431, 576, 493, 641, 551, 729

Offset: 20

Doug Bell, Jun 20 2015

			a(23) = 2 because there are 2 partitions of 23 into ten primes: [2,2,2,2,2,2,2,2,2,5] and [2,2,2,2,2,2,2,3,3,3].

Alois P. Heinz, Table of n, a(n) for n = 20..10000
Index entries for sequences related to partitions

Column k=10 of A117278.
Number of partitions of n into r primes for r = 1-9: A010051, A061358, A068307, A259194, A259195, A259196, A259197, A259198, A259200.
Cf. A000040, A010051.

[#RestrictedPartitions(k,10,Set(PrimesUpTo(1000))):k in [20..80]] ; // Marius A. Burtea, Jul 13 2019

a(n) = [x^n y^10] Product_{k>=1} 1/(1 - y*x^prime(k)). - Ilya Gutkovskiy, Apr 18 2019

a(n) = Sum_{r=1..floor(n/10)} Sum_{q=r..floor((n-r)/9)} Sum_{p=q..floor((n-q-r)/8)} Sum_{o=p..floor((n-p-q-r)/7)} Sum_{m=o..floor((n-o-p-q-r)/6)} Sum_{l=m..floor((n-m-o-p-q-r)/5)} Sum_{k=l..floor((n-l-m-o-p-q-r)/4)} Sum_{j=k..floor((n-k-l-m-o-p-q-r)/3)} Sum_{i=j..floor((n-j-k-l-m-o-p-q-r)/2)} A010051(r) * A010051(q) * A010051(p) * A010051(o) * A010051(m) * A010051(l) * A010051(k) * A010051(j) * A010051(i) * A010051(n-i-j-k-l-m-o-p-q-r). - Wesley Ivan Hurt, Jul 13 2019

cp's OEIS Frontend

A141340 Positive integers n such that A061358(n) = #{primes p | n/2 <= p < n-1}.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A162197 a(n) = A161912(n) - A061358(n).

Views

Author

Keywords

Crossrefs

Formula

A002375 From Goldbach conjecture: number of decompositions of 2n into an unordered sum of two odd primes.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

MuPAD

PARI

PARI

Python

Sage

Formula

Extensions

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

A014092 Numbers that are not the sum of 2 primes.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

A047845 a(n) = (m-1)/2, where m is the n-th odd nonprime (A014076(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica