login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A339309
Numbers k such that A228553(k) is twice a prime.
2
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,1
LINKS
EXAMPLE
a(3) = 7 is a term because A228553(7) = 82 = 2*41 and 41 is prime.
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]);
CROSSREFS
Sequence in context: A080639 A186306 A047483 * A167408 A047388 A284529
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Nov 29 2020
STATUS
approved