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