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