/-
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 FormalConjecturesUtilFactor bounds for Fermat numbers
$a(n)$ is the minimum integer $k$ such that the smallest prime factor of the $n$-th Fermat number exceeds $2^{2^n - k}$.
References:
namespace OeisA358684open NatA358684: $a(n)$ is the minimum integer $k$ such that the smallest prime factor of the $n$-th Fermat number exceeds $2^{2^n - k}$. Let $F_n = 2^{2^n} + 1$ be the $n$-th Fermat number, and $P_n$ be its smallest prime factor. The definition of $a(n)$ is equivalent to the closed form: $$a(n) = 2^n - \lfloor \log_2(P_n) \rfloor$$ where $P_n = \operatorname{minFac}(\operatorname{fermatNumber} n)$. The subtraction is defined in $\mathbb{N}$ and is safe since $P_n \le F_n$, implying $\log_2 P_n < 2^n$.
def a (n : ℕ) : ℕ :=
letI pn := minFac (fermatNumber n)
(2 ^ n) - (log2 pn)
The "original" definition: $a'(n)$ is the minimum $k$ such that $P_n > 2^{2^n - k}$.
We use Nat.find which returns the smallest natural number satisfying a predicate.
noncomputable def a' (n : ℕ) : ℕ :=
letI Pn := minFac (fermatNumber n)
Nat.find (show ∃ k, Pn > 2 ^ (2 ^ n - k) from n:ℕPn:ℕ := n.fermatNumber.minFac⊢ ∃ k, Pn > 2 ^ (2 ^ n - k)
n:ℕPn:ℕ := n.fermatNumber.minFac⊢ Pn > 2 ^ (2 ^ n - 2 ^ n)
n:ℕPn:ℕ := n.fermatNumber.minFac⊢ Pn > 1
n:ℕPn:ℕ := n.fermatNumber.minFac⊢ 1 < n.fermatNumber.minFac
n:ℕPn:ℕ := n.fermatNumber.minFac⊢ 1 < (2 ^ 2 ^ n + 1).minFac
exact (Nat.minFac_prime (n:ℕPn:ℕ := n.fermatNumber.minFac⊢ 2 ^ 2 ^ n + 1 ≠ 1 All goals completed! 🐙)).one_lt
)The log2 of the smallest prime factor of $F_n$ is at most $2^n$.
n:ℕ⊢ 2 ^ 2 ^ n + 1 < 2 ^ 2 ^ n + 2 ^ 2 ^ n
gcongr bc n:ℕ⊢ 1 < 2 ^ 2 ^ n
exact one_lt_pow (by n:ℕ⊢ 2 ^ n ≠ 0 norm_num All goals completed! 🐙) (by n:ℕ⊢ 1 < 2 norm_num All goals completed! 🐙)The minimization definition is equivalent to the closed form.
@[category API, AMS 11]
theorem a_equiv_a' (n : ℕ) : a n = a' n := by n:ℕ⊢ a n = a' n
unfold a a' n:ℕ⊢ 2 ^ n - n.fermatNumber.minFac.log2 = Nat.find ⋯; set Pn := (fermatNumber n).minFac n:ℕPn:ℕ := n.fermatNumber.minFac⊢ 2 ^ n - n.fermatNumber.minFac.log2 = Nat.find ⋯
rw [Eq.comm, n:ℕPn:ℕ := n.fermatNumber.minFac⊢ Nat.find ⋯ = 2 ^ n - n.fermatNumber.minFac.log2 n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ (2 ^ n - (2 ^ n - n.fermatNumber.minFac.log2)) ∧
∀ n_1 < 2 ^ n - n.fermatNumber.minFac.log2, ¬n.fermatNumber.minFac > 2 ^ (2 ^ n - n_1) Nat.find_eq_iff n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ (2 ^ n - (2 ^ n - n.fermatNumber.minFac.log2)) ∧
∀ n_1 < 2 ^ n - n.fermatNumber.minFac.log2, ¬n.fermatNumber.minFac > 2 ^ (2 ^ n - n_1) n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ (2 ^ n - (2 ^ n - n.fermatNumber.minFac.log2)) ∧
∀ n_1 < 2 ^ n - n.fermatNumber.minFac.log2, ¬n.fermatNumber.minFac > 2 ^ (2 ^ n - n_1)] n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ (2 ^ n - (2 ^ n - n.fermatNumber.minFac.log2)) ∧
∀ n_1 < 2 ^ n - n.fermatNumber.minFac.log2, ¬n.fermatNumber.minFac > 2 ^ (2 ^ n - n_1)
constructor left n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ (2 ^ n - (2 ^ n - n.fermatNumber.minFac.log2))right n:ℕPn:ℕ := n.fermatNumber.minFac⊢ ∀ n_1 < 2 ^ n - n.fermatNumber.minFac.log2, ¬n.fermatNumber.minFac > 2 ^ (2 ^ n - n_1)
· left n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ (2 ^ n - (2 ^ n - n.fermatNumber.minFac.log2)) rw [tsub_tsub_cancel_of_le (log2_minFac_le n), left n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ n.fermatNumber.minFac.log2 left n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ log 2 n.fermatNumber.minFac log2_eq_log_two left n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ log 2 n.fermatNumber.minFacleft n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ log 2 n.fermatNumber.minFac]left n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber.minFac > 2 ^ log 2 n.fermatNumber.minFac
refine lt_of_le_of_ne (pow_log_le_self 2 (Nat.Prime.ne_zero (minFac_prime ?_))) fun h => ?_ left.refine_1 n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber ≠ 1left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFac⊢ False
· left.refine_1 n:ℕPn:ℕ := n.fermatNumber.minFac⊢ n.fermatNumber ≠ 1 rw [fermatNumber left.refine_1 n:ℕPn:ℕ := n.fermatNumber.minFac⊢ 2 ^ 2 ^ n + 1 ≠ 1 left.refine_1 n:ℕPn:ℕ := n.fermatNumber.minFac⊢ 2 ^ 2 ^ n + 1 ≠ 1]left.refine_1 n:ℕPn:ℕ := n.fermatNumber.minFac⊢ 2 ^ 2 ^ n + 1 ≠ 1
norm_num All goals completed! 🐙
· left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFac⊢ False have hPn : Pn.Prime := minFac_prime (by n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFac⊢ n.fermatNumber ≠ 1 left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pn⊢ False rw [fermatNumber n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFac⊢ 2 ^ 2 ^ n + 1 ≠ 1 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFac⊢ 2 ^ 2 ^ n + 1 ≠ 1left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pn⊢ False] n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFac⊢ 2 ^ 2 ^ n + 1 ≠ 1left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pn⊢ False; norm_num All goals completed! 🐙left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pn⊢ False)left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pn⊢ False
have hnz := log2_eq_log_two ▸ (Nat.log_pos one_lt_two hPn.two_le).ne' left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pnhnz:Pn.log2 ≠ 0⊢ False
refine (Nat.not_even_iff_odd.mpr <| (odd_pow_iff hnz).mp ?_) even_two left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pnhnz:Pn.log2 ≠ 0⊢ Odd (2 ^ Pn.log2)
rw [log2_eq_log_two, left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pnhnz:Pn.log2 ≠ 0⊢ Odd (2 ^ log 2 Pn) left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pnhnz:Pn.log2 ≠ 0⊢ Odd n.fermatNumber.minFac h left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pnhnz:Pn.log2 ≠ 0⊢ Odd n.fermatNumber.minFacleft.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pnhnz:Pn.log2 ≠ 0⊢ Odd n.fermatNumber.minFac]left.refine_2 n:ℕPn:ℕ := n.fermatNumber.minFach:2 ^ log 2 n.fermatNumber.minFac = n.fermatNumber.minFachPn:Nat.Prime Pnhnz:Pn.log2 ≠ 0⊢ Odd n.fermatNumber.minFac
exact Odd.of_dvd_nat (odd_fermatNumber n) (minFac_dvd _) All goals completed! 🐙
· right n:ℕPn:ℕ := n.fermatNumber.minFac⊢ ∀ n_1 < 2 ^ n - n.fermatNumber.minFac.log2, ¬n.fermatNumber.minFac > 2 ^ (2 ^ n - n_1) intro m hm right n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:m < 2 ^ n - n.fermatNumber.minFac.log2⊢ ¬n.fermatNumber.minFac > 2 ^ (2 ^ n - m); simp only [not_lt] right n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:m < 2 ^ n - n.fermatNumber.minFac.log2⊢ n.fermatNumber.minFac ≤ 2 ^ (2 ^ n - m)
rw [lt_tsub_iff_left right n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:n.fermatNumber.minFac.log2 + m < 2 ^ n⊢ n.fermatNumber.minFac ≤ 2 ^ (2 ^ n - m) right n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:n.fermatNumber.minFac.log2 + m < 2 ^ n⊢ n.fermatNumber.minFac ≤ 2 ^ (2 ^ n - m)] at hmright n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:n.fermatNumber.minFac.log2 + m < 2 ^ n⊢ n.fermatNumber.minFac ≤ 2 ^ (2 ^ n - m)
rw [log2_eq_log_two right n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:log 2 n.fermatNumber.minFac + m < 2 ^ n⊢ n.fermatNumber.minFac ≤ 2 ^ (2 ^ n - m) right n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:log 2 n.fermatNumber.minFac + m < 2 ^ n⊢ n.fermatNumber.minFac ≤ 2 ^ (2 ^ n - m)] at hmright n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:log 2 n.fermatNumber.minFac + m < 2 ^ n⊢ n.fermatNumber.minFac ≤ 2 ^ (2 ^ n - m)
refine (lt_pow_succ_log_self (b:=2) (by n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:log 2 n.fermatNumber.minFac + m < 2 ^ n⊢ 1 < 2 norm_num All goals completed! 🐙) _).le.trans ?_
apply Nat.pow_le_pow_right (by n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:log 2 n.fermatNumber.minFac + m < 2 ^ n⊢ 2 > 0 norm_num All goals completed! 🐙)
apply le_tsub_of_add_le_right right n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:log 2 n.fermatNumber.minFac + m < 2 ^ n⊢ (log 2 n.fermatNumber.minFac).succ + m ≤ 2 ^ n
have := succ_le_of_lt hm right n:ℕPn:ℕ := n.fermatNumber.minFacm:ℕhm:log 2 n.fermatNumber.minFac + m < 2 ^ nthis:(log 2 n.fermatNumber.minFac + m).succ ≤ 2 ^ n⊢ (log 2 n.fermatNumber.minFac).succ + m ≤ 2 ^ n
omega All goals completed! 🐙@[category test, AMS 11]
theorem a_0 : a 0 = 0 := by ⊢ a 0 = 0 norm_num [a] ⊢ 1 - log2 3 = 0; simp [log2_def] All goals completed! 🐙@[category test, AMS 11]
theorem a_1 : a 1 = 0 := by ⊢ a 1 = 0 norm_num [a] ⊢ 2 - log2 5 = 0; simp [log2_def] All goals completed! 🐙@[category test, AMS 11]
theorem a_2 : a 2 = 0 := by ⊢ a 2 = 0 norm_num [a] ⊢ 4 - log2 17 = 0; simp [log2_def] All goals completed! 🐙@[category test, AMS 11]
theorem a_3 : a 3 = 0 := by ⊢ a 3 = 0
norm_num only [a, Nat.log2_eq_log_two,Nat.fermatNumber] All goals completed! 🐙@[category test, AMS 11]
theorem a_4 : a 4 = 0 := by ⊢ a 4 = 0
norm_num [a, fermatNumber, Nat.log2_eq_log_two] All goals completed! 🐙@[category test, AMS 11]
theorem a_5 : a 5 = 23 := by ⊢ a 5 = 23
norm_num [a, fermatNumber,Nat.log2_eq_log_two] All goals completed! 🐙@[category test, AMS 11]
theorem a_6 : a 6 = 46 := by ⊢ a 6 = 46
decide +native All goals completed! 🐙@[category test, AMS 11]
theorem a_7 : a 7 = 73 := by ⊢ a 7 = 73
sorry All goals completed! 🐙Conjecture: the dyadic valuation of A93179(n) - 1 does not exceed 2^n - a(n).
A93179(n) is minFac(fermatNumber n), the smallest prime factor of the n-th Fermat number. The conjecture states that if $P_n$ is the smallest prime factor of the $n$-th Fermat number, then $\nu_2(P_n - 1) \le 2^n - a(n)$. Substituting the definition of $a(n)$, this is equivalent to $\nu_2(P_n - 1) \le \lfloor \log_2(P_n) \rfloor$.
This is Conjecture 3.4 in [SA22].
@[category research solved, AMS 11]
theorem oeis_358684_conjecture_0 (n : ℕ) :
padicValNat 2 (minFac (fermatNumber n) - 1) ≤ 2 ^ n - a n := by n:ℕ⊢ padicValNat 2 (n.fermatNumber.minFac - 1) ≤ 2 ^ n - a n
delta fermatNumber and a n:ℕ⊢ padicValNat 2 ((2 ^ 2 ^ n + 1).minFac - 1) ≤ 2 ^ n - (2 ^ n - (2 ^ 2 ^ n + 1).minFac.log2)
rw [Nat.sub_sub_self n:ℕ⊢ padicValNat 2 ((2 ^ 2 ^ n + 1).minFac - 1) ≤ (2 ^ 2 ^ n + 1).minFac.log2n:ℕ⊢ (2 ^ 2 ^ n + 1).minFac.log2 ≤ 2 ^ n n:ℕ⊢ padicValNat 2 ((2 ^ 2 ^ n + 1).minFac - 1) ≤ (2 ^ 2 ^ n + 1).minFac.log2n:ℕ⊢ (2 ^ 2 ^ n + 1).minFac.log2 ≤ 2 ^ n] n:ℕ⊢ padicValNat 2 ((2 ^ 2 ^ n + 1).minFac - 1) ≤ (2 ^ 2 ^ n + 1).minFac.log2n:ℕ⊢ (2 ^ 2 ^ n + 1).minFac.log2 ≤ 2 ^ n
· n:ℕ⊢ padicValNat 2 ((2 ^ 2 ^ n + 1).minFac - 1) ≤ (2 ^ 2 ^ n + 1).minFac.log2 rw [Nat.log2_eq_log_two n:ℕ⊢ padicValNat 2 ((2 ^ 2 ^ n + 1).minFac - 1) ≤ log 2 (2 ^ 2 ^ n + 1).minFac n:ℕ⊢ padicValNat 2 ((2 ^ 2 ^ n + 1).minFac - 1) ≤ log 2 (2 ^ 2 ^ n + 1).minFac] n:ℕ⊢ padicValNat 2 ((2 ^ 2 ^ n + 1).minFac - 1) ≤ log 2 (2 ^ 2 ^ n + 1).minFac
apply Nat.le_log_of_pow_le (by n:ℕ⊢ 1 < 2 decide All goals completed! 🐙)
refine le_trans ?_ <| sub_le _ 1 n:ℕ⊢ 2 ^ padicValNat 2 ((2 ^ 2 ^ n + 1).minFac - 1) ≤ (2 ^ 2 ^ n + 1).minFac - 1
apply Nat.ordProj_le 2 n:ℕ⊢ (2 ^ 2 ^ n + 1).minFac - 1 ≠ 0
exact Nat.sub_ne_zero_of_lt (Nat.minFac_prime (by n:ℕ⊢ 2 ^ 2 ^ n + 1 ≠ 1 norm_num All goals completed! 🐙)).one_lt
· n:ℕ⊢ (2 ^ 2 ^ n + 1).minFac.log2 ≤ 2 ^ n rw [Nat.log2_eq_log_two n:ℕ⊢ log 2 (2 ^ 2 ^ n + 1).minFac ≤ 2 ^ n n:ℕ⊢ log 2 (2 ^ 2 ^ n + 1).minFac ≤ 2 ^ n] n:ℕ⊢ log 2 (2 ^ 2 ^ n + 1).minFac ≤ 2 ^ n
have : (2 ^ 2 ^ n) + 1 < 2 ^ ((2 ^ n) + 1) := by n:ℕ⊢ padicValNat 2 (n.fermatNumber.minFac - 1) ≤ 2 ^ n - a n n:ℕthis:2 ^ 2 ^ n + 1 < 2 ^ (2 ^ n + 1)⊢ log 2 (2 ^ 2 ^ n + 1).minFac ≤ 2 ^ n
simp [pow_add, mul_two] n:ℕthis:2 ^ 2 ^ n + 1 < 2 ^ (2 ^ n + 1)⊢ log 2 (2 ^ 2 ^ n + 1).minFac ≤ 2 ^ n n:ℕthis:2 ^ 2 ^ n + 1 < 2 ^ (2 ^ n + 1)⊢ log 2 (2 ^ 2 ^ n + 1).minFac ≤ 2 ^ n
refine Nat.le_of_lt_succ <| (2).log_lt_of_lt_pow ?_ ?_ refine_1 n:ℕthis:2 ^ 2 ^ n + 1 < 2 ^ (2 ^ n + 1)⊢ (2 ^ 2 ^ n + 1).minFac ≠ 0refine_2 n:ℕthis:2 ^ 2 ^ n + 1 < 2 ^ (2 ^ n + 1)⊢ (2 ^ 2 ^ n + 1).minFac < 2 ^ (2 ^ n).succ
· refine_1 n:ℕthis:2 ^ 2 ^ n + 1 < 2 ^ (2 ^ n + 1)⊢ (2 ^ 2 ^ n + 1).minFac ≠ 0 exact Nat.minFac_pos _|>.ne' All goals completed! 🐙
· refine_2 n:ℕthis:2 ^ 2 ^ n + 1 < 2 ^ (2 ^ n + 1)⊢ (2 ^ 2 ^ n + 1).minFac < 2 ^ (2 ^ n).succ exact (Nat.minFac_le (by n:ℕthis:2 ^ 2 ^ n + 1 < 2 ^ (2 ^ n + 1)⊢ 0 < 2 ^ 2 ^ n + 1 bound All goals completed! 🐙)).trans_lt thisend OeisA358684