OFFSET
1,2
MATHEMATICA
a[n_] := DivisorSum[n, 1 &, PowerMod[#, #, n] == Mod[-PowerMod[-#, n, n], n] &]; Array[a, 100] (* Amiram Eldar, Sep 20 2025 *)
PROG
(Magma) [1 + #[d: d in [1..n-1] | n mod d eq 0 and Modexp(d, d, n) eq -Modexp(-d, n, n) mod n]: n in [1..85]];
(PARI) a(n) = sumdiv(n, d, Mod(d, n)^d == - Mod(-d, n)^n); \\ Michel Marcus, Sep 16 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Sep 15 2025
STATUS
approved