A363752 Primes prime(k) such that prime(k) mod k is prime.

5, 7, 17, 19, 23, 41, 47, 53, 61, 71, 79, 89, 101, 107, 113, 127, 131, 137, 139, 151, 163, 167, 173, 181, 191, 193, 197, 211, 223, 227, 229, 233, 239, 241, 257, 269, 277, 281, 313, 317, 347, 359, 367, 373, 383, 397, 421, 433, 443, 457, 463, 479, 503, 521, 541

Offset

1,1

Links

Table of n, a(n) for n=1..55.

Formula

a(n) = A000040(A363751(n)).

Example

The 9th prime is 23 and 23 mod 9 = 5, which is prime, so 23 is a term.

Mathematica

Table[If[PrimeQ[Mod[Prime[k], k]], Prime[k], Nothing], {n, k, 100}]

Prog

(Python)
from sympy import prime, isprime
a363752=[]
for k in range(1, 101):
    if isprime(prime(k)%k):
        a363752.append(prime(k))

Crossrefs

Cf. A000040, A004648, A363751.

Sequence in context: A253085 A078764 A216731 * A045318 A040102 A100021

Adjacent sequences: A363749 A363750 A363751 * A363753 A363754 A363755

Keyword

nonn

Author

Nicholas Leonard, Jun 19 2023

Status

approved