OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
f:= n -> add(t^2 mod n, t = numtheory:-divisors(n)):
select(isprime, map(f, [$1..1000]));
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Nov 16 2020
STATUS
approved