OFFSET
1,2
MATHEMATICA
f[n_] := DivisorSum[n, 1 &, PowerMod[-#, #, n] == n - # &]; With[{t = Array[f, 10^6]}, TakeWhile[FirstPosition[t, #] & /@ Range[Max[t]] // Flatten, ! MissingQ[#] &]] (* Amiram Eldar, Jul 26 2025 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Juri-Stepan Gerasimov, Jul 25 2025
EXTENSIONS
a(25)-a(29) from Amiram Eldar, Jul 26 2025
STATUS
approved
