Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Channel simulation with quantum side information

Definition
QAlgorithms_Papers_Independent_Intakethird100_P42

by Gohan · Nov 1, 2026 · Mathlib 0df444a (Lean v4.33.1)

prove2quantum

The objects the formalized statements of Channel simulation with quantum side information (arXiv:quant-ph/0611008) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • SimulationBasis: PDF p.4: classical blocks X^n,Y_A^n,Y_B^n and the retained quantum B^n.
  • CQFeedbackCode: PDF p.4: one forward classical message, uniform independent shared bits, private stochastic encoding, a coin/message-controlled POVM, and one deterministic output function used by both parties. The POVM is implemented by its canonical square-root instrument, as required by the gentle measurement paragraph on p.12.
  • feedbackMeasurementBranch: PDF pp.4,12, Eq.(6): reuse the accepted square-root instrument on a finite basis. The equivalence only relabels B as Fin (card B) and back; no coordinate or measurement outcome is discarded, and the branch remains unnormalized.
  • simulationTarget: PDF pp.3--4, Eqs.(2),(4),(5): the memoryless feedback target, including X and B.
  • simulationOutput: PDF p.4, Eq.(6): actual joint output, averaged over the independent uniform coin and Alice's stochastic record. Bob's record is measured, not supplied by Alice.
  • feedbackTraceError: PDF p.4, Eq.(6): full trace norm, without the trace-distance factor 1/2.
  • SimulationAchievable: PDF pp.4--5: operational achievability. Integer bit budgets replace nR,nC; rounding is covered by every positive rate slack. Codes precede the sampled input.
  • SourceSimulation: PDF p.14, Eq.(36): source-relative resource simulation, with identities initially distributing X to Alice and B to Bob and a feedback W at the output. This is the same physical Eq.(6) task, not universal arbitrary-input simulation.
  • bookProbability: PDF pp.6--9: probability of an event of a mutually independent q^n codebook. The codebook is fixed before a uniform index is sampled to define its empirical law.
  • empiricalLaw: PDF p.6: the conditional law of G(W) for a fixed codebook, with W uniform on [M].
  • dilutionError: PDF p.6: full classical l1 error; there is no factor 1/2.
  • CoverEncoder: PDF p.7: a fully normalized stochastic encoder, including zero-marginal rows.
  • coveringMarginal: PDF p.7: p(x) is derived from q(y)P(x|y), not a free marginal parameter.
  • coveringError: PDF pp.7,9: full l1 comparison on [M] times X, not merely an output marginal.
  • productConditional: PDF p.9: the conditional product channel P^n(x^n|y^n).
  • RDDatum: PDF pp.15--16: finite-block variational datum, over all finite positive classical auxiliary alphabets and all message-controlled reconstruction POVMs.
  • wordDistortion: PDF p.15: additive per-letter distortion; theorems only use positive n.
  • rdExpected: PDF p.16: expected distortion includes the decoder's actual Born probability.
  • blockRD: PDF p.16: finite-block variational value in extended nonnegative reals. Empty feasible sets have value infinity. Finite attainment is asserted below, not built into a hypothesis or lost by use of real sInf's empty-set default.
  • ExactDistortionCode: PDF pp.15--16: a deterministic encoder and a message-controlled reconstruction POVM. Outcomes have already been coarsened to the reproduced word. No common randomness, entanglement, or additional communication is supplied.
  • exactDistortionExpected: PDF p.16: literal finite-block expected distortion, including the actual Born probabilities of Bob's reconstruction; the source is not replaced by its mean.
  • ExactDistortionAchievable: PDF p.16: the printed operational convention has rate slack but exact expected distortion at each sufficiently large block. It is kept here for the p.17 converse; no equality with Theorem 4.2's variational limit is asserted.
  • exactDistortionRate: PDF p.16: infimum of literal achievable nonnegative rates. The empty set has value infinity, so an impossible distortion target cannot become a zero rate.
  • QuantumRDChoice: PDF p.17: a finite-dimensional CQ channel X to C and a joint C B reconstruction POVM. The positive dimension is part of the datum and has NO upper bound. Every channel output is a density, even at source labels of weight zero.
  • quantumRDExpected: PDF p.17: expected one-letter distortion. Given X=x the quantum systems C and B are in the PRODUCT state W(x) tensor rho_x; the decoder is joint on CB.
  • quantumRDCost: PDF p.17: I(X;C)-I(C;B) in bits. The accepted label information and conditional-product average compute both information terms from the same source.
  • quantumRDValue: PDF p.17: optimize over ALL finite positive auxiliary dimensions and all joint reconstruction POVMs. Nonnegativity follows from the conditional-product source; ofReal does not clip valid costs. Finite attainment is a theorem conclusion.
  • redistributionReceiverInformation: PDF p.18: I(Bhat;B), computed in the same pure four-register source.
  • senderMutual: PDF p.18: I(Bhat;Ahat), computed by actual partial traces.
  • fqrsQuantum: PDF p.18: FQRS communication I(Bhat;RB)/2, without dropping Bob's register.
  • printedQuantumCost: PDF p.18, Proposition 5.1: printed I(Bhat;R|Ahat)/2. For a pure source this equals the accepted quantumCost I(Bhat;R|B)/2 by complementary entropy duality; retaining the printed conditioning avoids silently changing the result.
  • diagonalCQSource: PDF pp.17,19: classical side information as diagonal conditional densities. The channel Z(.|x) is normalized even on zero-probability source labels.
  • CQFeedbackCode.IsClassical: PDF p.19: diagonal POVMs are classical stochastic decoding of the retained Z^n. Their square-root instruments leave each classical basis label intact.
  • ClassicalSimulationAchievable: PDF p.19: fully classical version of the Eq.(6) task. Diagonal instruments preserve the original Z^n; X^n and both Y^n outputs are still in the error test.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_MatrixDistances
import Definitions.Def_QAlgorithms_Papers_Independent_Claimed50_P40
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P07
import Definitions.Def_QAlgorithms_Papers_Independent_Gap48_P20

/-!
Luo--Devetak, arXiv:quant-ph/0611008v1. Entropies and resource rates are in bits.
The accepted paper modules imported here are frozen, reviewed library interfaces.
Distinct source-specific definitions are paper-local; the accepted Luders branch
is reused through finite-basis relabelling. The simulation relation is computed from an
encoder, a square-root measurement instrument and two feedback outputs; it is
not defined by its proposed entropy region. The two source-ambiguity deferrals and the
exact Theorem 4.1 reuse are recorded in the unit artifacts. Compilation pending.
-/

noncomputable section

namespace QAlgorithms.Papers.Independent.Intakethird100.P42

open scoped BigOperators ComplexConjugate ENNReal
open Filter
open QAlgorithms.Papers.Independent.Dependency50.P07
open QAlgorithms.Papers.Independent.Gap48.P20
open QAlgorithms.Papers.Independent.Information100.P27 (IsProbability shannon)

attribute [local instance] Classical.decEq Classical.propDecidable

/-- PDF p.4: classical blocks X^n,Y_A^n,Y_B^n and the retained quantum B^n. -/
abbrev SimulationBasis (X B : Type) (y n : ℕ) :=
  Copies X n × (Copies (Fin y) n × (Copies (Fin y) n × Copies B n))

/-- PDF p.4: one forward classical message, uniform independent shared bits,
private stochastic encoding, a coin/message-controlled POVM, and one deterministic
output function used by both parties. The POVM is implemented by its canonical
square-root instrument, as required by the gentle measurement paragraph on p.12. -/
structure CQFeedbackCode (X B : Type) [Fintype B] (y n : ℕ) where
  coinBits : ℕ
  messageBits : ℕ
  privateBits : ℕ
  encode : Copies X n → Qubits coinBits → Qubits messageBits → Qubits privateBits → ℝ
  encode_nonnegative : ∀ x l m s, 0 ≤ encode x l m s
  encode_normalized : ∀ x l, (∑ m, ∑ s, encode x l m s) = 1
  measure : Qubits coinBits → Qubits messageBits → Qubits privateBits →
    Matrix (Copies B n) (Copies B n) ℂ
  measure_valid : ∀ l m, IsPOVM (measure l m)
  decode : Qubits coinBits → Qubits messageBits → Qubits privateBits → Copies (Fin y) n

/-- PDF pp.4,12, Eq.(6): reuse the accepted square-root instrument on a finite
basis. The equivalence only relabels B as Fin (card B) and back; no coordinate or
measurement outcome is discarded, and the branch remains unnormalized. -/
def feedbackMeasurementBranch {B : Type} [Fintype B]
    (effect rho : Matrix B B ℂ) : Matrix B B ℂ :=
  let e := Fintype.equivFin B
  (QAlgorithms.Papers.Independent.Claimed50.P40.ludersBranch
    (fun _ : Unit => effect.submatrix e.symm e.symm)
    (rho.submatrix e.symm e.symm) ()).submatrix e e

/-- PDF pp.3--4, Eqs.(2),(4),(5): the memoryless feedback target, including X and B. -/
def simulationTarget {X B : Type} (E : CQData X B) (V : Auxiliary X) (n : ℕ) :
    Matrix (SimulationBasis X B V.size n) (SimulationBasis X B V.size n) ℂ :=
  fun a b =>
    if a.1 = b.1 ∧ a.2.1 = b.2.1 ∧ a.2.2.1 = b.2.2.1 ∧ a.2.1 = a.2.2.1 then
      ((wordProb E.prob n a.1 * ∏ i, V.channel (a.1 i) (a.2.1 i)) : ℂ) *
        wordState E n a.1 a.2.2.2 b.2.2.2
    else 0

/-- PDF p.4, Eq.(6): actual joint output, averaged over the independent uniform
coin and Alice's stochastic record. Bob's record is measured, not supplied by Alice. -/
def simulationOutput {X B : Type} [Fintype B] (E : CQData X B) {y n : ℕ}
    (P : CQFeedbackCode X B y n) :
    Matrix (SimulationBasis X B y n) (SimulationBasis X B y n) ℂ :=
  fun a b => if a.1 = b.1 ∧ a.2.1 = b.2.1 ∧ a.2.2.1 = b.2.2.1 then
    ∑ l : Qubits P.coinBits, ∑ m : Qubits P.messageBits,
      ∑ s : Qubits P.privateBits, ∑ t : Qubits P.privateBits,
        if P.decode l m s = a.2.1 ∧ P.decode l m t = a.2.2.1 then
          ((wordProb E.prob n a.1 * ((2 : ℝ) ^ P.coinBits)⁻¹ *
              P.encode a.1 l m s) : ℂ) *
            feedbackMeasurementBranch (P.measure l m t) (wordState E n a.1) a.2.2.2 b.2.2.2
        else 0
    else 0

/-- PDF p.4, Eq.(6): full trace norm, without the trace-distance factor 1/2. -/
def feedbackTraceError {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (V : Auxiliary X) {n : ℕ} (P : CQFeedbackCode X B V.size n) : ℝ :=
  traceNorm (simulationTarget E V n - simulationOutput E P)

/-- PDF pp.4--5: operational achievability. Integer bit budgets replace nR,nC;
rounding is covered by every positive rate slack. Codes precede the sampled input. -/
def SimulationAchievable {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (V : Auxiliary X) (R C : ℝ) : Prop :=
  0 ≤ R ∧ 0 ≤ C ∧ ∀ eps > (0 : ℝ), ∀ delta > (0 : ℝ), ∃ N : ℕ,
    ∀ n ≥ N, 1 ≤ n → ∃ P : CQFeedbackCode X B V.size n,
      feedbackTraceError E V P ≤ eps ∧ (P.messageBits : ℝ) ≤ n * (R + delta) ∧
        (P.coinBits : ℝ) ≤ n * (C + delta)

/-- PDF p.14, Eq.(36): source-relative resource simulation, with identities
initially distributing X to Alice and B to Bob and a feedback W at the output.
This is the same physical Eq.(6) task, not universal arbitrary-input simulation. -/
def SourceSimulation {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (V : Auxiliary X) (R C : ℝ) : Prop :=
  SimulationAchievable E V R C

example {B : Type} [Fintype B] (rho : Matrix B B ℂ) :
    feedbackMeasurementBranch 0 rho = 0 := by
  simp [feedbackMeasurementBranch,
    QAlgorithms.Papers.Independent.Claimed50.P40.ludersBranch, cfc_zero]

example {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (V : Auxiliary X) (R C : ℝ) :
    SourceSimulation E V R C ↔ SimulationAchievable E V R C := Iff.rfl

/-- PDF pp.6--9: probability of an event of a mutually independent q^n codebook.
The codebook is fixed before a uniform index is sampled to define its empirical law. -/
def bookProbability {Y : Type} [Fintype Y] (q : Y → ℝ) (n M : ℕ)
    (event : (Fin M → Copies Y n) → Prop) : ℝ :=
  ∑ book : Fin M → Copies Y n,
    if event book then ∏ i, wordProb q n (book i) else 0

/-- PDF p.6: the conditional law of G(W) for a fixed codebook, with W uniform on [M]. -/
def empiricalLaw {Y : Type} [Fintype Y] {M : ℕ} (book : Fin M → Y) (y : Y) : ℝ :=
  (∑ i, if book i = y then (1 : ℝ) else 0) / M

/-- PDF p.6: full classical l1 error; there is no factor 1/2. -/
def dilutionError {Y : Type} [Fintype Y] (q : Y → ℝ) (n M : ℕ)
    (book : Fin M → Copies Y n) : ℝ :=
  ∑ y, |wordProb q n y - empiricalLaw book y|

/-- PDF p.7: a fully normalized stochastic encoder, including zero-marginal rows. -/
abbrev CoverEncoder (X : Type) [Fintype X] (M : ℕ) :=
  {e : X → Fin M → ℝ // ∀ x, IsProbability (e x)}

/-- PDF p.7: p(x) is derived from q(y)P(x|y), not a free marginal parameter. -/
def coveringMarginal {X Y : Type} [Fintype Y] (q : Y → ℝ) (P : Y → X → ℝ) (x : X) : ℝ :=
  ∑ y, q y * P y x

/-- PDF pp.7,9: full l1 comparison on [M] times X, not merely an output marginal. -/
def coveringError {X Y : Type} [Fintype X] [Fintype Y]
    (q : Y → ℝ) (P : Y → X → ℝ) (M : ℕ) (book : Fin M → Y)
    (e : CoverEncoder X M) : ℝ :=
  ∑ i, ∑ x, |P (book i) x / M - coveringMarginal q P x * e.val x i|

/-- PDF p.9: the conditional product channel P^n(x^n|y^n). -/
def productConditional {X Y : Type} (P : Y → X → ℝ) (n : ℕ)
    (y : Copies Y n) (x : Copies X n) : ℝ := ∏ i, P (y i) (x i)

example {Y : Type} [Fintype Y] (q : Y → ℝ) (n M : ℕ) :
    bookProbability q n M (fun _ => False) = 0 := by
  simp [bookProbability]

example {X Y : Type} (P : Y → X → ℝ) (x : Copies X 0) (y : Copies Y 0) :
    productConditional P 0 y x = 1 := by
  simp [productConditional]

example {Y : Type} [Fintype Y] [Subsingleton Y] [Nonempty Y]
    (M : ℕ) (hM : 0 < M) (book : Fin M → Y) (y : Y) :
    empiricalLaw book y = 1 := by
  have hbook : ∀ i, book i = y := fun _ => Subsingleton.elim _ _
  simp [empiricalLaw, hbook, Nat.ne_of_gt hM]

example {X : Type} [Fintype X] : Nonempty (CoverEncoder X 1) := by
  refine ⟨⟨fun _ _ => 1, ?_⟩⟩
  intro x
  exact ⟨fun _ => by norm_num, by simp⟩

/-- PDF pp.15--16: finite-block variational datum, over all finite positive
classical auxiliary alphabets and all message-controlled reconstruction POVMs. -/
structure RDDatum (X B : Type) [Fintype X] [Fintype B] (n : ℕ) where
  auxiliary : Auxiliary (Copies X n)
  decoder : Fin auxiliary.size → Copies X n → Matrix (Copies B n) (Copies B n) ℂ
  decoder_valid : ∀ y, IsPOVM (decoder y)

/-- PDF p.15: additive per-letter distortion; theorems only use positive n. -/
def wordDistortion {X : Type} (d : X → X → ℝ) (n : ℕ)
    (x z : Copies X n) : ℝ := (∑ i, d (x i) (z i)) / n

/-- PDF p.16: expected distortion includes the decoder's actual Born probability. -/
def rdExpected {X B : Type} [Fintype X] [Fintype B] (E : CQData X B)
    (d : X → X → ℝ) {n : ℕ} (P : RDDatum X B n) : ℝ :=
  ∑ x, ∑ y, ∑ z, wordProb E.prob n x * P.auxiliary.channel x y *
    effectProb (P.decoder y z) (wordState E n x) * wordDistortion d n x z

/-- PDF p.16: finite-block variational value in extended nonnegative reals.
Empty feasible sets have value infinity. Finite attainment is asserted below,
not built into a hypothesis or lost by use of real sInf's empty-set default. -/
def blockRD {X B : Type} [Fintype X] [Fintype B] (E : CQData X B)
    (d : X → X → ℝ) (D : ℝ) (n : ℕ) : ℝ≥0∞ :=
  sInf {r | ∃ P : RDDatum X B n, rdExpected E d P ≤ D ∧
    r = ENNReal.ofReal (auxCost (iidCQ E n) P.auxiliary / n)}

example {X : Type} (d : X → X → ℝ) (x z : Copies X 0) :
    wordDistortion d 0 x z = 0 := by simp [wordDistortion]

/-- PDF pp.15--16: a deterministic encoder and a message-controlled reconstruction
POVM. Outcomes have already been coarsened to the reproduced word. No common
randomness, entanglement, or additional communication is supplied. -/
structure ExactDistortionCode (X B : Type) [Fintype X] [Fintype B] (n : ℕ) where
  messageBits : ℕ
  encode : Copies X n → Qubits messageBits
  decoder : Qubits messageBits → Copies X n → Matrix (Copies B n) (Copies B n) ℂ
  decoder_valid : ∀ m, IsPOVM (decoder m)

/-- PDF p.16: literal finite-block expected distortion, including the actual
Born probabilities of Bob's reconstruction; the source is not replaced by its mean. -/
def exactDistortionExpected {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (d : X → X → ℝ) {n : ℕ} (P : ExactDistortionCode X B n) : ℝ :=
  ∑ x, ∑ z, wordProb E.prob n x *
    effectProb (P.decoder (P.encode x) z) (wordState E n x) * wordDistortion d n x z

/-- PDF p.16: the printed operational convention has rate slack but exact
expected distortion at each sufficiently large block. It is kept here for the
p.17 converse; no equality with Theorem 4.2's variational limit is asserted. -/
def ExactDistortionAchievable {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (d : X → X → ℝ) (D R : ℝ) : Prop :=
  0 ≤ R ∧ ∀ delta > (0 : ℝ), ∃ N : ℕ, ∀ n ≥ N, 1 ≤ n →
    ∃ P : ExactDistortionCode X B n,
      exactDistortionExpected E d P ≤ D ∧ (P.messageBits : ℝ) ≤ n * (R + delta)

/-- PDF p.16: infimum of literal achievable nonnegative rates. The empty set
has value infinity, so an impossible distortion target cannot become a zero rate. -/
def exactDistortionRate {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (d : X → X → ℝ) (D : ℝ) : ℝ≥0∞ :=
  sInf {r | ∃ R : ℝ, ExactDistortionAchievable E d D R ∧ r = ENNReal.ofReal R}

/-- PDF p.17: a finite-dimensional CQ channel X to C and a joint C B
reconstruction POVM. The positive dimension is part of the datum and has NO
upper bound. Every channel output is a density, even at source labels of weight zero. -/
structure QuantumRDChoice (X B : Type) [Fintype X] [Fintype B] where
  dimension : ℕ
  positive : 0 < dimension
  output : X → Matrix (Fin dimension) (Fin dimension) ℂ
  output_valid : ∀ x, IsDensityMatrix (output x)
  decoder : X → Matrix (Fin dimension × B) (Fin dimension × B) ℂ
  decoder_valid : IsPOVM decoder

/-- PDF p.17: expected one-letter distortion. Given X=x the quantum systems
C and B are in the PRODUCT state W(x) tensor rho_x; the decoder is joint on CB. -/
def quantumRDExpected {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (d : X → X → ℝ) (P : QuantumRDChoice X B) : ℝ :=
  ∑ x, ∑ z, E.prob x *
    effectProb (P.decoder z) (Matrix.kronecker (P.output x) (E.state x)) * d x z

/-- PDF p.17: I(X;C)-I(C;B) in bits. The accepted label information and
conditional-product average compute both information terms from the same source. -/
def quantumRDCost {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (P : QuantumRDChoice X B) : ℝ :=
  QAlgorithms.Papers.Independent.Information100.P27.labelInformation E.prob P.output -
    (vnEntropy (QAlgorithms.Papers.Independent.Information100.P27.average E.prob P.output) +
      vnEntropy (ensembleMean E) -
      vnEntropy (QAlgorithms.Papers.Independent.Information100.P27.memoryEnvironment
        E.prob E.state P.output))

/-- PDF p.17: optimize over ALL finite positive auxiliary dimensions and all
joint reconstruction POVMs. Nonnegativity follows from the conditional-product
source; ofReal does not clip valid costs. Finite attainment is a theorem conclusion. -/
def quantumRDValue {X B : Type} [Fintype X] [Fintype B]
    (E : CQData X B) (d : X → X → ℝ) (D : ℝ) : ℝ≥0∞ :=
  sInf {r | ∃ P : QuantumRDChoice X B, quantumRDExpected E d P ≤ D ∧
    r = ENNReal.ofReal (quantumRDCost E P)}

example {X B : Type} [Fintype X] [Fintype B] (E : CQData X B)
    (P : ExactDistortionCode X B 0) (d : X → X → ℝ) :
    exactDistortionExpected E d P = 0 := by
  simp [exactDistortionExpected, wordDistortion]

example {X B : Type} [Fintype X] [Fintype B] (E : CQData X B)
    (P : QuantumRDChoice X B) : quantumRDExpected E (fun _ _ => 0) P = 0 := by
  simp [quantumRDExpected]

/-- PDF p.18: I(Bhat;B), computed in the same pure four-register source. -/
def redistributionReceiverInformation {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
    (psi : EuclideanSpace ℂ (RedistributionSource A C B R)) : ℝ :=
  subsystemEntropy (pureDensity psi) (fun p => p.2.2.1)
      (fun p => (p.1, p.2.1, p.2.2.2)) +
    subsystemEntropy (pureDensity psi) (fun p => p.2.2.2)
      (fun p => (p.1, p.2.1, p.2.2.1)) -
    subsystemEntropy (pureDensity psi) (fun p => p.2.2) (fun p => (p.1, p.2.1))

/-- PDF p.18: I(Bhat;Ahat), computed by actual partial traces. -/
def senderMutual {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
    (psi : EuclideanSpace ℂ (RedistributionSource A C B R)) : ℝ :=
  subsystemEntropy (pureDensity psi) (fun p => p.2.2.1)
      (fun p => (p.1, p.2.1, p.2.2.2)) +
    subsystemEntropy (pureDensity psi) (fun p => p.2.1) (fun p => (p.1, p.2.2)) -
    subsystemEntropy (pureDensity psi) (fun p => (p.2.1, p.2.2.1))
      (fun p => (p.1, p.2.2.2))

/-- PDF p.18: FQRS communication I(Bhat;RB)/2, without dropping Bob's register. -/
def fqrsQuantum {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
    (psi : EuclideanSpace ℂ (RedistributionSource A C B R)) : ℝ :=
  (subsystemEntropy (pureDensity psi) (fun p => p.2.2.1)
      (fun p => (p.1, p.2.1, p.2.2.2)) +
    subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.2.2))
      (fun p => (p.2.1, p.2.2.1)) -
    subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.2))
      (fun p => p.2.1)) / 2

/-- PDF p.18, Proposition 5.1: printed I(Bhat;R|Ahat)/2. For a pure source
this equals the accepted quantumCost I(Bhat;R|B)/2 by complementary entropy
duality; retaining the printed conditioning avoids silently changing the result. -/
def printedQuantumCost {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
    (psi : EuclideanSpace ℂ (RedistributionSource A C B R)) : ℝ :=
  (subsystemEntropy (pureDensity psi) (fun p => (p.2.1, p.2.2.1))
      (fun p => (p.1, p.2.2.2)) +
    subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.1))
      (fun p => p.2.2) -
    subsystemEntropy (pureDensity psi) (fun p => p.2.1)
      (fun p => (p.1, p.2.2)) -
    subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.1, p.2.2.1))
      (fun p => p.2.2.2)) / 2

example {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
    (psi : EuclideanSpace ℂ (RedistributionSource A C B R)) :
    fqrsQuantum psi = (subsystemEntropy (pureDensity psi) (fun p => p.2.2.1)
      (fun p => (p.1, p.2.1, p.2.2.2)) +
      subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.2.2))
        (fun p => (p.2.1, p.2.2.1)) -
      subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.2))
        (fun p => p.2.1)) / 2 := rfl

/-- PDF pp.17,19: classical side information as diagonal conditional densities.
The channel Z(.|x) is normalized even on zero-probability source labels. -/
def diagonalCQSource {X : Type} (p : X → ℝ) (Z : Auxiliary X) : CQData X (Fin Z.size) where
  prob := p
  state := fun x => Matrix.diagonal (fun z => (Z.channel x z : ℂ))

/-- PDF p.19: diagonal POVMs are classical stochastic decoding of the retained
Z^n. Their square-root instruments leave each classical basis label intact. -/
def CQFeedbackCode.IsClassical {X B : Type} [Fintype B] {y n : ℕ}
    (P : CQFeedbackCode X B y n) : Prop :=
  ∀ l m s u v, u ≠ v → P.measure l m s u v = 0

/-- PDF p.19: fully classical version of the Eq.(6) task. Diagonal instruments
preserve the original Z^n; X^n and both Y^n outputs are still in the error test. -/
def ClassicalSimulationAchievable {X : Type} [Fintype X] (p : X → ℝ)
    (Z V : Auxiliary X) (R C : ℝ) : Prop :=
  0 ≤ R ∧ 0 ≤ C ∧ ∀ eps > (0 : ℝ), ∀ delta > (0 : ℝ), ∃ N : ℕ,
    ∀ n ≥ N, 1 ≤ n → ∃ P : CQFeedbackCode X (Fin Z.size) V.size n,
      P.IsClassical ∧ feedbackTraceError (diagonalCQSource p Z) V P ≤ eps ∧
        (P.messageBits : ℝ) ≤ n * (R + delta) ∧ (P.coinBits : ℝ) ≤ n * (C + delta)

example {X : Type} (p : X → ℝ) (Z : Auxiliary X) (x : X) (z : Fin Z.size) :
    (diagonalCQSource p Z).state x z z = (Z.channel x z : ℂ) := by
  simp [diagonalCQSource]

end QAlgorithms.Papers.Independent.Intakethird100.P42
Source
Channel simulation with quantum side information, arXiv:quant-ph/0611008v1, https://arxiv.org/abs/quant-ph/0611008v1
Read-back

What the Lean code literally says, in plain math · gpt-6.1-sol

QAlgorithms.Papers.Independent.Intakethird100.P42.CQFeedbackCode.IsClassical

Let XXX be any type, let BBB be a finite type, and let y,ny,ny,n be arbitrary natural numbers, including zero. Let PPP consist of natural numbers c,m,hc,m,hc,m,h, a function e(x,l,a,s)≥0e(x,l,a,s)\ge0e(x,l,a,s)≥0 for x∈Xnx\in X^nx∈Xn, l∈{0,1}cl\in\{0,1\}^cl∈{0,1}c, a∈{0,1}ma\in\{0,1\}^ma∈{0,1}m, and s∈{0,1}hs\in\{0,1\}^hs∈{0,1}h, satisfying ∑a,se(x,l,a,s)=1\sum_{a,s}e(x,l,a,s)=1∑a,s​e(x,l,a,s)=1 for every x,lx,lx,l; matrices Ml,a,sM_{l,a,s}Ml,a,s​ on CBn\mathbb C^{B^n}CBn, positive semidefinite and satisfying ∑sMl,a,s=I\sum_sM_{l,a,s}=I∑s​Ml,a,s​=I for every l,al,al,a; and a decoding function f(l,a,s)∈{0,…,y−1}nf(l,a,s)\in\{0,\ldots,y-1\}^nf(l,a,s)∈{0,…,y−1}n. The predicate defined here is

∀l,a,s, ∀u,v∈Bn,u≠v ⟹ Ml,a,s(u,v)=0.\forall l,a,s,\ \forall u,v\in B^n,\qquad u\ne v\ \Longrightarrow\ M_{l,a,s}(u,v)=0.∀l,a,s, ∀u,v∈Bn,u=v ⟹ Ml,a,s​(u,v)=0.

Thus it requires every one of the specified measurement matrices to be diagonal. When n=0n=0n=0, BnB^nBn has one element and this off-diagonal condition is vacuous. When y=0<ny=0<ny=0<n, no such PPP exists because its decoding function would have a nonempty domain and an empty codomain.

QAlgorithms.Papers.Independent.Intakethird100.P42.CQFeedbackCode

For an arbitrary type XXX, a finite type BBB with a specified finite enumeration, and natural numbers y,ny,ny,n, a code consists of natural numbers c,h,pc,h,pc,h,p, called its coin, message, and private bit counts; a real function e(x,ℓ,m,s)e(x,\ell,m,s)e(x,ℓ,m,s), where x∈X[n]x\in X^{[n]}x∈X[n], ℓ∈{0,1}c\ell\in\{0,1\}^cℓ∈{0,1}c, m∈{0,1}hm\in\{0,1\}^hm∈{0,1}h, and s∈{0,1}ps\in\{0,1\}^ps∈{0,1}p, satisfying

e(x,ℓ,m,s)≥0,∑m∑se(x,ℓ,m,s)=1for every x,ℓ;e(x,\ell,m,s)\ge0,\qquad \sum_m\sum_s e(x,\ell,m,s)=1 \quad\text{for every }x,\ell;e(x,ℓ,m,s)≥0,m∑​s∑​e(x,ℓ,m,s)=1for every x,ℓ;

matrices Fℓ,m,sF_{\ell,m,s}Fℓ,m,s​ on B[n]B^{[n]}B[n], such that for every ℓ,m\ell,mℓ,m, every Fℓ,m,sF_{\ell,m,s}Fℓ,m,s​ is positive semidefinite and ∑sFℓ,m,s=I\sum_sF_{\ell,m,s}=I∑s​Fℓ,m,s​=I; and a function d(ℓ,m,s)∈[y][n]d(\ell,m,s)\in[y]^{[n]}d(ℓ,m,s)∈[y][n]. Here [k]={0,…,k−1}[k]=\{0,\ldots,k-1\}[k]={0,…,k−1}, and words are functions on [n][n][n]. The type XXX need not be finite, and BBB may be empty. Zero bit counts are allowed and give singleton bitstring sets. Length n=0n=0n=0 is allowed and gives singleton word sets; if y=0<ny=0<ny=0<n, no such code exists because the required decoder would map a nonempty bitstring domain into an empty word type.

QAlgorithms.Papers.Independent.Intakethird100.P42.ClassicalSimulationAchievable

Let XXX be a finite type, possibly empty, let p:X→Rp:X\to\mathbb Rp:X→R be arbitrary, and let Z,VZ,VZ,V be auxiliary data with positive finite output sizes z,vz,vz,v, respectively, and conditional probabilities Z(b∣x)Z(b\mid x)Z(b∣x), V(u∣x)V(u\mid x)V(u∣x), each nonnegative and summing to 111 over its output for every xxx. For real numbers R,CR,CR,C, this predicate means R≥0R\ge0R≥0, C≥0C\ge0C≥0, and that for every ε,δ>0\varepsilon,\delta>0ε,δ>0 there is a natural NNN such that every n≥Nn\ge Nn≥N with n≥1n\ge1n≥1 admits the following data: natural numbers c,m,hc,m,hc,m,h; encoder weights e(x,l,a,s)≥0e(x,l,a,s)\ge0e(x,l,a,s)≥0, for x∈Xnx\in X^nx∈Xn, l∈{0,1}cl\in\{0,1\}^cl∈{0,1}c, a∈{0,1}ma\in\{0,1\}^ma∈{0,1}m, s∈{0,1}hs\in\{0,1\}^hs∈{0,1}h, satisfying ∑a,se(x,l,a,s)=1\sum_{a,s}e(x,l,a,s)=1∑a,s​e(x,l,a,s)=1; diagonal positive semidefinite matrices Ml,a,tM_{l,a,t}Ml,a,t​ on C{0,…,z−1}n\mathbb C^{\{0,\ldots,z-1\}^n}C{0,…,z−1}n, satisfying ∑tMl,a,t=I\sum_tM_{l,a,t}=I∑t​Ml,a,t​=I for every l,al,al,a; and a decoding function f(l,a,s)∈{0,…,v−1}nf(l,a,s)\in\{0,\ldots,v-1\}^nf(l,a,s)∈{0,…,v−1}n. To specify the required error, put pn(x)=∏ip(xi)p_n(x)=\prod_ip(x_i)pn​(x)=∏i​p(xi​), Zn(b∣x)=∏iZ(bi∣xi)Z_n(b\mid x)=\prod_iZ(b_i\mid x_i)Zn​(b∣x)=∏i​Z(bi​∣xi​), and Vn(u∣x)=∏iV(ui∣xi)V_n(u\mid x)=\prod_iV(u_i\mid x_i)Vn​(u∣x)=∏i​V(ui​∣xi​). On the basis indexed by (x,u,w,b)∈Xn×{0,…,v−1}n×{0,…,v−1}n×{0,…,z−1}n(x,u,w,b)\in X^n\times\{0,\ldots,v-1\}^n\times\{0,\ldots,v-1\}^n\times\{0,\ldots,z-1\}^n(x,u,w,b)∈Xn×{0,…,v−1}n×{0,…,v−1}n×{0,…,z−1}n, let TTT and OOO be diagonal matrices with entries

Tx,u,w,b=1{u=w} pn(x)Vn(u∣x)Zn(b∣x),T_{x,u,w,b}=\mathbf1_{\{u=w\}}\,p_n(x)V_n(u\mid x)Z_n(b\mid x),Tx,u,w,b​=1{u=w}​pn​(x)Vn​(u∣x)Zn​(b∣x), Ox,u,w,b=pn(x)2−cZn(b∣x)∑l,a,s,te(x,l,a,s) 1{f(l,a,s)=u}1{f(l,a,t)=w} Ml,a,t(b,b).O_{x,u,w,b} =p_n(x)2^{-c}Z_n(b\mid x) \sum_{l,a,s,t}e(x,l,a,s)\, \mathbf1_{\{f(l,a,s)=u\}}\mathbf1_{\{f(l,a,t)=w\}}\, M_{l,a,t}(b,b).Ox,u,w,b​=pn​(x)2−cZn​(b∣x)l,a,s,t∑​e(x,l,a,s)1{f(l,a,s)=u}​1{f(l,a,t)=w}​Ml,a,t​(b,b).

The required bounds are ∥T−O∥1≤ε\|T-O\|_1\le\varepsilon∥T−O∥1​≤ε, m≤n(R+δ)m\le n(R+\delta)m≤n(R+δ), and c≤n(C+δ)c\le n(C+\delta)c≤n(C+δ), where ∥M∥1=Re⁡tr⁡M†M\|M\|_1=\operatorname{Re}\operatorname{tr}\sqrt{M^\dagger M}∥M∥1​=RetrM†M​. No bound is imposed on hhh, and c,m,hc,m,hc,m,h may be zero. No sign or normalization condition is imposed on ppp; the displayed matrices and error are used even when ppp is not a probability law.

QAlgorithms.Papers.Independent.Intakethird100.P42.CoverEncoder

For a finite type XXX and a natural number MMM, a cover encoder is a function e:X×{0,…,M−1}→Re:X\times\{0,\ldots,M-1\}\to\mathbb Re:X×{0,…,M−1}→R satisfying, for every x∈Xx\in Xx∈X, e(x,i)≥0e(x,i)\ge0e(x,i)≥0 for every i<Mi<Mi<M and ∑i=0M−1e(x,i)=1\sum_{i=0}^{M-1}e(x,i)=1∑i=0M−1​e(x,i)=1. No nonemptiness of XXX or positivity of MMM is required in the definition: if XXX is empty these conditions are vacuous, whereas if XXX is nonempty and M=0M=0M=0, no cover encoder exists.

QAlgorithms.Papers.Independent.Intakethird100.P42.ExactDistortionAchievable

For finite, finitely enumerated types X,BX,BX,B, arbitrary data E=(p,ρ)E=(p,\rho)E=(p,ρ) with p:X→Rp:X\to\mathbb Rp:X→R and complex matrices ρx\rho_xρx​ on CB\mathbb C^BCB, any real-valued function d:X×X→Rd:X\times X\to\mathbb Rd:X×X→R, and real numbers D,RD,RD,R, this property means that R≥0R\ge0R≥0 and that, for every real δ>0\delta>0δ>0, there exists a natural number NNN such that every natural number n≥Nn\ge Nn≥N with n≥1n\ge1n≥1 admits a natural number bbb, an encoder e:Xn→{0,1}be:X^n\to\{0,1\}^be:Xn→{0,1}b, and matrices Mm,zM_{m,z}Mm,z​ on CBn\mathbb C^{B^n}CBn, positive semidefinite and summing to the identity over z∈Xnz\in X^nz∈Xn for each message mmm, satisfying

∑x,z∈Xn(∏i=0n−1p(xi))Re⁡tr⁡ ⁣(Me(x),zρx(n))∑i=0n−1d(xi,zi)n≤D,b≤n(R+δ),\sum_{x,z\in X^n} \left(\prod_{i=0}^{n-1}p(x_i)\right) \operatorname{Re}\operatorname{tr}\!\left(M_{e(x),z}\rho_x^{(n)}\right) \frac{\sum_{i=0}^{n-1}d(x_i,z_i)}{n}\le D, \qquad b\le n(R+\delta),x,z∈Xn∑​(i=0∏n−1​p(xi​))Retr(Me(x),z​ρx(n)​)n∑i=0n−1​d(xi​,zi​)​≤D,b≤n(R+δ),

where (ρx(n))u,v=∏i(ρxi)ui,vi(\rho_x^{(n)})_{u,v}=\prod_i(\rho_{x_i})_{u_i,v_i}(ρx(n)​)u,v​=∏i​(ρxi​​)ui​,vi​​. The code may depend on nnn and δ\deltaδ; the distortion bound is exactly DDD. There are no validity assumptions on EEE, no sign assumptions on d,Dd,Dd,D, and no nonemptiness assumptions on X,BX,BX,B. The block length n=0n=0n=0 is excluded from the code requirement, while zero message bits are allowed.

QAlgorithms.Papers.Independent.Intakethird100.P42.ExactDistortionCode

For finite types XXX and BBB, each equipped with a finite enumeration, and any natural number nnn, an exact distortion code consists of a natural number bbb of message bits, a function e:Xn→{0,1}be:X^n\to\{0,1\}^be:Xn→{0,1}b, and matrices Mm,zM_{m,z}Mm,z​ on CBn\mathbb C^{B^n}CBn for every message m∈{0,1}bm\in\{0,1\}^bm∈{0,1}b and reconstruction z∈Xnz\in X^nz∈Xn, where XnX^nXn and BnB^nBn mean functions from {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} into XXX and BBB, respectively. For every message mmm, the matrices satisfy Mm,z⪰0M_{m,z}\succeq0Mm,z​⪰0 for all zzz and ∑z∈XnMm,z=I\sum_{z\in X^n}M_{m,z}=I∑z∈Xn​Mm,z​=I; this requirement includes messages outside the image of eee. Neither nnn nor bbb must be positive, and neither type must be nonempty. At n=0n=0n=0, each word type has one element, including when its alphabet is empty.

QAlgorithms.Papers.Independent.Intakethird100.P42.QuantumRDChoice

For finite types X,BX,BX,B, each equipped with a finite enumeration, a quantum rate-distortion choice consists of a positive natural number aaa, matrices τx\tau_xτx​ on Ca\mathbb C^aCa for every x∈Xx\in Xx∈X, and matrices MzM_zMz​ on C{0,…,a−1}×B\mathbb C^{\{0,\ldots,a-1\}\times B}C{0,…,a−1}×B for every z∈Xz\in Xz∈X, satisfying

τx⪰0,tr⁡τx=1(x∈X),Mz⪰0(z∈X),∑z∈XMz=I.\tau_x\succeq0,\qquad \operatorname{tr}\tau_x=1 \quad(x\in X),\qquad M_z\succeq0\quad(z\in X),\qquad \sum_{z\in X}M_z=I.τx​⪰0,trτx​=1(x∈X),Mz​⪰0(z∈X),z∈X∑​Mz​=I.

Neither XXX nor BBB is required to be nonempty. If XXX is empty, the output-state conditions are vacuous, but the decoder condition still requires its empty sum to equal the identity; consequently, when BBB is nonempty, no such choice exists for empty XXX.

QAlgorithms.Papers.Independent.Intakethird100.P42.RDDatum

For finite types X,BX,BX,B and a natural number nnn, a rate-distortion datum consists of a positive natural number LLL, a real function V:Xn×{0,…,L−1}→RV:X^n\times\{0,\ldots,L-1\}\to\mathbb RV:Xn×{0,…,L−1}→R satisfying V(y∣x)≥0V(y\mid x)\ge0V(y∣x)≥0 and ∑y<LV(y∣x)=1\sum_{y<L}V(y\mid x)=1∑y<L​V(y∣x)=1 for every word x∈Xnx\in X^nx∈Xn, and matrices Dy,zD_{y,z}Dy,z​ on CBn\mathbb C^{B^n}CBn, indexed by y<Ly<Ly<L and z∈Xnz\in X^nz∈Xn, such that for each yyy, every Dy,zD_{y,z}Dy,z​ is positive semidefinite and ∑z∈XnDy,z=I\sum_{z\in X^n}D_{y,z}=I∑z∈Xn​Dy,z​=I. Here XnX^nXn and BnB^nBn are the functions from {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} into XXX and BBB, respectively. The definition allows n=0n=0n=0 and empty types X,BX,BX,B, subject to these normalization and matrix conditions; for example, when XnX^nXn is empty and BnB^nBn is nonempty, the required empty sum cannot equal the identity, so no such datum exists.

QAlgorithms.Papers.Independent.Intakethird100.P42.SimulationAchievable

For finite types X,BX,BX,B with specified finite enumerations, arbitrary data px∈R,ρx∈CB×Bp_x\in\mathbb R,\rho_x\in\mathbb C^{B\times B}px​∈R,ρx​∈CB×B, an auxiliary system of positive size yyy with V(u∣x)≥0V(u\mid x)\ge0V(u∣x)≥0 and ∑uV(u∣x)=1\sum_uV(u\mid x)=1∑u​V(u∣x)=1, and real numbers R,CR,CR,C, this property means that R,C≥0R,C\ge0R,C≥0 and, for every ε,δ>0\varepsilon,\delta>0ε,δ>0, there is N∈NN\in\mathbb NN∈N such that every n≥Nn\ge Nn≥N with n≥1n\ge1n≥1 admits a code with coin, message, and private bit counts c,h,pc,h,pc,h,p satisfying

∥T−O∥1≤ε,h≤n(R+δ),c≤n(C+δ).\|T-O\|_1\le\varepsilon,\qquad h\le n(R+\delta),\qquad c\le n(C+\delta).∥T−O∥1​≤ε,h≤n(R+δ),c≤n(C+δ).

The code has nonnegative weights e(x,ℓ,m,s)e(x,\ell,m,s)e(x,ℓ,m,s) summing to 111 over m,sm,sm,s for every x,ℓx,\ellx,ℓ, decoded words d(ℓ,m,s)∈[y][n]d(\ell,m,s)\in[y]^{[n]}d(ℓ,m,s)∈[y][n], and matrices Fℓ,m,t⪰0F_{\ell,m,t}\succeq0Fℓ,m,t​⪰0 on B[n]B^{[n]}B[n] summing to III over ttt for every ℓ,m\ell,mℓ,m; ℓ,m,s,t\ell,m,s,tℓ,m,s,t are bitstrings of lengths c,h,p,pc,h,p,pc,h,p,p. With pn(x)=∏ipxip_n(x)=\prod_i p_{x_i}pn​(x)=∏i​pxi​​ and (ρx(n))b,b′=∏i(ρxi)bi,bi′(\rho_x^{(n)})_{b,b'}=\prod_i(\rho_{x_i})_{b_i,b'_i}(ρx(n)​)b,b′​=∏i​(ρxi​​)bi​,bi′​​, the matrices on tuples (x,u,v,b)(x,u,v,b)(x,u,v,b) have entries

T(x,u,v,b),(x′,u′,v′,b′)=1{x=x′, u=u′, v=v′, u=v} pn(x)∏iV(ui∣xi)(ρx(n))b,b′,O(x,u,v,b),(x′,u′,v′,b′)=1{x=x′, u=u′, v=v′}∑ℓ,m,s,t1{d(ℓ,m,s)=u, d(ℓ,m,t)=v} pn(x)2−ce(x,ℓ,m,s)[LFℓ,m,t(ρx(n))]b,b′.\begin{aligned} T_{(x,u,v,b),(x',u',v',b')} &=\mathbf1_{\{x=x',\,u=u',\,v=v',\,u=v\}}\, p_n(x)\prod_iV(u_i\mid x_i)(\rho_x^{(n)})_{b,b'},\\ O_{(x,u,v,b),(x',u',v',b')} &=\mathbf1_{\{x=x',\,u=u',\,v=v'\}} \sum_{\ell,m,s,t}\mathbf1_{\{d(\ell,m,s)=u,\,d(\ell,m,t)=v\}}\, p_n(x)2^{-c}e(x,\ell,m,s) [\mathcal L_{F_{\ell,m,t}}(\rho_x^{(n)})]_{b,b'}. \end{aligned}T(x,u,v,b),(x′,u′,v′,b′)​O(x,u,v,b),(x′,u′,v′,b′)​​=1{x=x′,u=u′,v=v′,u=v}​pn​(x)i∏​V(ui​∣xi​)(ρx(n)​)b,b′​,=1{x=x′,u=u′,v=v′}​ℓ,m,s,t∑​1{d(ℓ,m,s)=u,d(ℓ,m,t)=v}​pn​(x)2−ce(x,ℓ,m,s)[LFℓ,m,t​​(ρx(n)​)]b,b′​.​

Here LF(σ)\mathcal L_F(\sigma)LF​(σ) relabels onto a finite ordinal, forms cfc⁡(⋅,F)σcfc⁡(⋅,F)\operatorname{cfc}(\sqrt{\cdot},F)\sigma\operatorname{cfc}(\sqrt{\cdot},F)cfc(⋅​,F)σcfc(⋅​,F), and relabels back; ∥Z∥1=Re⁡tr⁡cfc⁡(⋅,Z†Z)\|Z\|_1=\operatorname{Re}\operatorname{tr}\operatorname{cfc}(\sqrt{\cdot},Z^\dagger Z)∥Z∥1​=Retrcfc(⋅​,Z†Z). The private bit count is unrestricted. No validity assumption on the ensemble data or nonemptiness assumption on either finite type is part of this definition; only positive block lengths are tested.

QAlgorithms.Papers.Independent.Intakethird100.P42.SimulationBasis

For arbitrary types X,BX,BX,B and natural numbers y,ny,ny,n, this is the type

X[n]×([y][n]×([y][n]×B[n])),X^{[n]}\times\bigl([y]^{[n]}\times([y]^{[n]}\times B^{[n]})\bigr),X[n]×([y][n]×([y][n]×B[n])),

where [k]={0,…,k−1}[k]=\{0,\ldots,k-1\}[k]={0,…,k−1} and Z[n]Z^{[n]}Z[n] denotes functions from [n][n][n] to ZZZ. Thus an element consists of an XXX-word, two separate [y][y][y]-words, and a BBB-word, all of length nnn. No finiteness or nonemptiness assumptions on XXX or BBB are imposed. When n=0n=0n=0, every word type here is a singleton, even if its alphabet is empty; when n>0n>0n>0, an empty alphabet makes the corresponding word type empty.

QAlgorithms.Papers.Independent.Intakethird100.P42.SourceSimulation

For finite types X,BX,BX,B with specified finite enumerations, arbitrary real weights pxp_xpx​, complex matrices ρx\rho_xρx​ on BBB, an auxiliary system of positive size yyy with probabilities V(u∣x)V(u\mid x)V(u∣x), and real numbers R,CR,CR,C, this is the following property: R,C≥0R,C\ge0R,C≥0, and for every ε,δ>0\varepsilon,\delta>0ε,δ>0, all sufficiently large positive integers nnn admit natural bit counts c,h,pc,h,pc,h,p, encoding weights e(x,ℓ,m,s)≥0e(x,\ell,m,s)\ge0e(x,ℓ,m,s)≥0 summing to 111 over m,sm,sm,s, a decoder d(ℓ,m,s)∈[y][n]d(\ell,m,s)\in[y]^{[n]}d(ℓ,m,s)∈[y][n], and matrices Fℓ,m,t⪰0F_{\ell,m,t}\succeq0Fℓ,m,t​⪰0 on B[n]B^{[n]}B[n] summing to III over ttt, such that h≤n(R+δ)h\le n(R+\delta)h≤n(R+δ), c≤n(C+δ)c\le n(C+\delta)c≤n(C+δ), and ∥T−O∥1≤ε\|T-O\|_1\le\varepsilon∥T−O∥1​≤ε. Here x∈X[n]x\in X^{[n]}x∈X[n], the bitstrings ℓ,m,s,t\ell,m,s,tℓ,m,s,t have lengths c,h,p,pc,h,p,pc,h,p,p, pn(x)=∏ipxip_n(x)=\prod_i p_{x_i}pn​(x)=∏i​pxi​​, (ρx(n))b,b′=∏i(ρxi)bi,bi′(\rho_x^{(n)})_{b,b'}=\prod_i(\rho_{x_i})_{b_i,b'_i}(ρx(n)​)b,b′​=∏i​(ρxi​​)bi​,bi′​​, and

T(x,u,v,b),(x′,u′,v′,b′)=1{x=x′, u=u′, v=v′, u=v} pn(x)∏iV(ui∣xi)(ρx(n))b,b′,O(x,u,v,b),(x′,u′,v′,b′)=1{x=x′, u=u′, v=v′}∑ℓ,m,s,t1{d(ℓ,m,s)=u, d(ℓ,m,t)=v} pn(x)2−ce(x,ℓ,m,s)[LFℓ,m,t(ρx(n))]b,b′.\begin{aligned} T_{(x,u,v,b),(x',u',v',b')} &=\mathbf1_{\{x=x',\,u=u',\,v=v',\,u=v\}}\, p_n(x)\prod_iV(u_i\mid x_i)(\rho_x^{(n)})_{b,b'},\\ O_{(x,u,v,b),(x',u',v',b')} &=\mathbf1_{\{x=x',\,u=u',\,v=v'\}} \sum_{\ell,m,s,t}\mathbf1_{\{d(\ell,m,s)=u,\,d(\ell,m,t)=v\}}\, p_n(x)2^{-c}e(x,\ell,m,s) [\mathcal L_{F_{\ell,m,t}}(\rho_x^{(n)})]_{b,b'}. \end{aligned}T(x,u,v,b),(x′,u′,v′,b′)​O(x,u,v,b),(x′,u′,v′,b′)​​=1{x=x′,u=u′,v=v′,u=v}​pn​(x)i∏​V(ui​∣xi​)(ρx(n)​)b,b′​,=1{x=x′,u=u′,v=v′}​ℓ,m,s,t∑​1{d(ℓ,m,s)=u,d(ℓ,m,t)=v}​pn​(x)2−ce(x,ℓ,m,s)[LFℓ,m,t​​(ρx(n)​)]b,b′​.​

The matrix indices range over X[n]×[y][n]×[y][n]×B[n]X^{[n]}\times[y]^{[n]}\times[y]^{[n]}\times B^{[n]}X[n]×[y][n]×[y][n]×B[n]; LF(σ)\mathcal L_F(\sigma)LF​(σ) relabels onto a finite ordinal, forms cfc⁡(⋅,F)σcfc⁡(⋅,F)\operatorname{cfc}(\sqrt{\cdot},F)\sigma\operatorname{cfc}(\sqrt{\cdot},F)cfc(⋅​,F)σcfc(⋅​,F), and relabels back, while ∥Z∥1=Re⁡tr⁡cfc⁡(⋅,Z†Z)\|Z\|_1=\operatorname{Re}\operatorname{tr}\operatorname{cfc}(\sqrt{\cdot},Z^\dagger Z)∥Z∥1​=Retrcfc(⋅​,Z†Z). There is no bound on ppp, no ensemble-validity hypothesis, and no nonemptiness hypothesis on XXX or BBB.

QAlgorithms.Papers.Independent.Intakethird100.P42.blockRD

For finite types X,BX,BX,B, arbitrary data E=(p,ρ)E=(p,\rho)E=(p,ρ) with p:X→Rp:X\to\mathbb Rp:X→R and complex matrices ρa\rho_aρa​ on BBB, an arbitrary real function d:X×X→Rd:X\times X\to\mathbb Rd:X×X→R, a real threshold DDD, and a natural number nnn, this definition takes an infimum in [0,∞][0,\infty][0,∞] over the following quantities. Choose any L>0L>0L>0, nonnegative conditional weights V(u∣x)V(u\mid x)V(u∣x) with ∑u<LV(u∣x)=1\sum_{u<L}V(u\mid x)=1∑u<L​V(u∣x)=1 for every x∈Xnx\in X^nx∈Xn, and positive semidefinite matrices Du,zD_{u,z}Du,z​ on BnB^nBn with ∑z∈XnDu,z=I\sum_{z\in X^n}D_{u,z}=I∑z∈Xn​Du,z​=I, subject to

∑x∈Xn∑u<L∑z∈Xnpn(x)V(u∣x)Re⁡tr⁡(Du,zρx(n))∑i<nd(xi,zi)n≤D,\sum_{x\in X^n}\sum_{u<L}\sum_{z\in X^n} p_n(x)V(u\mid x)\operatorname{Re}\operatorname{tr}(D_{u,z}\rho_x^{(n)}) \frac{\sum_{i<n}d(x_i,z_i)}{n}\le D,x∈Xn∑​u<L∑​z∈Xn∑​pn​(x)V(u∣x)Retr(Du,z​ρx(n)​)n∑i<n​d(xi​,zi​)​≤D,

where pn(x)=∏i<np(xi)p_n(x)=\prod_{i<n}p(x_i)pn​(x)=∏i<n​p(xi​) and (ρx(n))a,b=∏i<n(ρxi)ai,bi(\rho_x^{(n)})_{a,b}=\prod_{i<n}(\rho_{x_i})_{a_i,b_i}(ρx(n)​)a,b​=∏i<n​(ρxi​​)ai​,bi​​. For each such choice set J(u,x)=pn(x)V(u∣x)J(u,x)=p_n(x)V(u\mid x)J(u,x)=pn​(x)V(u∣x), r(u)=∑xJ(u,x)r(u)=\sum_xJ(u,x)r(u)=∑x​J(u,x), ρˉ=∑xpn(x)ρx(n)\bar\rho=\sum_xp_n(x)\rho_x^{(n)}ρˉ​=∑x​pn​(x)ρx(n)​, and σu=r(u)−1∑xJ(u,x)ρx(n)\sigma_u=r(u)^{-1}\sum_xJ(u,x)\rho_x^{(n)}σu​=r(u)−1∑x​J(u,x)ρx(n)​. The quantity entered into the infimum is

max⁡ ⁣{0, H(pn)+H(r)−H(J)−S(ρˉ)+∑u<Lr(u)S(σu)n},\max\!\left\{0,\, \frac{H(p_n)+H(r)-H(J)-S(\bar\rho)+\sum_{u<L}r(u)S(\sigma_u)}{n}\right\},max{0,nH(pn​)+H(r)−H(J)−S(ρˉ​)+∑u<L​r(u)S(σu​)​},

embedded as a finite extended nonnegative real, where H(a)=−∑za(z)log⁡2a(z)H(a)=-\sum_za(z)\log_2a(z)H(a)=−∑z​a(z)log2​a(z) and S(A)=−Re⁡tr⁡(f(A))S(A)=-\operatorname{Re}\operatorname{tr}(f(A))S(A)=−Retr(f(A)) for f(t)=tlog⁡2tf(t)=t\log_2tf(t)=tlog2​t, evaluated by continuous functional calculus. No probability or density assumptions on EEE, or sign assumptions on d,Dd,Dd,D, are imposed. All divisions are total, so σu=0\sigma_u=0σu​=0 if r(u)=0r(u)=0r(u)=0, and the displayed quotient is 000 when n=0n=0n=0. The infimum of an empty feasible set is +∞+\infty+∞; at n=0n=0n=0, every feasible datum contributes 000.

QAlgorithms.Papers.Independent.Intakethird100.P42.bookProbability

For a finite type YYY with a specified finite enumeration, an arbitrary real function q:Y→Rq:Y\to\mathbb Rq:Y→R, natural numbers n,Mn,Mn,M, and any predicate E\mathcal EE on functions b:[M]→Y[n]b:[M]\to Y^{[n]}b:[M]→Y[n], this defines the real number

∑b:[M]→Y[n]{∏i∈[M]∏j∈[n]q(b(i)(j)),E(b),0,¬E(b).\sum_{b:[M]\to Y^{[n]}} \begin{cases} \displaystyle\prod_{i\in[M]}\prod_{j\in[n]}q\bigl(b(i)(j)\bigr),&\mathcal E(b),\\ 0,&\neg\mathcal E(b). \end{cases}b:[M]→Y[n]∑​⎩⎨⎧​i∈[M]∏​j∈[n]∏​q(b(i)(j)),0,​E(b),¬E(b).​

There is no assumption that qqq is nonnegative or sums to 111. Both n=0n=0n=0 and M=0M=0M=0 are included: in either case there is exactly one book and its product weight is 111, so the result is 111 or 000 according to the predicate on that book. If YYY is empty and both n,Mn,Mn,M are positive, there are no books and the sum is 000.

QAlgorithms.Papers.Independent.Intakethird100.P42.coveringError

For finite types X,YX,YX,Y, arbitrary real functions q:Y→Rq:Y\to\mathbb Rq:Y→R and P:Y×X→RP:Y\times X\to\mathbb RP:Y×X→R, a natural number MMM, an indexed book b:{0,…,M−1}→Yb:\{0,\ldots,M-1\}\to Yb:{0,…,M−1}→Y, and a function e(x,i)≥0e(x,i)\ge0e(x,i)≥0 with ∑i<Me(x,i)=1\sum_{i<M}e(x,i)=1∑i<M​e(x,i)=1 for every x∈Xx\in Xx∈X, the covering error is

∑i=0M−1∑x∈X∣P(bi,x)M−(∑y∈Yq(y)P(y,x))e(x,i)∣.\sum_{i=0}^{M-1}\sum_{x\in X} \left|\frac{P(b_i,x)}{M} -\left(\sum_{y\in Y}q(y)P(y,x)\right)e(x,i)\right|.i=0∑M−1​x∈X∑​​MP(bi​,x)​−​y∈Y∑​q(y)P(y,x)​e(x,i)​.

No probability assumptions on qqq or PPP are imposed. Division by 000 has value 000, and the outer sum is empty when M=0M=0M=0, giving error 000 whenever the supplied encoder exists; such an encoder cannot exist for M=0M=0M=0 and nonempty XXX.

QAlgorithms.Papers.Independent.Intakethird100.P42.coveringMarginal

For arbitrary types X,YX,YX,Y with YYY finite, functions q:Y→Rq:Y\to\mathbb Rq:Y→R and P:Y×X→RP:Y\times X\to\mathbb RP:Y×X→R, and x∈Xx\in Xx∈X, the covering marginal at xxx is ∑y∈Yq(y)P(y,x)\sum_{y\in Y}q(y)P(y,x)∑y∈Y​q(y)P(y,x). Neither function is required to be nonnegative or normalized, and XXX need not be finite. If YYY is empty, the value is 000.

QAlgorithms.Papers.Independent.Intakethird100.P42.diagonalCQSource

For any type XXX, any function p:X→Rp:X\to\mathbb Rp:X→R, and auxiliary data ZZZ consisting of a positive natural number zzz and numbers Z(u∣x)≥0Z(u\mid x)\ge0Z(u∣x)≥0, indexed by x∈Xx\in Xx∈X and 0≤u<z0\le u<z0≤u<z, with ∑u<zZ(u∣x)=1\sum_{u<z}Z(u\mid x)=1∑u<z​Z(u∣x)=1 for every xxx, this definition constructs data with label type XXX, matrix index type {0,…,z−1}\{0,\ldots,z-1\}{0,…,z−1}, weight function ppp, and matrices

ρx(u,v)={Z(u∣x),u=v,0,u≠v,\rho_x(u,v)= \begin{cases} Z(u\mid x),&u=v,\\ 0,&u\ne v, \end{cases}ρx​(u,v)={Z(u∣x),0,​u=v,u=v,​

viewing the real entries as complex numbers. There is no finiteness or nonemptiness assumption on XXX, and no nonnegativity or normalization assumption on ppp.

QAlgorithms.Papers.Independent.Intakethird100.P42.dilutionError

For any finite type YYY, any function q:Y→Rq:Y\to\mathbb Rq:Y→R, natural numbers n,Mn,Mn,M, and list of indexed words b:{0,…,M−1}→Ynb:\{0,\ldots,M-1\}\to Y^nb:{0,…,M−1}→Yn, where YnY^nYn denotes the functions {0,…,n−1}→Y\{0,\ldots,n-1\}\to Y{0,…,n−1}→Y, the dilution error is

∑y∈Yn∣∏j=0n−1q(yj)−#{i<M:bi=y}M∣.\sum_{y\in Y^n}\left|\prod_{j=0}^{n-1}q(y_j)-\frac{\#\{i<M:b_i=y\}}{M}\right|.y∈Yn∑​​j=0∏n−1​q(yj​)−M#{i<M:bi​=y}​​.

No nonnegativity or normalization of qqq is assumed. Empty products are 111, empty sums are 000, and division by 000 is defined as 000; thus the empirical term is 000 when M=0M=0M=0. When n=0n=0n=0, there is exactly one empty word, and the error is 000 for M>0M>0M>0 and 111 for M=0M=0M=0.

QAlgorithms.Papers.Independent.Intakethird100.P42.empiricalLaw

For a finite type YYY with a specified finite enumeration, any natural number MMM, any function b:[M]→Yb:[M]\to Yb:[M]→Y, and any y∈Yy\in Yy∈Y, this defines

∣{i∈[M]:b(i)=y}∣M,\frac{|\{i\in[M]:b(i)=y\}|}{M},M∣{i∈[M]:b(i)=y}∣​,

with numerator and denominator regarded as real numbers. The definition includes M=0M=0M=0: the numerator is then zero, and real division by zero is defined to be zero, so the value is 000 for every yyy. No positive-size assumption is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P42.exactDistortionExpected

For finite, finitely enumerated types X,BX,BX,B, arbitrary data E=(p,ρ)E=(p,\rho)E=(p,ρ) consisting of a real function p:X→Rp:X\to\mathbb Rp:X→R and complex matrices ρx\rho_xρx​ on CB\mathbb C^BCB, an arbitrary function d:X×X→Rd:X\times X\to\mathbb Rd:X×X→R, a natural number nnn, and an exact distortion code PPP, the defined real number is

∑x,z∈Xn(∏i=0n−1p(xi))Re⁡tr⁡ ⁣(Me(x),zρx(n))∑i=0n−1d(xi,zi)n,(ρx(n))u,v=∏i=0n−1(ρxi)ui,vi.\sum_{x,z\in X^n} \left(\prod_{i=0}^{n-1}p(x_i)\right) \operatorname{Re}\operatorname{tr}\!\left(M_{e(x),z}\rho_x^{(n)}\right) \frac{\sum_{i=0}^{n-1}d(x_i,z_i)}{n}, \qquad (\rho_x^{(n)})_{u,v} =\prod_{i=0}^{n-1}(\rho_{x_i})_{u_i,v_i}.x,z∈Xn∑​(i=0∏n−1​p(xi​))Retr(Me(x),z​ρx(n)​)n∑i=0n−1​d(xi​,zi​)​,(ρx(n)​)u,v​=i=0∏n−1​(ρxi​​)ui​,vi​​.

Here Xn,BnX^n,B^nXn,Bn are length-nnn word types, PPP has some number b≥0b\ge0b≥0 of message bits, e:Xn→{0,1}be:X^n\to\{0,1\}^be:Xn→{0,1}b, and, for every message mmm, Mm,z⪰0M_{m,z}\succeq0Mm,z​⪰0 and ∑zMm,z=I\sum_zM_{m,z}=I∑z​Mm,z​=I on CBn\mathbb C^{B^n}CBn. No normalization or positivity of ppp or ρx\rho_xρx​, and no sign condition on ddd, is assumed. Empty alphabets and n=0n=0n=0 are allowed; at n=0n=0n=0, empty products equal 111, and the distortion factor is 0/0=00/0=00/0=0, so the defined value is 000.

QAlgorithms.Papers.Independent.Intakethird100.P42.exactDistortionRate

For finite, finitely enumerated types X,BX,BX,B, arbitrary data E=(p,ρ)E=(p,\rho)E=(p,ρ) with real weights p(x)p(x)p(x) and complex matrices ρx\rho_xρx​ on CB\mathbb C^BCB, any d:X×X→Rd:X\times X\to\mathbb Rd:X×X→R, and any real DDD, the defined value is the infimum in [0,∞][0,\infty][0,∞] of all nonnegative real rates RRR having the following property: for every δ>0\delta>0δ>0, all sufficiently large positive integers nnn admit an encoder e:Xn→{0,1}be:X^n\to\{0,1\}^be:Xn→{0,1}b, for some integer b≥0b\ge0b≥0, and matrices Mm,zM_{m,z}Mm,z​ on CBn\mathbb C^{B^n}CBn with Mm,z⪰0M_{m,z}\succeq0Mm,z​⪰0 and ∑z∈XnMm,z=I\sum_{z\in X^n}M_{m,z}=I∑z∈Xn​Mm,z​=I for every message mmm, such that

b≤n(R+δ),∑x,z∈Xn(∏ip(xi))Re⁡tr⁡ ⁣(Me(x),zρx(n))∑id(xi,zi)n≤D,(ρx(n))u,v=∏i(ρxi)ui,vi.b\le n(R+\delta),\qquad \sum_{x,z\in X^n} \left(\prod_i p(x_i)\right) \operatorname{Re}\operatorname{tr}\!\left(M_{e(x),z}\rho_x^{(n)}\right) \frac{\sum_i d(x_i,z_i)}{n}\le D, \quad (\rho_x^{(n)})_{u,v}=\prod_i(\rho_{x_i})_{u_i,v_i}.b≤n(R+δ),x,z∈Xn∑​(i∏​p(xi​))Retr(Me(x),z​ρx(n)​)n∑i​d(xi​,zi​)​≤D,(ρx(n)​)u,v​=i∏​(ρxi​​)ui​,vi​​.

More precisely, each such RRR enters the infimum through its nonnegative extended-real image, which equals RRR because R≥0R\ge0R≥0. If no such rate exists, the value is +∞+\infty+∞; attainment is not part of this definition. Empty alphabets are permitted, and no distribution, density-matrix, or nonnegativity assumptions on E,d,DE,d,DE,d,D are imposed.

QAlgorithms.Papers.Independent.Intakethird100.P42.feedbackMeasurementBranch

For a finite type BBB with a specified finite enumeration and arbitrary complex matrices F,ρF,\rhoF,ρ indexed by BBB, choose the enumeration equivalence e:B→[∣B∣]e:B\to[|B|]e:B→[∣B∣], and write F~ij=Fe−1(i),e−1(j)\widetilde F_{ij}=F_{e^{-1}(i),e^{-1}(j)}Fij​=Fe−1(i),e−1(j)​ and ρ~ij=ρe−1(i),e−1(j)\widetilde\rho_{ij}=\rho_{e^{-1}(i),e^{-1}(j)}ρ​ij​=ρe−1(i),e−1(j)​. The resulting matrix has entries

(cfc⁡(⋅,F~) ρ~ cfc⁡(⋅,F~))e(a),e(b)(a,b∈B),\bigl(\operatorname{cfc}(\sqrt{\cdot},\widetilde F)\, \widetilde\rho\, \operatorname{cfc}(\sqrt{\cdot},\widetilde F)\bigr)_{e(a),e(b)} \qquad(a,b\in B),(cfc(⋅​,F)ρ​cfc(⋅​,F))e(a),e(b)​(a,b∈B),

where cfc⁡\operatorname{cfc}cfc is the continuous functional calculus and the square root is the real square-root function, which is zero on negative inputs. No positivity, Hermiticity, trace, or normalization hypothesis is imposed on either matrix; the functional-calculus expression uses its total definition even outside its usual hypotheses. The empty finite type BBB is also allowed.

QAlgorithms.Papers.Independent.Intakethird100.P42.feedbackTraceError

For finite types X,BX,BX,B with specified finite enumerations, arbitrary ensemble data px∈R,ρx∈CB×Bp_x\in\mathbb R,\rho_x\in\mathbb C^{B\times B}px​∈R,ρx​∈CB×B, an auxiliary system with positive size yyy, probabilities V(u∣x)≥0V(u\mid x)\ge0V(u∣x)≥0 summing to 111 over u∈[y]u\in[y]u∈[y], a natural number nnn, and a code, this is

Re⁡tr⁡cfc⁡ ⁣(⋅,(T−O)†(T−O)).\operatorname{Re}\operatorname{tr} \operatorname{cfc}\!\left(\sqrt{\cdot},(T-O)^\dagger(T-O)\right).Retrcfc(⋅​,(T−O)†(T−O)).

Here both matrices are indexed by (x,u,v,b)∈X[n]×[y][n]×[y][n]×B[n](x,u,v,b)\in X^{[n]}\times[y]^{[n]}\times[y]^{[n]}\times B^{[n]}(x,u,v,b)∈X[n]×[y][n]×[y][n]×B[n]. Write pn(x)=∏ipxip_n(x)=\prod_i p_{x_i}pn​(x)=∏i​pxi​​ and (ρx(n))b,b′=∏i(ρxi)bi,bi′(\rho_x^{(n)})_{b,b'}=\prod_i(\rho_{x_i})_{b_i,b'_i}(ρx(n)​)b,b′​=∏i​(ρxi​​)bi​,bi′​​. The target entry is zero unless x=x′,u=u′,v=v′,u=vx=x',u=u',v=v',u=vx=x′,u=u′,v=v′,u=v, and otherwise is pn(x)∏iV(ui∣xi)(ρx(n))b,b′p_n(x)\prod_iV(u_i\mid x_i)(\rho_x^{(n)})_{b,b'}pn​(x)∏i​V(ui​∣xi​)(ρx(n)​)b,b′​. To define OOO, the code supplies bit counts c,h,pc,h,pc,h,p, nonnegative encoding weights e(x,ℓ,m,s)e(x,\ell,m,s)e(x,ℓ,m,s) with ∑m,se(x,ℓ,m,s)=1\sum_{m,s}e(x,\ell,m,s)=1∑m,s​e(x,ℓ,m,s)=1, decoded words d(ℓ,m,s)d(\ell,m,s)d(ℓ,m,s), and positive semidefinite matrices Fℓ,m,tF_{\ell,m,t}Fℓ,m,t​ on B[n]B^{[n]}B[n] with ∑tFℓ,m,t=I\sum_tF_{\ell,m,t}=I∑t​Fℓ,m,t​=I. Its entry is zero unless x=x′,u=u′,v=v′x=x',u=u',v=v'x=x′,u=u′,v=v′, and otherwise equals

∑ℓ,m,s,t1{d(ℓ,m,s)=u, d(ℓ,m,t)=v} pn(x)2−ce(x,ℓ,m,s)[LFℓ,m,t(ρx(n))]b,b′,\sum_{\ell,m,s,t} \mathbf1_{\{d(\ell,m,s)=u,\ d(\ell,m,t)=v\}}\, p_n(x)2^{-c}e(x,\ell,m,s) \bigl[\mathcal L_{F_{\ell,m,t}}(\rho_x^{(n)})\bigr]_{b,b'},ℓ,m,s,t∑​1{d(ℓ,m,s)=u, d(ℓ,m,t)=v}​pn​(x)2−ce(x,ℓ,m,s)[LFℓ,m,t​​(ρx(n)​)]b,b′​,

where ℓ,m,s,t\ell,m,s,tℓ,m,s,t are bitstrings of lengths c,h,p,pc,h,p,pc,h,p,p, and LF(σ)\mathcal L_F(\sigma)LF​(σ) is obtained by relabeling onto a finite ordinal, applying the functional-calculus expression cfc⁡(⋅,F)σcfc⁡(⋅,F)\operatorname{cfc}(\sqrt{\cdot},F)\sigma\operatorname{cfc}(\sqrt{\cdot},F)cfc(⋅​,F)σcfc(⋅​,F), and relabeling back. No probability or density-matrix assumptions on px,ρxp_x,\rho_xpx​,ρx​, and no nonemptiness assumptions on X,BX,BX,B, are imposed. The definition includes n=0n=0n=0, when products are empty products and word sets are singleton.

QAlgorithms.Papers.Independent.Intakethird100.P42.fqrsQuantum

For arbitrary finite types A,C,B,RA,C,B,RA,C,B,R, possibly empty, and any vector ψ∈CR×(A×(C×B))\psi\in\mathbb C^{R\times(A\times(C\times B))}ψ∈CR×(A×(C×B)), form ρuv=ψuψv‾\rho_{uv}=\psi_u\overline{\psi_v}ρuv​=ψu​ψv​​. For a register collection UUU, obtain ρU\rho_UρU​ by summing over equal complementary labels, so that (ρU)u,u′=∑wψ(u,w)ψ(u′,w)‾(\rho_U)_{u,u'}=\sum_w\psi_{(u,w)}\overline{\psi_{(u',w)}}(ρU​)u,u′​=∑w​ψ(u,w)​ψ(u′,w)​​, with coordinates in their original order. Define S(U)=−Re⁡tr⁡(f(ρU))S(U)=-\operatorname{Re}\operatorname{tr}(f(\rho_U))S(U)=−Retr(f(ρU​)), where f(t)=tlog⁡2tf(t)=t\log_2tf(t)=tlog2​t is applied by continuous functional calculus and 0log⁡20=00\log_20=00log2​0=0. The assigned real number is exactly

S(C)+S(RB)−S(RCB)2.\frac{S(C)+S(RB)-S(RCB)}{2}.2S(C)+S(RB)−S(RCB)​.

No normalization or nonzero condition is imposed on ψ\psiψ, and the definition also applies when a register type is empty.

QAlgorithms.Papers.Independent.Intakethird100.P42.printedQuantumCost

For arbitrary finite types A,C,B,RA,C,B,RA,C,B,R, possibly empty, and any vector ψ∈CR×(A×(C×B))\psi\in\mathbb C^{R\times(A\times(C\times B))}ψ∈CR×(A×(C×B)), let ρuv=ψuψv‾\rho_{uv}=\psi_u\overline{\psi_v}ρuv​=ψu​ψv​​. For any collection UUU of registers, let (ρU)u,u′=∑wψ(u,w)ψ(u′,w)‾(\rho_U)_{u,u'}=\sum_w\psi_{(u,w)}\overline{\psi_{(u',w)}}(ρU​)u,u′​=∑w​ψ(u,w)​ψ(u′,w)​​, summing over the complementary labels with coordinates restored to their original order, and put S(U)=−Re⁡tr⁡(f(ρU))S(U)=-\operatorname{Re}\operatorname{tr}(f(\rho_U))S(U)=−Retr(f(ρU​)), where f(t)=tlog⁡2tf(t)=t\log_2tf(t)=tlog2​t is applied by continuous functional calculus, with 0log⁡20=00\log_20=00log2​0=0. The definition assigns

S(AC)+S(RA)−S(A)−S(RAC)2.\frac{S(AC)+S(RA)-S(A)-S(RAC)}{2}.2S(AC)+S(RA)−S(A)−S(RAC)​.

This formula requires neither ∥ψ∥=1\|\psi\|=1∥ψ∥=1 nor ψ≠0\psi\ne0ψ=0, and includes empty register types.

QAlgorithms.Papers.Independent.Intakethird100.P42.productConditional

For arbitrary types X,YX,YX,Y, any real function P:Y×X→RP:Y\times X\to\mathbb RP:Y×X→R, any natural number nnn, and words y:{0,…,n−1}→Yy:\{0,\ldots,n-1\}\to Yy:{0,…,n−1}→Y and x:{0,…,n−1}→Xx:\{0,\ldots,n-1\}\to Xx:{0,…,n−1}→X, the product conditional is ∏i=0n−1P(yi,xi)\prod_{i=0}^{n-1}P(y_i,x_i)∏i=0n−1​P(yi​,xi​). No finiteness, nonemptiness, nonnegativity, or normalization assumptions are imposed on X,Y,PX,Y,PX,Y,P. For n=0n=0n=0, the product is 111.

QAlgorithms.Papers.Independent.Intakethird100.P42.quantumRDCost

For finite, finitely enumerated types X,BX,BX,B, arbitrary data E=(p,ρ)E=(p,\rho)E=(p,ρ) with p:X→Rp:X\to\mathbb Rp:X→R and complex matrices ρx\rho_xρx​ on CB\mathbb C^BCB, and a quantum rate-distortion choice PPP, let a>0a>0a>0 be its output dimension and let τx\tau_xτx​ be its positive semidefinite, trace-one output matrices on Ca\mathbb C^aCa. Write

τˉ=∑xp(x)τx,ρˉ=∑xp(x)ρx,Ω=∑xp(x)(τx⊗ρx),S(A)=−Re⁡tr⁡ ⁣(f(A)),f(t)=tlog⁡2t,\bar\tau=\sum_xp(x)\tau_x,\qquad \bar\rho=\sum_xp(x)\rho_x,\qquad \Omega=\sum_xp(x)(\tau_x\otimes\rho_x), \qquad S(A)=-\operatorname{Re}\operatorname{tr}\!\bigl(f(A)\bigr), \quad f(t)=t\log_2t,τˉ=x∑​p(x)τx​,ρˉ​=x∑​p(x)ρx​,Ω=x∑​p(x)(τx​⊗ρx​),S(A)=−Retr(f(A)),f(t)=tlog2​t,

where f(A)f(A)f(A) is evaluated by continuous functional calculus. The defined cost is

[S(τˉ)−∑xp(x)S(τx)]−[S(τˉ)+S(ρˉ)−S(Ω)].\left[S(\bar\tau)-\sum_xp(x)S(\tau_x)\right] -\left[S(\bar\tau)+S(\bar\rho)-S(\Omega)\right].[S(τˉ)−x∑​p(x)S(τx​)]−[S(τˉ)+S(ρˉ​)−S(Ω)].

The choice PPP also includes matrices Mz⪰0M_z\succeq0Mz​⪰0 on Ca×B\mathbb C^{a\times B}Ca×B with ∑zMz=I\sum_zM_z=I∑z​Mz​=I, although these matrices do not occur in the cost. No validity assumptions on EEE or nonemptiness assumptions on X,BX,BX,B are made, and no sign condition on the resulting cost is included.

QAlgorithms.Papers.Independent.Intakethird100.P42.quantumRDExpected

For finite, finitely enumerated types X,BX,BX,B, arbitrary data E=(p,ρ)E=(p,\rho)E=(p,ρ) consisting of p:X→Rp:X\to\mathbb Rp:X→R and complex matrices ρx\rho_xρx​ on CB\mathbb C^BCB, any function d:X×X→Rd:X\times X\to\mathbb Rd:X×X→R, and a quantum rate-distortion choice PPP, the defined real number is

∑x,z∈Xp(x) Re⁡tr⁡ ⁣(Mz(τx⊗ρx))d(x,z).\sum_{x,z\in X}p(x)\, \operatorname{Re}\operatorname{tr}\!\left(M_z(\tau_x\otimes\rho_x)\right)d(x,z).x,z∈X∑​p(x)Retr(Mz​(τx​⊗ρx​))d(x,z).

Here PPP supplies an integer a>0a>0a>0, positive semidefinite matrices τx\tau_xτx​ on Ca\mathbb C^aCa with trace 111, and positive semidefinite matrices MzM_zMz​ on Ca×B\mathbb C^{a\times B}Ca×B satisfying ∑zMz=I\sum_zM_z=I∑z​Mz​=I. The tensor product has entries (τx⊗ρx)(i,u),(j,v)=(τx)i,j(ρx)u,v(\tau_x\otimes\rho_x)_{(i,u),(j,v)}=(\tau_x)_{i,j}(\rho_x)_{u,v}(τx​⊗ρx​)(i,u),(j,v)​=(τx​)i,j​(ρx​)u,v​. No probability or density-matrix conditions are imposed on EEE, and ddd may take negative values. Empty XXX or BBB is permitted whenever such a choice PPP exists; if XXX is empty, the displayed sum is 000.

QAlgorithms.Papers.Independent.Intakethird100.P42.quantumRDValue

For finite, finitely enumerated types X,BX,BX,B, arbitrary data E=(p,ρ)E=(p,\rho)E=(p,ρ) with real weights p(x)p(x)p(x) and complex matrices ρx\rho_xρx​ on CB\mathbb C^BCB, any d:X×X→Rd:X\times X\to\mathbb Rd:X×X→R, and any real DDD, consider all choices of an integer a>0a>0a>0, matrices τx⪰0\tau_x\succeq0τx​⪰0 on Ca\mathbb C^aCa with tr⁡τx=1\operatorname{tr}\tau_x=1trτx​=1, and matrices Mz⪰0M_z\succeq0Mz​⪰0 on Ca×B\mathbb C^{a\times B}Ca×B with ∑zMz=I\sum_zM_z=I∑z​Mz​=I, for which

∑x,zp(x)Re⁡tr⁡ ⁣(Mz(τx⊗ρx))d(x,z)≤D.\sum_{x,z}p(x)\operatorname{Re}\operatorname{tr}\!\left(M_z(\tau_x\otimes\rho_x)\right)d(x,z)\le D.x,z∑​p(x)Retr(Mz​(τx​⊗ρx​))d(x,z)≤D.

For each choice, set τˉ=∑xp(x)τx\bar\tau=\sum_xp(x)\tau_xτˉ=∑x​p(x)τx​, ρˉ=∑xp(x)ρx\bar\rho=\sum_xp(x)\rho_xρˉ​=∑x​p(x)ρx​, Ω=∑xp(x)(τx⊗ρx)\Omega=\sum_xp(x)(\tau_x\otimes\rho_x)Ω=∑x​p(x)(τx​⊗ρx​), and

C=[S(τˉ)−∑xp(x)S(τx)]−[S(τˉ)+S(ρˉ)−S(Ω)],S(A)=−Re⁡tr⁡ ⁣((t↦tlog⁡2t)(A)),C=\left[S(\bar\tau)-\sum_xp(x)S(\tau_x)\right] -\left[S(\bar\tau)+S(\bar\rho)-S(\Omega)\right], \qquad S(A)=-\operatorname{Re}\operatorname{tr}\!\bigl((t\mapsto t\log_2t)(A)\bigr),C=[S(τˉ)−x∑​p(x)S(τx​)]−[S(τˉ)+S(ρˉ​)−S(Ω)],S(A)=−Retr((t↦tlog2​t)(A)),

using continuous functional calculus. The defined value is the infimum in [0,∞][0,\infty][0,∞] of the numbers max⁡(C,0)\max(C,0)max(C,0) from these choices. Thus negative costs are mapped to 000, and an empty set of admissible choices gives +∞+\infty+∞. Neither attainment, validity of EEE, nonnegativity of d,Dd,Dd,D, nor nonemptiness of X,BX,BX,B is assumed.

QAlgorithms.Papers.Independent.Intakethird100.P42.rdExpected

For finite types X,BX,BX,B, data EEE consisting of an arbitrary real function p:X→Rp:X\to\mathbb Rp:X→R and arbitrary complex matrices ρa\rho_aρa​ on BBB, an arbitrary function d:X×X→Rd:X\times X\to\mathbb Rd:X×X→R, a natural number nnn, and a rate-distortion datum consisting of L>0L>0L>0, nonnegative normalized conditional weights V(y∣x)V(y\mid x)V(y∣x) for y<Ly<Ly<L, and positive semidefinite matrices Dy,zD_{y,z}Dy,z​ on BnB^nBn satisfying ∑z∈XnDy,z=I\sum_{z\in X^n}D_{y,z}=I∑z∈Xn​Dy,z​=I, the expected-distortion expression is

∑x∈Xn∑y<L∑z∈Xn(∏i<np(xi))V(y∣x) Re⁡tr⁡(Dy,zρx(n)) ∑i<nd(xi,zi)n,(ρx(n))a,b=∏i<n(ρxi)ai,bi.\sum_{x\in X^n}\sum_{y<L}\sum_{z\in X^n} \left(\prod_{i<n}p(x_i)\right)V(y\mid x)\, \operatorname{Re}\operatorname{tr}(D_{y,z}\rho_x^{(n)})\, \frac{\sum_{i<n}d(x_i,z_i)}{n}, \qquad (\rho_x^{(n)})_{a,b}=\prod_{i<n}(\rho_{x_i})_{a_i,b_i}.x∈Xn∑​y<L∑​z∈Xn∑​(i<n∏​p(xi​))V(y∣x)Retr(Dy,z​ρx(n)​)n∑i<n​d(xi​,zi​)​,(ρx(n)​)a,b​=i<n∏​(ρxi​​)ai​,bi​​.

No validity assumption on EEE is imposed: ppp need not be a probability function and the matrices ρa\rho_aρa​ need not be density matrices. The function ddd may take arbitrary real values. Empty products are 111, empty sums are 000, and division by 000 is 000, so this expression is 000 at n=0n=0n=0.

QAlgorithms.Papers.Independent.Intakethird100.P42.redistributionReceiverInformation

For arbitrary finite types A,C,B,RA,C,B,RA,C,B,R, possibly empty, and any vector ψ∈CR×(A×(C×B))\psi\in\mathbb C^{R\times(A\times(C\times B))}ψ∈CR×(A×(C×B)), form the matrix ρuv=ψuψv‾\rho_{uv}=\psi_u\overline{\psi_v}ρuv​=ψu​ψv​​. For a collection UUU of these registers, let ρU\rho_UρU​ be the matrix obtained by summing ρ\rhoρ's entries over equal labels in all complementary registers; equivalently, (ρU)u,u′=∑wψ(u,w)ψ(u′,w)‾(\rho_U)_{u,u'}=\sum_w\psi_{(u,w)}\overline{\psi_{(u',w)}}(ρU​)u,u′​=∑w​ψ(u,w)​ψ(u′,w)​​, with the coordinates placed in their original order. Write S(U)=−Re⁡tr⁡(f(ρU))S(U)=-\operatorname{Re}\operatorname{tr}(f(\rho_U))S(U)=−Retr(f(ρU​)), where f(t)=tlog⁡2tf(t)=t\log_2tf(t)=tlog2​t is applied by continuous functional calculus, with 0log⁡20=00\log_20=00log2​0=0. The definition assigns the real number

S(C)+S(B)−S(CB).S(C)+S(B)-S(CB).S(C)+S(B)−S(CB).

There is no assumption that ∥ψ∥=1\|\psi\|=1∥ψ∥=1; the same formula applies to the zero vector and to empty register types.

QAlgorithms.Papers.Independent.Intakethird100.P42.senderMutual

For arbitrary finite types A,C,B,RA,C,B,RA,C,B,R, possibly empty, and any vector ψ∈CR×(A×(C×B))\psi\in\mathbb C^{R\times(A\times(C\times B))}ψ∈CR×(A×(C×B)), set ρuv=ψuψv‾\rho_{uv}=\psi_u\overline{\psi_v}ρuv​=ψu​ψv​​. For each collection UUU of registers, define its reduced matrix by (ρU)u,u′=∑wψ(u,w)ψ(u′,w)‾(\rho_U)_{u,u'}=\sum_w\psi_{(u,w)}\overline{\psi_{(u',w)}}(ρU​)u,u′​=∑w​ψ(u,w)​ψ(u′,w)​​, summing over all complementary register labels and restoring the original coordinate order. Let S(U)=−Re⁡tr⁡(f(ρU))S(U)=-\operatorname{Re}\operatorname{tr}(f(\rho_U))S(U)=−Retr(f(ρU​)), with f(t)=tlog⁡2tf(t)=t\log_2tf(t)=tlog2​t applied by continuous functional calculus and 0log⁡20=00\log_20=00log2​0=0. The definition assigns

S(C)+S(A)−S(AC).S(C)+S(A)-S(AC).S(C)+S(A)−S(AC).

The vector need not be normalized, and zero vectors and empty register types are included.

QAlgorithms.Papers.Independent.Intakethird100.P42.simulationOutput

For an arbitrary type XXX, a finite type BBB with a specified finite enumeration, arbitrary real weights pxp_xpx​ and complex matrices ρx\rho_xρx​ on BBB, natural numbers y,ny,ny,n, and a code, let its bit counts be c,h,pc,h,pc,h,p, its encoding numbers be e(x,ℓ,m,s)≥0e(x,\ell,m,s)\ge0e(x,ℓ,m,s)≥0 with ∑m,se(x,ℓ,m,s)=1\sum_{m,s}e(x,\ell,m,s)=1∑m,s​e(x,ℓ,m,s)=1, its decoder be d(ℓ,m,s)∈[y][n]d(\ell,m,s)\in[y]^{[n]}d(ℓ,m,s)∈[y][n], and its matrices Fℓ,m,tF_{\ell,m,t}Fℓ,m,t​ on B[n]B^{[n]}B[n] satisfy Fℓ,m,t⪰0F_{\ell,m,t}\succeq0Fℓ,m,t​⪰0 and ∑tFℓ,m,t=I\sum_tF_{\ell,m,t}=I∑t​Fℓ,m,t​=I for every ℓ,m\ell,mℓ,m; the indices ℓ,m,s,t\ell,m,s,tℓ,m,s,t range respectively over bitstrings of lengths c,h,p,pc,h,p,pc,h,p,p. Put pn(x)=∏ipxip_n(x)=\prod_i p_{x_i}pn​(x)=∏i​pxi​​ and (ρx(n))b,b′=∏i(ρxi)bi,bi′(\rho_x^{(n)})_{b,b'}=\prod_i(\rho_{x_i})_{b_i,b'_i}(ρx(n)​)b,b′​=∏i​(ρxi​​)bi​,bi′​​. For matrices F,σF,\sigmaF,σ on B[n]B^{[n]}B[n], let LF(σ)\mathcal L_F(\sigma)LF​(σ) mean: relabel both matrices along the specified equivalence B[n]≃[∣B[n]∣]B^{[n]}\simeq[|B^{[n]}|]B[n]≃[∣B[n]∣], form cfc⁡(⋅,F)σcfc⁡(⋅,F)\operatorname{cfc}(\sqrt{\cdot},F)\sigma\operatorname{cfc}(\sqrt{\cdot},F)cfc(⋅​,F)σcfc(⋅​,F) in those coordinates, and relabel back. The defined output matrix is

O(x,u,v,b),(x′,u′,v′,b′)={∑ℓ,m,s,t1{d(ℓ,m,s)=u, d(ℓ,m,t)=v} pn(x) 2−c e(x,ℓ,m,s) [LFℓ,m,t(ρx(n))]b,b′,x=x′, u=u′, v=v′,0,otherwise.O_{(x,u,v,b),(x',u',v',b')}= \begin{cases} \displaystyle \sum_{\ell,m,s,t} \mathbf1_{\{d(\ell,m,s)=u,\ d(\ell,m,t)=v\}}\, p_n(x)\,2^{-c}\,e(x,\ell,m,s)\, \bigl[\mathcal L_{F_{\ell,m,t}}(\rho_x^{(n)})\bigr]_{b,b'}, &x=x',\ u=u',\ v=v',\\ 0,&\text{otherwise}. \end{cases}O(x,u,v,b),(x′,u′,v′,b′)​=⎩⎨⎧​ℓ,m,s,t∑​1{d(ℓ,m,s)=u, d(ℓ,m,t)=v}​pn​(x)2−ce(x,ℓ,m,s)[LFℓ,m,t​​(ρx(n)​)]b,b′​,0,​x=x′, u=u′, v=v′,otherwise.​

Thus sss supplies the encoding weight and the first decoded word, while the independently summed ttt supplies the measurement matrix and the second decoded word; equality of those two words is not required. No validity hypothesis is imposed on the ensemble data. Zero lengths and zero bit counts are permitted, with empty products equal to 111 and zero-length bitstring sets singleton.

QAlgorithms.Papers.Independent.Intakethird100.P42.simulationTarget

For arbitrary types X,BX,BX,B, data consisting of a real function p:X→Rp:X\to\mathbb Rp:X→R and complex matrices ρx\rho_xρx​ on BBB, an auxiliary system consisting of a positive natural number yyy and numbers V(u∣x)≥0V(u\mid x)\ge0V(u∣x)≥0 with ∑u∈[y]V(u∣x)=1\sum_{u\in[y]}V(u\mid x)=1∑u∈[y]​V(u∣x)=1 for every xxx, and any natural number nnn, this defines a matrix indexed by tuples (x,u,v,b)∈X[n]×([y][n]×([y][n]×B[n]))(x,u,v,b)\in X^{[n]}\times([y]^{[n]}\times([y]^{[n]}\times B^{[n]}))(x,u,v,b)∈X[n]×([y][n]×([y][n]×B[n])). Its entries are

T(x,u,v,b),(x′,u′,v′,b′)={(∏i∈[n]p(xi)V(ui∣xi))(∏i∈[n](ρxi)bi,bi′),x=x′, u=u′, v=v′, u=v,0,otherwise.T_{(x,u,v,b),(x',u',v',b')}= \begin{cases} \displaystyle \left(\prod_{i\in[n]}p(x_i)V(u_i\mid x_i)\right) \left(\prod_{i\in[n]}(\rho_{x_i})_{b_i,b'_i}\right), &x=x',\ u=u',\ v=v',\ u=v,\\ 0,&\text{otherwise}. \end{cases}T(x,u,v,b),(x′,u′,v′,b′)​=⎩⎨⎧​​i∈[n]∏​p(xi​)V(ui​∣xi​)​​i∈[n]∏​(ρxi​​)bi​,bi′​​​,0,​x=x′, u=u′, v=v′, u=v,otherwise.​

Real coefficients in this expression are regarded as complex numbers. The functions p,ρxp,\rho_xp,ρx​ are not assumed to form valid ensemble data: neither probability normalization nor density-matrix conditions are required. Neither XXX nor BBB is assumed finite or nonempty. Empty products equal 111, so at n=0n=0n=0 this is the one-by-one matrix with entry 111.

QAlgorithms.Papers.Independent.Intakethird100.P42.wordDistortion

For any type XXX, any function d:X×X→Rd:X\times X\to\mathbb Rd:X×X→R, any natural number nnn, and words x,z:{0,…,n−1}→Xx,z:\{0,\ldots,n-1\}\to Xx,z:{0,…,n−1}→X, the word distortion is 1n∑i=0n−1d(xi,zi)\frac1n\sum_{i=0}^{n-1}d(x_i,z_i)n1​∑i=0n−1​d(xi​,zi​). No sign, symmetry, diagonal, or metric conditions are imposed on ddd, and XXX need not be finite. When n=0n=0n=0, both the empty sum and division by 000 yield 000.

Human review
  • Endorsed by Shuze Chen · Nov 1, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Gohan · Nov 1, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me