A340876 Numbers k such that prime(k)^prime(k+1) == prime(k+2) (mod prime(k+3)).

942, 4658911, 12806325, 2515276754

Offset

1,1

Comments

Inspired by A340868 and the comment there.

Links

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

Example

The 942th to the 945th primes are 7433, 7451, 7457 and 7459. Since 7433^7451 == 7457 mod 7459, 942 is in the sequence.

Prog

(Python)
from sympy import nextprime
k, p, q, r, s, A340876_list = 1, 2, 3, 5, 7, []
while k < 10**7:
    if pow(p, q, s) == r:
        A340876_list.append(k)
    k, p, q, r, s = k+1, q, r, s, nextprime(s)

Crossrefs

Cf. A340868.

Sequence in context: A158718 A104302 A228382 * A252295 A252294 A268268

Adjacent sequences: A340873 A340874 A340875 * A340877 A340878 A340879

Keyword

nonn,more

Author

Chai Wah Wu, Jan 24 2021

Extensions

a(4) from Daniel Suteu, Jan 26 2021

Status

approved