A386410 Smallest k for which A385541(k) = n.
1, 2, 66, 30, 1386, 858, 34650, 1287678, 6380946, 11411400, 24877650
Offset: 1
Programs
-
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 *)
Extensions
a(8)-a(11) from Amiram Eldar, Jul 21 2025