1 = i * i0 + j * j1, gcd(i, j) = 1, i | j * (j0 * i * i0 + j0 * j * j1) ⊦ i | j0 * i * i0 + j0 * j * j1