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