A225945 Numbers k such that prime(k) divides k^k + 1.

1, 6, 60, 136, 124796, 3919272, 18363918, 153037808, 965108649, 3140421892, 5961162423, 20437804784

Offset

1,2

Comments

a(12) > 2*10^10. - Giovanni Resta, May 23 2013

Links

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

Mathematica

Select[Range[10^6], (p = Prime[#]; PowerMod[#, #, p] == p - 1) &] (* Giovanni Resta, May 23 2013 *)

Prog

(Python)
from sympy import nextprime, prime
from itertools import count, islice
def agen(startn=1): # generator of terms
    pn = prime(startn)
    for n in count(startn):
        if pow(n, n, pn) == pn - 1:
            yield n
        pn = nextprime(pn)
print(list(islice(agen(), 5))) # Michael S. Branicky, May 25 2023

Crossrefs

Cf. A000040, A000312, A014566, A225944.

Sequence in context: A250070 A036283 A228838 * A337691 A126576 A232181

Adjacent sequences: A225942 A225943 A225944 * A225946 A225947 A225948

Keyword

nonn,more

Author

Alex Ratushnyak, May 21 2013

Extensions

a(6)-a(11) from Giovanni Resta, May 23 2013

a(12) from Michael S. Branicky, May 25 2023

Status

approved