Skip to content
Merged
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
155 changes: 155 additions & 0 deletions FormalConjectures/OEIS/358684.lean
Original file line number Diff line number Diff line change
@@ -0,0 +1,155 @@
/-
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 FormalConjectures.Util.ProblemImports

/-!
$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:*
- [A358684](https://oeis.org/A358684)
- [SA22](https://doi.org/10.26493/2590-9770.1473.ec5) Lorenzo Sauras-Altuzarra, *Some properties of the factors of Fermat numbers*, Art Discrete Appl. Math. (2022).
-/


open Nat

/--
A358684: $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 by
use 2^n
simp only [tsub_self, pow_zero]
simp [Pn]
unfold Nat.fermatNumber
exact (Nat.minFac_prime (by norm_num)).one_lt
)


/-
The log2 of the smallest prime factor of $F_n$ is at most $2^n$.
-/
@[category undergraduate, AMS 11]
private lemma log2_minFac_le (n : ℕ) : log2 (fermatNumber n).minFac ≤ 2^n := by
rw [log2_eq_log_two]
refine (log_mono_right (minFac_le (by rw [fermatNumber]; norm_num))).trans_eq ?_
rw [fermatNumber, log_eq_of_pow_le_of_lt_pow (le_add_right _ _)]
rw [pow_succ, mul_two]; apply add_lt_add_left (one_lt_pow (by norm_num) (by norm_num))


/--
The minimization definition is equivalent to the closed form.
-/
@[category API, AMS 11]
theorem a_equiv_a' (n : ℕ) : a n = a' n := by
unfold a a'; set Pn := (fermatNumber n).minFac
rw [Eq.comm, Nat.find_eq_iff]
constructor
· rw [tsub_tsub_cancel_of_le (log2_minFac_le n), log2_eq_log_two]
refine lt_of_le_of_ne (pow_log_le_self 2 (Nat.Prime.ne_zero (minFac_prime ?_))) fun h => ?_
· rw [fermatNumber]
norm_num
· have hPn : Pn.Prime := minFac_prime (by rw [fermatNumber]; norm_num)
have hnz := log2_eq_log_two ▸ (Nat.log_pos one_lt_two hPn.two_le).ne'
refine (Nat.not_even_iff_odd.mpr <| (odd_pow_iff hnz).mp ?_) even_two
rw [log2_eq_log_two, h]
exact Odd.of_dvd_nat (odd_fermatNumber n) (minFac_dvd _)
· intro m hm; simp only [not_lt]
rw [lt_tsub_iff_left] at hm
rw [log2_eq_log_two] at hm
refine (lt_pow_succ_log_self (b:=2) (by norm_num) _).le.trans ?_
apply Nat.pow_le_pow_right (by norm_num)
apply le_tsub_of_add_le_right
have := succ_le_of_lt hm
omega


@[category test, AMS 11]
theorem zero : a 0 = 0 := by
norm_num [a, Nat.log2]

@[category test, AMS 11]
theorem one : a 1 = 0 := by
norm_num [a, Nat.log2]

@[category test, AMS 11]
theorem two : a 2 = 0 := by
norm_num [a, Nat.log2]

@[category test, AMS 11]
theorem three : a 3 = 0 := by
norm_num only [a, Nat.log2_eq_log_two,Nat.fermatNumber]

@[category test, AMS 11]
theorem four : a 0 = 0 := by
norm_num[a, Nat.log2]

@[category test, AMS 11]
theorem five : a 5 = 23 := by
norm_num [a, fermatNumber,Nat.log2_eq_log_two]

@[category test, AMS 11]
theorem six : a 6 = 46 := by
-- AlphaProof failed to close this goal
sorry

@[category test, AMS 11]
theorem seven : a 7 = 73 := by
sorry

/--
Conjecture: the dyadic valuation of A093179(n) - 1 does not exceed 2^n - a(n).

A093179(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
delta fermatNumber and a
rw [Nat.sub_sub_self]
· rw [Nat.log2_eq_log_two]
apply Nat.le_log_of_pow_le (by decide)
refine le_trans ?_ <| sub_le _ 1
apply Nat.ordProj_le 2
exact Nat.sub_ne_zero_of_lt (Nat.minFac_prime (by norm_num)).one_lt
· rw [Nat.log2_eq_log_two]
have : (2 ^ 2 ^ n) + 1 < 2 ^ ((2 ^ n) + 1) := by
simp [pow_add, mul_two]
refine Nat.le_of_lt_succ <| (2).log_lt_of_lt_pow ?_ ?_
· exact Nat.minFac_pos _|>.ne'
· exact (Nat.minFac_le (by bound)).trans_lt this
Loading