OFFSET
1,2
MATHEMATICA
A386410[n_] := DivisorSum[n, 1 &, PowerMod[#, #, n] == PowerMod[-#, #, n] == Mod[#, n] &]; seq[lim_] := Module[{s = Array[A386410, lim]}, TakeWhile[Flatten[FirstPosition[s, #] & /@ Range[Max[s]]], !MissingQ[#] &]]; seq[10^5] (* Amiram Eldar, Jul 21 2025 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Juri-Stepan Gerasimov, Jul 21 2025
EXTENSIONS
a(8)-a(11) from Amiram Eldar, Jul 21 2025
STATUS
approved
