OFFSET
0,3
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..5000
FORMULA
a(n) <= A003573(n). - Michel Marcus, Jun 10 2019
MATHEMATICA
Table[Module[{m=1}, While[PowerMod[4, m, 4 n + 1]!=1 && PowerMod[4, m, 4 n + 1]!=4 n, m++]; m], {n, 0, 70}] (* Vincenzo Librandi, Jun 10 2019 *)
PROG
(PARI) a(n) = {my(m=1); while (((mm=Mod(4, 4*n+1)^m) != 1) && (mm != 4*n), m++); m; } \\ Michel Marcus, Jun 10 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved