OFFSET
1,1
COMMENTS
The union of [{5, 7}, this sequence, A282324 and A282326] is the greater of twin primes sequence A006512.
A181604 without the 3. [Proof: working mod 10 we see that each value here is in A181604. For the other direction: Except 3 all twin primes in A181604 are upper twin primes; they cannot be lower twin primes because the upper ones would be multiples of 5. The twin primes in A181604 could be == 3 (mod 30) or == 13 (mod 30) or == 23 (mod 30). The first case is excluded because they would be multiples of 3; the third case is excluded because the lower twin primes would be == 21 (mod 30) and also multiples of 3. So only the case == 13 (mod 30) remains.] - R. J. Mathar, Feb 14 2017
Number of terms < 10^k for k >= 1: 0, 3, 13, 67, 401, 2736, 19797, 146841, 1141217, 9137078, ..., . - Robert G. Wilson v, Jan 07 2018
LINKS
Muniru A Asiru, Table of n, a(n) for n = 1..20000
MAPLE
a:={}:
for i from 1 to 1229 do
if isprime(ithprime(i)-2) and ithprime(i) mod 30 = 13 then
a:={op(a), ithprime(i)}:
fi:
od:
a;
MATHEMATICA
Select[13 + 30 Range[0, 200], PrimeQ[# - 2] && PrimeQ[#] &] (* Robert G. Wilson v, Jan 07 2018 *)
PROG
(PARI) list(lim)=my(v=List(), p=2); forprime(q=3, lim, if(q-p==2 && q%30==13, listput(v, q)); p=q); Vec(v) \\ Charles R Greathouse IV, Feb 14 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Martin Renner, Feb 11 2017
STATUS
approved