OFFSET
1,1
COMMENTS
From Robert Israel, Apr 08 2025: (Start)
Numbers k such that k * (k + 1) = (10^d + 1) * m for some d and m where m has d digits.
Includes (10^(3*d)-1)/3 + (10^d-1)*10^d/3 and 2*(10^(3*d)-1)/3 - (10^d-1)*10^d/3 + 1 for all d >= 1. (End)
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
A161356(n) = a(n)*(a(n)+1). - Michael S. Branicky, Jul 11 2025
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:
seq(op(q(d, 1)), 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
