OFFSET
1,2
MATHEMATICA
a[n_] := DivisorSum[n, 1 &, PowerMod[#, #, n] == PowerMod[-#, #, n] == Mod[#, n] &]; Array[a, 100] (* Amiram Eldar, Jul 03 2025 *)
PROG
(Magma) [1 + #[d: d in [1..n-1] | n mod d eq 0 and Modexp(d, d, n) eq d and Modexp(-d, d, n) eq d]: n in [1..100]];
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jul 02 2025
STATUS
approved
