|
MAPLE
|
F:= proc(d) local t, g, Cands:
t:= 10^d+1;
if NumberTheory:-QuadraticResidue(85, t) <> 1 then return NULL fi;
Cands:= map(s -> rhs(op(s)), [msolve(x^2 + 9*x - 1, t)]);
g:= proc(r) local v; v:= r^2 + 9*r - 1; v >= t*(t-11)/10 and v < t*(t-2) end proc;
op(sort(select(g, Cands)));
end proc:
map
|