| i0 * j0 = max_prime_divisor(j0) * i0 * i0, i0 * j0 = i * i * i1, max_prime_divisor(j0) * i0 = j0, prime(max_prime_divisor(j0)), 0 < i, 0 < j0, 0 < max_prime_divisor(j0) * i1, ∀ j1, j, i. i→n(j1) < i→n(j0) → j1 * j1 = j * i * i ∧ 0 < j1 ∧ 0 < j ∧ 0 < i → i = 0 ∨ isSquare(j) ⊦ i = 0, max_prime_divisor(j0) = 0, max_prime_divisor(j0) | i, isSquare(max_prime_divisor(j0) * i1), prime(j0) |
max_prime_divisor(j0) * i0 * i0 = i * i * i1, max_prime_divisor(j0) * i0 = j0, prime(max_prime_divisor(j0)), 0 < i, 0 < j0, 0 < max_prime_divisor(j0) * i1, ∀ j1, j, i. i→n(j1) < i→n(j0) → j1 * j1 = j * i * i ∧ 0 < j1 ∧ 0 < j ∧ 0 < i → i = 0 ∨ isSquare(j) ⊦ i = 0, max_prime_divisor(j0) = 0, max_prime_divisor(j0) | i, isSquare(max_prime_divisor(j0) * i1), prime(j0)