OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(3) = 10 is a term because A318996(10) = 11 is prime.
MAPLE
f:= proc(n) local s, D, t;
D:= numtheory:-divisors(n);
s:= convert(D, `+`);
add(s mod t, t=D)
end proc:
select(t -> isprime(f(t)), [$1..1000]);
PROG
(PARI) f(n) = my(sn = sigma(n)); sumdiv(n, d, sn % d); \\ A318996
isok(k) = isprime(f(k)); \\ Michel Marcus, Feb 04 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Feb 03 2021
STATUS
approved