OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
PROG
(PARI) pnpk(n, m=2, k=2) = { local(x, v1, v2); for(x=1, n, v1 = prime(x)-k; v2 = prime(x+m)-k; if(isprime(v1)&isprime(v2), print1(x, ", ") ) ) ; } \\ corrected by Amiram Eldar, Oct 04 2024
(PARI) lista(pmax) = {my(k = 1, p = primes(4)); forprime(p1 = p[#p], pmax, k++; p[#p] = p1; if(p[2]- p[1] == 2 && p[4] - p[3] == 2, print1(k, ", ")); for(i = 1, #p-1, p[i] = p[i+1])); } \\ Amiram Eldar, Oct 04 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Cino Hilliard, May 01 2005
EXTENSIONS
Offset corrected by Amiram Eldar, Oct 04 2024
STATUS
approved