OFFSET
1,1
COMMENTS
From Robert Israel, Apr 08 2025: (Start)
Numbers k such that k * (k + 2) = (10^d + 1) * m for some d and m where m has d digits.
Contains 10^d - 1 for all d >= 1. (End)
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
q:= proc(d, m) local R, t, a, b, x, q;
t:= 10^d+1;
R:= NULL;
for a in numtheory:-divisors(t) do
b:= t/a;
if igcd(a, b) > 1 then next fi;
for x from chrem([0, -m], [a, b]) by t do
q:= x*(x+m)/t;
if q >= 10^d then break fi;
if q >= 10^(d-1) then R:= R, x fi;
od od;
sort(convert({R}, list));
end proc:
A:=[seq(op(q(d, 2)), d=1..10)]; # Robert Israel, Apr 08 2025
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Giovanni Resta, Feb 06 2006
EXTENSIONS
More terms from Robert Israel, Apr 08 2025
STATUS
approved
