| abs(i) = i0 * abs(j) + j0, gcd(abs(i), abs(j)) * i1 = abs(j), gcd(abs(i), abs(j)) * i2 = abs(i), j0 < abs(j), 0 ≤ j0, 0 ≤ i0, gcd(abs(j), j0) < gcd(abs(i), abs(j)) ⊦ j = 0, gcd(abs(i), abs(j)) | gcd(abs(i), abs(j)) * i2 - i0 * abs(j) |
abs(i) = i0 * abs(j) + j0, gcd(abs(i), abs(j)) * i1 = abs(j), gcd(abs(i), abs(j)) * i2 = abs(i), j0 < abs(j), 0 ≤ j0, 0 ≤ i0, gcd(abs(j), j0) < gcd(abs(i), abs(j)) ⊦ j = 0, gcd(abs(i), abs(j)) | gcd(abs(i), abs(j)) * i2 - i0 * gcd(abs(i), abs(j)) * i1