OFFSET
1,3
LINKS
Karl-Heinz Hofmann, Table of n, a(n) for n = 1..10000
MATHEMATICA
p = Select[Prime[Range[400]], PrimeQ[2*# + 1] &]; Mod[p, Range[Length[p]]] (* Amiram Eldar, Jan 11 2022 *)
PROG
(Python)
from sympy import isprime
n = 1
for p in range (2, 10000):
if isprime(p) and isprime(2*p+1):
print (p % n, end=", ")
n += 1
(PARI) lista(nn) = my(v=select(p->isprime(2*p+1), primes(nn))); vector(#v, k, v[k] % k); \\ Michel Marcus, Jan 11 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Karl-Heinz Hofmann, Jan 10 2022
STATUS
approved