Approximate Parametric Computation of Minimum-Cost Flows with Convex Costs: Linear Splines of the Marginal Costs with Bounded Step Sizes Give an (α, β)-Approximate Minimum-Cost Flow FunctionResearch Paper
Motivation
Flow networks model the movement of material, traffic, electricity, or gas through connected systems. In many such models, the required net flow at each vertex varies with a parameter, so the object of interest is a whole family of optimal flows rather than one solution. When edge costs are nonlinear, solving a new optimization problem for every parameter value gives no finite description of that family. Joachims, Klimm and Warode study two approximations of this parametric problem: one approximates each edge's marginal cost, and the other interpolates selected optimal flows. This mission concerns their marginal cost approximation result, Theorem 5.
The paper builds on an exact method for piecewise quadratic edge costs by Klimm and Warode. A piecewise linear approximation of each marginal cost has a piecewise quadratic integral, so it brings general convex costs into the class handled by that method. The question is how fine the linear interpolation mesh must be for the resulting flows to be close to optimal under the original costs. The paper gives local step size bounds involving the second derivative of each marginal cost. The target here is a corrected version of that guarantee, with the precise corrections stated below.
Setting
A network consists of finite sets of vertices and edges . Every edge has a tail and a head. The graph is strongly connected and has neither self loops nor parallel directed edges. For an edge flow , its net inflow at a vertex is the sum of flows entering that vertex minus the sum leaving it. Two balanced demand vectors determine a demand for . A flow is feasible when its net inflow equals that demand. In the directed case all edge flows are nonnegative; in the undirected case their signs are unrestricted.
Each edge has a marginal cost that is strictly increasing, differentiable, nonpositive on negative arguments, and nonnegative on positive arguments. Its edge cost is . The total cost is the sum of edge costs, and is the minimum total cost of a feasible flow at parameter . These sign conditions make the edge costs nonnegative. The paper calls them homogeneous marginal costs, although this is a sign condition rather than algebraic homogeneity. Section 5.2 gives the model and its approximation criterion.
For each edge, choose a strictly increasing mesh from to , where
The linear spline joins the values of at adjacent mesh points by straight chords. The minimum cost flow function for these spline marginal costs is tested under the original costs: . This convention is essential because a small optimal cost for the spline instance alone says nothing about the cost under .
Formalization targets
The milestones follow the paper's argument in order.
- Edge-flow bound (§5.2.1, p. 16): every optimal edge flow satisfies .
- Equation (17) (proof of Lemma 4, p. 16): the closeness condition (15) of marginal costs on gives
- The (16) case (proof of Lemma 4, p. 17): and (16) give .
- Lemma 4 (p. 16): uniform control of on in the sense (15) or (16) makes the minimum-cost flow function of the instance an -approximation.
- Equation (18) (p. 17): the interpolation estimate
- Spline domination under (ii) (proof of Theorem 5, p. 18): on .
- Endpoint maximum under (ii) (proof of Theorem 5, p. 18): .
The goal is the corrected Theorem 5. It states that either of the paper's two explicit step size regimes, involving , , , and , makes the spline instance's minimum cost flow function satisfy
The first regime bounds each step by a local maximum of and the magnitude of at a mesh endpoint. The second adds monotonicity and sign conditions on and uses its endpoint values. The formal theorem retains the constants , the factor in (18), and both distinct denominators in Theorem 5. Theorem 5, pp. 17–18 is the source for the exact formulas.
Significance
The guarantee applies to an entire interval of demand parameters. It converts local numerical accuracy of the marginal cost splines into a bound on the original total cost at every parameter value, while preserving feasibility. The edge flow bound explains why it suffices to approximate marginal costs on a finite interval even when edge flows are otherwise real valued.
The paper proves the intended result by a paper argument, but its printed Theorem 5 has a counterexample when a mesh crosses zero without including zero: the spline can lose the marginal cost sign condition, and at zero demand the spline optimum can have positive original cost while the true optimum has cost zero. The corrected goal records a valid target for machine checked proof. Formalizing it also yields reusable definitions of finite flow networks, convex separable edge costs, parametric approximation, and chord splines. No machine checked proof of this theorem is claimed by this mission.
Difficulty
Pointwise closeness of marginal costs does not directly compare the two optimized flows, which may differ on every edge. A direct comparison at the spline optimum must also account for the original optimum and for how the edge costs accumulate across the network. The interpolation estimate alone does not settle this: a mesh interval crossing zero can change the sign of its chord at zero, even when every step size bound in the printed theorem holds. The shape condition in Theorem 5(ii), read literally with the third differentiability assumption, also forces the second derivative to vanish on the negative half line. The mission keeps that condition as printed rather than silently replacing it with a more symmetric one.
Formalization scope
Lean represents vertices and edges by finite types, net inflow by finite sums, and each edge cost by a real interval integral based at zero. A Boolean flag selects the directed case with or the undirected case with no sign constraint. The general finite lower and upper bounds of problem (14) are outside this mission. The optimal cost is the real infimum of feasible costs; the graph and balanced demand assumptions give feasible flows, and the sign assumptions bound their costs below. The spline is computed from its mesh values, with its end chords extended linearly outside . Meshes are indexed from zero in Lean, regardless of the paper's two indexing styles.
The corrected Theorem 5 requires zero to be a mesh point on every edge. In regime (i), it replaces by . In regime (ii), both entries of the factor are square roots of ; this matches the equality used in the printed proof. Lemma 4 and Theorem 5 explicitly require , since both step-size bounds and (15), (16) divide by it (with zero demands, a one-point mesh and , Theorem 5 would otherwise fail), and its second marginal cost family is continuous, strictly increasing, and sign preserving. These conventions and corrections are recorded in the source analysis for Theorem 5.
A proof must use the chord spline determined by and the mesh; replacing it by an arbitrary close function would discard the numerical claim. Likewise, the bound is never a hypothesis of Lemma 4 or of the goal: it is milestone 1, and assuming it would remove the combinatorial half of the argument. Milestone 6 is stated on rather than on all of as printed, because the linearly extended end chord of a convex function lies below it. Contributions to the flow bound, Lemma 4, the spline estimate, or the final connection between them are all within scope. The flow and spline definitions can also support other approximation results beyond this theorem.
Selected references
- Per Joachims, Max Klimm and Philipp Warode, Approximate Parametric Computation of Minimum-Cost Flows with Convex Costs, arXiv preprint, 2022. arXiv:2203.13146v1.
- Max Klimm and Philipp Warode, Parametric Computation of Minimum-Cost Flows with Piecewise Quadratic Costs, Mathematics of Operations Research, articles in advance, 2021. DOI:10.1287/moor.2021.1151.