2,1

a(n) = Max{A000120(p)+A000120(q): p,q prime and p+q=2*n);

a(n) = A108423(n) + A108421(n).

Robert Israel, Table of n, a(n) for n = 2..10000

Index entries for sequences related to Goldbach conjecture

Index entries for sequences related to binary expansion of n

N:= 200: # to get a(2)..a(N)

Primes:= select(isprime, [seq(i, i=3..2*N-3, 2)]):

Ones:= map(t -> convert(convert(t, base, 2), `+`), Primes):

V:= Vector(N): V[2]:= 2:

for i from 1 to nops(Primes) do

p:= Primes[i];

for j from 1 to i do

k:= (p+Primes[j])/2;

if k > N then break fi;

t:= Ones[i]+Ones[j];

if t > V[k] then V[k]:= t fi

od

od:

convert(V[2..N], list); # Robert Israel, Mar 26 2018

Cf. A004676, A005843, A007088, A108421.

Sequence in context: A036443 A036437 A053306 * A276523 A244320 A084616

Adjacent sequences: A108419 A108420 A108421 * A108423 A108424 A108425

nonn,base

Reinhard Zumkeller, Jun 03 2005

approved