import Mathlib
import Optlib


open Set

theorem extracted_formal_statement_98.{u_1} {E : Type u_1} [inst : NormedAddCommGroup E] [inst_1 : InnerProductSpace ℝ E]
  [inst_2 : CompleteSpace E] {f : E → ℝ} {m : ℝ} {f' : E → E} {xm x₀ : E} {alg : Gradient_Descent_fix_stepsize f f' x₀}
  (hm : m > 0) (min : IsMinOn f univ xm)
  (step₂ : Gradient_Descent_fix_stepsize.a f f' x₀ ≤ 2 / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))
  (hsc : StrongConvexOn univ m f) (this_1 : LipschitzWith (Gradient_Descent_fix_stepsize.l f f' x₀) f')
  (this : ↑(Gradient_Descent_fix_stepsize.l f f' x₀) > 0)
  (reduction :
    ∀ (k : ℕ),
      ‖Gradient_Descent_fix_stepsize.x f f' x₀ (k + 1) - xm‖ ^ 2 ≤
        (1 -
            Gradient_Descent_fix_stepsize.a f f' x₀ *
              (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
          ‖Gradient_Descent_fix_stepsize.x f f' x₀ k - xm‖ ^ 2)
  (eq :
    0 ≤
      1 -
        Gradient_Descent_fix_stepsize.a f f' x₀ *
          (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀))))
  (q : ℕ)
  (IH1 :
    ‖Gradient_Descent_fix_stepsize.x f f' x₀ q - xm‖ ^ 2 ≤
      (1 -
            Gradient_Descent_fix_stepsize.a f f' x₀ *
              (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) ^
          q *
        ‖x₀ - xm‖ ^ 2)
  (h_original_goal :
    (1 -
          Gradient_Descent_fix_stepsize.a f f' x₀ *
            (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
        ‖Gradient_Descent_fix_stepsize.x f f' x₀ q - xm‖ ^ 2 ≤
      (1 -
            Gradient_Descent_fix_stepsize.a f f' x₀ *
              (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
          (1 -
              Gradient_Descent_fix_stepsize.a f f' x₀ *
                (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) /
                  (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) ^
            q *
        ‖x₀ - xm‖ ^ 2) :
  (1 -
        Gradient_Descent_fix_stepsize.a f f' x₀ *
          (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
      ‖Gradient_Descent_fix_stepsize.x f f' x₀ q - xm‖ ^ 2 ≤
    (1 -
        Gradient_Descent_fix_stepsize.a f f' x₀ *
          (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
      ((1 -
            Gradient_Descent_fix_stepsize.a f f' x₀ *
              (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) ^
          q *
        ‖x₀ - xm‖ ^ 2) := sorry

Tactic state before:
E : Type u_1
inst✝² : NormedAddCommGroup E
inst✝¹ : InnerProductSpace ℝ E
inst✝ : CompleteSpace E
f : E → ℝ
m : ℝ
f' : E → E
xm x₀ : E
alg : Gradient_Descent_fix_stepsize f f' x₀
hm : m > 0
min : IsMinOn f univ xm
step₂ : Gradient_Descent_fix_stepsize.a f f' x₀ ≤ 2 / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀))
hsc : StrongConvexOn univ m f
this✝ : LipschitzWith (Gradient_Descent_fix_stepsize.l f f' x₀) f'
this : ↑(Gradient_Descent_fix_stepsize.l f f' x₀) > 0
reduction :
  ∀ (k : ℕ),
    ‖Gradient_Descent_fix_stepsize.x f f' x₀ (k + 1) - xm‖ ^ 2 ≤
      (1 -
          Gradient_Descent_fix_stepsize.a f f' x₀ *
            (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
        ‖Gradient_Descent_fix_stepsize.x f f' x₀ k - xm‖ ^ 2
eq :
  0 ≤
    1 -
      Gradient_Descent_fix_stepsize.a f f' x₀ *
        (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))
q : ℕ
IH1 :
  ‖Gradient_Descent_fix_stepsize.x f f' x₀ q - xm‖ ^ 2 ≤
    (1 -
          Gradient_Descent_fix_stepsize.a f f' x₀ *
            (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) ^
        q *
      ‖x₀ - xm‖ ^ 2
⊢ (1 -
        Gradient_Descent_fix_stepsize.a f f' x₀ *
          (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
      ‖Gradient_Descent_fix_stepsize.x f f' x₀ q - xm‖ ^ 2 ≤
    (1 -
          Gradient_Descent_fix_stepsize.a f f' x₀ *
            (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
        (1 -
            Gradient_Descent_fix_stepsize.a f f' x₀ *
              (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) ^
          q *
      ‖x₀ - xm‖ ^ 2

Tactic:
rw [mul_assoc _ _ (‖x₀ - xm‖ ^ 2)]

Tactic state after:
[E : Type u_1
 inst✝² : NormedAddCommGroup E
 inst✝¹ : InnerProductSpace ℝ E
 inst✝ : CompleteSpace E
 f : E → ℝ
 m : ℝ
 f' : E → E
 xm x₀ : E
 alg : Gradient_Descent_fix_stepsize f f' x₀
 hm : m > 0
 min : IsMinOn f univ xm
 step₂ : Gradient_Descent_fix_stepsize.a f f' x₀ ≤ 2 / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀))
 hsc : StrongConvexOn univ m f
 this✝ : LipschitzWith (Gradient_Descent_fix_stepsize.l f f' x₀) f'
 this : ↑(Gradient_Descent_fix_stepsize.l f f' x₀) > 0
 reduction :
   ∀ (k : ℕ),
     ‖Gradient_Descent_fix_stepsize.x f f' x₀ (k + 1) - xm‖ ^ 2 ≤
       (1 -
           Gradient_Descent_fix_stepsize.a f f' x₀ *
             (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
         ‖Gradient_Descent_fix_stepsize.x f f' x₀ k - xm‖ ^ 2
 eq :
   0 ≤
     1 -
       Gradient_Descent_fix_stepsize.a f f' x₀ *
         (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))
 q : ℕ
 IH1 :
   ‖Gradient_Descent_fix_stepsize.x f f' x₀ q - xm‖ ^ 2 ≤
     (1 -
           Gradient_Descent_fix_stepsize.a f f' x₀ *
             (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) ^
         q *
       ‖x₀ - xm‖ ^ 2
 ⊢ (1 -
         Gradient_Descent_fix_stepsize.a f f' x₀ *
           (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
       ‖Gradient_Descent_fix_stepsize.x f f' x₀ q - xm‖ ^ 2 ≤
     (1 -
         Gradient_Descent_fix_stepsize.a f f' x₀ *
           (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) / (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) *
       ((1 -
             Gradient_Descent_fix_stepsize.a f f' x₀ *
               (2 * m * ↑(Gradient_Descent_fix_stepsize.l f f' x₀) /
                 (m + ↑(Gradient_Descent_fix_stepsize.l f f' x₀)))) ^
           q *
         ‖x₀ - xm‖ ^ 2)]