OAI/NumberTheory/JointDickman/Arithmetic/SieveInterval.lean
4 declarations, 0 contain sorry, 0 declared as axioms
The interval sieve with its actual arithmetic remainder
- def
OAI.JointDickman.rootEventsline 14Source documentation
Forbidden primes of an actual integer point.
noncomputable def rootEvents (P : Finset ℕ) (A : ∀ p : ℕ, Finset (ZMod p)) (n : ℕ) : Finset ℕ := P.filter (fun p => (n : ZMod p) ∈ A p) - theorem
OAI.JointDickman.root_sieve_remainderline 18Source documentation
CRT evaluates the actual remainder used by the upper sieve.
theorem root_sieve_remainder (P : Finset ℕ) (hP : ∀ p ∈ P, p.Prime) (A : ∀ p : ℕ, Finset (ZMod p)) {a b : ℕ} (hab : a ≤ b) {D : Finset ℕ} (hD : D ⊆ P) : |sieveRemainder (fun n : Ico a b => rootEvents P A n) (fun _ => 1) ((b : ℝ) - a) (fun p => ((A p).card : ℝ) / p) D| ≤ 2 * ∏ p ∈ D, ((A p).card : ℝ) := by …Mentions: OAI.JointDickman.sieveRemainder, OAI.JointDickman.rootEvents
- theorem
OAI.JointDickman.weighted_interval_sieveline 44Source documentation
Uniform weighted interval bound, derived from Ford and Mertens.
theorem weighted_interval_sieve (hFord : PublishedInputs.FordUpperSieveInput) (hM : PublishedInputs.PrimeReciprocalMertensInput) {k : ℝ} (hk : 0 < k) : ∃ C : ℝ, 0 < C ∧ ∀ (P : Finset ℕ) (A : ∀ p : ℕ, Finset (ZMod p)) (θ : ℕ → ℝ) (a b : ℕ) (z : ℝ), a ≤ b → 2 ≤ z → (∀ p ∈ P, p.Prime ∧ (p : ℝ) ≤ z ∧ 2 * k ≤ p) → (∀ p ∈ P, ((A p).card : ℝ) ≤ k) → (∀ p ∈ P, 0 ≤ θ p ∧ θ p ≤ 1) → (∑ n ∈ Ico a b, ∏ p ∈ rootEvents P A n, θ p) ≤ C * ((b : ℝ) - a) * (∏ p ∈ P, (1 - ((A p).card : ℝ) / p + ((A p).card : ℝ) / p * θ p)) + 2 * ∑ D ∈ P.powerset.filter (fun D => (∏ p ∈ D, p : ℕ) ≤ z), ∏ p ∈ D, ((A p).card : ℝ) := by …Mentions: OAI.JointDickman.PublishedInputs.FordUpperSieveInput, OAI.JointDickman.PublishedInputs.PrimeReciprocalMertensInput, OAI.JointDickman.rootEvents
- theorem
OAI.JointDickman.weighted_interval_sieve_polynomialline 86Source documentation
A completely explicit polynomial remainder suffices at a small fixed power cutoff. All interval positions and reducing weights are uniform.
theorem weighted_interval_sieve_polynomial (hFord : PublishedInputs.FordUpperSieveInput) (hM : PublishedInputs.PrimeReciprocalMertensInput) {k : ℕ} (hk : 0 < k) : ∃ C : ℝ, 0 < C ∧ ∀ (P : Finset ℕ) (A : ∀ p : ℕ, Finset (ZMod p)) (θ : ℕ → ℝ) (a b Z : ℕ), a ≤ b → 2 ≤ Z → (∀ p ∈ P, p.Prime ∧ p ≤ Z ∧ 2 * k ≤ p) → (∀ p ∈ P, (A p).card ≤ k) → (∀ p ∈ P, 0 ≤ θ p ∧ θ p ≤ 1) → (∑ n ∈ Ico a b, ∏ p ∈ rootEvents P A n, θ p) ≤ C * ((b : ℝ) - a) * (∏ p ∈ P, (1 - ((A p).card : ℝ) / p + ((A p).card : ℝ) / p * θ p)) + 2 * (Z + 1 : ℝ) * (Z : ℝ) ^ k := by …Mentions: OAI.JointDickman.PublishedInputs.FordUpperSieveInput, OAI.JointDickman.PublishedInputs.PrimeReciprocalMertensInput, OAI.JointDickman.rootEvents