|
%I #34 Jul 31 2018 07:22:54
%S 2,3,3,3,5,3,3,5,3,3,5,3,5,7,3,3,5,7,3,5,3,3,5,3,5,7,3,5,7,3,3,5,7,3,
%T 5,3,3,5,7,3,5,3,5,7,3,5,7,19,3,5,3,3,5,3,3,5,3,5,7,13,11,13,19,3,5,3,
%U 5,7,3,3,5,7,11,11,3,3,5,7,3,5,7,3,5,3,5,7,3,5,7,3,3,5,7,11,11,3,3,5,3,3,5,7
%N Least p with p, q both prime, p+q = 2n.
%C Essentially the same as A002373, which does not have the a(2) term. - _T. D. Noe_, Sep 24 2007
%C a(n) = A171637(n,1). - _Reinhard Zumkeller_, Mar 03 2014
%C Conjecture: a(n) ~ O(n^1/2). - _Jon Perry_, Apr 29 2014
%H H. J. Smith, <a href="/A020481/b020481.txt">Table of n, a(n) for n = 2..20000</a>
%H <a href="/index/Go#Goldbach">Index entries for sequences related to Goldbach conjecture</a>
%F a(n) = n - A047949(n). - _Jason Kimberley_, Oct 09 2012
%t a[n_] := For[p = 2, True, p = NextPrime[p], If[PrimeQ[2n-p], Return[p]]];
%t Table[a[n], {n, 2, 103}] (* _Jean-François Alcover_, Jul 31 2018 *)
%o (PARI) A020481(n) = {local(np);np=1;while(!isprime(2*n-prime(np)),np++);prime(np)} \\ _Michael B. Porter_, Dec 11 2009
%o (PARI) A020481(n)=forprime(p=1,n,isprime(2*n-p)&return(p)) \\ _M. F. Hasler_, Sep 18 2012
%o (Haskell)
%o a020481 n = head [p | p <- a000040_list, a010051' (2 * n - p) == 1]
%o -- _Reinhard Zumkeller_, Jul 07 2014, Mar 03 2014
%Y Cf. A020482.
%K nonn
%O 2,1
%A _David W. Wilson_
|
