Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

LP hardness of approximation for polynomial-time solvable problems

Contributed by Sebastian Pokutta
Published 2026-10-20 · last changed in 29ff25c
License: CC BY 4.0

Cite as
@misc{lemmata:2610.00004,
  title        = {LP hardness of approximation for polynomial-time solvable problems},
  author       = {Sebastian Pokutta},
  howpublished = {Lemmata, lemmata:2610.00004, \url{https://lemmata.iol.berlin/abs/2610.00004}},
  year         = {2026},
  url          = {https://lemmata.iol.berlin/abs/2610.00004},
  note         = {Version 29ff25c1aac1; formal statements checked in Lean 4}
}

Polynomial-time optimization does not imply the existence of a small LP formulation. This contribution gives approximation-preserving reductions for four problems: minimum-cost perfect matching, minimum-cost TT-join, minimum odd cut, and minimum-cost odd-red bipartite perfect matching. Each problem admits an exact polynomial-time algorithm. Nevertheless, shifted matching matrices yield strong lower bounds on a common LP formulation supporting all objectives.

We give the four applications together with the factorization argument that underlies them. The Lean development includes the actual graph and matching domains, objective sums, attaining solutions, and reductions. The published matching-matrix rank bound is an explicit theorem hypothesis. For parity matching, Lean also proves the new finite extension of the rectangle argument, from explicitly stated hypotheses describing the Jia–Svensson–Yuan partition framework. No external theorem is introduced as a Lean axiom.

For every fixed α>1\alpha>1, the implications are as follows. Here nn is the number of vertices in the first three rows and the number of vertices in each bipartition class in the last row.

ProblemRestricted objectivesLower bound on LP sizeNonnegative-rank lower bound used
Minimum-cost perfect matchingStrictly positive metric costs on KnK_n2Ωα(n)2^{\Omega_\alpha(n)}Shifted ordinary matching matrix
Minimum-cost TT-joinT=V(Kn)T=V(K_n), the same metric costs2Ωα(n)2^{\Omega_\alpha(n)}Shifted ordinary matching matrix
Minimum odd cutWeights in {0,1}\{0,1\} on KnK_n2Ωα(n)2^{\Omega_\alpha(n)}Transpose of the shifted ordinary matching matrix
Minimum-cost odd-red bipartite perfect matchingPositive costs on a fixed two-colored complete bipartite multigraph2Ωα(n)2^{\Omega_\alpha(n)}Shifted parity-matching matrix

The first three applications use the exponential nonnegative-rank bound for the shifted odd-cut versus perfect-matching matrix proved by Braun and Pokutta. The fourth uses a different matrix, from Jia, Svensson, and Yuan. Their published exact bound alone does not give the shift needed for a fixed-factor approximation. Below we prove the required extension; Lean checks its new construction, charging estimate, and finite rank consequence, while their partition and entropy assertions remain explicit hypotheses.

Formulations and source matrices

Fix a set II of instances, a set XX of solutions, and objective values v(i,x)v(i,x). An α\alpha-approximate formulation for this minimization problem has one polyhedron PP, with mm inequalities, a realization zx∈Pz_x\in P for every solution, and an affine objective ℓi\ell_i for each instance such that

ℓi(zx)=v(i,x),inf⁡z∈Pℓi(z)≥OPT⁡(i)α.\ell_i(z_x)=v(i,x),\qquad \inf_{z\in P}\ell_i(z)\ge \frac{\operatorname{OPT}(i)}{\alpha}.

Both PP and the realization x↦zxx\mapsto z_x are independent of the instance. The encoding is otherwise arbitrary. Equations can be removed by parametrizing their affine solution space and do not count toward mm. The lower bounds concern this model, which allows arbitrary extended formulations. They do not bound the running time of an algorithm. For the usual rational-input model, it suffices to use a fixed rational approximation factor: hardness for any rational factor larger than a given real α\alpha also rules out an α\alpha-approximate formulation. The resulting cost families have rational entries of polynomial encoding length.

Write r+(S)r_+(S) for the nonnegative rank of a real matrix SS: the least number of nonnegative rank-one summands, or ∞\infty if no finite factorization exists. The Lean definition is shared with the spanning-tree contribution.

For even nn, let Un\mathcal U_n be the odd subsets of [n][n] and let Mn\mathcal M_n be the perfect matchings of KnK_n. Set

qn(U,M)=∣δ(U)∩M∣,Qn,r(U,M)=qn(U,M)−r.q_n(U,M)=|\delta(U)\cap M|, \qquad Q_{n,r}(U,M)=q_n(U,M)-r.

The theorem of Braun and Pokutta used here states that

r+(Qn,r)=2Ωr(n)for every fixed 0<r<1.r_+(Q_{n,r})=2^{\Omega_r(n)} \quad\text{for every fixed }0<r<1.

In Section 2.2.3 and the proof of Theorem 3.1 of Braun and Pokutta, the relevant submatrix has entries (qn−1+β)/2(q_n-1+\beta)/2 for fixed 0<β<10<\beta<1. Taking β=1−r\beta=1-r and rescaling gives the stated nonnegative-rank bound. Rothvoss’s exact matching bound corresponds to r=1r=1; it is not by itself a bound for r<1r<1.

We use the formulation-complexity factorization theorem in its usual affine Farkas form, as in Braun, Pokutta, and Zink and Braun, Pokutta, and Roy. For precision, the formalization starts from the certificate data: if P={z:sl(z)≥0, 1≤l≤m}P=\{z:s_l(z)\ge0,\ 1\le l\le m\} and βi\beta_i is the claimed objective bound, it assumes nonnegative coefficients yil,κiy_{il},\kappa_i with

ℓi(z)−βi=∑l=1myilsl(z)+κifor all z.\ell_i(z)-\beta_i =\sum_{l=1}^m y_{il}s_l(z)+\kappa_i \quad\text{for all }z.

The Lean structure CertifiedLP contains the inequalities, feasible realizations, exact objective values, and this identity. The existence of these certificates from validity on a finite-dimensional polyhedron is the external affine Farkas theorem.

In each application below, LL is a nonnegative integer. A source bound L≤r+(S)L\le r_+(S) therefore gives L≤m+1L\le m+1. Substituting an exponential source bound gives the corresponding exponential lower bound on formulation size.

Minimum-cost perfect matching

Let n≥2n\ge2 be even, put H=n/2H=n/2, and fix α>1\alpha>1. For U∈UnU\in\mathcal U_n give the edges of KnK_n costs

ceU=1+λ 1[e∈δ(U)],λ=(α−1)n=2(α−1)H.c^U_e=1+\lambda\,\mathbf1[e\in\delta(U)], \qquad \lambda=(\alpha-1)n=2(\alpha-1)H.

In particular, cUc^U is a strictly positive metric. It is also a tree metric: attach each vertex to its side’s center by an edge of length 1/21/2, and join the centers by an edge of length λ\lambda.

Applying the bound r+(Qn,r)=2Ωr(n)r_+(Q_{n,r})=2^{\Omega_r(n)} of Braun and Pokutta with r=1/(2α)r=1/(2\alpha) gives m=2Ωα(n)m=2^{\Omega_\alpha(n)} for every fixed α>1\alpha>1, although weighted perfect matching is polynomial-time solvable. The costs are 1 and 1+(α−1)n1+(\alpha-1)n. Thus their aspect ratio grows with nn; this is not a bounded-aspect-ratio claim.

Minimum-cost TT-join

Again let n≥2n\ge2 be even and take T=V(Kn)T=V(K_n). A TT-join is an edge set JJ whose odd-degree vertex set is TT. Use the same costs cUc^U.

Thus minimum-cost TT-join has formulation size 2Ωα(n)2^{\Omega_\alpha(n)} already for T=V(Kn)T=V(K_n) and these positive metric costs. Minimum-cost TT-join is polynomial-time solvable, for example through shortest-path distances and minimum-weight perfect matching on TT.

Minimum odd cut

Let n≥2n\ge2 be even. The solutions are now odd subsets UU of [n][n]. Each perfect matching MM specifies an objective with edge weights weM=1[e∈M]w^M_e=\mathbf1[e\in M]. The cost of the cut is qn(U,M)q_n(U,M).

The resulting bound is 2Ωα(n)2^{\Omega_\alpha(n)} for fixed α>1\alpha>1. Polynomial-time solvability is supplied by the odd-cut algorithm of Padberg and Rao.

Minimum-cost odd-red bipartite perfect matching

Let GnG_n have left vertices u1,…,unu_1,\ldots,u_n and right vertices v1,…,vnv_1,\ldots,v_n, with one red and one blue edge between each opposite pair. A feasible solution is a perfect matching with an odd number of red edges. Equivalently, it is a permutation σ\sigma of [n][n] together with a Boolean color for each edge (ui,vσ(i))(u_i,v_{\sigma(i)}), with odd total red count. This is the concrete Lean solution type ParityMatching.

Rows are Boolean labelings aa of the 2n2n vertices such that the number of labels equal to one is congruent to nn modulo two. Call an edge a violation if it is blue with equal endpoint labels or red with unequal endpoint labels. Let tn(a,M)t_n(a,M) count the violations in MM and put

Pn,r(a,M)=tn(a,M)−r.P_{n,r}(a,M)=t_n(a,M)-r.

The nonnegative-rank bound needed for this application is

r+(Pn,r)=2Ωr(n)for fixed 0<r<1.r_+(P_{n,r})=2^{\Omega_r(n)} \quad\text{for fixed }0<r<1.

The mathematical proof of this separate shifted parity-matching assertion is given in Theorem 7. The original reduction accepts this rank bound as a hypothesis. The additional Lean theorem Corollary 4 instead derives a finite lower bound from explicit partition and probability-comparison hypotheses. The argument extends the rectangle estimates of Jia, Svensson, and Yuan; their exact bound at r=1r=1 alone would not suffice, and the ordinary matching bound of Braun and Pokutta concerns a different matrix.

Applying Theorem 7 gives 2Ωα(n)2^{\Omega_\alpha(n)} inequalities for every fixed α>1\alpha>1 on a host with 2n2n vertices and 2n22n^2 colored edges. All costs are positive and at most 1+2(α−1)n1+2(\alpha-1)n. No metric restriction is asserted here. The weighted odd-red problem is deterministically polynomial-time solvable; see the discussion of weighted parity matching in Jia, Svensson, and Yuan, Section 1. The condition is an odd red count, rather than an exactly prescribed red count.

The common reduction mechanism

The preceding arguments use simple scaling, restriction, and transposition. The accompanying library also proves the more general reduction used in Affine reductions for LPs and SDPs and Strong reductions for extended formulations.

The maps ff and gg are independent. Dependence on both indices may occur in AA and BB, and its nonnegative factorization size is accounted for in the inequality.

Bounded additive baselines

Objective reversal

The same contribution retains the near-exact objective-reversal transfer. This auxiliary result is not needed for the four fixed-factor applications above.

Interpreting oio_i as the source maximum requires the corresponding upper bound and an attaining solution. The library proves optimum transfer from those additional hypotheses; the corollary does not silently assume them.

Source bound for the shifted parity matrix

This section proves the additional rank bound for odd-red bipartite matching. The new extension is formalized in Lean, conditional on explicit assertions from the earlier partition framework. These assertions include the small- and bad-triple estimates, Equations (4.5)–(4.6), and the permutation lemma, Lemma 4.4, of Jia, Svensson, and Yuan, Sections 4.2–4.6. Their entropy proofs and the complete outside-partition construction are not formalized here. The one-violation construction, its admissibility, the new charging estimate, conversion to finite rank bounds, numerical parameter choice, and padding are proved in Lean. The precise conditional statements appear at the end of this section.

Write LL for a vertex labeling and ELE_L for its violating edges; thus t(L,M)=tn(L,M)=∣M∩EL∣t(L,M)=t_n(L,M)=|M\cap E_L| and Pn,r(L,M)=t(L,M)−rP_{n,r}(L,M)=t(L,M)-r. Fix a sufficiently large integer kk, let h=4k+3h=4k+3, and first take n=4k(2s+1)+3n=4k(2s+1)+3, where ss tends to infinity. A partition consists of ordered blocks A1,…,As,B1,…,BsA_1,\ldots,A_s,B_1,\ldots,B_s of 4k4k vertex pairs each, and an unordered core CC of hh pairs. On an ordered block DD, define MA(D)M_A(D) by the red edges (ui,v2k+i)(u_i,v_{2k+i}) and (u2k+i,vi)(u_{2k+i},v_i) for 1≤i≤2k1\le i\le2k. Let LB(D)L_B(D) assign one to both endpoints at positions 1,…,k,2k+1,…,3k1,\ldots,k,2k+1,\ldots,3k, and zero at the other positions. These matching and labeling choices have no violations.

Outside the core, fix MAM_A on every AA-block and choose any labeling consistent with it; fix LBL_B on every BB-block and choose any matching consistent with it. All such choices have the required even parities outside CC: consistency gives equal labels on each fixed red edge, and the red-count parity of a consistent matching is the number of pairs minus the number of labels equal to one. For a fixed rectangle R=L×M\mathcal R=\mathcal L\times\mathcal M, write pM(J)p_M(J) and pL(K)p_L(K) for the respective probabilities of membership in its factors when the core matching is JJ, or the core labeling is KK, and the permitted outside completions are sampled uniformly. The two completions are independent.

Formal statements for the extension

The finite statements below specify exactly which part of the preceding argument Lean proves. They use the same costs and violation matrices as the four applications.

For clarity, call the following the finite partition hypotheses. There is a finite set of partitions with nonnegative averaging weights and a finite set of pp triples. Each partition specifies its good triples, canonical and high-violation completions, and nonnegative row and column membership probabilities. The following assertions are supplied:

These are explicit fields of ConcreteJSYRectangleData. Its new matching and labeling are the concrete core constructions, not arbitrary replacement objects. Lean proves that they are permitted and derives their probability comparisons from the generic goodness condition. The sampling identities, permutation-family assertion, and small-/bad-triple estimates remain hypotheses; their proofs in the earlier framework are outside this formalization.

Scope of the formalization

The four final Lean theorems use actual combinatorial solutions. Perfect matchings are finite simple graphs of degree one, VV-joins are finite simple graphs with all degrees odd, odd cuts are odd finite vertex subsets, and colored bipartite matchings are permutations with edge colors. Edge counts and objective sums are defined on these objects.

The formal proofs include the one-crossing matching construction, degree-sum arguments, the singleton odd-cut optimum, the parity identity and one-violation matching construction, all four exact optimum calculations, and the transfers to L≤m+1L\le m+1. The metric property is proved for the actual cost function. The first three applications use a graph with nn vertices; the parity application uses a multigraph with 2n2n vertices.

The remaining external hypotheses and results are:

The original four transfer theorems remain available with an assumed source rank bound. The additional parity theorem instead derives that rank bound, in finite quantitative form, from the displayed partition hypotheses. The mathematical exponential statement also uses the exponential decay in the external small-/bad-triple estimates. The Lean development does not formalize the earlier paper in its entirety.

All four applications can be imported from one module:

Main.lean
/-
The four application theorems concern concrete combinatorial domains and
actual edge costs. Each concludes L ≤ m + 1 from an explicit shifted
nonnegative-rank lower bound and a formulation with affine certificates.
The parity extension below derives the rank bound from explicit JSY
partition and entropy hypotheses. General Farkas certificate existence
and the earlier source theorems remain external.
-/
export Matching (matching_lp_lower_bound oddCut_lp_lower_bound)
export Joins (vjoin_lp_lower_bound)
export Parity (parity_matching_lp_lower_bound)
export ParityExtension (parity_lp_lower_bound_from_jsy)
Formal verification
  • ✓ kernel-checked (Lean v4.35.0) Theorem 1 — Lemmata.LPHardnessInP.Formulation.certified_slack_factorization (source)

  • ✓ kernel-checked (Lean v4.35.0) Lemma 1 — Lemmata.LPHardnessInP.Matching.matching_minimum (source)

  • ✓ kernel-checked (Lean v4.35.0) Theorem 2 — Lemmata.LPHardnessInP.Matching.matching_lp_lower_bound (source)

  • ✓ kernel-checked (Lean v4.35.0) Proposition 1 — Lemmata.LPHardnessInP.CostTransfer.cut_metric_properties (source)

  • ✓ kernel-checked (Lean v4.35.0) Lemma 2 — Lemmata.LPHardnessInP.Joins.vjoin_minimum (source)

  • ✓ kernel-checked (Lean v4.35.0) Theorem 3 — Lemmata.LPHardnessInP.Joins.vjoin_lp_lower_bound (source)

  • ✓ kernel-checked (Lean v4.35.0) Lemma 3 — Lemmata.LPHardnessInP.Matching.oddCut_minimum (source)

  • ✓ kernel-checked (Lean v4.35.0) Theorem 4 — Lemmata.LPHardnessInP.Matching.oddCut_lp_lower_bound (source)

  • ✓ kernel-checked (Lean v4.35.0) Lemma 4 — Lemmata.LPHardnessInP.Parity.matchingCost_minimum (source)

  • ✓ kernel-checked (Lean v4.35.0) Theorem 5 — Lemmata.LPHardnessInP.Parity.parity_matching_lp_lower_bound (source)

  • ✓ kernel-checked (Lean v4.35.0) Theorem 6 — Lemmata.LPHardnessInP.RankReduction.factorization_reduction (source)

  • ✓ kernel-checked (Lean v4.35.0) Corollary 1 — Lemmata.LPHardnessInP.Formulation.certified_reduction_rank_lower_bound (source)

  • ✓ kernel-checked (Lean v4.35.0) Proposition 2 — Lemmata.LPHardnessInP.RankReduction.rank_transpose (source)

  • ✓ kernel-checked (Lean v4.35.0) Proposition 3 — Lemmata.LPHardnessInP.RankReduction.rank_rowScale_eq (source)

  • ✓ kernel-checked (Lean v4.35.0) Corollary 2 — Lemmata.LPHardnessInP.Formulation.bounded_baseline_lp_lower_bound (source)

  • ✓ kernel-checked (Lean v4.35.0) Corollary 3 — Lemmata.LPHardnessInP.Formulation.near_exact_reversal_lp_lower_bound (source)

  • ✓ kernel-checked (Lean v4.35.0) Proposition 4 — Lemmata.LPHardnessInP.ParityExtension.coreProbability_supported (source)

  • ✓ kernel-checked (Lean v4.35.0) Theorem 8 — Lemmata.LPHardnessInP.ParityExtension.ConcreteJSYRectangleData.corruption_bound (source)

  • ✓ kernel-checked (Lean v4.35.0) Theorem 9 — Lemmata.LPHardnessInP.ParityExtension.pairing_factorization_le (source)

  • ✓ kernel-checked (Lean v4.35.0) Proposition 5 — Lemmata.LPHardnessInP.ParityExtension.calibrated_witness_gap (source)

  • ✓ kernel-checked (Lean v4.35.0) Theorem 10 — Lemmata.LPHardnessInP.ParityExtension.calibrated_shifted_rank_from_jsy (source)

  • ✓ kernel-checked (Lean v4.35.0) Corollary 4 — Lemmata.LPHardnessInP.ParityExtension.parity_lp_lower_bound_from_jsy (source)

  • ✓ kernel-checked (Lean v4.35.0) Proposition 6 — Lemmata.LPHardnessInP.ParityExtension.shifted_rank_mono (source)

Checked at commit 29ff25c1aac1. The formal statement is machine-checked; that it matches the prose is checked by review. What “checked” means.

AI use. The exposition and Lean formalization were developed with OpenAI Codex. The published matching rank theorem, the earlier parity partition and entropy assertions, and affine Farkas certificate existence remain external. The new parity extension is formalized conditionally on precise hypotheses; no custom axioms are introduced.

Acknowledgements. Research reported in this paper was partially supported by the Berlin Mathematics Research Center MATH+^+ (EXC-2046/2, project ID 390685689), funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy. This work has been partially supported by OpenAI’s Academic Researchers program.

References
  1. Padberg, M. W., & Rao, M. R. (1982). Odd Minimum Cut-Sets and b-Matchings. Mathematics of Operations Research, 7(1), 67–80. 10.1287/moor.7.1.67