robinhouston

sage: digits = 3
....: 
....: for a in range(10^(digits-1), 10^digits - 10^(digits - 1) + 2):
....:     x = (10^digits + 1)*a - a^2
....:     for b in divisors(x):
....:         if 10^(digits - 1) <= b < 10^digits and 10^(digits - 1) <= x//b < 10^digits:
....:             print(a, b, "|", x//b, 10^digits + 1 - a)
....:             
100 100 | 901 901
100 106 | 850 901