Задача по теории чисел - программированием

from itertools import *
 
s = 11425
p = 22612
 
# генерируем все варианты трехзначных чисел, соответствующих условию
a = product('0123456789', repeat=3)
words = []
for w in a:
    if w[0] != '0' and w[1] == '0' and w[2] != '0':
        words.append(w)
 
# объединяем цифры в одну строку, а затем в число
words = list(map(''.join, words))
words = list(map(int, words))
 
 
### ПОМЕСТИТЬ В ЦИКЛ 2..81
 
# Выбираем двойки чисел и перебираем суммы двоек
a = permutations(words, r = 2)
a = list(a)
 
for i in range(len(a)):
    if a[i][0] + a[i][1] == s:
        print('2')
 
print('no')

print(27*804 + 904)
print(27*408 + 409)
print((22612 - 11425)%99 == 0)
 
print(901//109)
 
print(29997*(901//109))

def gcd_ext(a, b):
    if b == 0:
        return a, 1, 0
    d, x, y = gcd_ext(b, a % b)
    x, y = y, x - (a//b) * y
    return d, x, y
 
 
diff = [i*99 for i in range(8,0,-1)]
pi = [[900 + i - 101*j for j in range(i)] for i in range(1, 9)]
 
p = 22612
s = 11425
c = p - s
for i in range(7):
    a = diff[i]
    for j in range(i+1, 8):
        b = diff[j]
        d, x, y = gcd_ext(a, b)
 
        if c % d == 0:
            k = c//d
            t = k*x // (b//d)
            x = k*x - (b//d) * t
 
            while x <= 0 :
                x = k*x + b//d
            y = (c - a*x) // b
            if y > 0:
                for ei in pi[i]:
                    for ej in pi[j]:
                        if ei*x + ej*y == p:
                            print((ei, x), (ej, y), ei*x + ej*y)
                            ei = int(str(ei)[::-1])
                            ej = int(str(ej)[::-1])
                            print((ei, x), (ej, y), ei*x + ej*y)
                            print()