OFFSET
1,2
MATHEMATICA
Table[Length[Select[Divisors[n], PowerMod[#, #, n] == Mod[n - PowerMod[#, #, n], n] &]], {n, 1, 100}] (* Vaclav Kotesovec, Aug 23 2025 *)
PROG
(Magma) [1 + #[d: d in [1..n-1] | n mod d eq 0 and Modexp(d, d, n) eq -Modexp(d, d, n) mod n]: n in [1..100]];
(PARI) a(n) = sumdiv(n, d, 2*Mod(d, n)^d == 0); \\ Michel Marcus, Aug 30 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Aug 20 2025
STATUS
approved
