-- Lean 4 proof step from Mathlib.Algebra.Polynomial.PartialFractions (goal)
Goal:
R : Type u_1
inst : CommRing R
K : Type u_2
inst_1 : CommRing K
inst_2 : Algebra.{u_1, u_2} R[X] K
inst_3 : Nontrivial R
ι : Type u_3
s : Finset ι
f : R[X]
g : ι → R[X]
hg : ∀ i ∈ s, (g i).Monic
hgg : (↑s : Set ι).Pairwise fun (i j : ι) ↦ IsCoprime (g i) (g j)
n : ι → ℕ
gi : ι → K
hgi : ∀ i ∈ s, HMul.hMul (β := K) (gi i) ((algebraMap R[X] K : (a : R[X]) → K) (g i) : K) = 1
obtain :
  ∀ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]) (hr : ∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree)
    (hf : f = q * ∏ i ∈ s, g i ^ n i + ∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k),
    ∃ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]),
      (∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree) ∧
        HMul.hMul (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) f : K) (∏ i ∈ s, gi i ^ n i : K) =
          HAdd.hAdd (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) q : K)
            (∑ i ∈ s,
                ∑ j : Fin (n i),
                  HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) (r i j) : K) (gi i ^ ((↑j : ℕ) + 1)) :
              K)
obtain₁ :
  ∀ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]) (hr : ∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree)
    (hf : f = q * ∏ i ∈ s, g i ^ n i + ∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k),
    HMul.hMul (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) f : K) (∏ i ∈ s, gi i ^ n i : K) =
      HAdd.hAdd (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) q : K)
        (∑ i ∈ s,
            ∑ j : Fin (n i),
              HMul.hMul (α := K)
                ((algebraMap R[X] K : (a : R[X]) → K) ((fun (i : ι) (j : Fin (n i)) ↦ r i j.rev) i j) : K)
                (gi i ^ ((↑j : ℕ) + 1)) :
          K)
obtain₂ :
  ∀ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]) (hr : ∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree)
    (hf : f = q * ∏ i ∈ s, g i ^ n i + ∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k),
    HMul.hMul (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) q : K)
          (∏ x ∈ s, HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) (g x ^ n x) : K) (gi x ^ n x) : K) +
        HMul.hMul (α := K) (β := K)
          ((algebraMap R[X] K : (a : R[X]) → K)
              (∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k) :
            K)
          (∏ i ∈ s, gi i ^ n i : K) =
      HAdd.hAdd (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) q : K)
        (∑ i ∈ s,
            ∑ j : Fin (n i),
              HMul.hMul (α := K)
                ((algebraMap R[X] K : (a : R[X]) → K) ((fun (i : ι) (j : Fin (n i)) ↦ r i j.rev) i j) : K)
                (gi i ^ ((↑j : ℕ) + 1)) :
          K)
obtain₃ :
  ∀ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]) (hr : ∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree)
    (hf : f = q * ∏ i ∈ s, g i ^ n i + ∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k)
    (hc :
      ∀ (x : ι) (hx : x ∈ s),
        HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) (g x ^ n x) : K) (gi x ^ n x) = 1),
    HMul.hMul (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) q : K)
          (∏ x ∈ s, HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) (g x ^ n x) : K) (gi x ^ n x) : K) +
        HMul.hMul (α := K) (β := K)
          ((algebraMap R[X] K : (a : R[X]) → K)
              (∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k) :
            K)
          (∏ i ∈ s, gi i ^ n i : K) =
      HAdd.hAdd (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) q : K)
        (∑ i ∈ s,
            ∑ j : Fin (n i),
              HMul.hMul (α := K)
                ((algebraMap R[X] K : (a : R[X]) → K) ((fun (i : ι) (j : Fin (n i)) ↦ r i j.rev) i j) : K)
                (gi i ^ ((↑j : ℕ) + 1)) :
          K)
obtain₄ :
  ∀ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]) (hr : ∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree)
    (hf : f = q * ∏ i ∈ s, g i ^ n i + ∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k)
    (hc :
      ∀ (x : ι) (hx : x ∈ s),
        HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) (g x ^ n x) : K) (gi x ^ n x) = 1),
    Eq.{u_2 + 1} (α := K)
      (∑ i ∈ s,
        HMul.hMul (α := K) (β := K)
          ((algebraMap R[X] K : (a : R[X]) → K) (∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k) :
            K)
          (∏ i ∈ s, gi i ^ n i : K))
      (∑ i ∈ s,
        ∑ j : Fin (n i),
          HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) ((fun (i : ι) (j : Fin (n i)) ↦ r i j.rev) i j) : K)
            (gi i ^ ((↑j : ℕ) + 1)))
obtain₅ :
  ∀ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]) (hr : ∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree)
    (hf : f = q * ∏ i ∈ s, g i ^ n i + ∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k)
    (hc :
      ∀ (x : ι) (hx : x ∈ s),
        HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) (g x ^ n x) : K) (gi x ^ n x) = 1)
    (i : ι) (hi : i ∈ s),
    HMul.hMul (α := K) (β := K)
        ((algebraMap R[X] K : (a : R[X]) → K) (∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k) :
          K)
        (∏ i ∈ s, gi i ^ n i : K) =
      ∑ j : Fin (n i),
        HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) ((fun (i : ι) (j : Fin (n i)) ↦ r i j.rev) i j) : K)
          (gi i ^ ((↑j : ℕ) + 1))
obtain₆ :
  ∀ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]) (hr : ∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree)
    (hf : f = q * ∏ i ∈ s, g i ^ n i + ∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k)
    (hc :
      ∀ (x : ι) (hx : x ∈ s),
        HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) (g x ^ n x) : K) (gi x ^ n x) = 1)
    (i : ι) (hi : i ∈ s),
    Eq.{u_2 + 1} (α := K)
      (∑ i_1 : Fin (n i),
        HMul.hMul (α := K) (β := K)
          ((algebraMap R[X] K : (a : R[X]) → K)
              (r i ((Fin.revPerm : (a : Fin (n i)) → Fin (n i)) i_1) *
                  g i ^ (↑((Fin.revPerm : (a : Fin (n i)) → Fin (n i)) i_1) : ℕ) *
                ∏ k ∈ s.erase i, g k ^ n k) :
            K)
          (∏ i ∈ s, gi i ^ n i : K))
      (∑ j : Fin (n i),
        HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) ((fun (i : ι) (j : Fin (n i)) ↦ r i j.rev) i j) : K)
          (gi i ^ ((↑j : ℕ) + 1)))
⊢ ∃ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]),
    (∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree) ∧
      HMul.hMul (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) f : K) (∏ i ∈ s, gi i ^ n i : K) =
        HAdd.hAdd (α := K) (β := K) ((algebraMap R[X] K : (a : R[X]) → K) q : K)
          (∑ i ∈ s,
              ∑ j : Fin (n i),
                HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) (r i j) : K) (gi i ^ ((↑j : ℕ) + 1)) :
            K)
have obtain₇ : ∀ (q : R[X]) (r : (i : ι) → Fin (n i) → R[X]) (hr : ∀ i ∈ s, ∀ (j : Fin (n i)), (r i j).degree < (g i).degree)
  (hf : f = q * ∏ i ∈ s, g i ^ n i + ∑ i ∈ s, ∑ j : Fin (n i), r i j * g i ^ (↑j : ℕ) * ∏ k ∈ s.erase i, g k ^ n k)
  (hc :
    ∀ (x : ι) (hx : x ∈ s), HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) (g x ^ n x) : K) (gi x ^ n x) = 1)
  (i : ι) (hi : i ∈ s) (j : Fin (n i)),
  HMul.hMul (α := K) (β := K)
      ((algebraMap R[X] K : (a : R[X]) → K)
          (r i ((Fin.revPerm : (a : Fin (n i)) → Fin (n i)) j) *
              g i ^ (↑((Fin.revPerm : (a : Fin (n i)) → Fin (n i)) j) : ℕ) *
            ∏ k ∈ s.erase i, g k ^ n k) :
        K)
      (∏ i ∈ s, gi i ^ n i : K) =
    HMul.hMul (α := K) ((algebraMap R[X] K : (a : R[X]) → K) ((fun (i : ι) (j : Fin (n i)) ↦ r i j.rev) i j) : K)
      (gi i ^ ((↑j : ℕ) + 1))
Tactic: classical
  obtain ⟨q, r, hr, hf⟩ := eq_quo_mul_prod_pow_add_sum_rem_mul_prod_pow f hg hgg n
  refine ⟨q, fun i j => r i j.rev, fun i hi j => hr i hi j.rev, ?_⟩
  rw [hf, map_add, map_mul, map_prod, add_mul, mul_assoc, ← Finset.prod_mul_distrib]
  have hc (x : ι) (hx : x ∈ s) : (algebraMap R[X] K) (g x ^ n x) * gi x ^ n x = 1 := by
    rw [map_pow, ← mul_pow, mul_comm, hgi x hx, one_pow]
  rw [Finset.prod_congr rfl hc, Finset.prod_const_one,
    mul_one, add_right_inj, map_sum, Finset.sum_mul]
  refine Finset.sum_congr rfl fun i hi => ?_
  rw [map_sum, Finset.sum_mul, ← Equiv.sum_comp Fin.revPerm]
  refine Fintype.sum_congr _ _ fun j => ?_
  rw [Fin.revPerm_apply, map_mul, map_prod, ← Finset.prod_erase_mul s _ hi,
    ← mul_rotate', mul_assoc, ← Finset.prod_mul_distrib,
    Finset.prod_congr rfl fun x hx => hc x (Finset.mem_of_mem_erase hx),
    Finset.prod_const_one, mul_one, map_mul, map_pow, mul_left_comm]
  refine congrArg (_ * ·) ?_
  rw [← mul_one (gi i ^ (j.1 + 1)), ← @one_pow K _ j.rev, ← hgi i hi,
    mul_pow, ← mul_assoc, ← pow_add, Fin.val_rev, Nat.add_sub_cancel' (by lia)]