OFFSET
1,1
COMMENTS
The union of [{5, 7}, A282322, this sequence and A282326] is the greater of twin primes sequence A006512.
Number of terms less than 10^k, k=2,3,4,...: 1, 9, 64, 414, 2734, 19674, 146953, ... - Muniru A Asiru, Feb 09 2018
LINKS
Muniru A Asiru, Table of n, a(n) for n = 1..2000
MAPLE
a:={}:
for i from 1 to 1229 do
if isprime(ithprime(i)-2) and ithprime(i) mod 30 = 19 then
a:={op(a), ithprime(i)}:
fi:
od:
a;
# More efficient
select(n -> isprime(n-2) and isprime(n), [seq(30*k+19, k=0..220)]); # Muniru A Asiru, Jan 30 2018
MATHEMATICA
Select[Prime[Range[1000]], PrimeQ[# - 2] && Mod[#, 30] == 19 &] (* Vincenzo Librandi, Feb 13 2017 *)
PROG
(Magma) [p: p in PrimesUpTo(7000) | IsPrime(p-2) and p mod 30 eq 19 ]; // Vincenzo Librandi, Feb 13 2017
(PARI) list(lim)=my(v=List(), p=2); forprime(q=3, lim, if(q-p==2 && q%30==19, listput(v, q)); p=q); Vec(v) \\ Charles R Greathouse IV, Feb 14 2017
(GAP) Filtered(List([1..220], k -> 30*k-11), n -> IsPrime(n) and IsPrime(n-2)); # Muniru A Asiru, Feb 02 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Martin Renner, Feb 11 2017
STATUS
approved