diff --git a/Cslib.lean b/Cslib.lean
index 8c6984eb..2be39bd8 100644
--- a/Cslib.lean
+++ b/Cslib.lean
@@ -36,6 +36,7 @@ public import Cslib.Foundations.Data.OmegaSequence.Flatten
 public import Cslib.Foundations.Data.OmegaSequence.InfOcc
 public import Cslib.Foundations.Data.OmegaSequence.Init
 public import Cslib.Foundations.Data.OmegaSequence.Temporal
+public import Cslib.Foundations.Data.RelatesInSteps
 public import Cslib.Foundations.Data.Relation
 public import Cslib.Foundations.Lint.Basic
 public import Cslib.Foundations.Semantics.FLTS.Basic
diff --git a/Cslib/Foundations/Data/RelatesInSteps.lean b/Cslib/Foundations/Data/RelatesInSteps.lean
new file mode 100644
index 00000000..07c66d0f
--- /dev/null
+++ b/Cslib/Foundations/Data/RelatesInSteps.lean
@@ -0,0 +1,222 @@
+/-
+Copyright (c) 2025 Bolton Bailey. All rights reserved.
+Released under Apache 2.0 license as described in the file LICENSE.
+Authors: Bolton Bailey
+-/
+
+module
+
+public import Cslib.Init
+public import Mathlib.Logic.Relation
+
+@[expose] public section
+
+variable {α : Type*} {r : α → α → Prop} {a b c : α}
+
+/-! # Relations Across Steps
+
+This file defines `Relation.RelatesInSteps` (and `Relation.RelatesWithinSteps`).
+These are inductively defines propositions that communicate whether a relation forms a
+chain of length `n` (or at most `n`) between two elements.
+-/
+
+namespace Relation
+
+/--
+A relation `r` relates two elements of `α` in `n` steps
+if there is a chain of `n` pairs `(t_i, t_{i+1})` such that `r t_i t_{i+1}` for each `i`,
+starting from the first element and ending at the second.
+-/
+inductive RelatesInSteps (r : α → α → Prop) : α → α → ℕ → Prop
+  | refl (a : α) : RelatesInSteps r a a 0
+  | tail (t t' t'' : α) (n : ℕ) (h₁ : RelatesInSteps r t t' n) (h₂ : r t' t'') :
+      RelatesInSteps r t t'' (n + 1)
+
+theorem RelatesInSteps.reflTransGen (h : RelatesInSteps r a b n) : ReflTransGen r a b := by
+  induction h with
+  | refl => rfl
+  | tail _ _ _ _ h ih => exact .tail ih h
+
+theorem ReflTransGen.relatesInSteps (h : ReflTransGen r a b) : ∃ n, RelatesInSteps r a b n := by
+  induction h with
+  | refl => exact ⟨0, .refl a⟩
+  | tail _ _ ih =>
+    obtain ⟨n, _⟩ := ih
+    exact ⟨n + 1, by grind [RelatesInSteps]⟩
+
+lemma RelatesInSteps.single {a b : α} (h : r a b) : RelatesInSteps r a b 1 :=
+  tail a a b 0 (refl a) h
+
+theorem RelatesInSteps.head (t t' t'' : α) (n : ℕ) (h₁ : r t t')
+    (h₂ : RelatesInSteps r t' t'' n) : RelatesInSteps r t t'' (n+1) := by
+  induction h₂ with
+  | refl =>
+    exact single h₁
+  | tail _ _ n _ hcd hac =>
+    exact tail _ _ _ (n + 1) hac hcd
+
+@[elab_as_elim]
+theorem RelatesInSteps.head_induction_on {motive : ∀ (a : α) (n : ℕ), RelatesInSteps r a b n → Prop}
+    {a : α} {n : ℕ} (h : RelatesInSteps r a b n) (hrefl : motive b 0 (.refl b))
+    (hhead : ∀ {a c n} (h' : r a c) (h : RelatesInSteps r c b n),
+      motive c n h → motive a (n + 1) (h.head a c b n h')) :
+    motive a n h := by
+  induction h with
+  | refl => exact hrefl
+  | tail t' b' m hstep hrt'b' hstep_ih =>
+    apply hstep_ih
+    · exact hhead hrt'b' _ hrefl
+    · exact fun h1 h2 ↦ hhead h1 (.tail _ t' b' _ h2 hrt'b')
+
+lemma RelatesInSteps.zero {a b : α} (h : RelatesInSteps r a b 0) : a = b := by
+  cases h
+  rfl
+
+@[simp]
+lemma RelatesInSteps.zero_iff {a b : α} : RelatesInSteps r a b 0 ↔ a = b := by
+  constructor
+  · exact RelatesInSteps.zero
+  · intro rfl
+    exact RelatesInSteps.refl a
+
+lemma RelatesInSteps.trans {a b c : α} {n m : ℕ}
+    (h₁ : RelatesInSteps r a b n) (h₂ : RelatesInSteps r b c m) :
+    RelatesInSteps r a c (n + m) := by
+  induction h₂ generalizing a n with
+  | refl => simp [h₁]
+  | tail t' t'' k hsteps hstep ih =>
+    rw [← Nat.add_assoc]
+    exact .tail _ t' t'' (n + k) (ih h₁) hstep
+
+lemma RelatesInSteps.succ {n : ℕ} (h : RelatesInSteps r a b (n + 1)) :
+    ∃ t', RelatesInSteps r a t' n ∧ r t' b := by
+  cases h with
+  | tail t' _ _ hsteps hstep => exact ⟨t', hsteps, hstep⟩
+
+lemma RelatesInSteps.succ_iff {a b : α} {n : ℕ} :
+    RelatesInSteps r a b (n + 1) ↔ ∃ t', RelatesInSteps r a t' n ∧ r t' b := by
+  constructor
+  · exact RelatesInSteps.succ
+  · rintro ⟨t', h_steps, h_red⟩
+    exact .tail _ t' b n h_steps h_red
+
+lemma RelatesInSteps.succ' {a b : α} : ∀ {n : ℕ}, RelatesInSteps r a b (n + 1) →
+    ∃ t', r a t' ∧ RelatesInSteps r t' b n := by
+  intro n h
+  obtain ⟨t', hsteps, hstep⟩ := succ h
+  cases n with
+  | zero =>
+    rw [zero_iff] at hsteps
+    subst hsteps
+    exact ⟨b, hstep, .refl _⟩
+  | succ k' =>
+    obtain ⟨t''', h_red''', h_steps'''⟩ := succ' hsteps
+    exact ⟨t''', h_red''', .tail _ _ b k' h_steps''' hstep⟩
+
+lemma RelatesInSteps.succ'_iff {a b : α} {n : ℕ} :
+    RelatesInSteps r a b (n + 1) ↔ ∃ t', r a t' ∧ RelatesInSteps r t' b n := by
+  constructor
+  · exact succ'
+  · rintro ⟨t', h_red, h_steps⟩
+    exact h_steps.head a t' b n h_red
+
+/--
+If `h : α → ℕ` increases by at most 1 on each step of `r`,
+then the value of `h` at the output is at most `h` at the input plus the number of steps.
+-/
+lemma RelatesInSteps.apply_le_apply_add {a b : α} (h : α → ℕ)
+    (h_step : ∀ a b, r a b → h b ≤ h a + 1)
+    (m : ℕ) (hevals : RelatesInSteps r a b m) :
+    h b ≤ h a + m := by
+  induction hevals with
+  | refl => simp
+  | tail t' t'' k _ hstep ih =>
+    have h_step' := h_step t' t'' hstep
+    lia
+
+/--
+If `g` is a homomorphism from `r` to `r'` (i.e., it preserves the reduction relation),
+then `RelatesInSteps` is preserved under `g`.
+-/
+lemma RelatesInSteps.map {α α' : Type*}
+    {r : α → α → Prop} {r' : α' → α' → Prop}
+    (g : α → α') (hg : ∀ a b, r a b → r' (g a) (g b))
+    {a b : α} {n : ℕ} (h : RelatesInSteps r a b n) :
+    RelatesInSteps r' (g a) (g b) n := by
+  induction h with
+  | refl => exact RelatesInSteps.refl (g _)
+  | tail t' t'' m _ hstep ih =>
+    exact .tail (g _) (g t') (g t'') m ih (hg t' t'' hstep)
+
+/--
+`RelatesWithinSteps` is a variant of `RelatesInSteps` that allows for a loose bound.
+It states that `a` relates to `b` in *at most* `n` steps.
+-/
+def RelatesWithinSteps (r : α → α → Prop) (a b : α) (n : ℕ) : Prop :=
+  ∃ m ≤ n, RelatesInSteps r a b m
+
+/-- `RelatesInSteps` implies `RelatesWithinSteps` with the same bound. -/
+lemma RelatesWithinSteps.of_relatesInSteps {a b : α} {n : ℕ} (h : RelatesInSteps r a b n) :
+    RelatesWithinSteps r a b n :=
+  ⟨n, Nat.le_refl n, h⟩
+
+lemma RelatesWithinSteps.refl (a : α) : RelatesWithinSteps r a a 0 :=
+  RelatesWithinSteps.of_relatesInSteps (RelatesInSteps.refl a)
+
+lemma RelatesWithinSteps.single {a b : α} (h : r a b) : RelatesWithinSteps r a b 1 :=
+  RelatesWithinSteps.of_relatesInSteps (RelatesInSteps.single h)
+
+lemma RelatesWithinSteps.zero {a b : α} (h : RelatesWithinSteps r a b 0) : a = b := by
+  obtain ⟨m, hm, hevals⟩ := h
+  have : m = 0 := Nat.le_zero.mp hm
+  subst this
+  exact RelatesInSteps.zero hevals
+
+@[simp]
+lemma RelatesWithinSteps.zero_iff {a b : α} : RelatesWithinSteps r a b 0 ↔ a = b := by
+  constructor
+  · exact RelatesWithinSteps.zero
+  · intro h
+    subst h
+    exact RelatesWithinSteps.refl a
+
+/-- Transitivity of `RelatesWithinSteps` in the sum of the step bounds. -/
+@[trans]
+lemma RelatesWithinSteps.trans {a b c : α} {n₁ n₂ : ℕ}
+    (h₁ : RelatesWithinSteps r a b n₁) (h₂ : RelatesWithinSteps r b c n₂) :
+    RelatesWithinSteps r a c (n₁ + n₂) := by
+  obtain ⟨m₁, hm₁, hevals₁⟩ := h₁
+  obtain ⟨m₂, hm₂, hevals₂⟩ := h₂
+  use m₁ + m₂
+  constructor
+  · lia
+  · exact RelatesInSteps.trans hevals₁ hevals₂
+
+lemma RelatesWithinSteps.of_le {a b : α} {n₁ n₂ : ℕ}
+    (h : RelatesWithinSteps r a b n₁) (hn : n₁ ≤ n₂) :
+    RelatesWithinSteps r a b n₂ := by
+  obtain ⟨m, hm, hevals⟩ := h
+  exact ⟨m, Nat.le_trans hm hn, hevals⟩
+
+/-- If `h : α → ℕ` increases by at most 1 on each step of `r`,
+then the value of `h` at the output is at most `h` at the input plus the step bound. -/
+lemma RelatesWithinSteps.apply_le_apply_add {a b : α} (h : α → ℕ)
+    (h_step : ∀ a b, r a b → h b ≤ h a + 1)
+    (n : ℕ) (hevals : RelatesWithinSteps r a b n) :
+    h b ≤ h a + n := by
+  obtain ⟨m, hm, hevals_m⟩ := hevals
+  have := RelatesInSteps.apply_le_apply_add h h_step m hevals_m
+  lia
+
+/--
+If `g` is a homomorphism from `r` to `r'` (i.e., it preserves the reduction relation),
+then `RelatesWithinSteps` is preserved under `g`.
+-/
+lemma RelatesWithinSteps.map {α α' : Type*} {r : α → α → Prop} {r' : α' → α' → Prop}
+    (g : α → α') (hg : ∀ a b, r a b → r' (g a) (g b))
+    {a b : α} {n : ℕ} (h : RelatesWithinSteps r a b n) :
+    RelatesWithinSteps r' (g a) (g b) n := by
+  obtain ⟨m, hm, hevals⟩ := h
+  exact ⟨m, hm, RelatesInSteps.map g hg hevals⟩
+
+end Relation