/-
Copyright 2025 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.
-/
module
public import FormalConjecturesForMathlib.Data.Sym.Sym2
public import Mathlib.Algebra.BigOperators.Ring.Finset
public import Mathlib.Data.Finset.Sym
public import Mathlib.Data.Nat.Choose.Basic
public import Mathlib.Data.Sym.Card
public import Mathlib.Data.ZMod.Basic
public import Mathlib.Order.Lattice.Nat
public import Mathlib.Topology.MetricSpace.Defs@[expose] public sectionopen scoped Finsetvariable {X : Type*} [MetricSpace X]
The number of pairs of points of a finite set s in a metric space that are distance 1 apart.
noncomputable def unitDistNum (s : Finset X) : ℕ := #{p ∈ s.sym2 | dist p.out.1 p.out.2 = 1}The set of distances determined by a finite set of points in a metric space.
noncomputable def distanceSet (points : Finset X) : Finset ℝ :=
points.offDiag.image fun (pair : X × X) => dist pair.1 pair.2Given a finite set of points in a metric space, we define the number of distinct distances between pairs of points.
noncomputable def distinctDistances (points : Finset X) : ℕ :=
#(distanceSet points)variable (X) in
The minimum number of distinct distances determined by a set of n points in X.
noncomputable def minimalDistinctDistances (n : ℕ) : ℕ :=
sInf {m : ℕ | ∃ points : Finset X, #points = n ∧ distinctDistances points = m}
The multiplicity of the distance d determined by points, that is, the number of unordered
pairs of distinct points at distance d apart.
noncomputable def distanceMultiplicity (points : Finset X) (d : ℝ) : ℕ :=
#(points.offDiag.filter fun (pair : X × X) => dist pair.1 pair.2 = d) / 2open Classical inGiven a finite set of points in a metric space, we define the number of distinct distances between a given point and all other points.
noncomputable def distinctDistancesFrom (points : Finset X) (pt : X) : ℕ :=
#((points.erase pt).image fun x => dist x pt)open Classical in
The number of unit-distance pairs of a finite set of n points is at most $\binom{n}{2}$,
the total number of unordered pairs of distinct points.
X:Type u_1inst✝:MetricSpace Xs:Finset Xp:Sym2 Xhp:p ∈ {p ∈ s.sym2 | dist (Quot.out p).1 (Quot.out p).2 = 1}hps:p ∈ Finset.image (Function.uncurry Sym2.mk) s.diag ∨ p ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiaghpd:dist (Quot.out p).1 (Quot.out p).2 = 1⊢ p ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiag
rcases hps with h | h inl X:Type u_1inst✝:MetricSpace Xs:Finset Xp:Sym2 Xhp:p ∈ {p ∈ s.sym2 | dist (Quot.out p).1 (Quot.out p).2 = 1}hpd:dist (Quot.out p).1 (Quot.out p).2 = 1h:p ∈ Finset.image (Function.uncurry Sym2.mk) s.diag⊢ p ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiaginr X:Type u_1inst✝:MetricSpace Xs:Finset Xp:Sym2 Xhp:p ∈ {p ∈ s.sym2 | dist (Quot.out p).1 (Quot.out p).2 = 1}hpd:dist (Quot.out p).1 (Quot.out p).2 = 1h:p ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiag⊢ p ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiag
· inl X:Type u_1inst✝:MetricSpace Xs:Finset Xp:Sym2 Xhp:p ∈ {p ∈ s.sym2 | dist (Quot.out p).1 (Quot.out p).2 = 1}hpd:dist (Quot.out p).1 (Quot.out p).2 = 1h:p ∈ Finset.image (Function.uncurry Sym2.mk) s.diag⊢ p ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiag obtain ⟨⟨x, y⟩, hxy, rfl⟩ := Finset.mem_image.mp h inl X:Type u_1inst✝:MetricSpace Xs:Finset Xx:Xy:Xhxy:(x, y) ∈ s.diaghp:Function.uncurry Sym2.mk (x, y) ∈ {p ∈ s.sym2 | dist (Quot.out p).1 (Quot.out p).2 = 1}hpd:dist (Quot.out (Function.uncurry Sym2.mk (x, y))).1 (Quot.out (Function.uncurry Sym2.mk (x, y))).2 = 1h:Function.uncurry Sym2.mk (x, y) ∈ Finset.image (Function.uncurry Sym2.mk) s.diag⊢ Function.uncurry Sym2.mk (x, y) ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiag
obtain ⟨-, rfl⟩ : _ ∧ x = y := Finset.mem_diag.mp hxy inl X:Type u_1inst✝:MetricSpace Xs:Finset Xx:Xhxy:(x, x) ∈ s.diaghp:Function.uncurry Sym2.mk (x, x) ∈ {p ∈ s.sym2 | dist (Quot.out p).1 (Quot.out p).2 = 1}hpd:dist (Quot.out (Function.uncurry Sym2.mk (x, x))).1 (Quot.out (Function.uncurry Sym2.mk (x, x))).2 = 1h:Function.uncurry Sym2.mk (x, x) ∈ Finset.image (Function.uncurry Sym2.mk) s.diag⊢ Function.uncurry Sym2.mk (x, x) ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiag
simp at hpd All goals completed! 🐙
· inr X:Type u_1inst✝:MetricSpace Xs:Finset Xp:Sym2 Xhp:p ∈ {p ∈ s.sym2 | dist (Quot.out p).1 (Quot.out p).2 = 1}hpd:dist (Quot.out p).1 (Quot.out p).2 = 1h:p ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiag⊢ p ∈ Finset.image (Function.uncurry Sym2.mk) s.offDiag exact h All goals completed! 🐙
The ordered pairs of distinct points of points at distance d come in swapped pairs, so
there are evenly many of them.
theorem even_card_offDiag_filter_dist_eq (points : Finset X) (d : ℝ) :
Even #(points.offDiag.filter fun (pair : X × X) => dist pair.1 pair.2 = d) := by X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝ⊢ Even #({pair ∈ points.offDiag | dist pair.1 pair.2 = d})
set F := points.offDiag.filter fun (pair : X × X) => dist pair.1 pair.2 = d with hF X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}⊢ Even #F
rw [← ZMod.natCast_eq_zero_iff_even X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}⊢ ↑(#F) = 0 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}⊢ ↑(#F) = 0] X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}⊢ ↑(#F) = 0
have h : ∑ _a ∈ F, (1 : ZMod 2) = 0 := by X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝ⊢ Even #({pair ∈ points.offDiag | dist pair.1 pair.2 = d}) X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0
refine Finset.sum_involution (fun a _ => a.swap) (fun a _ => by X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}a:X × Xx✝:a ∈ F⊢ 1 + 1 = 0 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0 decide All goals completed! 🐙 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0) ?_ ?_
(fun a _ => Prod.swap_swap a)
· refine_1 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}⊢ ∀ a ∈ F, 1 ≠ 0 → a.swap ≠ a X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0 intro a ha _ hsw refine_1 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}a:X × Xha:a ∈ Fa✝:1 ≠ 0hsw:a.swap = a⊢ False X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0
have := (Finset.mem_offDiag.1 (Finset.mem_filter.1 ha).1).2.2 refine_1 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}a:X × Xha:a ∈ Fa✝:1 ≠ 0hsw:a.swap = athis:a.1 ≠ a.2⊢ False X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0
exact this (Prod.ext_iff.1 hsw).1.symm All goals completed! 🐙 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0
· refine_2 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}⊢ ∀ a ∈ F, a.swap ∈ F X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0 intro a ha refine_2 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}a:X × Xha:a ∈ F⊢ a.swap ∈ F X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0
obtain ⟨hoff, hd⟩ := Finset.mem_filter.1 ha refine_2 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}a:X × Xha:a ∈ Fhoff:a ∈ points.offDiaghd:dist a.1 a.2 = d⊢ a.swap ∈ F X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0
obtain ⟨h1, h2, h3⟩ := Finset.mem_offDiag.1 hoff refine_2 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}a:X × Xha:a ∈ Fhoff:a ∈ points.offDiaghd:dist a.1 a.2 = dh1:a.1 ∈ pointsh2:a.2 ∈ pointsh3:a.1 ≠ a.2⊢ a.swap ∈ F X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0
exact Finset.mem_filter.2 ⟨Finset.mem_offDiag.2 ⟨h2, h1, Ne.symm h3⟩, by X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}a:X × Xha:a ∈ Fhoff:a ∈ points.offDiaghd:dist a.1 a.2 = dh1:a.1 ∈ pointsh2:a.2 ∈ pointsh3:a.1 ≠ a.2⊢ dist a.swap.1 a.swap.2 = d X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0
simpa [dist_comm] using hd All goals completed! 🐙 X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0⟩ X:Type u_1inst✝:MetricSpace Xpoints:Finset Xd:ℝF:Finset (X × X) := {pair ∈ points.offDiag | dist pair.1 pair.2 = d}hF:F = {pair ∈ points.offDiag | dist pair.1 pair.2 = d}h:∑ _a ∈ F, 1 = 0⊢ ↑(#F) = 0
simpa [Finset.sum_const, nsmul_eq_mul] using h All goals completed! 🐙
The multiplicities of the distances determined by points add up to the number of
unordered pairs of distinct points.
theorem sum_distanceMultiplicity (points : Finset X) :
∑ d ∈ distanceSet points, distanceMultiplicity points d = (#points).choose 2 := by X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ ∑ d ∈ distanceSet points, distanceMultiplicity points d = (#points).choose 2
unfold distanceSet distanceMultiplicity X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ ∑ d ∈ Finset.image (fun pair ↦ dist pair.1 pair.2) points.offDiag,
#({pair ∈ points.offDiag | dist pair.1 pair.2 = d}) / 2 =
(#points).choose 2
rw [← Nat.sum_div (fun d _ => even_iff_two_dvd.1 (even_card_offDiag_filter_dist_eq points d)), X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ (∑ i ∈ Finset.image (fun pair ↦ dist pair.1 pair.2) points.offDiag,
#({pair ∈ points.offDiag | dist pair.1 pair.2 = i})) /
2 =
(#points).choose 2 X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2
← Finset.card_eq_sum_card_image (fun (pair : X × X) => dist pair.1 pair.2) points.offDiag, X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ #points.offDiag / 2 = (#points).choose 2 X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2
Finset.offDiag_card, X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ (#points * #points - #points) / 2 = (#points).choose 2 X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2 Nat.choose_two_right X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2 X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2] X:Type u_1inst✝:MetricSpace Xpoints:Finset X⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2
rcases Nat.eq_zero_or_pos #points with h | h inl X:Type u_1inst✝:MetricSpace Xpoints:Finset Xh:#points = 0⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2inr X:Type u_1inst✝:MetricSpace Xpoints:Finset Xh:#points > 0⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2
· inl X:Type u_1inst✝:MetricSpace Xpoints:Finset Xh:#points = 0⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2 simp [h] All goals completed! 🐙
· inr X:Type u_1inst✝:MetricSpace Xpoints:Finset Xh:#points > 0⊢ (#points * #points - #points) / 2 = #points * (#points - 1) / 2 rw [Nat.mul_sub_one inr X:Type u_1inst✝:MetricSpace Xpoints:Finset Xh:#points > 0⊢ (#points * #points - #points) / 2 = (#points * #points - #points) / 2 All goals completed! 🐙] All goals completed! 🐙