Introduction

In a mean-payoff game, two players choose edges of a finite graph and compete over the long-run average weight of the resulting infinite walk. The objective is numerical, but the question of who can enforce a nonnegative average is a finite decision problem. A central algorithmic issue is the representation of the weights: a procedure polynomial in their absolute values may take exponential time in the length of the input. We give a deterministic algorithm whose running time is quasipolynomial in that binary length and which determines all winning starting vertices together.

The problem and the main theorem

Let G=(V,E)G=(V,E) be a finite directed graph with V={1,…,n}V=\{1,\ldots,n\} and n≥1n\ge1. Each vertex has at least one outgoing edge and belongs to exactly one of VMaxV_{\mathrm{Max}} and VMinV_{\mathrm{Min}}. The edge weights w:E→Zw:E\to\mathbb{Z} are arbitrary signed binary integers. Self-loops and parallel edges are allowed; parallel edges are separate entries of the explicit edge list. At a vertex its owner chooses the next edge. Both players’ strategies may depend on the entire finite history. For an infinite play π=(e0,e1,…)\pi=(e_0,e_1,\ldots), write

MP⁡w(π)=lim inf⁡T→∞1T∑t=0T−1w(et).\operatorname{MP}_{w}(\pi)=\liminf_{T\to\infty}\frac{1}{T}\sum_{t=0}^{T-1}w(e_t).

We count the complete input length LL in a fixed conventional binary encoding of the graph, endpoints, ownership, and weights. For a start vv and strategies σ,τ\sigma,\tau, denote their induced play by πv,σ,τ\pi_{v,\sigma,\tau}.

Theorem 1.1. A uniform deterministic algorithm returns exactly the set

U={v∈V:∃Max⁡ strategy σ ∀Min⁡ strategies τ, MP⁡w(πv,σ,τ)≥0}.U=\{v\in V:\exists\operatorname{Max}\text{ strategy }\sigma\ \forall\operatorname{Min}\text{ strategies }\tau,\ \operatorname{MP}_{w}(\pi_{v,\sigma,\tau})\ge0\}.

for every game specified above. For an absolute constant C>0C>0, the complete computation uses at most

2C(log⁡2(L+2))22^{C(\log_2(L+2))^2}

bit operations in a standard model polynomially equivalent to a uniform deterministic Turing machine. The bound includes preprocessing, exact arithmetic, recursion, and output of the entire set UU.

A positional strategy chooses its outgoing edge using only the current vertex. Section 7 derives the exact rational value at every starting vertex and globally optimal positional strategies for both players, within the same quasipolynomial bit bound. Section 8 decides whether a prescribed initial credit can keep every finite accumulated weight nonnegative in a one-resource energy game, and computes minimum winning credits. Section 9 applies classical reductions to tropical feasibility and fixed-corank tropical-rank tests. A polynomial bound for the complete binary input would be a further improvement.

Earlier algorithms and binary weights

Ehrenfeucht and Mycielski proved the existence of positional optimal strategies for mean-payoff games [13]. Positional determinacy places the threshold problem in NP ∩\cap coNP. Zwick and Paterson developed pseudopolynomial algorithms and complexity reductions [30]. Such bounds can be polynomial in the number of vertices and the largest absolute weight while remaining exponential in the binary encoding length. This distinction is central to Theorem 1.1.

The energy formulation asks whether a finite initial credit can keep all finite accumulated weights nonnegative. Brim, Chaloupka, Doyen, Gentilini, and Raskin developed bounded integer progress measures and monotone lifting, with an O(∣E∣nW0)O(|E|nW_0) threshold and winning-region bound when W0W_0 bounds the absolute weights [7] [Theorem 8]. Dorfman, Kaplan, and Zwick developed a deterministic exponential-time approach based on accelerated potential updates and scaling; Austin and Dell’Erba corrected its update procedure [11, 3]. Kozachinskiy gives a polyhedral interpretation and a weight-independent nO(1)2n/2n^{O(1)}2^{n/2} arithmetic bound for energy winners [19].

Strategy-improvement algorithms provide another important comparison. Björklund and Vorobyov proved an expected 2O(nlog⁡n)2^{O(\sqrt{n}\log n)} arithmetic bound for the threshold partition on graphs without parallel edges [5] [Theorem 7.1], following work with Sandberg [4]. Recent analyses concern Random-Action-Removal [29] and Switch-All and Random-Edge for energy games [12]. Iteration counts and arithmetic bounds must be distinguished from bit bounds: strategy evaluation and the sizes of intermediate numbers contribute to the latter.

Parity games use an extremal priority occurring infinitely often to determine the winner. They reduce in polynomial time to mean-payoff threshold games using weights of exponential magnitude but polynomial binary length [16] [Section 3]. Calude, Jain, Khoussainov, Li, and Stephan gave a quasipolynomial algorithm for parity games [9]. The reduction goes from parity to mean-payoff, so this result does not by itself give a quasipolynomial algorithm for arbitrary binary-weight mean-payoff games. Daviaud, Jurdziński, and Lazić combined succinct parity progress measures with energy bounds in a pseudo-quasipolynomial algorithm for mean-payoff parity games [10]; that bound still contains a numerical weight parameter. Smoothed polynomial results for independent Gaussian perturbations of payoffs on ergodic game graphs concern a different input model [22].

How the proof works

An integer potential assigns a number to each vertex. First replace w(e)w(e) by w^(e)=(n+1)w(e)+1\widehat{w}(e)=(n+1)w(e)+1, which makes every simple cycle sum nonzero while separating negative original sums from nonnegative ones. For z∈ZVz\in\mathbb{Z}^{V}, let Fi(z)F_i(z) be the maximum of w^(e)+zj\widehat{w}(e)+z_j over edges e:i→je:i\to j at a Max vertex, and the minimum at a Min vertex. For integer vectors A≤uA\leq u, the box [A,u][A,u] consists of all integer vectors between them, coordinatewise. In this box, a subsolution is a vector x∈[A,u]x\in[A,u] satisfying xi≤Fi(x)x_i\leq F_i(x) wherever xi>Aix_i>A_i; a supersolution is a vector y∈[A,u]y\in[A,u] satisfying yi≥Fi(y)y_i\geq F_i(y) wherever yi<uiy_i<u_i. Clipping a vector to the box means replacing coordinate ziz_i by min⁡{ui,max⁡{Ai,zi}}\min\{u_i,\max\{A_i,z_i\}\}. The comparison theorem in Randomized quasipolynomial-time mean-payoff games, our companion article, gives a unique fixed point hh of z↦min⁡{u,max⁡{A,F(z)}}z\mapsto\min\{u,\max\{A,F(z)\}\}, and the inequalities x≤h≤yx\leq h\leq y for every such pair. These finite-box facts and their basic iteration procedure are the only results imported from the companion in the proof of Theorem 1.1 [24] [Definition 2.2 and Lemmas 2.3–2.4]. Section 2 states that interface in full.

The same section places each coordinate of hh near one of the two boundaries of a sufficiently wide box [0,D]V[0,D]^V. The upper band is winning for Max and the lower band for Min. Iterating the clipped operator from the bottom can require a number of steps proportional to the numerical width DD. We instead ask which boundary labels are forced, without computing hh.

For a box [A,A+D][A,A+D], give the vertices positive rational masses aia_i and bib_i. For a vertex set SS, its masses are a(S)=∑i∈Saia(S)=\sum_{i\in S}a_i and b(S)=∑i∈Sbib(S)=\sum_{i\in S}b_i. A Max witness is a subsolution whose support {i:xi>Ai}\{i:x_i>A_i\} has aa-mass at most one; it requires a positive label at coordinates near the upper boundary. A Min witness is the dual supersolution with support {i:yi<Ai+D}\{i:y_i<A_i+D\} of bb-mass at most one; it requires a negative label near the lower boundary. The task is to meet every requirement of every witness at once. With initial masses 1/n1/n, the box solution qualifies on both sides, so these requirements determine every vertex’s winner.

To solve that task recursively, we move an integer pivot vector kk near the middle of the box. For a Max witness, a common translation aligns the largest deviation xi−kix_i-k_i with the top of a narrower box centered at kk. Clipping at the bottom leaves only coordinates whose deviations are close to that largest value. A dual operation handles Min witnesses. Section 3 proves that these operations preserve the required inequalities and restrict supports to subsets of the original supports.

The algorithm in Section 4 first moves kk downward, then upward. Narrow recursive calls choose which coordinates move. Each pass protects the extremal deviation of all witnesses of one type. When a pass stops, any still-needed witness of the other type must have a set of coordinates close to an extremum with aa-mass greater than 3/43/4 for a Max witness or bb-mass greater than 3/43/4 for a Min witness. Otherwise the last recursive call would have forced another movement. This use of the final, stationary call is the key deterministic certificate.

One call on a half-width box labels those large-mass sets. Its labels need not settle the original boundary requirements. They instead specify new masses: coordinates labelled in a witness’s favor become lighter, and the others heavier. Section 5 shows that every unsettled witness still has support mass at most one. Meanwhile, every product aibia_i b_i increases by a fixed factor. Once every product exceeds one, each coordinate is excluded from at least one type of witness support, and the labels can be chosen directly. Starting with masses 1/n1/n, only O(log⁡(n+1))O(\log(n+1)) such product increases can therefore occur along a branch. The final call on the original box uses the reweighted masses to settle the remaining requirements. Every child except the half-width call increases the products; that exceptional child instead reduces the width. Section 6 counts the possible positions of the product increases along a branch and bounds the exact arithmetic in each call, giving the stated bit complexity.

Potential methods and recursive precision

Gurvich, Karzanov, and Khachiyan developed potential transformations in their algorithm for cyclic games [15, 18]: adding differences of endpoint potentials to edge weights preserves cycle sums while changing the local inequalities. Paired potentials appear in the work of Lishfits and Pavlov [21]. Cadilhac, Casares, and Ohlmann give a general framework for accelerated value iteration [8]. These methods and the progress-measure approach described above provide the setting for our finite-box inequalities.

The closest recursive ancestry is Parys’s reduced-precision parity algorithm [25] (Algorithm 2, Section 5, and Lemma 6.1) and its development by Lehtinen, Parys, Schewe, and Wojtczak [20]. Their guarantees protect all dominions—sets within which one player can keep play and win—within specified size bounds, and their recursion separates one child retaining full precision from children of reduced precision. We likewise use simultaneous obligations for small objects, but smallness here means rational mass of a potential’s support. Translation near an extremal deviation, the stationary-pivot certificates, and reweighting establish the needed progress for that notion of smallness.

Ohlmann’s symmetric mean-payoff recursion supplies closely related potential-based context [23]; its subexponential time analysis is left open in the cited version. Jurdziński, Morvan, Ohlmann, and Thejaswini use two simultaneous labellings for symmetric acceleration of parity algorithms [17]; a related development by Thejaswini, Ohlmann, and Jurdziński retains both decompositions to accelerate symmetric recursion [26]. Arnold, Niwiński, and Parys study recursive interval restriction in a quasipolynomial black-box algorithm for nested fixed-point evaluation [2]. Their bound is expressed in Boolean-lattice dimension and fixed-point nesting depth. The integer box here has chains of length nD+1nD + 1, so that result does not directly give the dependence on log⁡D\log D established by our translation and mass arguments.

Fijalkow, Gawrychowski, and Ohlmann studied mean-payoff value iteration using universal graphs and lower bounds for those representations [14]. Those lower bounds constrain universal-graph representations. The present recursion constructs and reweights boxes for the given instance without using such a representation.

The companion algorithm [24] returns potential vectors with probabilistic comparison guarantees stated separately for each fixed admissible vector. Here the recursive guarantee is simultaneous over all qualifying witnesses and concerns their boundary labels. The two algorithms share the finite-box theory stated in Section 2; the companion’s randomized proof is independent of the deterministic recursion. Section 7 also explains how its winning-set interface yields randomized exact values and positional strategies.

From potential comparisons to a labelling problem

This section specifies the finite-box input used from the companion, then proves the game reduction at the width required by our recursion. The outcome is a labelling problem whose obligations are expressed by potential vectors. Later sections solve that problem without finding the potential that certifies the winning partition.

The common finite-box interface

Fix the input graph and ownership throughout. Following the integer cycle perturbation of Björklund and Vorobyov [5], put

w^(e)=(n+1)w(e)+1,W=max⁡e∈E∣w^(e)∣.(1)\widehat{w}(e) = (n+1)w(e) + 1,\qquad W = \max_{e\in E} \lvert\widehat{w}(e)\rvert. \tag*{(1)}

The edge set is nonempty. Each transformed weight is a nonzero integer, so W≥1W \ge1. If a simple directed cycle has length ℓ\ell and original sum cc, then

w^(C)=(n+1)c+ℓ,1≤ℓ≤n.(2)\widehat{w}(C) = (n+1)c + \ell,\qquad1 \le\ell\le n. \tag*{(2)}

For c≥0c \ge0 the transformed sum is positive; for c≤−1c \le-1 it is negative. Thus no simple cycle has zero transformed sum, and positive transformed sums correspond exactly to nonnegative original sums. A self-loop has ℓ=1\ell= 1, and the calculation applies to every actual edge sequence when parallel edges are present.

For all z∈ZVz \in\mathbb{Z}^{V}, define

Fi(z)={max⁡e:i→j(w^(e)+zj),i∈VMax,min⁡e:i→j(w^(e)+zj),i∈VMin.(3)F_i(z) = \begin{cases} \max_{e:i\to j}\bigl(\widehat{w}(e)+z_j\bigr), & i\in V_{\mathrm{Max}},\\ \min_{e:i\to j}\bigl(\widehat{w}(e)+z_j\bigr), & i\in V_{\mathrm{Min}}. \end{cases} \tag*{(3)}

The finite, nonempty outgoing-edge sets ensure that these extrema are attained. With coordinate-wise vector order and common scalar shifts,

z≤z′⟹F(z)≤F(z′),F(z+t)=F(z)+t(t∈Z).(4)z \le z' \Longrightarrow F(z) \le F(z'),\qquad F(z+t)=F(z)+t\quad(t\in\mathbb{Z}). \tag*{(4)}

For integer vectors A≤uA \leq u, the notation [A,u][A,u] denotes all integer vectors between them. A subsolution is a full vector x∈[A,u]x \in[A,u] satisfying xi≤Fi(x)x_i \leq F_i(x) at every coordinate where xi>Aix_i > A_i. A supersolution is a full vector y∈[A,u]y \in[A,u] satisfying yi≥Fi(y)y_i \geq F_i(y) wherever yi<uiy_i < u_i. No inequality is required on the respective excluded boundary. These conventions agree with the companion’s Definition 2.2 [24].

Lemma 2.1 (Finite-box comparison and basic iteration). For the operator in Equation (3) and arbitrary finite integer bounds A≤uA \leq u, every subsolution xx and supersolution yy in [A,u][A,u] satisfy xi≤yix_i \leq y_i for all i∈Vi \in V. There is a unique h∈[A,u]h \in[A,u] that is both a subsolution and a supersolution. It is precisely the unique fixed point of

TA,u(z)=min⁡{u,max⁡{A,F(z)}}.T_{A,u}(z) = \min\{u,\max\{A,F(z)\}\}.

Consequently every such pair satisfies x≤h≤yx \leq h \leq y.

The following synchronous iteration finds hh. Start with p=Ap = A; evaluate q=TA,u(p)q = T_{A,u}(p); return pp if q=pq = p, and otherwise set p=qp = q and repeat. It uses at most

1+∑i∈V(ui−Ai)1 + \sum_{i \in V}(u_i-A_i)

evaluations of TA,uT_{A,u}, including the final stability test.

We use this statement from the complete September 25 companion, Lemmas 2.3–2.4 [24]. It permits negative endpoints and collapsed coordinates Ai=uiA_i = u_i. In particular its comparison conclusion concerns all vertices, even though each hypothesis is imposed only on its stated support. The operator and integer hypotheses are exactly those above, and Equation (2) verifies the required absence of zero simple-cycle sums. The sandwich inequality follows by comparing with hh on each side. We call hh the box solution and the iteration Basic. The algorithm will invoke Basic only when the box width is bounded by an absolute constant.

The two bands of a wide box

Choose the initial dyadic width

D0=min⁡{2d:d∈Z≥0, 2d≥64(n+1)(W+1)}.(5)D_0 = \min\{2^d : d \in\mathbb{Z}_{\geq0},\ 2^d \geq64(n+1)(W+1)\}. \tag*{(5)}

The next proposition establishes the game reduction for this width. The companion uses width 8nW8nW in its own reduction; that different choice is not an input to the proof here.

Proposition 2.2 (Winning partition from the wide box). Let hh be the box solution in [0,D0]V[0,D_0]^V. Then

Vlo={i:hi≤nW},Vhi={i:hi≥D0−nW}V_{\mathrm{lo}} = \{i : h_i \leq nW\}, \qquad V_{\mathrm{hi}} = \{i : h_i \geq D_0-nW\}

partition VV. There is a positional Max strategy under which every play starting in VhiV_{\mathrm{hi}} has original mean payoff at least zero. There is a positional Min strategy under which every play starting in VloV_{\mathrm{lo}} has original upper limiting average at most −1/n-1/n. Both assertions allow arbitrary history-dependent opposition.

Proof. We first locate the coordinates of hh, then construct the two strategies. If 0<hi<D00 < h_i < D_0, both box inequalities apply and give hi=Fi(h)h_i = F_i(h). Follow an edge attaining that equality and continue from its head as long as the current coordinate is interior. A first repeated vertex would yield a simple cycle of tight equalities hi=w^(e)+hjh_i = \widehat{w}(e) + h_j, with transformed sum zero. This is impossible by (2). A boundary is therefore reached within nn edges. Each edge changes the potential by at most WW. Including the length-zero path from a boundary coordinate, every coordinate is within nWnW of either 00 or D0D_0. Since D0>2nWD_0 > 2nW, the two bands are disjoint and exhaust VV.

Every high coordinate is above the lower boundary. At a high Max vertex choose one outgoing edge satisfying

hi≤w^(e)+hj;(6)h_i \le\widehat{w}(e) + h_j; \tag*{(6)}

such an edge exists by the subsolution inequality. At a high Min vertex every outgoing edge satisfies this inequality. Each of these permitted edges has

hj≥D0−(n+1)W>nW,h_j \ge D_0 - (n+1)W > nW,

so its head is high by the partition just proved. Fixing the Max choices and extending them arbitrarily outside the high band gives a positional strategy keeping all consistent plays from high vertices in that band.

On a simple cycle of permitted edges, Equation (6) telescopes to a nonnegative transformed sum. Its original sum is therefore nonnegative. Let W0=max⁡e∣w(e)∣W_0 = \max_e |w(e)|. Delete simple cycles successively from any finite play prefix until fewer than nn edges remain. The removed cycles have nonnegative original sums, and the remaining walk has weight at least −(n−1)W0-(n-1)W_0. Hence every prefix of every consistent play has at least that weight. Dividing by its length and taking a lower limit proves Max’s assertion.

For the other band, every low coordinate is below the upper boundary. At a low Min vertex fix an edge attaining the minimum in Fi(h)F_i(h), and at a low Max vertex permit all outgoing edges. The supersolution inequality gives

hi≥w^(e)+hj(7)h_i \ge\widehat{w}(e) + h_j \tag*{(7)}

for each permitted edge. Its head satisfies hj≤(n+1)W<D0−nWh_j \le(n+1)W < D_0 - nW, so it is low. These choices, extended arbitrarily elsewhere, give a positional Min strategy preserving the low band.

Every simple cycle of permitted edges now has nonpositive transformed sum and hence original integer sum at most −1-1. Deleting simple cycles from a prefix of length TT leaves at most n−1n-1 edges. Each deleted cycle has at most nn edges, so the number removed is at least (T−(n−1))/n(T-(n-1))/n. The original prefix weight is at most

−T−(n−1)n+(n−1)W0.-\frac{T-(n-1)}{n} + (n-1)W_0.

Its upper limiting average is at most −1/n-1/n. The prefix estimates on both sides apply to every consistent edge sequence, independently of how the opposing player uses the history.

Simultaneous boundary obligations

The recursion works on inputs

A∈ZV,D=2d (d∈Z≥0),a,b∈Q>0V.(8)A \in\mathbb{Z}^{V}, \qquad D = 2^d\ (d \in\mathbb{Z}_{\ge0}), \qquad a,b \in\mathbb{Q}_{>0}^{V}. \tag*{(8)}

Its box is [A,u][A,u] with u=A+Du = A + D, and its claim margin is P=D/32P = D/32. For S⊆VS \subseteq V let a(S)=∑i∈Saia(S) = \sum_{i \in S} a_i and b(S)=∑i∈Sbib(S) = \sum_{i \in S} b_i. The positive rational vectors a,ba,b are called masses.

Definition 2.3 (Witnesses and their claims). A Max witness for (A,D,a,b)(A,D,a,b) is a subsolution x∈[A,u]x \in[A,u] with support Sx={i:xi>Ai}S_x = \{i : x_i > A_i\} and a(Sx)≤1a(S_x) \le1. It claims each ii with xi≥ui−Px_i \ge u_i - P. A Min witness is a supersolution y∈[A,u]y \in[A,u] with support Sy={i:yi<ui}S_y = \{i : y_i < u_i\} and b(Sy)≤1b(S_y) \le1. It claims each ii with yi≤Ai+Py_i \le A_i + P. A witness is relevant when it claims at least one index.

Since P<DP < D, every claim belongs to its witness’s support. Claims of opposite types cannot coincide: otherwise box comparison would give ui−P≤xi≤yi≤Ai+Pu_i - P \le x_i \le y_i \le A_i + P, contradicting 2P<D2P < D. Masses limit the witnesses covered by an obligation; the graph itself always retains every vertex and edge.

Theorem 2.4 (Simultaneous labelling). For every input in (8), the procedure Label⁡(A,D,a,b)\operatorname{Label}(A,D,a,b) of Section 4 terminates and returns a vector in {+,−}V\{+,-\}^{V}. It assigns ++ to every index claimed by any Max witness and −- to every index claimed by any Min witness.

All witnesses of a call are covered by this one output. They are neither enumerated nor supplied to the procedure. On vertices with no claim the theorem permits either label.

Apply this guarantee to

A=0,D=D0,ai=bi=1/n.A = 0,\qquad D = D_0,\qquad a_i = b_i = 1/n.

Every support then has mass at most one. The box solution hh is therefore a witness of both types. Moreover P=D0/32≥nWP = D_0/32 \ge nW, so its Max claims include all high vertices and its Min claims include all low vertices. By Proposition 2.2, a labelling satisfying Theorem 2.4 returns the exact winning set by taking its positively labelled vertices. This uses the existence and boundary location of hh, without executing Basic on the wide initial box.

Moving witnesses to a narrower box

A recursive child will receive a box centered at a pivot kk near the middle of its parent box. To transfer a parent’s obligations to this child, we translate a witness by one scalar and clip it at a boundary. This section establishes the exact support and claims of the translated vector. It is the geometric part of the argument and uses only monotonicity, common-shift invariance, and the witness definition.

Consider an input (A,D,a,b)(A,D,a,b) with dyadic D≥4096D \ge4096. In addition to u=A+Du = A + D and P=D/32P = D/32, fix

c=A+D/2,M=D/8,q=D/128,p=q/32=D/4096.c = A + D/2,\qquad M = D/8,\qquad q = D/128,\qquad p = q/32 = D/4096.

The vector cc is the center, MM bounds the pivot’s displacement, qq is the narrow child width, and pp will be one pivot step. All these quantities, as well as q/2q/2 and D/4D/4, are integral. We consider integer pivots c−M≤k≤c+Mc - M \le k \le c + M.

For a relevant Max witness xx, let

sx(k)=max⁡i∈V(xi−ki),gix(k)=sx(k)−(xi−ki).s_x(k) = \max_{i \in V}(x_i-k_i),\qquad g_i^x(k) = s_x(k) - (x_i-k_i).

We call sx(k)s_x(k) its peak and gix(k)g_i^x(k) its gap at ii. For a relevant Min witness yy the dual definitions are

sy(k)=max⁡i∈V(ki−yi),giy(k)=sy(k)−(ki−yi).s_y(k) = \max_{i \in V}(k_i-y_i),\qquad g_i^y(k) = s_y(k) - (k_i-y_i).

Every gap is a nonnegative integer. The following lemma treats the two widths that the recursion will use.

Lemma 3.1 (Translated supports and claims). Let (A,D,a,b)(A,D,a,b) be an input with D≥4096D \ge4096, let c−M≤k≤c+Mc - M \le k \le c + M be integral, and let r∈{q,D/2}r \in\{q,D/2\}. For any relevant Max witness xx, the vector

xi′=max⁡{ki−r/2, xi−sx(k)+r/2}(9)x'_i = \max\{k_i-r/2,\ x_i-s_x(k)+r/2\} \tag*{(9)}

is an integer subsolution in [k−r/2,k+r/2][k-r/2,k+r/2]. Its support is

Sx′={i:gix(k)<r}⊆Sx.S_{x'}=\{i:g_i^x(k)<r\}\subseteq S_x.

For any positive new Max masses a′a' with a′(Sx′)≤1a'(S_{x'})\le1, it is a Max witness and claims every ii with gix(k)≤r/32g_i^x(k)\le r/32.

For any relevant Min witness yy, the vector

yi′=min⁡{ki+r/2,yi+sy(k)−r/2}(10)y_i'=\min\{k_i+r/2,y_i+s_y(k)-r/2\} \tag*{(10)}

is an integer supersolution in the same box, with support

Sy′={i:giy(k)<r}⊆Sy.S_{y'}=\{i:g_i^y(k)<r\}\subseteq S_y.

For any positive new Min masses b′b' with b′(Sy′)≤1b'(S_{y'})\le1, it is a Min witness and claims every ii with giy(k)≤r/32g_i^y(k)\le r/32. In particular, the corresponding unchanged masses aa and bb always make these translated vectors witnesses.

Proof. The support assertion requires knowing which coordinates could survive the clipping. Relevance supplies that information. A Max claim at jj gives xj≥Aj+D−Px_j\ge A_j+D-P, while kj≤Aj+D/2+Mk_j\le A_j+D/2+M; hence

sx(k)≥D/2−M−P.s_x(k)\ge D/2-M-P.

If i∉Sxi\notin S_x, then xi=Aix_i=A_i and xi−ki≤−D/2+Mx_i-k_i\le-D/2+M. Combining the two inequalities gives

gix(k)≥D−2M−P=23D/32>D/2.g_i^x(k)\ge D-2M-P=23D/32>D/2.

The untruncated term in (9) is

ki+r/2−gix(k).k_i+r/2-g_i^x(k).

It lies below or at the new upper boundary. Taking its maximum with the lower boundary puts x′x' inside the new box, with support exactly where gix(k)<rg_i^x(k)<r. Since r≤D/2r\le D/2, (3.7) shows that this support lies in SxS_x. Integrality follows from the integral half-widths and peaks.

Write t=−sx(k)+r/2t=-s_x(k)+r/2. Globally x′≥x+tx'\ge x+t; on its new support, xi′=xi+tx_i'=x_i+t and the parent’s subsolution inequality applies. (4) therefore yields

Fi(x′)≥Fi(x+t)=Fi(x)+t≥xi+t=xi′.F_i(x')\ge F_i(x+t)=F_i(x)+t\ge x_i+t=x_i'.

Thus x′x' is a subsolution. Its witness condition is precisely the stated mass bound. If gix(k)≤r/32g_i^x(k)\le r/32, its untruncated entry is at least ki+r/2−r/32k_i+r/2-r/32, the new Max claim threshold; taking a maximum with the lower boundary cannot destroy that claim.

For Min, relevance likewise gives sy(k)≥D/2−M−Ps_y(k)\ge D/2-M-P. Outside SyS_y, one has yi=Ai+Dy_i=A_i+D, so ki−yi≤−D/2+Mk_i-y_i\le-D/2+M. The gap is again at least 23D/32>D/223D/32>D/2. Now the translated term before clipping is

ki−r/2+giy(k).k_i-r/2+g_i^y(k).

It is above or at the lower boundary, and taking its minimum with the upper boundary leaves support exactly at gaps <r<r. The support is contained in SyS_y by the same outside-support estimate. Set t=sy(k)−r/2t=s_y(k)-r/2. The vector inequality y′≤y+ty'\le y+t and equality yi′=yi+ty_i'=y_i+t on the new support imply

Fi(y′)≤Fi(y+t)=Fi(y)+t≤yi+t=yi′.F_i(y')\le F_i(y+t)=F_i(y)+t\le y_i+t=y_i'.

This is the required supersolution inequality. The mass condition makes y′y' a witness, and a gap at most r/32r/32 puts it at most r/32r/32 above the new lower boundary, giving the asserted claim. Finally, support containment and positivity imply a(Sx′)≤a(Sx)≤1a(S_{x'}) \leq a(S_x) \leq1 and b(Sy′)≤b(Sy)≤1b(S_{y'}) \leq b(S_y) \leq1 for unchanged masses.

The two thresholds serve different purposes: gaps strictly below rr describe the surviving support, whereas gaps at most r/32r/32 give claims.

A deterministic recursion on width and mass

We now give the procedure whose output meets the simultaneous obligations of Theorem 2.42.4. Every step uses exact integers or rationals and label vectors from children. The witnesses from the preceding section are proof objects; no execution searches for or tests them.

For an integer pivot kk and an even width rr, a centered call with masses (a′,b′)(a',b') means Label⁡(k−r/2,r,a′,b′)\operatorname{Label}(k-r/2,r,a',b'), on the box [k−r/2,k+r/2][k-r/2,k+r/2]. This notation will only be used with integer half-widths.

In each pivot pass, one side’s masses remain unchanged. Its translated witnesses force labels at all gaps at most p=q/32p=q/32; those labels will keep these coordinates fixed while the pivot moves. The later width-D/2D/2 call has unchanged masses on both sides and makes all gaps below qq into claims, supplying the labels for reweighting.

Procedure 4.1 (Label⁡(A,D,a,b)\operatorname{Label}(A,D,a,b)). Set u=A+Du=A+D and P=D/32P=D/32. Execute the first applicable base case; if neither applies, proceed through Steps 3–6.

  1. Mass base case. If aibi>1a_i b_i>1 for all i∈Vi\in V, return ++ where bi>1b_i>1 and return −- at the remaining vertices.

  1. Width base case. If D<4096D<4096, compute the box solution zz by Basic. Return ++ at ii if zi>Ai+D/2z_i>A_i+D/2, and −- otherwise.

  1. Downward pass. Use cc, MM, qq, pp from (3.1) and initialize k=ck=c. Repeatedly make a centered width-qq call with masses (a,43b)(a,\frac{4}{3}b) and, from its labels, set simultaneously

kinew={max⁡{ci−M,ki−p},if the label at i is −,ki,if the label at i is +.k_i^{\mathrm{new}} = \begin{cases} \max\{c_i-M,k_i-p\}, & \text{if the label at }i\text{ is }-,\\ k_i, & \text{if the label at }i\text{ is }+. \end{cases}

Replace kk by knewk^{\mathrm{new}}. Stop when the update leaves kk unchanged, denote this last pivot by k↓k^\downarrow, and record

I−={i:ki↓<ci−P}.I_-=\{i:k_i^\downarrow<c_i-P\}.
  1. Upward pass. Start from k=k↓k=k^\downarrow. Repeatedly make a centered width-qq call with masses (43a,b)(\frac{4}{3}a,b) and set simultaneously

kinew={min⁡{ci+M,ki+p},if the label at i is +,ki,if the label at i is −.k_i^{\mathrm{new}} = \begin{cases} \min\{c_i+M,k_i+p\}, & \text{if the label at }i\text{ is }+,\\ k_i, & \text{if the label at }i\text{ is }-. \end{cases}

Replace kk by knewk^{\mathrm{new}}. Stop when it is unchanged, write kfk^f for the final pivot, and record

I+={i:kif>ci+P}.I_+=\{i:k_i^f>c_i+P\}.
  1. Provisional labels and new masses. Make one centered width-D/2D/2 call at kfk^f with masses (a,b)(a,b). For its returned vector λ0∈{+,−}V\lambda^0 \in\{+,-\}^{V}, define

(aˉi,bˉi)={(45ai,85bi),λi0=+,(85ai,45bi),λi0=−.(11)(\bar a_i,\bar b_i)= \begin{cases} \left(\frac{4}{5}a_i,\frac{8}{5}b_i\right), & \lambda_i^0=+,\\ \left(\frac{8}{5}a_i,\frac{4}{5}b_i\right), & \lambda_i^0=-. \end{cases} \tag*{(11)}
  1. Continuation and overrides. Call Label(A,D,aˉ,bˉ)(A,D,\bar a,\bar b). In that child’s output put −- on I−I_- and ++ on I+∖I−I_+\setminus I_-, leaving every other label unchanged. Return the resulting vector.

Each pass includes a final child call whose update changes no coordinate. Correctness will use that call’s guarantee together with the stationary update. The two passes run in the stated order; within one update every coordinate uses the same returned label vector and the same old pivot.

Integral nested boxes and bounded passes

The next facts depend only on the mechanics of the procedure, not on the meaning of the returned labels.

Lemma 4.2. Every recursive child is a valid input of the form Equation (8) , and its box is contained in its parent’s box. In a nonbase call the pivot satisfies

c−M≤k≤cduring the downward pass,c−M≤k≤c+Mduring the upward pass.\begin{aligned} c-M &\le k \le c && \text{during the downward pass},\\ c-M &\le k \le c+M && \text{during the upward pass}. \end{aligned}

If child calls terminate, each pass makes at most 1024n+11024n+1 calls.

Proof. In the recursive case, D≥4096D \ge4096 is dyadic, so all parameters in Equation (3.1), both child half-widths, and every pivot coordinate are integral. The child widths D/128D/128 and D/2D/2 are positive powers of two. The update caps imply the stated pivot bounds. A centered child has width r≤D/2r \le D/2, and its endpoints satisfy

c−M−r/2≥A+D/8≥A,c+M+r/2≤A+7D/8≤A+D.c-M-r/2 \ge A+D/8 \ge A,\qquad c+M+r/2 \le A+7D/8 \le A+D.

The continuation has the unchanged parent box. All mass multipliers preserve positivity and rationality.

The pivot lies on the grid of spacing pp through cc, since M/p=512M/p=512 is integral and both caps belong to that grid. In a changing update at least one coordinate moves by at least pp. Each pass is coordinatewise monotone over a range of at most 2M2M, so a coordinate moves at most 2M/p=10242M/p=1024 times. At most 1024n1024n updates can change the pivot. One further child call certifies the unchanged update at which the pass stops. □

A rank that decreases for every child

For positive masses define

μ(a,b)=min⁡i∈Vaibi,ρ=3225,B(μ)={0,μ>1,1+⌊log⁡ρ(1/μ)⌋,0<μ≤1.\mu(a,b)=\min_{i\in V}a_i b_i,\qquad\rho=\frac{32}{25},\qquad B(\mu)= \begin{cases} 0, & \mu>1,\\ 1+\lfloor\log_{\rho}(1/\mu)\rfloor, & 0<\mu\le1. \end{cases}

The integer BB measures how many guaranteed increases of the minimum mass product can precede the mass base case. It is used for analysis only; the executable procedure does not evaluate logarithms.

Lemma 4.3 (Well-founded recursion). The procedure terminates on every input in Equation (8). Writing D=2dD = 2^d, every child of a nonbase call has smaller rank

R=d+B(μ(a,b)).R = d + B(\mu(a,b)).

The provisional child in Step 5 is the only child that preserves BB. Every other child decreases BB by at least one.

Proof. If μ>1\mu> 1, the mass base case applies. Otherwise multiplying μ\mu by a factor at least ρ\rho decreases BB by at least one: either the new value exceeds one, or its logarithm in Equation (4.2) falls by at least one before rounding down.

In either pivot pass, each product aibia_i b_i is multiplied by 4/3≥ρ4/3 \ge\rho. In the continuation, each product is multiplied by (4/5)(8/5)=ρ(4/5)(8/5) = \rho, independent of the provisional label. Thus these children decrease BB and do not increase dd. The provisional child keeps the masses and halves the width, decreasing dd by one. Every child has strictly smaller nonnegative rank RR.

Induction on RR now proves termination. Base cases require only finite arithmetic and, in the width case, the finite iteration of Lemma 2.1. At a nonbase node, every child terminates by induction. Lemma 4.2 then bounds the two passes; the remaining two children and local operations are finite as well. This reasoning makes no correctness assumption about a label.

The distinction between width reduction and budget reduction is shown in Figure 1. It permits induction for correctness and a sharper count than a general branching bound: only one child can preserve the mass budget, irrespective of the width.

Children of a nonbase call, showing the parent input and three types of child

Figure 1. Children of a nonbase call. The two passes together produce at most 2(1024n+1)2(1024n + 1) narrow children, each decreasing BB. The provisional child preserves BB and halves the width; the continuation preserves the box and decreases BB. Edges show the recursion tree, not the execution order.

Why the labels satisfy every witness

Termination supplies a well-founded order for the correctness proof. At a nonbase call, assume that all children satisfy the simultaneous labelling guarantee. We will show first that the marks I−I_- and I+I_+ can safely receive their prescribed labels. We will then show that every witness with a remaining claim survives the mass change. These two statements allow the continuation to finish the call.

In the next three lemmas, fix a nonbase input (A,D,a,b)(A,D,a,b) and assume Theorem 2.4 for all its children. Witnesses, supports, and claims without a prime refer to this parent input. Peaks and gaps are those of Equations (3.2)–(3.3).

Protection during motion and mass at stopping

Lemma 5.1 (The downward pass). For every relevant Max witness xx, its peak sx(k)s_x(k) is constant throughout Step 3 and at most D/2D/2. In particular, no Max witness claims an index in I−I_-. For every Min witness yy that has a claim outside I−I_-,

b({i:giy(k↓)<q})>3/4.(12)b(\{i : g_i^y(k^\downarrow) < q\}) > 3/4. \tag*{(12)}

Proof. Fix xx and consider one downward update from pivot kk. Its width-qq translation is a witness for the child, since that child keeps Max masses aa. By Lemma [3], every coordinate of gap at most q/32=pq/32 = p is a Max claim of the translated witness. The child’s guarantee labels these coordinates ++, and they do not move. In particular at least one gap-zero coordinate continues to attain the old peak.

At every other coordinate, the downward movement δi\delta_i lies in [0,p][0,p]. For its gap gix(k)>pg_i^x(k) > p,

xi−kinew=sx(k)−gix(k)+δi<sx(k).x_i-k_i^{\mathrm{new}}=s_x(k)-g_i^x(k)+\delta_i<s_x(k).

Thus no coordinate can overtake the old peak, even though many may move in the same update. The peak is unchanged. Initially k=ck=c and x≤A+Dx\le A+D, so sx(c)≤D/2s_x(c)\le D/2. If xx claims ii,

ki↓≥xi−D/2≥Ai+D/2−P=ci−P,k_i^\downarrow\ge x_i-D/2\ge A_i+D/2-P=c_i-P,

which excludes ii from I−I_-.

For the mass assertion, let yy claim j∉I−j\notin I_-. The stopped pivot satisfies kj↓≥cj−Pk_j^\downarrow\ge c_j-P, and the claim gives yj≤Aj+Py_j\le A_j+P. Hence

sy(k↓)≥D/2−2P>D/2−M.(13)s_y(k^\downarrow)\ge D/2-2P>D/2-M. \tag*{(13)}

A coordinate at the lower cap would have ki↓−yi≤D/2−Mk_i^\downarrow-y_i\le D/2-M, since yi≥Aiy_i\ge A_i. Every maximizer of ki↓−yik_i^\downarrow-y_i is therefore strictly above its cap.

Suppose the mass in (12) were at most 3/43/4. The width-qq Min translation has exactly that support, so it would have support mass at most one for the last child’s masses 43b\frac{4}{3}b. Its gap-zero maximizers would all have to receive −-. At least one of those coordinates would then move down by a positive amount, because it is above the cap and p>0p>0. This contradicts the final unchanged update. Equality at 3/43/4 also leads to this contradiction, proving the strict mass bound.

The downward marks are now safe, but the pivot will still move upward. The next lemma both protects the other type of claim and shows why the first mass certificate survives that later movement.

Lemma 5.2 (The upward pass). For every relevant Min witness yy, its peak sy(k)s_y(k) is constant throughout Step 4 and at most D/2D/2. Thus no Min witness claims an index in I+I_+. At the final pivot,

b({i:giy(kf)<q})>34if y has a claim outside I−;(14)b\left(\{i : g_i^y(k^f) < q\}\right) > \frac{3}{4} \qquad\text{if } y \text{ has a claim outside } I_-; \tag*{(14)}
a({i:gix(kf)<q})>34if x has a claim outside I+.(15)a\left(\{i : g_i^x(k^f) < q\}\right) > \frac{3}{4} \qquad\text{if } x \text{ has a claim outside } I_+. \tag*{(15)}

Proof. In every upward child the Min masses are unchanged. Translate a fixed relevant Min witness yy to that width-qq box. All gaps at most pp become Min claims, so those coordinates are labelled −- and stay fixed. Every other coordinate has gap greater than pp and moves upward by at most pp. Its new deviation kinew−yik_i^{\mathrm{new}}-y_i is therefore strictly below the old peak. A gap-zero coordinate does not move, proving exact peak preservation. At the start, k↓≤ck^\downarrow\le c and y≥Ay \ge A imply sy≤D/2s_y \le D/2. A Min claim at ii consequently gives

kif≤yi+D/2≤ci+P,k_i^f \le y_i + D/2 \le c_i + P,

so it is outside I+I_+.

Each coordinate kik_i is nondecreasing during this pass while sys_y is fixed. It follows that

giy(kf)≤giy(k↓).(16)g_i^y(k^f) \le g_i^y(k^\downarrow). \tag*{(16)}

The set of Min gaps <q<q can only grow. Its bb-mass was greater than 3/43/4 at k↓k^\downarrow whenever yy had a claim outside I−I_-, by Lemma 5.1. Positivity of bb proves (14).

It remains to obtain the Max mass certificate at this second stopping point. If xx claims j∉I+j \notin I_+, then xj≥Aj+D−Px_j \ge A_j+D-P and kjf≤cj+Pk_j^f \le c_j+P, giving

sx(kf)≥D/2−2P>D/2−M.s_x(k^f) \ge D/2-2P > D/2-M.

At an upper-capped coordinate, the inequalities kif=ci+Mk_i^f=c_i+M and xi≤Ai+Dx_i \le A_i+D would give a deviation at most D/2−MD/2-M. Thus every maximizer lies strictly below the upper cap. If the aa-mass of the gaps <q<q were at most 3/43/4, the width-qq translation would be a Max witness for the final child’s masses 43a\frac{4}{3}a. Its maximizers would receive ++ and force a positive upward movement, contradicting the unchanged update. This proves (15).

The mass bounds are absolute: each qualifying set has mass greater than 3/43/4, while the entire witness support has mass at most one. They are not merely fractional lower bounds relative to a support that might have much smaller mass. This strength is exactly what the reweighting calculation needs.

Preserving the unsettled witnesses

The provisional child does not directly label the parent’s claims. Its purpose is to assign the favorable label on each large-mass set identified above. The larger width D/2D/2 makes these sets claims of the translated witnesses.

Lemma 5.3 (Witness survival). Every parent Max witness with a claim outside I+I_+ remains a Max witness for (A,D,aˉ,bˉ)(A,D,\bar a,\bar b). Every parent Min witness with a claim outside I−I_- remains a Min witness for that input.

Proof. Write r=D/2r=D/2 for the provisional width. Since

q=D/128<r/32=D/64,q=D/128<r/32=D/64,

the translated width-rr witness claims every coordinate whose parent gap is <q<q.

For a Max witness xx as in the statement, let Nx={i:gix(kf)<q}N_x=\{i:g_i^x(k^f)<q\} and Q+={i:λi0=+}Q_+=\{i:\lambda_i^0=+\}. Lemma 5.2 gives a(Nx)>3/4a(N_x)>3/4. By Lemma 3.1, Nx⊆SxN_x\subseteq S_x; its coordinates are claims of the width-rr translation, which qualifies for the unchanged masses aa. Correctness of that child therefore implies Nx⊆Q+N_x\subseteq Q_+. Reweighting the original support gives

aˉ(Sx)=85a(Sx)−45a(Sx∩Q+)<85⋅1−45⋅34=1.\begin{aligned} \bar a(S_x)=\frac{8}{5}a(S_x)-\frac{4}{5}a(S_x\cap Q_+) \\ &<\frac{8}{5}\cdot1-\frac{4}{5}\cdot\frac{3}{4}=1. \end{aligned}

The continuation has the same box and operator, so the subsolution inequalities are unchanged. The displayed support bound is all that is needed to keep xx a witness.

For Min set Ny={i:giy(kf)<q}N_y = \{i : g_i^y(k^f) < q\} and Q−={i:λi0=−}Q_- = \{i : \lambda_i^0 = -\}. The same two lemmas show Ny⊆Sy∩Q−N_y \subseteq S_y \cap Q_- and b(Ny)>3/4b(N_y) > 3/4. Hence

bˉ(Sy)=85b(Sy)−45b(Sy∩Q−)<85−45⋅34=1.\bar{b}(S_y) = \frac{8}{5}b(S_y) - \frac{4}{5}b(S_y \cap Q_-) < \frac{8}{5} - \frac{4}{5} \cdot\frac{3}{4} = 1.

The supersolution inequalities still hold in the unchanged box, so yy also survives.

Completing the simultaneous induction

Proof of Theorem 2.4. Use induction on R=d+B(μ)R = d + B(\mu) from Lemma 4.3, which has already established termination. We prove that a single returned vector meets every parent obligation.

In the mass base case, aibi>1a_i b_i > 1 at every vertex. If bi>1b_i > 1, no Min witness can contain ii in its support, and hence no Min witness can claim it; returning ++ is permitted. If bi≤1b_i \le1, then ai>1a_i > 1, so no Max witness can claim that vertex; returning −- is permitted. This proves all simultaneous obligations in Step 1.

In the width base case, comparison with the computed box solution zz gives x≤z≤yx \le z \le y for every Max witness xx and Min witness yy. At a Max claim,

zi≥Ai+D−P>Ai+D/2,z_i \ge A_i + D - P > A_i + D/2,

and Step 2 returns ++. At a Min claim, zi≤Ai+P<Ai+D/2z_i \le A_i + P < A_i + D/2, and it returns −-. The argument also applies to D=1D = 1, for which the midpoint is a rational threshold rather than an integer coordinate.

In a nonbase call, all children have smaller rank and satisfy the simultaneous guarantee by induction. The hypotheses of Lemmas 5.1–5.3 therefore hold. Fix any Max witness and any index ii it claims. Lemma 5.1 excludes ii from I−I_-. If i∈I+i \in I_+, the positive override settles it. Otherwise that witness has a claim outside I+I_+ and survives reweighting by Lemma 5.3. The continuation must label ii positively, and no override changes it.

For any Min claim, Lemma 5.2 excludes I+I_+. A claim in I−I_- receives the negative override. A claim outside I−I_- belongs to a surviving Min witness and receives −- from the continuation, with no later change. The marks I−I_- refer to the earlier pivot, but their safety was proved for the fixed parent witnesses and remains valid after the upward pass. The two marked sets may intersect; that intersection contains no claim of either type, so the stipulated negative precedence causes no conflict.

Every vertex receives exactly one label. Since the witnesses and claimed indices in the preceding argument were arbitrary, all their obligations hold together. This completes the induction.

Counting calls and bit operations

Correctness now gives the complete winning set. To finish the theorem, we must bound the work in terms of binary input length rather than the numerical width of the box. There are two relevant depths: the width can be halved polynomially many times in LL, but the mass budget can decrease only logarithmically many times in nn. Because only one child preserves that budget, these depths give a quasipolynomial call count. We then bound the exact arithmetic within each call.

The recursion tree

For the initial data of Equation (2.9), write

d0=log⁡2D0,B0=B(1/n2),H=d0+B0,s=2048n+3.d_0 = \log_2 D_0,\qquad B_0 = B(1/n^2),\qquad H = d_0 + B_0,\qquad s = 2048n + 3.

The minimum initial mass product is 1/n21/n^2. Every path in the call tree has at most HH edges by the rank decrease of Lemma 4.3. At most B0B_0 of its edges can decrease BB, even if an individual edge decreases it by more than one. At any node there is at most one child preserving BB, namely the provisional child. There are at most 2(1024n+1)2(1024n+1) loop children and one continuation, so at most ss children decrease BB.

Lemma 6.1 (Number of calls). The complete main computation makes at most

∑ℓ=0H∑j=0min⁡(ℓ,B0)(ℓj)sj≤(H+1)(B0+1)(1+Hs)B0(17)\sum_{\ell=0}^{H}\sum_{j=0}^{\min(\ell,B_0)}\binom{\ell}{j}s^j \le(H+1)(B_0+1)(1+Hs)^{B_0} \tag*{(17)}

calls, including the initial call.

Proof. Consider a node at depth ℓ\ell whose path uses jj budget-decreasing edges. Their positions can be chosen in (ℓj)\binom{\ell}{j} ways. Order each node’s budget-decreasing children by execution order. At each chosen position there are at most ss choices of that ordinal; at each other position there is at most one. This specifies at most (ℓj)sj\binom{\ell}{j}s^j possible nodes. Sum over all depths 0≤ℓ≤H0 \le\ell\le H and 0≤j≤min⁡(ℓ,B0)0 \le j \le\min(\ell,B_0). For the second inequality, each summand is at most (Hs)j≤(1+Hs)B0(Hs)^j \le(1+Hs)^{B_0} and there are at most (H+1)(B0+1)(H+1)(B_0+1) summands. □\square

This count uses the same recursion principle as the reduced-precision analysis of Parys [25], Section 5, further developed in [20]: one child retains the bounded parameter while all others improve it. Here the two parameters are the logarithm of an integer width and the budget derived from rational mass products.

The explicit encoding bounds both nn and m=∣E∣m = \lvert E\rvert polynomially in L+1L+1, and each original weight has O(L+1)O(L+1) bits. The perturbation and Equation (5) imply

log⁡2(W+1)=O(L+log⁡(n+1)),d0=O(L+log⁡(n+1)).\log_2(W+1) = O(L+\log(n+1)),\qquad d_0 = O(L+\log(n+1)).

Since ρ=32/25\rho= 32/25 is a fixed constant greater than one,

B0=1+⌊log⁡ρ(n2)⌋=O(1+log⁡(n+1)).B_0 = 1+\left\lfloor\log_\rho(n^2)\right\rfloor= O(1+\log(n+1)).

Thus HH and ss are polynomial in L+1L+1, and B0=O(log⁡(L+2))B_0 = O(\log(L+2)). Taking logarithms in the call bound gives

log⁡2(number of calls)≤log⁡2(H+1)+log⁡2(B0+1)+B0log⁡2(1+Hs)=O((log⁡2(L+2))2).\begin{aligned} \log_2(\text{number of calls}) &\le\log_2(H+1)+\log_2(B_0+1)+B_0\log_2(1+Hs)\\ &= O\left((\log_2(L+2))^2\right). \end{aligned}

This includes n=1n=1, for which B0=1B_0=1.

Exact representations

The nested-box property keeps every endpoint and pivot of the main computation inside [0,D0]V[0,D_0]^V. Each such coordinate and each width therefore uses O(d0+1)O(d_0+1) bits. In a basic iteration, an intermediate edge expression is a coordinate plus a transformed weight. It may be negative, but its signed binary length is O(d0+log⁡(W+1)+1)O(d_0+\log(W+1)+1), polynomial in L+1L+1.

Represent a mass by a positive numerator and a positive denominator; reducing fractions is unnecessary. Starting from 1/n1/n, a mass is multiplied along any path only by numbers in

{1,4/3,4/5,8/5}.\{1,4/3,4/5,8/5\}.

Only budget-decreasing edges change the masses, and a path contains at most B0B_0 such edges. Thus each numerator and denominator has O(B0+log⁡(n+1))O(B_0+\log(n+1)) bits. Children receive their own copies of the input masses. In particular, completing one sibling does not multiply the parent masses before another sibling starts. Testing aibi>1a_i b_i>1 uses integer cross-multiplication, and testing bi>1b_i>1 or applying a constant multiplier also uses numbers of polynomial length. All mass arithmetic is exact.

The proof quantities BB and μ\mu, and the peaks, gaps, supports, and witness vectors, need not be stored or computed. The algorithm uses neither a logarithm evaluation for its rank nor an oracle for the existence of a witness. The procedure evaluates FF only inside Basic, by explicit max/min scans of the edge list.

Local cost and uniform implementation

Each call has polynomial nonrecursive work. The mass and width tests are finite scans with the arithmetic just bounded. If the width base case applies, D<4096D<4096, and the imported iteration bound gives

1+nD≤1+4095n1+nD\le1+4095n

evaluations, including the final stability test. Every evaluation scans the mm explicit edges, adds their transformed weights, takes the relevant extrema, and clips at nn endpoints. It acts on integers of polynomial length. The midpoint comparison is implemented as 2zi>2Ai+D2z_i>2A_i+D, which also handles D=1D=1.

In a nonbase call, the constants defining the pivot are obtained by exact divisions by powers of two. Each pass has O(n+1)O(n+1) iterations apart from its children’s work. A scan of the nn coordinates performs one simultaneous update and tests whether it changed anything. Preparing child inputs and receiving their labels, forming I−I_- and I+I_+, reweighting, and applying the final overrides likewise use polynomially many operations on polynomial-length integers. Schoolbook arithmetic suffices throughout.

A depth-first implementation runs children sequentially and retains one record per active call. The record holds the parent masses, endpoints, pivot, labels, and current loop position. There are at most H+1H+1 active records, each of polynomial size. Stack management is therefore included in the same polynomial local bound. Preprocessing the transformed weights, finding the least power of two in (5), and creating masses 1/n1/n take polynomial time. Reading the input and producing the nn final labels also do so.

Multiply this polynomial per-call cost by the call count bounded in (6.3), and include preprocessing and output. For an absolute constant C>0C>0, the complete bound is

2C(log⁡2(L+2))2.2^{C(\log_2(L+2))^2}.

The constants in the procedure are fixed independently of the game, and every arithmetic and control operation has been specified. This gives a uniform deterministic bit algorithm.

Proof of Theorem 1.1. Compute (2.9), run Label, and output the positively labelled vertices. Theorem 2.4 guarantees the labels required by every witness. In the initial box, the box solution is a witness of both types, and the argument after (2.9) identifies all its required labels with the bands of Proposition 2.2. The positive band is exactly the zero-threshold winning set, including against every history-dependent opponent. (6.4) bounds the entire computation. Signed binary weights, self-loops, and separately listed parallel edges are all covered by the construction and accounting above.

Exact values and optimal positional strategies

We now use the winning-set algorithm to recover the numerical values and strategies of the game. For the input model of Theorem 1.1, write

val⁡w(v)=sup⁡σinf⁡τMP⁡w(πv,σ,τ),\operatorname{val}_{w}(v)=\sup_{\sigma}\inf_{\tau}\operatorname{MP}_{w}(\pi_v,\sigma,\tau),

where the strategies may depend on the full finite history.

Corollary 7.1 (Full game solution). A uniform deterministic algorithm returns the exact value val⁡w(v)\operatorname{val}_{w}(v) as a reduced rational for every v∈Vv \in V and one positional strategy σ⋆\sigma^\star for Max and one positional strategy τ⋆\tau^\star for Min such that, simultaneously for every v∈Vv \in V and every pair of history-dependent strategies σ,τ\sigma,\tau,

MP⁡w(πv,σ⋆,τ)≥val⁡w(v),MP⁡w(πv,σ,τ⋆)≤val⁡w(v).\operatorname{MP}_{w}(\pi_v,\sigma^\star,\tau) \ge\operatorname{val}_{w}(v), \qquad\operatorname{MP}_{w}(\pi_v,\sigma,\tau^\star) \le\operatorname{val}_{w}(v).

The complete computation takes 2O((log⁡2(L+2))2)2^{O((\log_2(L+2))^2)} bit operations.

Proof. Positional values. The uniform form of Ehrenfeucht–Mycielski positional determinacy [13], stated explicitly in [7], provides one positional strategy for each player attaining the value from all vertices at once, against arbitrary history-dependent opponents. The latter source uses an edge relation without parallel entries. For our convention, replace each explicit edge e:u→ve:u \to v by u→xe→vu \to x_e \to v through a fresh vertex with a single outgoing edge, with weights 2w(e)2w(e) and 00. Keep ownership at the original vertices; the owner of xex_e is immaterial. The resulting graph has no parallel edges or loops. Expanded histories retain the identity of every original edge. Their averages at even prefixes equal the original averages, and the odd-prefix discrepancy tends to zero. Thus payoffs and the uniform positional guarantees transfer back to our graph, including loops and parallel entries.

Set Win=max⁡e∈E∣w(e)∣W_{\mathrm{in}}=\max_{e\in E}|w(e)|. Fix both uniform optimal positional strategies. From any vertex the resulting play eventually repeats a simple directed cycle of at most nn edges in the original graph, and its mean is the value at that vertex. Consequently every value has a reduced representation a/ba/b with

1≤b≤n,∣a∣≤bWin,−Win≤a/b≤Win.1 \le b \le n,\qquad|a| \le bW_{\mathrm{in}},\qquad-W_{\mathrm{in}} \le a/b \le W_{\mathrm{in}}.

The same bounds hold after restricting outgoing choices, provided at least one outgoing edge remains at every vertex.

Threshold queries and exact reconstruction. We use the classical reduction from values to threshold queries [7], keeping track of binary lengths. For integers pp and q>0q > 0, replace each weight by wp/q(e)=qw(e)−pw_{p/q}(e) = qw(e) - p. For every infinite play,

MP⁡wp/q(π)=qMP⁡w(π)−p.\operatorname{MP}_{w_{p/q}}(\pi) = q\operatorname{MP}_{w}(\pi) - p.

Positional attainment therefore makes the winning set returned by Theorem 1.1 exactly {v:val⁡w(v)≥p/q}\{v : \operatorname{val}_{w}(v) \ge p/q\}.

Choose K=1+⌈log⁡2(2(Win+1)n2)⌉K = 1 + \lceil\log_{2}(2(W_{\mathrm{in}} + 1)n^{2})\rceil, computable from integer bit lengths. For each vertex, start with the closed interval [−Win,Win][-W_{\mathrm{in}}, W_{\mathrm{in}}] and make KK bisections. Query its midpoint and move the lower endpoint to the midpoint if the vertex belongs to the returned set; otherwise move the upper endpoint. The final interval contains its value and has width 2Win2−K<1/n22W_{\mathrm{in}}2^{-K} < 1/n^{2}. Distinct reduced rationals with positive denominators at most nn differ by at least 1/n21/n^{2}. To reconstruct the value from an interval [ℓ,h][\ell, h], for each 1≤b≤n1 \le b \le n test the integer a=⌈bℓ⌉a = \lceil b\ell\rceil for a/b≤ha/b \le h, reduce the retained fractions, and remove repetitions. Since b(h−ℓ)<1b(h-\ell) < 1, this tests every possible numerator. The value bound above guarantees exactly one resulting rational. Call this value-vector procedure Values⁡(H)\operatorname{Values}(H) for any restricted game HH; the original nn, WinW_{\mathrm{in}}, KK remain valid for all such calls.

Separate strategy self-reductions. An edge joining vertices of the same value need not be optimal. For example, let a Max vertex uu have a zero-weight self-loop and a zero-weight edge to a vertex vv whose only edge is a self-loop of weight 1. Both vertices have value 1, but staying at uu yields payoff 0. We therefore test restrictions of the current game and preserve its entire value vector.

First compute v=Values⁡(G)\mathbf{v} = \operatorname{Values}(G). For Max, start from a fresh copy H=GH = G and process the Max vertices in any fixed order. At a vertex uu, form a trial HeH_e for each outgoing edge ee by keeping only ee at uu. Keep a trial for which Values⁡(He)=v\operatorname{Values}(H_e) = \mathbf{v}, and continue from that trial. Such an edge always exists. Indeed, take the edge chosen at uu by a globally optimal positional Max strategy of the current HH. Restricting Max’s choices cannot increase any value, while that strategy remains available and secures the current entire value vector, so no value decreases either. This proves the invariant Values⁡(H)=v\operatorname{Values}(H) = \mathbf{v}.

At the end, Max has one allowed edge at every owned vertex, whereas all of Min’s original choices remain. The resulting single Max strategy therefore secures v\mathbf{v} from every vertex against every Min strategy. To recover Min’s strategy, start again from a fresh copy of the original GG, not from the Max-restricted game. At each Min vertex keep an outgoing edge whose trial preserves v\mathbf{v}. A globally optimal positional Min strategy of the current game supplies such an edge: restricting Min cannot decrease values, and that strategy prevents any increase. The final Min strategy is optimal against every original Max strategy. Thus the two separately recovered strategies have the simultaneous guarantees in the statement.

Call and bit bounds. Put m=∣E∣m = |E|. At a vertex, at most its original outdegree many trials are tested, so the two self-reductions together use at most mm trial value-vector computations. Including the original computation, there are at most

Q=nK(1+m)Q = nK(1+m)

threshold calls. Every midpoint has denominator dividing 2K2^{K} and a numerator of O(log⁡(Win+1)+K)O(\log(W_{\mathrm{in}}+1)+K) bits. Each transformed weight has the same order of bit length. Since K=O(L+log⁡(n+1))K = O(L+\log(n+1)), the full encoding length of every queried game is O(L+mK)O(L+mK), polynomial in L+1L+1. Reconstruction, reduced-rational comparison, and copying restricted edge lists also have polynomial bit cost. The explicit encoding bounds nn, mm polynomially in L+1L+1, so QQ is polynomial. Applying Theorem 1.1 to these polynomial-length games and multiplying by QQ gives 2O((log⁡2(L+2))2)2^{O((\log_2(L+2))^2)} bit operations, including the polynomial-size value and strategy output.

The preceding reduction also applies to the companion algorithm. Here the additional issue is that threshold queries depend on earlier answers, so the winning-set error must be controlled at each query.

Corollary 7.2 (Randomized full game solution). For the game model of Theorem 1.1, a uniform randomized algorithm returns all exact values and both globally optimal positional strategies specified in Corollary 7.1 with probability at least 7/8 for the entire output. Its running time is 2O((log⁡2(L+2))2)2^{O((\log_2(L+2))^2)} bit operations on every random tape.

Proof. The companion’s Theorem 1.1 [24] returns the whole winning set correctly with probability at least 7/8 and has the quasipolynomial bit bound on every tape. Use the reduction above, whose number of threshold queries is bounded by the polynomial QQ in (7.1). For each query, take r=2⌈log⁡2(8Q)⌉+1r=2\lceil\log_2(8Q)\rceil+1 independent fresh runs and use the coordinatewise majority of their winning sets. Conditional on the previous history, the current query is fixed. The amplified set is correct whenever more than half of the runs return the whole correct set, and its conditional failure probability is at most

∑j=⌈r/2⌉r(rj)(1/8)j≤2r(1/8)r/2=2−r/2<18Q.\sum_{j=\lceil r/2\rceil}^{r} \binom{r}{j}(1/8)^j \le2^r(1/8)^{r/2}=2^{-r/2}<\frac{1}{8Q}.

A first-error union bound over the adaptive queries makes every threshold answer correct with probability at least 7/8. On this event the preceding deterministic argument gives all exact values and both globally optimal strategies. On other tapes, if reconstruction does not yield a unique rational or no value-preserving trial is found, the reduction stops with zero at every value coordinate and the first listed outgoing edge at every owned vertex. All loops have the stated caps even then. Since r=O(log⁡(L+2))r=O(\log(L+2)) and all queried encodings retain their polynomial bound on every tape, the total bit bound remains quasipolynomial on every tape, including random-bit generation. This transfer uses only the winning-set interface.

Fixed- and minimum-credit energy games

The preceding algorithms also compute the initial credits needed in one-resource energy games. Retain the game model of Theorem 1.1. For a start v∈Vv\in V and an integer credit c≥0c\ge0, call (v,c)(v,c) energy-winning if Max has a strategy such that, against every history-dependent Min strategy, the induced play satisfies

c+∑t=0j−1w(et)≥0for every j≥0.c+\sum_{t=0}^{j-1}w(e_t)\ge0 \qquad\text{for every }j\ge0.

There is one accumulated resource and no upper bound or saturation. Write cmin⁡(v)c_{\min}(v) for the least nonnegative integer credit for which (v,c)(v,c) is energy-winning, and put cmin⁡(v)=∞c_{\min}(v)=\infty if there is none. A positional Max policy specifies one labelled outgoing edge at each Max vertex.

The reductions connecting these problems to mean-payoff games are classical. Bouyer, Fahrenberg, Larsen, Markey, and Srba formulate the fixed-credit lower-bound problem and reduce it to mean payoff using adversary-controlled returns to the start [6].

Our construction places the initial credit on each return edge. Brim, Chaloupka, Doyen, Gentilini, and Raskin distinguish the unknown-credit and minimum-credit problems, record the unknown-credit/mean-payoff equivalence, and characterize minimum credits by an energy progress measure [7]. The reset construction below lets us apply Corollary 7.1 to a prescribed-credit instance, recovering a policy when that instance is winning. A bound on finite minimum credits then permits exact binary search.

Corollary 8.1 (Fixed and minimum energy credits). Uniform deterministic algorithms have the following guarantees.

  1. Given vv and a nonnegative binary integer cc, decide whether (v,c)(v,c) is energy-winning. On a positive answer, return a positional Max policy satisfying (8.1) from that start and credit against every history-dependent Min strategy.

  1. Given vv, return the exact cmin⁡(v)c_{\min}(v), including the possible answer ∞\infty. On a finite answer, return a positional Max policy winning from vv at credit cmin⁡(v)c_{\min}(v).

The fixed-credit computation takes 2O((log⁡2(Lc+2))2)2^{O((\log_2(L_c+2))^2)} bit operations, where LcL_c is the full explicit input length including the start and credit. The minimum-credit computation takes 2O((log⁡2(L+2))2)2^{O((\log_2(L+2))^2)} bit operations for the original game length LL. Handling every vertex by separate per-start searches, including a separate policy for every finite answer, has the same form of total bit bound. Each policy may depend on the queried start and credit; it also wins from that same start at every larger credit.

Proof. Fixed credit by resetting. For the fixed pair (v,c)(v,c) construct a mean-payoff game Rv,cR_{v,c}. Replace every u∈Vu \in V by an entry uinu_{\mathrm{in}} owned by Min and an exit uoutu_{\mathrm{out}} owned by the original owner of uu. Add the two edges

uin→cvin,uin→0uout.u_{\mathrm{in}} \xrightarrow{c} v_{\mathrm{in}}, \qquad u_{\mathrm{in}} \xrightarrow{0} u_{\mathrm{out}}.

and replace each labelled original edge e:u→ze:u \to z by its own edge

uout→w(e)zin.u_{\mathrm{out}} \xrightarrow{w(e)} z_{\mathrm{in}}.

Start at vinv_{\mathrm{in}}. The game is total and has 2n2n vertices and m+2nm+2n edges, where m=∣E∣m = |E|. The reset at u=vu=v is a nonnegative self-loop, including when c=0c=0, and distinct original parallel edges remain distinct. If U(Rv,c)U(R_{v,c}) denotes its zero-threshold winning set, then

(v,c) is energy-winning in G⟺vin∈U(Rv,c).(v,c)\text{ is energy-winning in }G \Longleftrightarrow v_{\mathrm{in}} \in U(R_{v,c}).

For the forward implication, use an energy-winning strategy from (v,c)(v,c) and restart its history after each reset. In any transformed prefix, every completed segment has original-edge sum at least −c-c and then receives the reset weight +c+c. Its total is therefore nonnegative. The unfinished segment has sum at least −c-c; the intervening zero edges do not change this bound. Every transformed prefix sum is consequently at least −c-c. Dividing by its length shows that the lower limiting average is nonnegative. This argument also covers no resets and infinitely many resets, and does not assume that the original energy strategy is positional.

Conversely, if vin∈U(Rv,c)v_{\mathrm{in}} \in U(R_{v,c}), Corollary 7.1 returns a positional Max strategy for Rv,cR_{v,c} guaranteeing nonnegative mean payoff from that start. Restrict its choices at Max exits to obtain a positional policy in GG. Suppose a consistent original prefix of length jj from vv had weight q<−cq < -c. In the reset game Min can choose continuation at its entries, follow that prefix at its owned exits, and then reset. This is a closed walk from vinv_{\mathrm{in}} with weight q+c<0q+c < 0 and length 2j+12j+1. Positionality of the Max strategy makes the same walk consistent on every repetition. Min could therefore repeat it forever, giving mean payoff (q+c)/(2j+1)<0(q+c)/(2j+1)<0, a contradiction. Thus all original prefixes have weight at least −c-c. This proves the equivalence and the positive-answer policy guarantee; only the Max strategy from the full-solution call is used.

A finite credit bound and exact search. For the original weights put

Win=max⁡e∈E∣w(e)∣,H=(n−1)Win.W_{\mathrm{in}}=\max_{e\in E}|w(e)|,\qquad H=(n-1)W_{\mathrm{in}}.

In particular WinW_{\mathrm{in}} and HH may be zero. An energy winner with any finite credit is an original zero-threshold mean-payoff winner, since its prefix sums are bounded below. Conversely, the cycle-deletion estimate in the proof of Proposition 2.2 gives a positional Max policy on the original winning region for which every consistent prefix has weight at least −H-H. Indeed, its permitted simple cycles have nonnegative original sum, and deleting such cycles from a prefix leaves at most n−1n-1 edges, each of weight at least −Win-W_{\mathrm{in}}. Consequently, for the original winning set U(G)U(G),

cmin⁡(v)<∞⟺v∈U(G),v∈U(G)⟹0≤cmin⁡(v)≤H.c_{\min}(v)<\infty\quad\Longleftrightarrow\quad v\in U(G),\qquad v\in U(G)\quad\Longrightarrow\quad0\le c_{\min}(v)\le H.

The algorithm uses this existence bound as a known winning endpoint; it need not compute the potential used in that proposition.

Compute U(G)U(G) once by Theorem 1.1 and output ∞\infty for starts outside it. For a winning start, the fixed-credit predicate is monotone: a strategy that wins at cc still wins at every larger credit. Binary search with a formal losing lower endpoint ℓ=−1\ell=-1 and the known winning upper endpoint h=Hh=H. While h−ℓ>1h-\ell>1, query ⌊(ℓ+h)/2⌋\lfloor(\ell+h)/2\rfloor by (8.2) and replace the endpoint having the same answer. The sentinel itself is never queried. This takes at most ⌈log⁡2(H+1)⌉\lceil\log_2(H+1)\rceil fixed-credit decisions and returns the exact least winning integer. When H=0H=0 the initial gap is one, so no search query is made and the answer is zero. This includes both n=1n=1 and all-zero weights; any start without a finite credit has already been excluded by U(G)U(G). To return a policy at the minimum, apply Corollary 7.1 to Rv,cmin⁡(v)R_{v,c_{\min}(v)} and restrict its Max choices as above.

Bit complexity and output size. Let bcb_c be the binary length of the supplied nonnegative credit, so bc=O(1+log⁡2(c+1))b_c=O(1+\log_2(c+1)). The explicit reset game has encoding length

O(L+n(bc+log⁡2(n+1))),O\bigl(L+n(b_c+\log_2(n+1))\bigr),

which is polynomial in LcL_c. The credit is copied into only nn edge records; no array indexed by its numerical value is constructed. Theorem 1.1 decides membership in its winning set, and a single application of Corollary 7.1 on a positive answer supplies the policy, within the fixed-credit bound claimed.

For minimum credits, HH and every queried credit have O(1+log⁡2(n+1)+log⁡2(Win+1))=O(L+1)O(1+\log_2(n+1)+\log_2(W_{\mathrm{in}}+1))=O(L+1) bits. All queried encodings therefore have polynomial length in L+1L+1. For all vertices, the outer reduction makes at most 1+n⌈log⁡2(H+1)⌉1+n\lceil\log_2(H+1)\rceil direct winning-set calls, including the original one, and at most one full-solution call for each requested finite-credit policy. These counts are polynomial in L+1L+1. Polynomially many calls on polynomial-length inputs preserve 2O((log⁡2(L+2))2)2^{O((\log_2(L+2))^2)} bit complexity. Computing HH, the binary search arithmetic, and copying the policies have polynomial bit cost. The output also has polynomial length: each finite credit has the stated bit bound, and each per-start policy lists one original edge at each Max vertex. No energy-indexed array is used in the search.

Tropical feasibility and fixed corank

Akian, Gaubert, and Guterman identify the coordinates that can be finite in a solution of tropical inequalities with the winning states of an associated mean-payoff game [1] (Theorem 3.2 and Corollary 3.7). Together with Theorem 1.1, this gives binary-length bounds for tropical feasibility and for rank within a fixed distance of the number of columns. We use the conventions a+(−∞)=−∞a + (-\infty) = -\infty and max⁡∅=−∞\max\varnothing= -\infty.

Corollary 9.1 (Tropical feasibility and fixed corank). For explicitly listed m×nm \times n matrices A,BA,B and mm-dimensional vectors c,dc,d over Q∪{−∞}\mathbb{Q} \cup\{-\infty\}, feasibility of

max⁡{ci,max⁡j(Aij+xj)}≤max⁡{di,max⁡j(Bij+xj)}for every row i\max\{c_i,\max_j(A_{ij}+x_j)\} \le\max\{d_i,\max_j(B_{ij}+x_j)\} \qquad\text{for every row }i

is decidable deterministically in 2O((log⁡(L+2))2)2^{O((\log(L+2))^2)} bit operations, where LL is the complete explicit binary input length. Variables range over R∪{−∞}\mathbb{R} \cup\{-\infty\}; feasibility with all coordinates finite is decidable within the same bound. Equalities between the two displayed maxima are allowed, interpreted as pairs of inequalities. For a homogeneous system (c=d=−∞c=d=-\infty), one can also compute the maximal finite-coordinate support of a solution, including the empty support. In particular, the same bound decides feasibility of max-atom systems z≤max⁡(x,y)+az \le\max(x,y)+a, with rational constants aa and finite real, equivalently rational, variables.

For each fixed integer k≥0k \ge0, a deterministic algorithm with bit bound 2Ok((log⁡(L+2))2)2^{O_k((\log(L+2))^2)} decides whether an explicitly listed m×nm \times n matrix over Q∪{−∞}\mathbb{Q} \cup\{-\infty\}, with m≥n≥1m \ge n \ge1, has tropical rank at least n−kn-k. Here tropical rank is the largest order of a square submatrix whose permutation sums have a finite, unique maximum. The rank is zero if no such nonempty submatrix exists.

Proof. Integer data and homogeneous supports. Let QQ be the product of the positive denominators of the finite input coefficients, taking Q=1Q=1 when there are none. Multiplying every finite coefficient and variable by QQ preserves the inequalities and the sets of finite coordinates. It also preserves the comparisons between permutation sums that define tropical rank. The bit length of QQ is at most one plus the sum of the denominator bit lengths, so the scaled integer data have polynomial binary length.

Consider first a homogeneous system

max⁡j(Cij+zj)≤max⁡j(Dij+zj)for every i.(18)\max_j(C_{ij}+z_j) \le\max_j(D_{ij}+z_j) \qquad\text{for every }i. \tag*{(18)}

The support of zz is {j:zj≠−∞}\{j:z_j \ne-\infty\}. To obtain a total game, apply the preprocessing of [1] (Section 2.1). Discard a row whose left side is identically −∞-\infty. If a right side is identically −∞-\infty, every variable with a finite coefficient on that row’s left side must equal −∞-\infty. Record these coordinates, substitute −∞-\infty for them throughout the system, and repeat. Each step removes a row or a variable, so this takes polynomial time. If no variables remain, the maximal support is empty. Otherwise add the tautologies zj≤zjz_j \le z_j for the remaining variables.

The resulting bipartite game has a Min vertex for each remaining column and a Max vertex for each row. A finite CijC_{ij} gives an edge from column jj to row ii of weight −Cij-C_{ij}; a finite DijD_{ij} gives an edge from row ii to column jj of weight DijD_{ij}. Every vertex has an outgoing edge. Theorem 3.2 of [1] identifies its winning column vertices with the maximal support of a solution. The recorded deleted coordinates belong to no support. The union of all supports is attained, because coordinatewise maxima of finitely many solutions are again solutions; the identically −∞-\infty solution realizes the empty case. The cited game counts a column–row–column pair as one turn. Counting the two edges separately halves the mean payoff; the contribution of an unmatched edge tends to zero. The source uses upper limiting averages; positional optimality makes their nonnegative winning set coincide with that for our lower limiting averages. The game has polynomially many vertices and edges, with weights copied or negated from the scaled coefficients or equal to zero. Its full binary encoding is polynomial in L+1L+1, so Theorem 1.1 computes the maximal support within the asserted bound.

Affine and finite feasibility. After scaling, append cc to AA and dd to BB as a new column, and introduce one coordinate tt. The homogeneous system is

max⁡{ci+t,max⁡j(Aij+zj)}≤max⁡{di+t,max⁡j(Bij+zj)}.\max\{c_i+t,\max_j(A_{ij}+z_j)\}\leq\max\{d_i+t,\max_j(B_{ij}+z_j)\}.

Every affine solution xx gives a homogeneous solution (x,0)(x,0). Conversely, a homogeneous solution with finite tt gives the affine solution xj=zj−tx_j=z_j-t. Thus affine feasibility asks whether the new coordinate belongs to the maximal support. Feasibility with all original coordinates finite asks whether that support contains all n+1n+1 coordinates. A coordinate eliminated during preprocessing is absent from the computed support, so any test requiring it fails; this includes tt if it is forced to −∞-\infty. This is the homogenization underlying [1] (Theorem 3.5 and Corollaries 3.4 and 3.7). Equalities are handled by imposing both inequalities.

A max-atom z≤max⁡(x,y)+az\leq\max(x,y)+a is a homogeneous row with one term on the left and two on the right. Apply the all-finite test. With integer coefficients, taking the floor of each finite coordinate preserves every inequality: flooring commutes with integer translation and with finite maxima. Scaling back therefore gives a rational solution whenever a real one exists. This also explains the integral-solution interface in [1] (Proposition 3.9, Corollary 3.10 and Remark 3.15).

Fixed corank. If n≤kn\leq k, the rank threshold is nonpositive and the answer is yes. Otherwise set r=n−kr=n-k. The rank criterion of [1] (Corollaries 4.16–4.17) says that the rank is at least rr exactly when some set of rr columns is tropically independent. For each such set let MM be the resulting m×rm\times r matrix. Its columns are dependent precisely when there is a vector z≢−∞z\not\equiv-\infty such that every row maximum is attained at least twice or equals −∞-\infty. Equivalently, zz satisfies the homogeneous system

Mij+zj≤max⁡ℓ≠j(Miℓ+zℓ)(1≤i≤m,  1≤j≤r).(19)M_{ij}+z_j\leq\max_{\ell\ne j}(M_{i\ell}+z_\ell)\qquad(1\leq i\leq m,\;1\leq j\leq r). \tag*{(19)}

Its maximal support is therefore empty exactly when the selected columns are independent. This is the inequality form of the game reduction in [1] (Theorem 4.9). The row index (i,j)(i,j) retains the column excluded from Max’s next choice. Empty right sides and all-−∞-\infty columns are covered by the support computation above, including when r=1r=1.

There are at most mrmr inequalities, with rr coefficients on each side, so each test has polynomial encoding length. There are (nr)=(nk)≤nk\binom{n}{r}=\binom{n}{k}\leq n^k tests. For fixed kk, their total bit cost is 2O(k(log⁡(L+2))2)2^{O(k(\log(L+2))^2)}, as claimed.

Strongly polynomial algorithms have been claimed using max-atoms and tropical optimization [27, 28]. Appendix A of the randomized companion gives a two-variable counterexample to the printed variable-selection rule in arXiv:2603.26423v4 [24]. That comparison is specific to the version and rule examined; it neither assesses the distinct 2025 max-atom proof nor transfers to a revised procedure. It has no role in the proofs of this article.

References

  1. [1]Marianne Akian, Stéphane Gaubert, and Alexander Guterman. Tropical polyhedra are equivalent to mean payoff games. International Journal of Algebra and Computation, 22(1):1250001, 2012.
  2. [2]André Arnold, Damian Niwiński, and Paweł Parys. A quasi-polynomial black-box algorithm for fixed point evaluation. In 29th EACSL Annual Conference on Computer Science Logic (CSL 2021), volume 183 of Leibniz International Proceedings in Informatics (LIPIcs), pages 9:1–9:23. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2021.DOI
  3. [3]Peter Austin and Daniele Dell’Erba. Errata to: “Faster Deterministic Exponential Time Algorithm for Energy Games and Mean Payoff Games”. arXiv:2310.04130v1, October 6, 2023.arxiv.org/abs/2310.04130
  4. [4]Henrik Björklund, Sven Sandberg, and Sergei Vorobyov. A combinatorial strongly subexponential strategy improvement algorithm for mean payoff games. In Mathematical Foundations of Computer Science 2004, volume 3153 of Lecture Notes in Computer Science, pages 673–685. Springer, 2004. DOI: 10.1007/978-3-540-28629-5_52.DOI
  5. [5]Henrik Björklund and Sergei Vorobyov. A combinatorial strongly subexponential strategy improvement algorithm for mean payoff games. Discrete Applied Mathematics, 155(2):210–229, 2007. DOI: 10.1016/j.dam.2006.04.029.DOI
  6. [6]Patricia Bouyer, Uli Fahrenberg, Kim G. Larsen, Nicolas Markey, and Jiří Srba. Infinite runs in weighted timed automata with energy constraints. In Franck Cassez and Claude Jard, editors, Formal Modeling and Analysis of Timed Systems (FORMATS 2008), volume 5215 of Lecture Notes in Computer Science, pages 33–47. Springer, 2008. DOI: 10.1007/978-3-540-85778-5_4.DOI
  7. [7]Ľuboš Brim, Jakub Chaloupka, Laurent Doyen, Raffaella Gentilini, and Jean-François Raskin. Faster algorithms for mean-payoff games. Formal Methods in System Design, 38(2):97–118, 2011. DOI: 10.1007/s10703-010-0105-x.DOI
  8. [8]Michaël Cadilhac, Antonio Casares, and Pierre Ohlmann. Fast value iteration: A uniform approach to efficient algorithms for energy games. In Arie Gurfin̈kel and Marijn Heule, editors, Tools and Algorithms for the Construction and Analysis of Systems (TACAS 2025), Part II, volume 15697 of Lecture Notes in Computer Science, pages 323–342. Springer, 2025. DOI: 10.1007/978-3-031-90653-4_16.DOI
  9. [9]Cristian S. Calude, Sanjay Jain, Bakhadyr Khoussainov, Wei Li, and Frank Stephan. Deciding parity games in quasi-polynomial time. SIAM Journal on Computing, 51(2):STOC17–152–STOC17–188, 2022. DOI: 10.1137/17M1145288.
  10. [10]Laure Daviaud, Marcin Jurdziński, and Ranko Lazić. A pseudo-quasi-polynomial algorithm for mean-payoff parity games. In Proceedings of the 33rd Annual ACM/IEEE Symposium on Logic in Computer Science, LICS ’18, pages 325–334. ACM, 2018. DOI: 10.1145/3209108.3209162.DOI
  11. [11]Dani Dorfman, Haim Kaplan, and Uri Zwick. A faster deterministic exponential time algorithm for energy games and mean payoff games. In 46th International Colloquium on Automata, Languages, and Programming (ICALP 2019), volume 132 of Leibniz International Proceedings in Informatics (LIPIcs), pages 114:1–114:14. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2019. An updated author version contains the weight-independent bound in Theorem 3.10: https://danidorfman.com/publication/energy-games/energy-games.pdf.
  12. [12]Dani Dorfman, Haim Kaplan, and Uri Zwick. Improved bounds for strategy improvement algorithms for energy games. In Philip Bille, Seth Pettie, and Sabine Storandt, editors, 34th Annual European Symposium on Algorithms (ESA 2026), volume 388 of Leibniz International Proceedings in Informatics (LIPIcs), pages 140:1–140:18. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2026. DOI: 10.4230/LIPIcs.ESA.2026.140.DOI
  13. [13]Andrzej Ehrenfeucht and Jan Mycielski. Positional strategies for mean payoff games. International Journal of Game Theory, 8(2):109–113, 1979. DOI: 10.1007/BF01768705.DOI
  14. [14]Nathanaël Fijalkow, Pawel Gawrychowski, and Pierre Ohlmann. Value iteration using universal graphs and the complexity of mean payoff games. In 45th International Symposium on Mathematical Foundations of Computer Science (MFCS 2020), volume 170 of Leibniz International Proceedings in Informatics (LIPIcs), pages 34:1–34:15. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2020. DOI: 10.4230/LIPIcs.MFCS.2020.34.arxiv.org/abs/1812.07072
  15. [15]Vladimir A. Gurvich, Alexander V. Karzanov, and Leonid G. Khachiyan. Cyclic games and an algorithm to find minimax cycle means in directed graphs. USSR Computational Mathematics and Mathematical Physics, 28(5):85–91, 1988. DOI: 10.1016/0041-5553(88)90012-2.DOI
  16. [16]Marcin Jurdziński. Deciding the winner in parity games is in UP ∩ co-UP. Information Processing Letters, 68(3):119–124, 1998. DOI: 10.1016/S0020-0190(98)00150-1.DOI
  17. [17]Marcin Jurdziński, Rémi Morvan, Pierre Ohlmann, and K. S. Thejaswini. A symmetric attractor-decomposition lifting algorithm for parity games. arXiv:2010.08288v1, 2020.arxiv.org/abs/2010.08288
  18. [18]Alexander V. Karzanov and Vasilij N. Lebedev. Cyclical games with prohibitions. Mathematical Programming, 60:277–293, 1993. DOI: 10.1007/BF01580616.DOI
  19. [19]Alexander Kozachinskiy. Polyhedral value iteration for discounted games and energy games. In Dániel Marx, editor, Proceedings of the 2021 ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 600–616. SIAM, 2021. DOI: 10.1137/1.9781611976465.37.DOI
  20. [20]Karoliina Lehtinen, Pawel Parys, Sven Schewe, and Dominik Wojtczak. A recursive approach to solving parity games in quasipolynomial time. Logical Methods in Computer Science, 18(1):8:1–8:18, 2022. DOI: 10.46298/LMCS-18(1:8)2022.DOI
  21. [21]Yury M. Lifshits and Dmitri S. Pavlov. Potential theory for mean payoff games. Journal of Mathematical Sciences, 145(3):4967–4974, 2007. DOI: 10.1007/s10958-007-0331-y.DOI
  22. [22]Bruno Loff and Mateusz Skomra. Smoothed analysis of deterministic discounted and mean-payoff games. In Karl Bringmann, Martin Grohe, Gabriele Puppis, and Ola Svensson, editors, 51st International Colloquium on Automata, Languages, and Programming (ICALP 2024), volume 297 of Leibniz International Proceedings in Informatics (LIPIcs), pages 147:1–147:16. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2024. DOI: 10.4230/LIPIcs.ICALP.2024.147.DOI
  23. [23]Pierre Ohlmann. A symmetric recursive algorithm for mean-payoff games. arXiv:2603.07555v1, 8 March 2026. DOI: 10.48550/arXiv.2603.07555.DOI
  24. [24]OpenAI. Randomized quasipolynomial-time mean-payoff games. OpenAI Math Release preprint OAI:Randomized-quasipolynomial-time-mean-payoff-games-September-25-2026, 2026.
  25. [25]Pawel Parys. Parity games: Zielonka’s algorithm in quasi-polynomial time. In Peter Rossmanith, Pinar Heggernes, and Joost-Pieter Katoen, editors, 44th International Symposium on Mathematical Foundations of Computer Science (MFCS 2019), volume 138 of Leibniz International Proceedings in Informatics (LIPIcs), pages 10:1–10:13. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2019. DOI: 10.4230/LIPIcs.MFCS.2019.10.DOI
  26. [26]K. S. Thejaswini, Pierre Ohlmann, and Marcin Jurdziński. A technique to speed up symmetric attractor-based algorithms for parity games. In Anuj Dawar and Venkatesan Guruswami, editors, 42nd IARCS Annual Conference on Foundations of Software Technology and Theoretical Computer Science (FSTTCS 2022), volume 250 of Leibniz International Proceedings in Informatics (LIPIcs), pages 44:1–44:20. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2022. DOI: 10.4230/LIPIcs.FSTTCS.2022.44.
  27. [27]Laurent Truffet. Looking for all solutions of a set of max-atoms solves the max atom problem in strongly polynomial time. Modélisation des Systèmes Réactifs (MSR 2025), 2025. HAL: hal-05491586.
  28. [28]Laurent Truffet. Substitution for minimizing/maximizing a tropical linear (fractional) programming. arXiv:2603.26423v4, August 29, 2026. arXiv:2603.26423v4.arxiv.org/abs/2603.26423
  29. [29]Uri Zwick. Improved subexponential analysis of the Random-Action-Removal algorithm for 2-player turn-based games and non-binary AUSOs. arXiv:2607.06334v1, 2026. DOI: 10.48550/arXiv.2607.06334.DOI
  30. [30]Uri Zwick and Mike Paterson. The complexity of mean payoff games on graphs. Theoretical Computer Science, 158(1–2):343–359, 1996. DOI: 10.1016/0304-3975(95)00188-3.DOI

Paper details

Contents