login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A059998
Number of different primes occurring when n is expressed as p1+q1+r1 = ... = pk+qk+rk where pk,qk,rk are primes with pk <= qk <= rk.
3
0, 0, 0, 0, 0, 1, 2, 2, 3, 3, 4, 4, 3, 3, 5, 4, 6, 5, 5, 5, 7, 5, 8, 6, 7, 7, 9, 6, 8, 5, 8, 7, 10, 5, 11, 8, 10, 9, 10, 4, 12, 7, 11, 9, 13, 7, 14, 8, 13, 11, 15, 9, 14, 7, 14, 11, 16, 7, 15, 8, 15, 13, 17, 6, 18, 11, 17, 13, 17, 5, 19, 11, 18, 13, 20, 10, 21, 11, 20, 15, 20, 9, 22, 10, 21
OFFSET
1,7
COMMENTS
Goldbach conjectured that every integer >5 is the sum of three primes. 6=2+2+2, 7=2+2+3, 8=2+3+3, 9=3+3+3=2+2+5,......
The largest possible value of a(n) is PrimePi(n)-1, which is frequently attained. - T. D. Noe, May 05 2008
FORMULA
If n is of the form 2*(prime + 1) then a(n) is an even number.
EXAMPLE
For n=9: 9 = 3+3+3 = 2+2+5, we can see 3 different primes. so a(9) = 3.
MATHEMATICA
a[n_] := Select[ Reverse /@ IntegerPartitions[n, {3}] , LessEqual @@ # && PrimeQ[#[[1]]] && PrimeQ[#[[2]]] && PrimeQ[#[[3]]] &] // Flatten // Union // Length; Table[a[n], {n, 1, 85}] (* Jean-François Alcover, Oct 03 2012 *)
PROG
(PARI) a(n)=my(v=List()); forprime(r=(n+2)\3, n-4, forprime(q=(n-r+1)\2, n-r-2, if(isprime(n-r-q), listput(v, r); listput(v, q); listput(v, n-r-q)))); #vecsort(Vec(v), , 8) \\ Charles R Greathouse IV, Jul 14 2013
CROSSREFS
Sequence in context: A332250 A340321 A340320 * A361384 A339731 A234475
KEYWORD
easy,nonn,nice
AUTHOR
Naohiro Nomoto, Mar 10 2001
STATUS
approved