OFFSET
1,1
COMMENTS
For the first 30000 terms a(n) > A121762(n), see plot A121764(n) - A121762(n). But is it so for all n? - Zak Seidov, Apr 25 2015
Subsequence of A002476. - Michel Marcus, Apr 26 2015
LINKS
Zak Seidov, Table of n, a(n) for n = 1..10000
Zak Seidov, Plot of A121764(n) - A121762(n)
MATHEMATICA
Select[Table[6n + 1, {n, 200}], PrimeQ[#] && !PrimeQ[#-2] &] (* Ray Chandler, Aug 22 2006 *)
Select[Prime[Range[200]], Mod[#, 6]==1&&NoneTrue[#+{2, -2}, PrimeQ]&] (* Harvey P. Dale, Jul 16 2021 *)
PROG
(PARI) {is(n)=n%6==1 && isprime(n) && !isprime(n-2)}; \\ G. C. Greubel, Feb 26 2019
(Magma) [n: n in [1..1150] | (n mod 6 eq 1) and not IsPrime(n-2) and IsPrime(n)]; // G. C. Greubel, Feb 26 2019
(Sage) [n for n in (1..1150) if mod(n, 6)==1 and not is_prime(n-2) and is_prime(n)] # G. C. Greubel, Feb 26 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Aug 20 2006
EXTENSIONS
Extended by Ray Chandler, Aug 22 2006
STATUS
approved