/-
Copyright 2026 The Formal Conjectures Authors.
Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at
https://www.apache.org/licenses/LICENSE-2.0
Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
-/
import FormalConjecturesUtilBounds on fractional parts $n r^2 - a(n)$ for $r = (2+\sqrt{5})/2$
$a(n) = \lfloor r \cdot \lfloor r \cdot n \rfloor \rfloor$, where $r = (2 + \sqrt{5})/2$.
References:
arxiv/2605.22763 Advancing Mathematics Research with AI-Driven Formal Proof Search by George Tsoukalas et al.
namespace OeisA341254open RealThe constant $r = (2 + \sqrt{5})/2$.
noncomputable def rConst : ℝ := (2 + sqrt 5) / 2The constant $r^2$.
noncomputable def rSq : ℝ := rConst * rConst$a(n) = \lfloor r \cdot \lfloor r \cdot n \rfloor \rfloor$, where $r = (2 + \sqrt{5})/2$. Note: The original OEIS definition has $n$ starting at 1. We define $a(n)$ for all $\mathbb{N}$.
noncomputable def a (n : ℕ) : ℕ :=
let r := rConst
let inner_floor : ℤ := Int.floor (r * n)
(Int.floor (r * inner_floor.cast)).toNath1:⌊(2 + √5) / 2⌋ = 2h2:⌊(2 + √5) / 2 * 2⌋ = 4⊢ Int.toNat 4 = 4
rfl All goals completed! 🐙
@[category test, AMS 11]
lemma a_2 : a 2 = 8 := by ⊢ a 2 = 8
unfold a rConst ⊢ (have r := (2 + √5) / 2;
have inner_floor := ⌊r * ↑2⌋;
⌊r * ↑inner_floor⌋.toNat) =
8
dsimp only ⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * ↑2⌋⌋.toNat = 8
push_cast ⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8
have h1 : ⌊(2 + Real.sqrt 5) / 2 * 2⌋ = 4 := by ⊢ a 2 = 8 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8
rw [Int.floor_eq_iff ⊢ ↑4 ≤ (2 + √5) / 2 * 2 ∧ (2 + √5) / 2 * 2 < ↑4 + 1 ⊢ ↑4 ≤ (2 + √5) / 2 * 2 ∧ (2 + √5) / 2 * 2 < ↑4 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8] ⊢ ↑4 ≤ (2 + √5) / 2 * 2 ∧ (2 + √5) / 2 * 2 < ↑4 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8
have h5 : 2236 / 1000 < Real.sqrt 5 ∧ Real.sqrt 5 < 2237 / 1000 := by ⊢ a 2 = 8 h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑4 ≤ (2 + √5) / 2 * 2 ∧ (2 + √5) / 2 * 2 < ↑4 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8
norm_num [Real.sqrt_lt, Real.lt_sqrt] h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑4 ≤ (2 + √5) / 2 * 2 ∧ (2 + √5) / 2 * 2 < ↑4 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8 h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑4 ≤ (2 + √5) / 2 * 2 ∧ (2 + √5) / 2 * 2 < ↑4 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8
constructor left h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑4 ≤ (2 + √5) / 2 * 2right h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 2 < ↑4 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8 <;> left h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑4 ≤ (2 + √5) / 2 * 2right h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 2 < ↑4 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8 linarith h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 2⌋⌋.toNat = 8
rw [h1 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑4⌋.toNat = 8 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑4⌋.toNat = 8] h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * ↑4⌋.toNat = 8
push_cast h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8
have h2 : ⌊(2 + Real.sqrt 5) / 2 * 4⌋ = 8 := by ⊢ a 2 = 8 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8
rw [Int.floor_eq_iff h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8] h1:⌊(2 + √5) / 2 * 2⌋ = 4⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8
have h5 : 2236 / 1000 < Real.sqrt 5 ∧ Real.sqrt 5 < 2237 / 1000 := by ⊢ a 2 = 8 h1:⌊(2 + √5) / 2 * 2⌋ = 4h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8
norm_num [Real.sqrt_lt, Real.lt_sqrt] h1:⌊(2 + √5) / 2 * 2⌋ = 4h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8 h1:⌊(2 + √5) / 2 * 2⌋ = 4h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8
constructor left h1:⌊(2 + √5) / 2 * 2⌋ = 4h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4right h1:⌊(2 + √5) / 2 * 2⌋ = 4h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8 <;> left h1:⌊(2 + √5) / 2 * 2⌋ = 4h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4right h1:⌊(2 + √5) / 2 * 2⌋ = 4h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8 linarith h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 4⌋.toNat = 8
rw [h2 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ Int.toNat 8 = 8 h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ Int.toNat 8 = 8] h1:⌊(2 + √5) / 2 * 2⌋ = 4h2:⌊(2 + √5) / 2 * 4⌋ = 8⊢ Int.toNat 8 = 8
rfl All goals completed! 🐙
@[category test, AMS 11]
lemma a_3 : a 3 = 12 := by ⊢ a 3 = 12
unfold a rConst ⊢ (have r := (2 + √5) / 2;
have inner_floor := ⌊r * ↑3⌋;
⌊r * ↑inner_floor⌋.toNat) =
12
dsimp only ⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * ↑3⌋⌋.toNat = 12
push_cast ⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12
have h1 : ⌊(2 + Real.sqrt 5) / 2 * 3⌋ = 6 := by ⊢ a 3 = 12 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12
rw [Int.floor_eq_iff ⊢ ↑6 ≤ (2 + √5) / 2 * 3 ∧ (2 + √5) / 2 * 3 < ↑6 + 1 ⊢ ↑6 ≤ (2 + √5) / 2 * 3 ∧ (2 + √5) / 2 * 3 < ↑6 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12] ⊢ ↑6 ≤ (2 + √5) / 2 * 3 ∧ (2 + √5) / 2 * 3 < ↑6 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12
have h5 : 2236 / 1000 < Real.sqrt 5 ∧ Real.sqrt 5 < 2237 / 1000 := by ⊢ a 3 = 12 h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑6 ≤ (2 + √5) / 2 * 3 ∧ (2 + √5) / 2 * 3 < ↑6 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12
norm_num [Real.sqrt_lt, Real.lt_sqrt] h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑6 ≤ (2 + √5) / 2 * 3 ∧ (2 + √5) / 2 * 3 < ↑6 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12 h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑6 ≤ (2 + √5) / 2 * 3 ∧ (2 + √5) / 2 * 3 < ↑6 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12
constructor left h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑6 ≤ (2 + √5) / 2 * 3right h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 3 < ↑6 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12 <;> left h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑6 ≤ (2 + √5) / 2 * 3right h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 3 < ↑6 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12 linarith h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 3⌋⌋.toNat = 12
rw [h1 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑6⌋.toNat = 12 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑6⌋.toNat = 12] h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * ↑6⌋.toNat = 12
push_cast h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12
have h2 : ⌊(2 + Real.sqrt 5) / 2 * 6⌋ = 12 := by ⊢ a 3 = 12 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12
rw [Int.floor_eq_iff h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ↑12 ≤ (2 + √5) / 2 * 6 ∧ (2 + √5) / 2 * 6 < ↑12 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ↑12 ≤ (2 + √5) / 2 * 6 ∧ (2 + √5) / 2 * 6 < ↑12 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12] h1:⌊(2 + √5) / 2 * 3⌋ = 6⊢ ↑12 ≤ (2 + √5) / 2 * 6 ∧ (2 + √5) / 2 * 6 < ↑12 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12
have h5 : 2236 / 1000 < Real.sqrt 5 ∧ Real.sqrt 5 < 2237 / 1000 := by ⊢ a 3 = 12 h1:⌊(2 + √5) / 2 * 3⌋ = 6h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑12 ≤ (2 + √5) / 2 * 6 ∧ (2 + √5) / 2 * 6 < ↑12 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12
norm_num [Real.sqrt_lt, Real.lt_sqrt] h1:⌊(2 + √5) / 2 * 3⌋ = 6h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑12 ≤ (2 + √5) / 2 * 6 ∧ (2 + √5) / 2 * 6 < ↑12 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12 h1:⌊(2 + √5) / 2 * 3⌋ = 6h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑12 ≤ (2 + √5) / 2 * 6 ∧ (2 + √5) / 2 * 6 < ↑12 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12
constructor left h1:⌊(2 + √5) / 2 * 3⌋ = 6h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑12 ≤ (2 + √5) / 2 * 6right h1:⌊(2 + √5) / 2 * 3⌋ = 6h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 6 < ↑12 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12 <;> left h1:⌊(2 + √5) / 2 * 3⌋ = 6h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑12 ≤ (2 + √5) / 2 * 6right h1:⌊(2 + √5) / 2 * 3⌋ = 6h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 6 < ↑12 + 1 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12 linarith h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ ⌊(2 + √5) / 2 * 6⌋.toNat = 12
rw [h2 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ Int.toNat 12 = 12 h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ Int.toNat 12 = 12] h1:⌊(2 + √5) / 2 * 3⌋ = 6h2:⌊(2 + √5) / 2 * 6⌋ = 12⊢ Int.toNat 12 = 12
rfl All goals completed! 🐙
@[category test, AMS 11]
lemma a_4 : a 4 = 16 := by ⊢ a 4 = 16
unfold a rConst ⊢ (have r := (2 + √5) / 2;
have inner_floor := ⌊r * ↑4⌋;
⌊r * ↑inner_floor⌋.toNat) =
16
dsimp only ⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * ↑4⌋⌋.toNat = 16
push_cast ⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16
have h1 : ⌊(2 + Real.sqrt 5) / 2 * 4⌋ = 8 := by ⊢ a 4 = 16 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16
rw [Int.floor_eq_iff ⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 ⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16] ⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16
have h5 : 2236 / 1000 < Real.sqrt 5 ∧ Real.sqrt 5 < 2237 / 1000 := by ⊢ a 4 = 16 h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16
norm_num [Real.sqrt_lt, Real.lt_sqrt] h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16 h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4 ∧ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16
constructor left h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4right h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16 <;> left h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑8 ≤ (2 + √5) / 2 * 4right h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 4 < ↑8 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16 linarith h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 4⌋⌋.toNat = 16
rw [h1 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑8⌋.toNat = 16 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑8⌋.toNat = 16] h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * ↑8⌋.toNat = 16
push_cast h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16
have h2 : ⌊(2 + Real.sqrt 5) / 2 * 8⌋ = 16 := by ⊢ a 4 = 16 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16
rw [Int.floor_eq_iff h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ↑16 ≤ (2 + √5) / 2 * 8 ∧ (2 + √5) / 2 * 8 < ↑16 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ↑16 ≤ (2 + √5) / 2 * 8 ∧ (2 + √5) / 2 * 8 < ↑16 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16] h1:⌊(2 + √5) / 2 * 4⌋ = 8⊢ ↑16 ≤ (2 + √5) / 2 * 8 ∧ (2 + √5) / 2 * 8 < ↑16 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16
have h5 : 2236 / 1000 < Real.sqrt 5 ∧ Real.sqrt 5 < 2237 / 1000 := by ⊢ a 4 = 16 h1:⌊(2 + √5) / 2 * 4⌋ = 8h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑16 ≤ (2 + √5) / 2 * 8 ∧ (2 + √5) / 2 * 8 < ↑16 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16
norm_num [Real.sqrt_lt, Real.lt_sqrt] h1:⌊(2 + √5) / 2 * 4⌋ = 8h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑16 ≤ (2 + √5) / 2 * 8 ∧ (2 + √5) / 2 * 8 < ↑16 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16 h1:⌊(2 + √5) / 2 * 4⌋ = 8h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑16 ≤ (2 + √5) / 2 * 8 ∧ (2 + √5) / 2 * 8 < ↑16 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16
constructor left h1:⌊(2 + √5) / 2 * 4⌋ = 8h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑16 ≤ (2 + √5) / 2 * 8right h1:⌊(2 + √5) / 2 * 4⌋ = 8h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 8 < ↑16 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16 <;> left h1:⌊(2 + √5) / 2 * 4⌋ = 8h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑16 ≤ (2 + √5) / 2 * 8right h1:⌊(2 + √5) / 2 * 4⌋ = 8h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 8 < ↑16 + 1 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16 linarith h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ ⌊(2 + √5) / 2 * 8⌋.toNat = 16
rw [h2 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ Int.toNat 16 = 16 h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ Int.toNat 16 = 16] h1:⌊(2 + √5) / 2 * 4⌋ = 8h2:⌊(2 + √5) / 2 * 8⌋ = 16⊢ Int.toNat 16 = 16
rfl All goals completed! 🐙
@[category test, AMS 11]
lemma a_5 : a 5 = 21 := by ⊢ a 5 = 21
unfold a rConst ⊢ (have r := (2 + √5) / 2;
have inner_floor := ⌊r * ↑5⌋;
⌊r * ↑inner_floor⌋.toNat) =
21
dsimp only ⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * ↑5⌋⌋.toNat = 21
push_cast ⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21
have h1 : ⌊(2 + Real.sqrt 5) / 2 * 5⌋ = 10 := by ⊢ a 5 = 21 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21
rw [Int.floor_eq_iff ⊢ ↑10 ≤ (2 + √5) / 2 * 5 ∧ (2 + √5) / 2 * 5 < ↑10 + 1 ⊢ ↑10 ≤ (2 + √5) / 2 * 5 ∧ (2 + √5) / 2 * 5 < ↑10 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21] ⊢ ↑10 ≤ (2 + √5) / 2 * 5 ∧ (2 + √5) / 2 * 5 < ↑10 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21
have h5 : 2236 / 1000 < Real.sqrt 5 ∧ Real.sqrt 5 < 2237 / 1000 := by ⊢ a 5 = 21 h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑10 ≤ (2 + √5) / 2 * 5 ∧ (2 + √5) / 2 * 5 < ↑10 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21
norm_num [Real.sqrt_lt, Real.lt_sqrt] h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑10 ≤ (2 + √5) / 2 * 5 ∧ (2 + √5) / 2 * 5 < ↑10 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21 h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑10 ≤ (2 + √5) / 2 * 5 ∧ (2 + √5) / 2 * 5 < ↑10 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21
constructor left h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑10 ≤ (2 + √5) / 2 * 5right h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 5 < ↑10 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21 <;> left h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑10 ≤ (2 + √5) / 2 * 5right h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 5 < ↑10 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21 linarith h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑⌊(2 + √5) / 2 * 5⌋⌋.toNat = 21
rw [h1 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑10⌋.toNat = 21 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑10⌋.toNat = 21] h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * ↑10⌋.toNat = 21
push_cast h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21
have h2 : ⌊(2 + Real.sqrt 5) / 2 * 10⌋ = 21 := by ⊢ a 5 = 21 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21
rw [Int.floor_eq_iff h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ↑21 ≤ (2 + √5) / 2 * 10 ∧ (2 + √5) / 2 * 10 < ↑21 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ↑21 ≤ (2 + √5) / 2 * 10 ∧ (2 + √5) / 2 * 10 < ↑21 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21] h1:⌊(2 + √5) / 2 * 5⌋ = 10⊢ ↑21 ≤ (2 + √5) / 2 * 10 ∧ (2 + √5) / 2 * 10 < ↑21 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21
have h5 : 2236 / 1000 < Real.sqrt 5 ∧ Real.sqrt 5 < 2237 / 1000 := by ⊢ a 5 = 21 h1:⌊(2 + √5) / 2 * 5⌋ = 10h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑21 ≤ (2 + √5) / 2 * 10 ∧ (2 + √5) / 2 * 10 < ↑21 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21
norm_num [Real.sqrt_lt, Real.lt_sqrt] h1:⌊(2 + √5) / 2 * 5⌋ = 10h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑21 ≤ (2 + √5) / 2 * 10 ∧ (2 + √5) / 2 * 10 < ↑21 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21 h1:⌊(2 + √5) / 2 * 5⌋ = 10h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑21 ≤ (2 + √5) / 2 * 10 ∧ (2 + √5) / 2 * 10 < ↑21 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21
constructor left h1:⌊(2 + √5) / 2 * 5⌋ = 10h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑21 ≤ (2 + √5) / 2 * 10right h1:⌊(2 + √5) / 2 * 5⌋ = 10h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 10 < ↑21 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21 <;> left h1:⌊(2 + √5) / 2 * 5⌋ = 10h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ ↑21 ≤ (2 + √5) / 2 * 10right h1:⌊(2 + √5) / 2 * 5⌋ = 10h5:2236 / 1000 < √5 ∧ √5 < 2237 / 1000⊢ (2 + √5) / 2 * 10 < ↑21 + 1 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21 linarith h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ ⌊(2 + √5) / 2 * 10⌋.toNat = 21
rw [h2 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ Int.toNat 21 = 21 h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ Int.toNat 21 = 21] h1:⌊(2 + √5) / 2 * 5⌋ = 10h2:⌊(2 + √5) / 2 * 10⌋ = 21⊢ Int.toNat 21 = 21
rfl All goals completed! 🐙Conjecture: $1/4 < n r^2 - a(n) < 3$ for $n \ge 1$, where $r = (2 + \sqrt{5})/2$.
A formal proof has been found with the methods described in arxiv/2605.22763.
@[category research solved, AMS 11, formal_proof using formal_conjectures at
"https://github.com/mo271/formal-conjectures/blob/a32396489dcb8f86c3549b93aa358ac6a10a3a1f/FormalConjectures/OEIS/341254.wip.lean#L188"]
theorem a_bounds (n : ℕ) (hn : 1 ≤ n) : (1 / 4 : ℝ) < (n : ℝ) * rSq - (a n : ℝ) ∧
(n : ℝ) * rSq - (a n : ℝ) < 3 := by n:ℕhn:1 ≤ n⊢ 1 / 4 < ↑n * rSq - ↑(a n) ∧ ↑n * rSq - ↑(a n) < 3
sorry All goals completed! 🐙end OeisA341254