A338701 Sums of two odd-indexed primes.

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

Links

Table of n, a(n) for n=1..64.

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

Adjacent sequences: A338698 A338699 A338700 * A338702 A338703 A338704

Keyword

nonn

Author

Wesley Ivan Hurt, Apr 24 2021

Status

approved