j1 + j2 = 2 * i0 + 1, 2 * i = j1 ^2 + j2 ^2 ⊦ (j1 - j2) % 2 = 0, i = ((j1 - j2) / 2)^2 + ((j1 + j2) / 2)^2