OFFSET
1,1
COMMENTS
a(n) is the smallest 10th root of unity (mod 11^n) that is prime. - Robert Israel, Jan 14 2024
LINKS
Robert Israel, Table of n, a(n) for n = 1..954
W. Keller and J. Richstein, Fermat quotients that are divisible by p.
MAPLE
f:= proc(n) local R, r, i;
R:= sort(map(rhs@op, [msolve(x^10=1, 11^n)]));
for i from 0 do
for r in R do
if isprime(11^n * i + r) then return 11^n * i + r fi
od od;
end proc:
map(f, [$1..20]); # Robert Israel, Jan 14 2024
MATHEMATICA
spp[n_]:=Module[{p=2, c=11^n}, While[PowerMod[p, 10, c]!=1, p=NextPrime[p]]; p]; Array[spp, 16] (* Harvey P. Dale, Aug 08 2019 *)
PROG
(PARI) \\ See A125609
(Python)
from itertools import count
from sympy import nthroot_mod, isprime
def A125612(n):
m = 11**n
r = sorted(nthroot_mod(1, 10, m, all_roots=True))
for i in count(0, m):
for p in r:
if isprime(i+p): return i+p # Chai Wah Wu, May 02 2024
CROSSREFS
Cf. A125609 = Smallest prime p such that 3^n divides p^2 - 1. Cf. A125610 = Smallest prime p such that 5^n divides p^4 - 1. Cf. A125611 = Smallest prime p such that 7^n divides p^6 - 1. Cf. A125632 = Smallest prime p such that 13^n divides p^12 - 1. Cf. A125633 = Smallest prime p such that 17^n divides p^16 - 1. Cf. A125634 = Smallest prime p such that 19^n divides p^18 - 1. Cf. A125635 = Smallest prime p such that 257^n divides p^256 - 1.
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Nov 28 2006
EXTENSIONS
More terms from Ryan Propper, Jan 03 2007
More terms from Martin Fuller, Jan 11 2007
More terms from Robert Israel, Jan 14 2024
STATUS
approved