cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A339309 Numbers k such that A228553(k) is twice a prime.

Original entry on oeis.org

2, 5, 7, 8, 9, 10, 14, 16, 34, 53, 55, 83, 91, 140, 172, 180, 182, 187, 189, 242, 246, 254, 265, 288, 361, 398, 404, 421, 435, 490, 508, 559, 574, 623, 706, 727, 736, 747, 770, 773, 777, 984, 985, 1065, 1093, 1106, 1125, 1187, 1245, 1256, 1287, 1341, 1346, 1361, 1372, 1381, 1449, 1467, 1573, 1596
Offset: 1

Views

Author

J. M. Bergot and Robert Israel, Nov 29 2020

Keywords

Examples

			a(3) = 7 is a term because A228553(7) = 82 = 2*41 and 41 is prime.
		

Crossrefs

Programs

  • Maple
    f:= proc(n)
      local S;
      S:= select(t -> isprime(t) and isprime(2*n-t), [seq(i, i=3..n, 2)]);
      add(t*(2*n-t), t=S)
    end proc:
    f(2):= 4:
    filter:= proc(t) local v; v:= f(t); v::even and isprime(v) end proc:
    select(filter, [$1..2000]);