|
%I #18 Feb 15 2021 22:41:04
%S 0,0,0,0,0,0,1,1,0,1,2,1,1,1,1,1,2,2,1,1,1,2,2,3,2,2,3,2,2,2,1,3,3,1,
%T 3,3,3,3,5,3,2,4,4,2,4,3,3,4,1,3,4,2,4,4,3,4,5,4,4,6,2,3,5,2,4,5,3,3,
%U 4,3,4,5,2,2,5,2,4,5,3,4,5,3,3,8,5,3,6,4,4,8,4,4,7,3,4,6,5,6,7,5
%N Number of Goldbach partitions (p,q) of 2n such that |p-q| > n.
%H Robert Israel, <a href="/A279794/b279794.txt">Table of n, a(n) for n = 1..10000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/GoldbachPartition.html">Goldbach Partition</a>
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Goldbach%27s_conjecture">Goldbach's conjecture</a>
%H <a href="/index/Go#Goldbach">Index entries for sequences related to Goldbach conjecture</a>
%H <a href="/index/Par#part">Index entries for sequences related to partitions</a>
%F a(n) = Sum_{i=3..n-1} A010051(i) * A010051(2n-i) * (1-sign(floor(n/(2*(n-i))))).
%p with(numtheory): A279794:=n->add( (pi(i)-pi(i-1)) * (pi(2*n-i)-pi(2*n-i-1)) * (1-signum(floor(n/(2*(n-i))))), i=3..n-1): seq(A279794(n), n=1..100);
%p # Alternative:
%p f:= proc(n) local p;
%p nops(select(t -> isprime(t) and isprime(2*n-t), [seq(p,p=3..(n-1)/2,2)]))
%p end proc:
%p map(f, [$1..100]); # _Robert Israel_, Feb 15 2021
%t Table[Sum[Boole[PrimeQ@ i] Boole[PrimeQ[2 n - i]] (1 - Sign@ Floor[n/(2 (n - i))]), {i, 3, n - 1}], {n, 100}] (* _Michael De Vlieger_, Dec 21 2016 *)
%Y Cf. A010051, A002375, A279792.
%K nonn,easy
%O 1,11
%A _Wesley Ivan Hurt_, Dec 18 2016
|
