OAI/NumberTheory/JointDickman/Arithmetic/SieveInterval.lean

GitHub source file

4 declarations, 0 contain sorry, 0 declared as axioms

The interval sieve with its actual arithmetic remainder

  1. def OAI.JointDickman.rootEvents line 14

    Source 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)
  2. theorem OAI.JointDickman.root_sieve_remainder line 18

    Source 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

  3. theorem OAI.JointDickman.weighted_interval_sieve line 44

    Source 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

  4. theorem OAI.JointDickman.weighted_interval_sieve_polynomial line 86

    Source 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