Если можно с объяснением

uinput = input('трёхзначное натуральное число: ') # 123

numbers = [int(i) for i in uinput] # [1, 2, 3]

highest = max(numbers) # 3

lowest = min(numbers) # 1

numbers.pop(numbers.index(highest)) # [1, 2, 3] = [1, 2]

numbers.pop(numbers.index(lowest)) # [1, 2] = [2]

middle = numbers[0] # 2

# print(highest, lowest, middle) 3 1 2

num1 = str(highest) + str(middle) # '32'

if lowest == 0:

~~~~if middle == 0:

~~~~~~~~num2 = str(highest) + str(lowest)

~~~~else:

~~~~~~~~num2 = str(middle) + str(lowest)

else:

~~~~num2 = str(lowest) + str(middle) # '12'

print(num1 + ' - ' + num2 + ' = ' + str(int(num1) - int(num2))) # 32 - 12 = 20

Выполняя алгоритм, получаем следующий результат (15 итераций)

1. 0..65534 -> 32767
2. 0..32766 -> 16383
3. 0..16382 -> 8191
4. 0..8190  -> 4095
5. 0..4094  -> 2047
6. 2048..4094 -> 3071
7. 2048..3070 -> 2559
8. 2560..3070 -> 2815
9. 2816..3070 -> 2943
10. 2944..3070 -> 3007
11. 2944..3006 -> 2975
12. 2976..3006 -> 2991
13. 2992..3006 -> 2999
14. 3000..3006 -> 3003
15. 3000..3002 -> 3001

Если лень перебирать вручную, можно воспользоваться программой

var k,l,r,x,f:integer;
begin
f := 3001;
l := 0;
r := 65534;
x := (l + r) div 2;
k := 1;
while (x <> f) and (l < r) do
  begin
  writeln(k,' ',l,' ',r,' ',x);
  k := k + 1;
  if f < x then r := x - 1
    else l := x + 1;
  x := (l + r) div 2
  end;
writeln(k,' ',l,' ',r,' ',x);
end.