OFFSET
1,1
COMMENTS
A prime p belongs to A068209 if and only if p = 5 mod 6 and there are integers x with (x+1)^p - x^p - 1 = 0 mod p^2 and gcd(x^2+x,p) = 1.
This sequence is the subsequence of A068209 of primes p for which no such x solves x^x + (x+1)^(x+1) = 0 mod p^2.
For all other primes p < 1486577 in A068209, simultaneous solutions have been found by computing discrete logarithms.
LINKS
David Broadhurst, On roots of n^n + (n+1)^(n+1) = 0 mod p^2
Kevin Brown, On the Density of Some Exceptional Primes
EXAMPLE
To prove that a(3) = 58073, we first show that (x+1)^p - x^p - 1 mod p^2, with gcd(x^2+x,p) = 1, has solutions when p = 58073 only for the residues x = r, -r/(1+r), 1/r, -(1+r), -1/(1+r), -(1+1/r) mod p, with r = 1281. By examining the orders of 1+1/r, 1+r, -r mod p, we prove that no x in this equivalence class can satisfy x^x + (x+1)^(x+1) = 0 mod p^2.
Similarly, we prove the absence of simultaneous roots for p = 37493, with r = 3730, and for p = 51941, with r = 15579.
By computing discrete logarithms, we provide simultaneous solutions for all other primes in A068209 with p < 58073.
CROSSREFS
KEYWORD
nonn
AUTHOR
David Broadhurst, Sep 13 2009
STATUS
approved