login
A338701
Sums of two odd-indexed primes.
1
4, 7, 10, 13, 16, 19, 22, 25, 28, 33, 34, 36, 40, 42, 43, 46, 48, 49, 52, 54, 58, 61, 62, 64, 69, 70, 72, 75, 76, 78, 82, 84, 85, 88, 90, 94, 96, 98, 99, 100, 102, 104, 105, 106, 108, 111, 114, 118, 120, 124, 126, 128, 129, 130, 132, 134, 138, 139, 140, 142, 144, 146, 148, 150
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
Cf. A031368 (odd-indexed primes), A338703 (sums of two even-indexed primes).
Sequence in context: A016777 A308014 A143460 * A310679 A143459 A143458
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Apr 24 2021
STATUS
approved