OFFSET
1,1
MATHEMATICA
a[n_] := Module[{p = Prime[n], b = 2}, While[PowerMod[b, p - 1, p^5] != 1, b++]; b]; Array[a, 12] (* Amiram Eldar, May 12 2022 *)
PROG
(PARI) a(n) = my(p=prime(n)); for(b=2, oo, if(Mod(b, p^5)^(p-1)==1, return(b)))
(Python)
from sympy import prime
from sympy.ntheory.residue_ntheory import nthroot_mod
def A353938(n): return 2**5+1 if n == 1 else int(nthroot_mod(1, (p:= prime(n))-1, p**5, True)[1]) # Chai Wah Wu, May 17 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, May 12 2022
STATUS
approved