OFFSET
1,1
MATHEMATICA
Take[Union[Total /@ Tuples[Prime[2 Range[70]], 2]], 70]
PROG
(Python)
from sympy import primerange
def aupto(limit):
ep = [p for p in primerange(1, limit)][1::2]
ss = sorted(set(p1 + p2 for i, p1 in enumerate(ep) for p2 in ep[i:]))
return [s for s in ss if s <= limit]
print(aupto(162)) # Michael S. Branicky, Apr 24 2021
(PARI) my(N=18); select(x->(x<prime(2)+prime(2*N+2)), setbinop((x, y)->(x+y), vector(N, n, prime(2*n)))) \\ Michel Marcus, Apr 24 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Apr 24 2021
STATUS
approved