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