Documentation

AnalyticNumberTheory.LargeSieve.GeomSum

AnalyticNumberTheory.LargeSieve.GeomSum #

Geometric-sum bound #

The first analytic step in Montgomery's additive large-sieve argument: for e(x) = exp(2πix) and ‖x‖, the distance to the nearest integer (expressed using Int.fract), the bound for nonintegral x is

|Σ_{n<N} e(nx)| ≤ min(N, 1/(2‖x‖)).

This module proves the trivial bound ≤ N for every x (each summand has norm 1) and the nontrivial bound ≤ 1/(2‖x‖) for nonintegral x, using the geometric-series formula, |e(x)−1| = 2|sin(πx)|, and |sin(πt)| ≥ 2·min(t,1−t) for t ∈ [0,1].

References: Iwaniec--Kowalski, "Analytic Number Theory" (2004), Ch. 7 Lemma 7.6; Davenport, "Multiplicative Number Theory", Ch. 27.

1. Distance to the nearest integer #

The fractional part x − ⌊x⌋ ∈ [0,1) is nonnegative.

The fractional part is less than 1.

Distance to the nearest integer: ‖x‖ = min(x − ⌊x⌋, 1 − (x − ⌊x⌋)).

Equations
Instances For

    2. Additive characters and geometric series #

    e(x) = exp(2πix): the additive character parametrized by .

    Equations
    Instances For

      e(nx) = e(x)^n (n : ℕ), by Complex.exp_nsmul.

      e is 1-periodic: e(x+1) = e(x), since exp(2πi) = 1.

      e(a + b) = e(a)·e(b): the character maps addition in to multiplication in .

      Rearrangement of the argument: (2π·c)·I = c·(2π·I), a polynomial identity for every complex c (in particular, for every real c).

      The character depends only on the class modulo 1: e(x) = e(fract x).

      theorem AnalyticNumberTheory.LargeSieve.geomSum_exp_eq_geomSeries (N : ) {x : } (hz : charReal x 1) :
      nFinset.range N, charReal (n * x) = (charReal (N * x) - 1) / (charReal x - 1)

      Geometric-series identity: Σ_{n<N} e(nx) = (e(Nx)−1)/(e(x)−1) when e(x) ≠ 1.

      Trivial bound: |Σ_{n<N} e(nx)| ≤ N, since every term has norm 1.

      3. Trigonometric identities #

      |e(x) − 1| = 2|sin(πx)|, combining e(x) = exp(πix)² and exp(πix) − exp(−πix) = 2i·sin(πx).

      theorem AnalyticNumberTheory.LargeSieve.abs_sin_pi_mul_ge_two_min (t : ) (ht0 : 0 t) (ht1 : t 1) :
      2 * min t (1 - t) |Real.sin (Real.pi * t)|

      |sin(πt)| ≥ 2·min(t, 1−t) for t ∈ [0,1]: use sin(πt) ≥ 2t on [0,1/2], and reflection sin(πt) = sin(π(1−t)) on [1/2,1].

      4. Nontrivial bound #

      e(x) ≠ 1 for x ∉ ℤ, since |e(x)−1| = 2|sin(π·fract x)| > 0.

      theorem AnalyticNumberTheory.LargeSieve.geomSum_exp_bound_far (N : ) {x : } (hx : ¬∃ (k : ), k = x) :
      nFinset.range N, charReal (n * x) 1 / (2 * distToInt x)

      Geometric-sum bound (nontrivial part): for x ∉ ℤ, |Σ_{n<N} e(nx)| ≤ 1/(2‖x‖). Combine the geometric-series formula, |e(x)−1| = 2|sin(π·fract x)|, and |sin(πt)| ≥ 2·min(t,1−t).

      theorem AnalyticNumberTheory.LargeSieve.geomSum_exp_bound (N : ) (x : ) :
      nFinset.range N, charReal (n * x) N ((¬∃ (k : ), k = x) → nFinset.range N, charReal (n * x) 1 / (2 * distToInt x))

      Geometric-sum bound: the two components of |Σ_{n<N} e(nx)| ≤ min(N, 1/(2‖x‖)). The trivial component holds for all N, x; the nontrivial component requires x ∉ ℤ.

      5. Interval geometric-sum bound #

      The character at a negative argument: e(−x) = star(e(x)), by exp(conj z) = conj(exp z).

      Character identity: e(a − b) = e(a)·star(e(b)), using the homomorphism property and conjugation.

      noncomputable def AnalyticNumberTheory.LargeSieve.charRealSubIcc (M : ) (N : ) (x : ) :

      The exponential sum Σ_{M<n≤M+N} e(nx) over the shifted interval (M, M+N] (n : ℤ).

      Equations
      Instances For
        theorem AnalyticNumberTheory.LargeSieve.charRealSubIcc_succ (M : ) (N : ) (x : ) :
        charRealSubIcc M (N + 1) x = charRealSubIcc M N x + charReal (↑(M + N + 1) * x)
        theorem AnalyticNumberTheory.LargeSieve.charRealSubIcc_eq_shift (M : ) (N : ) (x : ) :
        charRealSubIcc M N x = charReal (↑(M + 1) * x) * kFinset.range N, charReal (k * x)

        Interval exponential-sum factorization: Σ_{M<n≤M+N} e(nx) = e((M+1)x)·Σ_{k<N} e(kx).

        Interval exponential sum (trivial bound): |Σ_{M<n≤M+N} e(nx)| ≤ N.

        theorem AnalyticNumberTheory.LargeSieve.geomSum_exp_bound_Icc_far (M : ) (N : ) {x : } (hx : ¬∃ (k : ), k = x) :

        Interval exponential sum (nontrivial bound): for x ∉ ℤ, |Σ_{M<n≤M+N} e(nx)| ≤ 1/(2‖x‖).

        theorem AnalyticNumberTheory.LargeSieve.geomSum_exp_bound_Icc (M : ) (N : ) (x : ) :
        charRealSubIcc M N x N ((¬∃ (k : ), k = x) → charRealSubIcc M N x 1 / (2 * distToInt x))

        Interval geometric-sum bound: the two components of |Σ_{M<n≤M+N} e(nx)| ≤ min(N, 1/(2‖x‖)), with the nontrivial component conditional on x ∉ ℤ.