OFFSET
1,2
MATHEMATICA
s[n_] := DivisorSum[n, 1 &, PowerMod[#, #, n] == n - # &]; With[{v = Array[s, 45000]}, TakeWhile[Flatten[FirstPosition[v, #] & /@ Range[Max[v]]], NumberQ]] (* Amiram Eldar, Jul 03 2025 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Juri-Stepan Gerasimov, Jun 30 2025
EXTENSIONS
a(8)-a(10) from Amiram Eldar, Jul 03 2025
STATUS
approved
