OFFSET
1,6
MATHEMATICA
a[n_] := Count[Range[0, n-1], _?(PowerMod[#, #, n] == PowerMod[-#, #, n] == # &)]; Array[a, 100] (* Amiram Eldar, Jul 03 2025 *)
PROG
(Magma) [#[s: s in [0..n-1] | Modexp(s, s, n) eq s and Modexp(-s, s, n) eq s]: n in [1..100]];
(PARI) a(n) = sum(s=0, n-1, (s == Mod(s, n)^s) && (s == Mod(-s, n)^s)); \\ Michel Marcus, Jul 09 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Jul 02 2025
STATUS
approved
