OFFSET
1,6
MATHEMATICA
a[n_] := DivisorSum[n, 1 &, Mod[-PowerMod[#, n, n], n] == # &]; Array[a, 100] (* Amiram Eldar, Feb 04 2025 *)
PROG
(Magma) [#[d: d in Divisors(n) | -d^n mod n eq d]: n in [1..100]];
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Feb 03 2025
STATUS
approved
