OFFSET
1,1
COMMENTS
Numbers k such that k + A000203(k) is the product of a twin prime pair.
EXAMPLE
a(3) = 50 is a term because 50 + A000203(50) = 143 = 11*13 where (11,13) is a twin prime pair.
MAPLE
filter:= proc(n) local t, s; t:= n + numtheory:-sigma(n);
if t::even or not issqr(t+1) then return false fi;
s:= sqrt(t+1);
isprime(s+1) and isprime(s-1)
end proc:
select(filter, [$1..3000000]);
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Mar 23 2022
STATUS
approved