---
title: Linear Matroid Parity Problem
url: https://www.emergentmind.com/topics/linear-matroid-parity-problem
type: topic
---

# Linear Matroid Parity Problem

Searching arXiv for the cited papers and closely related work on linear matroid parity.
The Linear Matroid Parity Problem is the linearly represented form of matroid parity, also called linear matroid matching in the line-based formulation. It asks for a feasible collection of paired columns, or equivalently 2-dimensional subspaces, such that each chosen pair contributes full dimension and the total cardinality is maximized; in weighted forms, the objective is to optimize line weights; in counting forms, the objective is to enumerate feasible or minimum-weight feasible solutions. In the linear setting, the problem is a common generalization of matching and linear matroid intersection, admits algebraic formulations through skew-symmetric matrices and Pfaffians, and has deterministic polynomial-time algorithms for optimization. Exact counting, however, is only known to be tractable on special cancellation-free subclasses such as Pfaffian parities [1912.00620].

## 1. Problem statement and equivalent formulations

A standard matrix formulation uses a matrix \(A \in \mathbb{K}^{2r \times 2n}\) whose \(2n\) columns are partitioned into \(n\) disjoint pairs called lines. Writing \(L\) for the set of lines, the pair \((A,L)\) is a linear matroid parity instance. For a subset \(J \subseteq L\), the submatrix \(A[J]\) consists of all columns in the lines of \(J\). A subset \(B \subseteq L\) is a parity base if \(|B|=r\) and \(A[B]\) is nonsingular; equivalently, \(A[B]\) is a \(2r \times 2r\) invertible matrix made of whole lines. The set of parity bases is denoted \(B(A,L)\). This yields three canonical tasks: the decision/search problem of determining whether \(B(A,L)\neq\varnothing\) and finding \(B\in B(A,L)\), the optimization problem of minimizing \(w(B)=\sum_{l\in B}w(l)\) over parity bases, and the counting problem of computing \(|B(A,L)|\) or counting minimum-weight parity bases [1912.00620].

An equivalent line-based formulation represents the input as a family of 2-dimensional subspaces
\[
L=\{l_1,\dots,l_m\},\qquad l_i=\langle a_i,b_i\rangle \le \mathbb F^n.
\]
A subset \(M\) of lines is a matroid matching if
\[
\dim\Big(\sum_{l\in M} l\Big)=2|M|.
\]
A perfect matroid matching has size \(n/2\). In the pair-partition specialization, each pair \((a_i,b_i)\) spans a line \(l_i=\operatorname{span}(a_i,b_i)\), and selecting a set of pairs is feasible exactly when all selected pair-elements are independent, equivalently when the chosen lines contribute full dimension \(2|M|\). In this sense, linear matroid matching is “also known as the linear matroid parity problem” [2402.18276].

This problem strictly generalizes ordinary matching. For a graph \(G=(V,E)\) with \(|V|=2n\), one can construct a totally unimodular matrix \(A\) with two columns per edge, viewed as a line, such that perfect matchings of \(G\) are exactly parity bases \(B(A,E)\). It also generalizes linear matroid intersection: given an \(r\times n\) matrix pair \((A_1,A_2)\), Lawler’s reduction produces a \(2r\times 2n\) parity instance \((A,E)\) whose lines correspond to the original columns and such that \(B(A_1,A_2)\leftrightarrow B(A,E)\) [1912.00620].

## 2. Algebraic and Pfaffian formulations

The algebraic core of linear matroid parity is a skew-symmetric encoding. For \(z=(z_l)_{l\in L}\), define the \(2n\times 2n\) skew-symmetric block-diagonal matrix \(\Delta_L(z)\) whose \(2\times 2\) block for a line \(l\) is
\[
\begin{pmatrix}
0 & z_l\\
-z_l & 0
\end{pmatrix},
\]
and define
\[
\Phi_{A,L}(z):=
\begin{pmatrix}
O & A\\
-A^\top & \Delta_L(z)
\end{pmatrix}.
\]
A standard characterization is
\[
(A,L)\text{ has a parity base} \iff A\Delta(z)A^\top \text{ is nonsingular} \iff \Phi(z)\text{ is nonsingular}
\]
over the rational function field \(\mathbb{K}(z)\). The skew-symmetric matrix \(A\Delta(z)A^\top\) can be written as
\[
A\Delta(z)A^\top
=
A_1D(z)A_2^\top-A_2D(z)A_1^\top
=
\sum_{l=(v,\bar v)\in L} z_l\bigl(a_v a_{\bar v}^\top-a_{\bar v} a_v^\top\bigr),
\]
which already exhibits the rank-two skew-symmetric coefficient structure characteristic of the problem [1912.00620].

Using the Ishikawa–Wakayama Pfaffian Cauchy–Binet formula, the parity bases appear as monomials in a Pfaffian expansion:
\[
\pf\!\bigl(A\Delta(z)A^\top\bigr)
=
\sum_{B\in B(A,L)} \det A[B]\prod_{l\in B} z_l,
\]
and
\[
\pf \Phi(z)
=
\sum_{B\in B(A,L)} \det A[B]\prod_{l\in L\setminus B} z_l.
\]
Thus the Pfaffian is a generating polynomial over parity bases, with coefficients \(\det A[B]\). In the weighted setting, if \(w:L\to \mathbb R\) and \(\theta\) is an indeterminate, then
\[
\pf\!\bigl(A\Delta(\theta^w)A^\top\bigr)
=
\sum_{B\in B(A,L)} \det A[B]\theta^{w(B)}.
\]
The coefficient of \(\theta^x\) is therefore
\[
\sum_{B\in B_x}\det A[B],\qquad
B_x:=\{B\in B(A,L):w(B)=x\}.
\]
This formulation is exact, but in general it does not by itself yield counts, because the coefficients \(\det A[B]\) can cancel algebraically [1912.00620].

A complementary weighted formulation uses the skew-symmetric polynomial matrix
\[
\Phi_A(\theta)=
\begin{pmatrix}
O & A\\
-A^\top & D(\theta)
\end{pmatrix},
\]
where \(D(\theta)\) is block diagonal with one \(2\times 2\) block per line \(\ell\):
\[
D_\ell(\theta)=
\begin{pmatrix}
0 & -\tau_\ell \theta^{w_\ell}\\
\tau_\ell \theta^{w_\ell} & 0
\end{pmatrix}.
\]
Its Pfaffian encodes the optimum weight through the identity
\[
\zeta(A,L,w)=\sum_{\ell\in L}w_\ell-\deg_\theta Pf\Phi_A(\theta),
\]
and \(Pf\Phi_A(\theta)=0\) if there is no parity base. This degree-of-Pfaffian viewpoint underlies the deterministic weighted algorithm described below [1905.13371].

## 3. Weighted optimization and deterministic polynomial-time algorithms

The weighted linear matroid parity problem remained open for decades in the represented setting, even though the unweighted linear case was already polynomial-time solvable. A deterministic polynomial-time algorithm is now known for the minimum-weight parity base problem. In the formulation with a row-full-rank matrix \(A\), row set \(U\), column set \(V\), and partition of \(V\) into lines \(L\), a parity base is a base \(B\) of \(M(A)\) such that \(B\) is a union of whole lines, and the objective is
\[
\zeta(A,L,w)
=
\min\left\{
\sum_{\ell\subseteq B} w_\ell
\;\middle|\;
B \text{ is a parity base}
\right\}.
\]
The main theorem is that Algorithm \({\sf Minimum\mbox{-}Weight\ Parity\ Base}\) finds a minimum-weight parity base, or proves infeasibility, in
\[
O(n^3m)
\]
arithmetic operations over \(K\). Over a fixed finite field this is strongly polynomial; over \(\mathbb Q\), polynomial-time solvability follows by reduction through a sequence of finite fields [1905.13371].

The algorithm is combinatorial and primal-dual in spirit, but it is formulated through a cocircuit matrix and Pfaffian algebra rather than through an explicit LP for matroid parity. For a current base \(B\), the fundamental cocircuit matrix is
\[
C=A[U,B]^{-1}A[U,V\setminus B].
\]
A standard fact is
\[
B\triangle S \text{ is a base } \iff C[S] \text{ is nonsingular}.
\]
Thus exchanges of elements are converted into nonsingularity tests on principal submatrices. The algorithm allows the current base to violate parity; lines containing exactly one chosen endpoint are source lines. The goal is to eliminate source lines by augmenting-path operations until the base becomes a parity base [1905.13371].

To support matching-style blossom manipulations, the method introduces an enlarged vertex set \(V^*\), a laminar family of blossoms \(\Lambda\), dummy lines formed by bud-tip pairs, a current augmented base \(B^*\), and an augmented matrix \(C^*\) obtained by pivoting. Dual variables consist of endpoint potentials \(p\) and blossom weights \(q\), constrained by
\[
p(v)+p(\bar v)=w_\ell
\]
for every line \(\ell=\{v,\bar v\}\), together with inequalities of the form
\[
p(v)-p(u)\ge Q_{uv},\qquad
Q_{uv}:=\sum_{i\in I_{uv}} q(H_i),
\]
for exchange edges \((u,v)\). Tight edges satisfy equality. Search is performed in the graph of tight edges, and augmenting paths are required to induce a unique tight perfect matching on their support, ensuring that the corresponding pivot is nonsingular [1905.13371].

Optimality is certified indirectly. Feasible dual variables \(p,q\) are embedded into the dual LP of a maximum-weight perfect matching problem on an auxiliary graph built from a transformed skew-symmetric matrix \(\Phi_A^*(\theta)\). Weak duality in that matching LP yields the required upper bound on \(\deg_\theta Pf\Phi_A(\theta)\), and hence on \(\zeta(A,L,w)\). This is why the algorithm behaves like a weighted blossom algorithm while remaining fundamentally algebraic [1905.13371].

## 4. Pfaffian parities and exact counting

For exact counting, the decisive obstruction is algebraic cancellation. In a general parity instance,
\[
\pf\!\bigl(A\Delta(z)A^\top\bigr)
=
\sum_{B\in B(A,L)} \det A[B]\prod_{l\in B} z_l
\]
has mixed coefficients \(\det A[B]\), so extracting the number of solutions is not straightforward. The subclass introduced to remove this obstruction is the Pfaffian parity: a matroid parity \((A,L)\) is Pfaffian if there exists \(c\in\mathbb K\setminus\{0\}\) such that
\[
\det A[B]=c \qquad \text{for all } B\in B(A,L).
\]
The value \(c\) is called the constant of \((A,L)\). In that case every feasible parity base contributes the same nonzero coefficient in the Pfaffian expansion, and the identities simplify to
\[
\pf\!\bigl(A\Delta(z)A^\top\bigr)
=
c\sum_{B\in B(A,L)}\prod_{l\in B} z_l,
\]
\[
\pf\Phi(z)
=
c\sum_{B\in B(A,L)}\prod_{l\in L\setminus B} z_l.
\]
At \(z=\mathbf 1\),
\[
|B(A,L)|=c^{-1}\pf\!\bigl(A\Delta(\mathbf 1)A^\top\bigr)=c^{-1}\pf\Phi(\mathbf 1).
\]
This is the parity analogue of Pfaffian orientations in matching and of Webb’s Pfaffian pairs in matroid intersection [1912.00620].

The resulting algorithmic consequences are strong. If \(\operatorname{char}(\mathbb K)=0\), one can count the number of parity bases of a \(2r\times 2n\) Pfaffian parity deterministically in
\[
O(nr^{\omega-1})
\]
time, and also construct one parity base in the same time. If the Pfaffian constant \(c\) is given, the number of parity bases can be computed modulo \(\operatorname{char}(\mathbb K)\) in deterministic
\[
O(nr^{\omega-1}+r^3)
\]
time, improving to \(O(nr^{\omega-1})\) in characteristic \(0\). For weighted Pfaffian parities with line weights \(w:L\to\mathbb R\), the number of minimum-weight parity bases modulo \(\operatorname{char}(\mathbb K)\) can be computed deterministically in
\[
O(n^3r).
\]
For \(\mathbb K=\mathbb Q\), the number of minimum-weight parity bases can be computed deterministically in time polynomial in the binary encoding length of \(A\) [1912.00620].

This tractable subclass includes several structurally important examples. If \(\vec G\) is a Pfaffian orientation of a graph \(G\), then the associated parity instance \((\vec A,E)\) is a Pfaffian parity, embedding classical Pfaffian perfect-matching counting into linear matroid parity. If a 3-uniform hypergraph has a 3-Pfaffian orientation, then Lovász’s reduction to graphic matroid parity yields a parity instance whose parity bases correspond to spanning hypertrees. The paper also proves that, in the \(|\mathcal S|=3\) case, weighted shortest disjoint \(S\)-\(T\)-\(U\) paths in an LGV position correspond bijectively to minimum-weight parity bases of a weighted Pfaffian parity [1912.00620].

A common misconception is that these results make exact counting easy for linear matroid parity in general. They do not. The tractable class is explicitly the class of Pfaffian parities, namely those instances for which feasible determinants are constant and algebraic cancellation disappears [1912.00620].

## 5. Fractional linear matroid parity, non-commutative rank, and quasi-parallelism

The fractional relaxation replaces integral line selection by a half-integral polytope. For lines \(\ell_i=\langle a_i,b_i\rangle\le K^n\), a vector \(y\in\mathbb R_{\ge 0}^m\) is a fractional matroid matching if for every subspace \(S\le K^n\),
\[
\sum_{i=1}^m \dim(S\cap \ell_i)\, y_i \le \dim S.
\]
The feasible region \(P\) satisfies \(P\subseteq [0,1]^m\), and every vertex of \(P\) is half-integral:
\[
y_i\in\{0,\tfrac12,1\}.
\]
Its integer vertices are exactly the integral matroid matching vertices. The dual problem is the minimum nested 2-cover problem, and the min-max theorem is
\[
2\max_{y\in P}|y|
=
\min_{(S,T):\ \text{nested 2-cover}} \bigl(\dim S+\dim T\bigr)
\]
[2207.07946].

The algebraic matrix used here is Lovász’s skew-symmetric symbolic matrix
\[
A=\sum_{i=1}^m x_i(a_i\wedge b_i),\qquad
a_i\wedge b_i:=a_i b_i^\top-b_i a_i^\top.
\]
Ordinary symbolic rank encodes the integral optimum. The key result is that non-commutative rank encodes the fractional optimum:
\[
\ncrank A = 2\max_{y\in P}|y|.
\]
Equivalently, with the second blow-up,
\[
\max_{y\in P}|y|
=
\frac14 \rank A^{\{2\}},
\]
so blow-up size \(2\) already suffices for this problem family. This matches the half-integrality of the fractional polytope: the rank \(0,1,2\) of a \(2\times2\) block corresponds to the values \(0,\tfrac12,1\) [2207.07946].

These identities yield algebraic algorithms for the fractional problem. A simple search-to-decision reduction finds the lexicographically minimum maximum fractional matroid matching in
\[
O(n^\omega + mn^2)
\]
field operations with probability at least \(1/2\) if \(|R|\ge 16mn\). A divide-and-conquer sparse-representation algorithm finds a maximum fractional matroid matching in
\[
O(mn^{\omega-1}),
\]
and a randomized algebraic algorithm finds the dominant 2-cover in
\[
O(mn^{\omega+1})
\]
time with probability at least \(1/2\). Combined with Gijswijt’s framework, these routines give a weighted fractional algorithm running in
\[
O(mn^{\omega+4})
\]
[2207.07946].

A later development establishes that fractional linear matroid matching is in quasi-NC. In that formulation, the same line-based problem is treated as the fractional relaxation of linear matroid parity, and the optimum is characterized by
\[
\max_{y\in P}|y| = \frac12 \nr(A)=\frac14 \rank(A^{\{2\}).
\]
The algorithm is deterministic and runs in \(\operatorname{polylog}(m)\) time using \(m^{O(\log m)}\) parallel processors. It also yields a deterministic black-box hitting set for non-commutative Edmonds’ problem with rank-two skew-symmetric coefficients. However, the exact integral linear matroid parity problem is not solved in quasi-NC; the conclusion stated is that “The parallel complexity of linear matroid matching is still open” [2402.18276].

## 6. Generalizations, methodological extensions, and current boundaries

Linear matroid parity also appears as a special case of a broader algebraic framework over linear and projected linear delta-matroids. In this setting, a key representation shift is the contraction representation \(D=(A)/T\), where \(A\) is skew-symmetric over \(V\cup T\) and a set \(F\subseteq V\) is feasible exactly when \(A[F\cup T]\) is nonsingular. This representation is equivalent to the standard twist representation up to \(O(n^\omega)\)-time transformations, but is more convenient for algorithmic tasks. For a linear matroid \(M\) represented by a matrix \(A\), its basis family becomes a linear delta-matroid, and ordinary matroid parity becomes a special case of DM Parity [2402.11596].

Within this framework, the union and delta-sum of linear or projected linear delta-matroids remain linear or projected linear, and several decision problems reduce to a single rank computation. In particular, the decision versions of DM Covering, DM Delta-Covering, DM Intersection, DM Partition, and DM Parity reduce to computing the rank of an \(O(n)\times O(n)\) skew-symmetric matrix. Search versions can then be solved in
\[
O(n^\omega)
\]
field operations with high probability by lifting Harvey’s inverse-update technique, the same algebraic style previously used for fast linear matroid parity algorithms. This improves over the \(O(n^4)\)-time augmenting path algorithm of Geelen, Iwata, and Murota, and over the \(O(n^{\omega+1})\) bound obtainable by self-reducibility [2402.11596].

These extensions do not supersede the classical linear matroid parity algorithms; rather, they subsume them into a larger skew-symmetric framework. At the same time, they clarify several current boundaries. Weighted perfect parity in the delta-matroid setting remains harder: the paper states that the more interesting weighted perfect parity version is equivalent to weighted DM intersection and only gets
\[
O(n^{\omega+1})
\]
search via self-reduction after computing the optimum weight [2402.11596]. Likewise, the determinant-degree method for weighted linear matroid intersection gives a linear-algebraic interpretation of Frank’s weight-splitting algorithm and is conceptually close to parity, but it is not itself a parity algorithm and does not provide a direct parity-specific complexity theorem [1908.11529].

The contemporary boundary is therefore sharply delineated. Exact integral optimization for linearly represented matroid parity is polynomial-time solvable, including the weighted case [1905.13371]. Exact counting is polynomial-time only for special cancellation-free subclasses such as Pfaffian parities [1912.00620]. Fractional parity is algebraically characterized by non-commutative rank and lies in quasi-NC [2207.07946], [2402.18276]. Exact integral quasi-NC for linear matroid parity remains open [2402.18276].

Source: https://www.emergentmind.com/topics/linear-matroid-parity-problem