OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..2000
MAPLE
filter:= proc(n) local F, t;
F:= ifactors(n)[2];
n*mul((t[1]-1)*t[1]^(t[2]-1), t=F)/2 mod add(t[1]*t[2], t=F) = 1;
end proc:
select(filter, [$2..50000]);
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Jan 03 2021
STATUS
approved