Abstract 1 Introduction 2 Preliminaries 3 Defining Problem Parameters 4 Problem Parameters Are Computable Efficiently 5 Problem Parameters Determine Distributed Complexity 6 Putting Things Together References

Classification of Local Optimization Problems in Directed Cycles

Thomas Boudier ORCID Gran Sasso Science Institute, L’Aquila, Italy    Fabian Kuhn ORCID University of Freiburg, Germany    Augusto Modanese ORCID CISPA Helmholtz Center for Information Security, Saarbrücken, Germany    Ronja Stimpert ORCID Aalto University, Espoo, Finland    Jukka Suomela ORCID Aalto University, Espoo, Finland
Abstract

We present a complete classification of the distributed computational complexity of local optimization problems in directed cycles for both the deterministic and the randomized LOCAL model. We show that for any local optimization problem Π (that can be of the form min-sum, max-sum, min-max, or max-min, for any local cost or utility function over some finite alphabet), and for any constant approximation ratio α, the task of finding an α-approximation of Π in directed cycles has one of the following complexities:

  1. 1.

    O(1) rounds in deterministic LOCAL, O(1) rounds in randomized LOCAL,

  2. 2.

    Θ(logn) rounds in deterministic LOCAL, O(1) rounds in randomized LOCAL,

  3. 3.

    Θ(logn) rounds in deterministic LOCAL, Θ(logn) rounds in randomized LOCAL,

  4. 4.

    Θ(n) rounds in deterministic LOCAL, Θ(n) rounds in randomized LOCAL.

Moreover, for any given Π and α, we can determine the complexity class automatically, with an efficient (centralized, sequential) meta-algorithm, and we can also efficiently synthesize an asymptotically optimal distributed algorithm.

Before this work, similar results were only known for local search problems (e.g., locally checkable labeling problems). The family of local optimization problems is a strict generalization of local search problems, and it contains numerous commonly studied distributed tasks, such as the problems of finding approximations of the maximum independent set, minimum vertex cover, minimum dominating set, and minimum vertex coloring.

Keywords and phrases:
LOCAL model, optimization, cycles, meta-algorithms
Category:
Track A: Algorithms, Complexity and Games
Copyright and License:
[Uncaptioned image] © Thomas Boudier, Fabian Kuhn, Augusto Modanese, Ronja Stimpert, and
Jukka Suomela; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Distributed algorithms
Related Version:
Full Version: https://arxiv.org/abs/2602.13046 [6]
Acknowledgements:
This project was initiated at the Research Workshop on Distributed Algorithms (RW-DIST 2025) in Freiburg, Germany; we would like to thank all workshop participants and organizers for inspiring discussions, and the anonymous reviewers for their helpful feedback on this paper. Most of this work was done while A.M. was affiliated with Aalto University.
Funding:
This work was supported in part by the Research Council of Finland, Grant 359104, and by the Quantum Doctoral Education Pilot, the Ministry of Education and Culture, decision VN/3137/2024-OKM-4.
Editors:
Sayan Bhattacharya, Danupon Nanongkai, Michael Benedikt, and Gabriele Puppis

1 Introduction

In this work we study distributed algorithms for local optimization problems – these include a wide range of familiar graph optimization problems such as maximum independent set, minimum vertex cover, minimum dominating set, and minimum vertex coloring. We show that in directed cycles, essentially all questions about the existence of distributed algorithms for local optimization problems can be answered in a mechanical and efficient manner. We give a meta-algorithm that answers the following questions, for any given local optimization problem Π and any given constant approximation factor α:

  1. 1.

    How fast can we find an α-approximation of Π in the deterministic 𝖫𝖮𝖢𝖠𝖫 model?

  2. 2.

    How fast can we find an α-approximation of Π in the randomized 𝖫𝖮𝖢𝖠𝖫 model?

1.1 Broader Context

There is a long line of prior work that has studied local search problems, and in particular locally checkable labeling problems or LCLs; these are problems where the task is merely to satisfy all local constraints, i.e., the task is to find a feasible solution that does not need to be optimal in any sense. Familiar examples include the problems of finding a maximal independent set, or a valid coloring with Δ+1 colors in a graph of maximum degree Δ.

While a lot of early work focused on specific concrete examples of such tasks, e.g., exactly what is the distributed complexity of finding a 3-coloring of a cycle [10, 20], Naor and Stockmeyer [22] initiated the study of meta-computational questions in this context: for any given local search problem Π, can we systematically determine what the distributed computational complexity of Π is?

It turns out that such questions are in general undecidable, already in the case of seemingly simple graph classes such as 2-dimensional grids [22, 7]. However, in graph families such as cycles, paths, and trees many such questions are decidable [9, 3, 4, 2]. Furthermore, while such questions are computationally hard if we have input labels [2], many of these questions can be solved efficiently (in time that is polynomial in the size of the description of problem Π) in the case of unlabeled graphs [9, 3].

In particular, if we look at LCL problems on unlabeled directed cycles, they are very well understood: the round complexity of any such problem is O(1), Θ(logn), or Θ(n) in all of our usual models of distributed computing (deterministic 𝖫𝖮𝖢𝖠𝖫, randomized 𝖫𝖮𝖢𝖠𝖫, deterministic 𝖢𝖮𝖭𝖦𝖤𝖲𝖳, and randomized 𝖢𝖮𝖭𝖦𝖤𝖲𝖳), and there is a computer program that automatically finds out the complexity class for any given problem [24].

However, many problems that have been studied in the field of distributed graph algorithms are local optimization problems, e.g., [18, 11, 14, 13, 19, 17, 16], and much less is currently known about them. There is a large collection of isolated examples that apply to specific problems, but there is no complete classification of possible distributed round complexities of local optimization problems, and there is no meta-algorithm that is able to automatically determine the complexity of a given problem.

Furthermore, what makes the question particularly intriguing is that it is known that the complexity landscape has to be fundamentally different, already in the simplest case of unlabeled directed cycles. To see this, consider the problem of finding a 1.1-approximation of a minimum dominating set in cycles. Thanks to [11], it is known that in the deterministic 𝖫𝖮𝖢𝖠𝖫 model, the complexity of this problem is Θ(logn) rounds, while in the randomized 𝖫𝖮𝖢𝖠𝖫 model, the complexity is O(1) rounds. The same never happens with any LCL problem: for them, O(1) rounds in randomized 𝖫𝖮𝖢𝖠𝖫 always implies O(1) rounds in deterministic 𝖫𝖮𝖢𝖠𝖫 [8, 22].

Hence we have three questions that we aim to answer in this work:

  1. 1.

    What is the right framework for formalizing and representing local optimization problems (analogous to LCL problems in the context of search problems)?

  2. 2.

    What does the complexity landscape look like for local optimization, especially when we compare deterministic 𝖫𝖮𝖢𝖠𝖫 with randomized 𝖫𝖮𝖢𝖠𝖫?

  3. 3.

    Can we design a meta-algorithm that automatically determines the complexity class for any given local optimization problem and any given approximation ratio?

1.2 Formalism

In this work, we look at the simplest nontrivial case: local optimization problems in directed cycles. Let us first recall how we can specify LCL problems in this setting: There is some finite set of output labels Γ and verification distance r, and we specify which (r+1)-tuples of output labels are valid. For example, for the task of 3-coloring, we have Γ={1,2,3}, and we can select r=1, and the set of valid tuples of output labels is {12,13,21,23,31,32}. A feasible solution is a labeling s:VΓ that is everywhere locally valid: for each node v, the sequence of output labels s(Nr(v)) formed by v and its r successors is a valid tuple.

 Remark.

In cycles, it does not really matter whether we look at v and its r successors, or v and its r predecessors, or v and its radius-r/2 symmetric neighborhood. In this work we follow the convention that a problem is defined in terms of r successors.

In local optimization problems we keep the same basic setup, but now we associate a value c(x) with each valid (r+1)-tuple x. This way, each labeling s associates a local value c(s(Nr(v))) with each node (we will use the terms cost or utility as appropriate to refer to the value, depending on the problem). The optimization task is one of the following, combining an objective function obj{max,min} and an aggregation function aggr{max,min,}:

  1. 1.

    maxmin: maximize the minimum of local values,

  2. 2.

    max: maximize the sum of local values,

  3. 3.

    minmax: minimize the maximum of local values,

  4. 4.

    min: minimize the sum of local values.

We will use the value to indicate forbidden label combinations (which can also be interpreted as e.g. infinite cost). We are interested in approximations, and we follow the convention that approximation ratio α is always at least 1, so e.g. a 2-approximation for a minimization problem is a feasible solution with total cost at most twice the optimum, and a 2-approximation for a maximization problem is a feasible solution with total utility at least half the optimum.

1.3 Examples of Problems

The following examples will hopefully clarify the formalism, and also demonstrate that it is broadly applicable. We will start with familiar, classic graph optimization problems.

Example 1.1 (maximum independent set).

The maximum independent set problem can be encoded as follows, with 1 indicating nodes that form the independent set:

  • Γ={0,1}, r=1, obj=max, aggr=,

  • c(00)=0, c(01)=0, c(10)=1, c(11)=.

Here any valid independent set has a nonnegative total value, and the value is equal to the size of the independent set.

Example 1.2 (minimum dominating set).

The minimum dominating set problem can be encoded as follows:

  • Γ={0,1}, r=2, obj=min, aggr=,

  • c(000)=, c(001)=0, c(010)=0, c(011)=0,
    c(100)=1, c(101)=1, c(110)=1, c(111)=1.

Example 1.3 (minimum vertex coloring).

To encode the minimum vertex coloring problem, we can use e.g. the following formulation:

  • Γ={1,2,3}, r=1, obj=min, aggr=max,

  • c(ij)=i if ij, c(ii)=.

Here for example a valid 2-coloring (using only colors 1 and 2) has total cost 2, and a valid 3-coloring has total cost 3.

Example 1.4 (maximum domatic partition).

Let us now look at a slightly more interesting example: domatic partition. In this problem, the task is to partition the set of nodes into as many disjoint dominating sets as possible; in a cycle we can have at most 3 such sets. While more clever ad-hoc encodings are possible, we will use the following approach that generalizes to a wide range of similar problems. We use labels of the form xi, where i indicates that we claim to have a solution with exactly i dominating sets, and x indicates that this node claims to be in a dominating set with label x. We then need to ensure that nodes also agree on the number of dominating sets, and each dominating set is also valid. This will suffice:

  • Γ={a1,a2,b2,a3,b3,c3}, r=2, obj=max, aggr=min,

  • c(a1a1a1)=1,

  • c(a2a2b2)=c(a2b2a2)=c(a2b2b2)=c(b2a2a2)=c(b2a2b2)=c(b2b2a2)=2,

  • c(a3b3c3)=c(b3c3a3)=c(c3a3b3)=c(a3c3b3)=c(c3b3a3)=c(b3a3c3)=3,

  • c(x)= otherwise.

1.4 A Classification Example

So far we have seen examples of natural max, min, maxmin, and minmax problems. We will introduce a somewhat artificial problem, which is engineered to illustrate many key phenomena in a single, concise package.

Example 1.5 (sloppy coloring).

In this problem, we can produce a proper 2-coloring with colors {,}, a proper 3-coloring with colors {1,2,3}, or a somewhat sloppy 3-coloring with colors {a,b,c}. The problem is defined as follows:

  • Γ={,,1,2,3,a,b,c}, r=1, obj=min, aggr=,

  • c()=c()=1,

  • c(12)=c(21)=c(23)=c(32)=c(13)=c(31)=2,

  • c(ab)=c(ba)=c(bc)=c(cb)=c(ac)=c(ca)=3,

  • c(aa)=100,

  • c(x)= otherwise.

Let us try to gain a bit more intuition on the problem first. If we manage to produce a proper 2-coloring, the total cost in an n-cycle is n, and if we manage to produce a proper 3-coloring using labels {1,2,3}, the total cost is 2n. If we resort to labels {a,b,c} and produce a proper 3-coloring, the total cost is 3n. If a sublinear number of edges uses the color pair aa, the total cost will be 3n+o(n), and if we produce a trivial solution with all nodes colored a, the total cost will be 100n.

Now the key question is this: given some constant α1, what is the distributed round complexity of finding an α-approximation of a sloppy coloring? It turns out that the complete answer in this case looks like this:

  • For 1α<2, the complexity is Θ(n) in both the deterministic 𝖫𝖮𝖢𝖠𝖫 and randomized 𝖫𝖮𝖢𝖠𝖫 models. The optimal algorithm is a trivial brute-force algorithm that finds a 2-coloring if it exists, and otherwise produces a 3-coloring. There is no o(n)-round algorithm for finding a 1.999-approximation or better.

  • For 2α3, the complexity is Θ(logn) in both the deterministic 𝖫𝖮𝖢𝖠𝖫 and randomized 𝖫𝖮𝖢𝖠𝖫 models. We can always produce a 3-coloring, using the Cole–Vishkin algorithm [10], and this is always a 2-approximation. There is no o(logn)-round deterministic or randomized algorithm for finding a 3-approximation or better.

  • For 3<α<100, the complexity is Θ(logn) in the deterministic 𝖫𝖮𝖢𝖠𝖫 model but O(1) in the randomized 𝖫𝖮𝖢𝖠𝖫 model. The randomized algorithm essentially finds a random ruling set to divide the cycle in segments. The algorithm uses color a at the segment boundaries and fills in all sufficiently short segments with b,c,b,c,. Super-constant-length segments are rare, and there we can use the more costly solution a,a,a,a,. By adjusting the parameters, we can achieve approximation ratio 3+ϵ for any constant ϵ>0 with high probability (the worst case being that the input is an even cycle that could be 2-colored). There is no o(logn)-round deterministic algorithm for finding a 99.999-approximation or better.

  • For α100, the complexity is O(1) in both the deterministic 𝖫𝖮𝖢𝖠𝖫 and randomized 𝖫𝖮𝖢𝖠𝖫 models. The algorithm simply outputs a everywhere.

As we will see, this turns out to be a well-representative problem: for any local optimization problem in our problem family, there are at most these four complexity classes, and the optimal algorithms are also essentially generalizations of the above four algorithms. We can automatically deduce all relevant α thresholds between the classes.

1.5 Contributions

We show that for any local optimization problem Π (in the formalism of Section 1.2), and for any constant α, the task of finding an α-approximation of Π in directed cycles has one of these complexities:

  1. 1.

    O(1) rounds in deterministic 𝖫𝖮𝖢𝖠𝖫, O(1) rounds in randomized 𝖫𝖮𝖢𝖠𝖫,

  2. 2.

    Θ(logn) rounds in deterministic 𝖫𝖮𝖢𝖠𝖫, O(1) rounds in randomized 𝖫𝖮𝖢𝖠𝖫,

  3. 3.

    Θ(logn) rounds in deterministic 𝖫𝖮𝖢𝖠𝖫, Θ(logn) rounds in randomized 𝖫𝖮𝖢𝖠𝖫,

  4. 4.

    Θ(n) rounds in deterministic 𝖫𝖮𝖢𝖠𝖫, Θ(n) rounds in randomized 𝖫𝖮𝖢𝖠𝖫.

Moreover, for any given Π and α, we can determine the complexity class automatically, with an efficient (centralized, sequential) meta-algorithm, and we can also efficiently synthesize an asymptotically optimal distributed algorithm. In particular, all local optimization problems similar to those in Sections 1.3 and 1.4 are now fully understood in the case of directed cycles.

One key consequence of our work is that increasing the running time from Θ(logn) to, say, Θ(logn) or Θ(n) does not give any constant advantage in the approximation ratio, for any local optimization problem in directed cycles. If we are interested in finding a (1+1/n)-approximation, for example, then the situation is more complex, and indeed we might still obtain some sub-constant advantage by doing so. However, as we show, we will never be able to improve our approximation ratio from, say, 1.2345 to 1.2344 by such means. You will need to jump from Θ(logn) rounds all the way to Θ(n) rounds for such an improvement.

We emphasize that the case of max-min and min-max problems can be handled by known techniques, as each fixed α essentially defines a sub-problem that is an LCL. The main novel part is min and max problems, which require genuinely new techniques and ideas. However, we wanted to give a unified formalization of all of these problems, and show that the same definitions and concepts indeed largely apply to both problem families.

1.6 Techniques

We take a three-step approach:

Step 1: Defining the problem parameters (Section 3).

We first define seven parameters βopt, βflex, δflex, βcoprime, βgap, δgap, and βconst that capture essentially everything that we need to know about a given optimization problem Π.

To do this, we define a node-weighted de Bruijn graph GdB that captures the entire structure of Π, including the values. The intuition is that for any input cycle G, a feasible solution of Π in G corresponds to a closed walk in GdB, with the same value. Then we proceed to define four subgraphs of GdB:

  1. 1.

    We let Gopt=GdB. We define βopt based on the cheapest closed walk in Gopt.

  2. 2.

    We zoom into flexible components GflexGopt, in the spirit of [9]. In brief, a node v of GdB is K-flexible if there are self-returning walks from v back to v of all possible lengths K,K+1,K+2,. A node is flexible if it is K-flexible for some K. It follows that in each strongly connected component either all or none of the nodes are flexible, and we say that a component is flexible if the nodes are flexible. We define parameters βflex and δflex based on the cheapest closed walks in Gflex.

  3. 3.

    We zoom into flexible components with self-loops GgapGflex; we simply take those flexible components that also contain nodes with a self-loop. We define parameters βgap and δgap based on the cheapest closed walks in Ggap.

  4. 4.

    We form the subgraph that consists only of self-loops GconstGgap. We define parameter βconst based on the cheapest closed walks (that is, cheapest self-loops) in Gconst.

The above six parameters suffice for min and max problems. For max-min and min-max problems we also need to introduce a parameter βcoprime that, in brief, captures the cheapest pair of closed walks that have coprime lengths.

Table 1: All β values defined in Section 3 for the examples in Section 1.3. Missing values are not essential to capture the problem’s classification, regardless of whether the value can be obtained or not.
Example βopt βflex δflex βcoprime βgap δgap βconst
Example 1.1: max. independent set 12 12 true 12 true 0
Example 1.2: min. dominating set 13 13 true 13 true 1
Example 1.3: min. vertex coloring 2 false 3 false
Example 1.4: max. domatic partition 3 false 2 1 false 1
Example 1.5: sloppy coloring 1 2 false 3 true 100

Table 1 shows concrete examples of these parameter values for the examples from Section 1.3. We emphasize that these parameters have a purely graph-theoretic definition; their definitions do not make any references to e.g. a particular model of distributed computing. However, in step 3 we will see how these parameters connect with the task at hand.

Step 2: The problem parameters are efficiently computable (Section 4).

We show that given a description of any local optimization problem Π, we can efficiently compute all necessary parameters out of βopt, βflex, δflex, βcoprime, βgap, δgap, and βconst. In each case, we first construct the de Bruijn graph GdB, which has size polynomial in the size of the problem description. Then we show that we can efficiently construct the relevant subgraph (Gopt,Gflex,Ggap,Gconst). Finally, we show that we can design an efficient graph algorithm for determining the corresponding parameter values based on the structure of the subgraph.

For example, to determine the value of βopt for a min problem, it is sufficient to find a cycle in Gopt=GdB that minimizes the average cost; we show that all of this can indeed be done efficiently, largely relying on standard graph algorithms.

Determining e.g. βflex and δflex is more interesting, and there we will build on and extend prior work related to flexibility [9]. In brief, it is no longer sufficient to find a single cheap cycle, but we also need to find a pair of low-average-cost cycles with coprime lengths. We refer to Section 4 for the details of how all seven parameters can be determined efficiently.

Step 3: The problem parameters determine distributed complexity (Section 5).

Finally, we show that the complexity of any given problem for any approximation ratio α can be determined using parameters βopt, βflex, δflex, βcoprime, βgap, δgap, and βconst. Furthermore, we can also automatically construct an asymptotically-fastest algorithm for any given α.

The intuition here is that βopt provides a baseline against which we can compare. For example, for min problems, for infinitely many values of n, there exists a solution of total cost βoptn, and for no value of n we have cheaper solutions. So an average cost of βopt per node represents the best-case scenario for the optimal solution. In particular, if we can find a solution of cost at most αβoptn for all n, it will clearly be an α-approximation.

For a full classification, we will now need both upper bounds and matching lower bounds. Let us first consider upper bounds. It is perhaps easier to consider what we can do with a given time budget, using min problems as an example:

  1. 1.

    If we can use O(n) rounds, we can brute-force an optimal solution, and hence trivially find a 1-approximation, regardless of the input size.

  2. 2.

    If we can use only O(logn) rounds, we can use the following idea: find a sufficiently-sparse ruling set that divides the cycle into segments, so that the length of each segment is at least A and at most B for some appropriately-chosen constants A and B. Then we resort to the notion of flexibility: we put some fixed “flexible output sequence” x at segment boundaries. Then we can follow a cheap closed walk from x to x in the de Bruijn graph Gflex to fill the gap between two such markers. There are two cases:

    1. (a)

      There are two closed walks from x to x in Gflex that have coprime lengths and that have average cost exactly βflex. Then by using a combination of such walks we can fill in gaps so that the total cost of the solution is exactly βflexn, and hence we can achieve an approximation ratio of α=βflex/βopt.

    2. (b)

      Otherwise we will occasionally use slightly more expensive fragments, but we can still achieve a total cost of βflexn+o(n), and hence approximation ratio α=βflex/βopt+ϵ for any constant ϵ>0.

  3. 3.

    If we further restrict the running time to O(1) rounds, but allow the use of randomness, we can use a faster but sloppier approach for dividing the cycle into segments by choosing a random ruling set. The key difference is that now we can no longer guarantee that all segments have length at most some constant B, yet we need to be able to solve the problem in constant time. This is where we resort to self-loops in Ggap: in those (rare) segments that are too long, we simply fill in the middle part of the segment by following a self-loop, that is, we produce in those areas a constant solution. For example, in the maximum independent set problem this means that we may occasionally have long stretches of 0s in our output. By choosing the parameters appropriately, we can achieve a total cost of βgapn+o(n), and hence an approximation ratio of α=βgap/βopt+ϵ.

  4. 4.

    Finally, if we must use an O(1)-round deterministic algorithm, we will produce a constant output that corresponds to following the cheapest self-loop of Gconst repeatedly, achieving a total cost of βconstn and hence approximation ratio of α=βconst/βopt.

Then we need to show that the above strategy is optimal: there are no other complexity classes, and for each class the above approximation ratios are the best possible. To give some spirit of the lower bounds, we for example reason as follows. Assume that we have an algorithm that runs in o(n) rounds. Then if the algorithm succeeds in producing a valid solution for every instance, it cannot use output values outside Gflex. Hence the output of the algorithm forms a closed walk in Gflex. But such walks have an average cost at least βflex, while the optimum might have an average cost of βopt, and it follows that the algorithm cannot find an approximation better than α=βflex/βopt.

The lower bounds get more technical in the O(logn) region. There we use Naor’s randomized version [21] of Linial’s lower bound [20] as well as the Ramsey-theoretic lower bound technique familiar from Czygrinow, Hanckowiak, and Wawrzyniak [11]. We refer to Section 5 for the details.

1.7 Full Classification

Tables 2 and 3 summarize how we can determine the exact complexity class for any given problem Π and any given approximation ratio α. The classification is complete; there are matching upper and lower bounds for each case. Table 4 shows what is obtained from applying our method to the previous examples from Section 1.3.

Table 2: Scheme for inferring the complexity class from the problem parameters of -LCLs. To determine the optimal complexity for a problem and a given α, one scans the “sufficient condition” column until a line is matched. Notice that the bounds given are tight in the sense that, if a condition on the right is satisfied and none of those above it apply, then the complexity on the left is also a lower bound for the corresponding approximation problem. For example, if we have a min problem where αβoptβgap and δgap is false but αβopt<βconst holds, then there is an Ω(logn) lower bound for deterministic 𝖫𝖮𝖢𝖠𝖫.
Complexity Sufficient condition
det. 𝖫𝖮𝖢𝖠𝖫 rand. 𝖫𝖮𝖢𝖠𝖫 min max
O(1) O(1) αβoptβconst 1αβoptβconst
Θ(logn) O(1) αβoptβgap and δgap is false 1αβoptβgap and δgap is false
αβopt>βgap and δgap is true 1αβopt<βgap and δgap is true
Θ(logn) Θ(logn) αβoptβflex and δflex is false 1αβoptβflex and δflex is false
αβopt>βflex and δflex is true 1αβopt<βflex and δflex is true
Θ(n) Θ(n) αβopt<βflex 1αβopt>βflex
unsolvable unsolvable βflex= βflex=
Table 3: Scheme for deriving the complexity class from the problem parameters of min-max and max-min-LCLs. The same remarks regarding determining the complexity of a problem as in Table 2 apply. Notice that, unlike the case of -LCLs, there is not the possibility of a deterministic Θ(logn) and randomized O(1) complexity class here.
Complexity Sufficient condition
det. 𝖫𝖮𝖢𝖠𝖫 rand. 𝖫𝖮𝖢𝖠𝖫 minmax maxmin
O(1) O(1) αβoptβconst 1αβoptβconst
Θ(logn) Θ(logn) αβoptβcoprime 1αβoptβcoprime
Θ(n) Θ(n) αβopt<βcoprime 1αβopt>βcoprime
unsolvable unsolvable βcoprime= βcoprime=
Table 4: Complexity of computing an α-approximation in the deterministic and randomized 𝖫𝖮𝖢𝖠𝖫 models for each of our examples from Section 1.3.
Complexity
Example Range of α Strategy det. 𝖫𝖮𝖢𝖠𝖫 rand. 𝖫𝖮𝖢𝖠𝖫
Example 1.1: 1 Optimal Θ(n) Θ(n)
max independent set (1,) Constant fragment Θ(logn) O(1)
Example 1.2: 1 Optimal Θ(n) Θ(n)
min. dominating set (1,3) Constant fragment Θ(logn) O(1)
[3,) Constant solution O(1) O(1)
Example 1.3: [1,32) Optimal Θ(n) Θ(n)
min. vertex coloring [32,) Flexible Θ(logn) Θ(logn)
Example 1.4: [1,32) Optimal Θ(n) Θ(n)
max. domatic partition [32,3) Flexible Θ(logn) Θ(logn)
[3,) Constant solution O(1) O(1)
Example 1.5: [1,2) Optimal Θ(n) Θ(n)
sloppy coloring [2,3] Flexible Θ(logn) Θ(logn)
(3,100) Constant fragment Θ(logn) O(1)
[100,) Constant solution O(1) O(1)

1.8 Future Directions and Open Questions

Our results are easy to adapt to weaker models: the results generalize essentially verbatim also to the 𝖢𝖮𝖭𝖦𝖤𝖲𝖳 model [23], and they are easy to adapt to the port-numbering model, where we essentially lose the deterministic Θ(logn) class. However, generalizing the results to stronger classes leads to interesting open questions. One obvious candidate is the quantum-𝖫𝖮𝖢𝖠𝖫 model; however, whether symmetry-breaking problems require Ω(logn) rounds in quantum-𝖫𝖮𝖢𝖠𝖫 is a major open question [1], and until this is resolved, we cannot fully characterize local optimization problems in the quantum-𝖫𝖮𝖢𝖠𝖫 model, either.

An orthogonal direction is generalizations of our questions to more diverse families of inputs. The first steps would be to generalize our results from directed cycles to undirected cycles, directed paths, and rooted trees. Already the case of undirected cycles leads to new technical challenges: constant-round deterministic algorithms for optimization problems in undirected cycles may achieve non-trivial approximations (consider e.g. the task of orienting the edges, and assigning a higher cost to sink nodes). We believe that by introducing notions analogous to mirror-flexibility from [9] to our setting, it is also possible to eventually classify all these cases.

We have not considered inputs, and as soon as we have inputs, the classification will be 𝖯𝖲𝖯𝖠𝖢𝖤-hard [2]. Whether it is nevertheless decidable for local optimization problems with inputs is yet another open question for future work.

2 Preliminaries

We denote the set of positive integers by + and +{0} by 0+. Similarly, 0+ is the set of non-negative reals. We also write [n]={1,,n} for the set of the first n elements in +.

Our input graph is always a directed cycle G=(V,E). Without restriction, let us name the nodes of a directed cycle of length n such that V={v0,v1,,vn1} and (vi,v(i+1)modn)E for i{0,,n1}. The (one-way) neighborhood of radius r0+ of a node is a function Nr:VVr+1 such that

Nr(vi)=(vi,v(i+1)modn,,v(i+r)modn).

Hence the neighborhood of a node extends itself along a single direction of the cycle, namely following the orientation of the edges. We will use the one-way neighborhood to define the local costs incurred by every node; our computational model, however, allows communication in either direction (i.e., disregarding the orientation of the edges).

2.1 Locally Checkable Optimization Problems

A locally checkable optimization problem (opt-LCL) Π on directed cycles is a tuple Π=(Γ,r,c,aggr,obj) where:

  • Γ is a finite alphabet,

  • r0+ is the (constant) radius,

  • c:Γr+10+{} assigns every neighborhood a certain cost,

  • aggr{,min,max} is an aggregation function,

  • obj{min,max} is an objective function.

The label is used to indicate neighborhoods that are not valid. In the following, such neighborhoods and their cost can be immediately disregarded. In order to distinguish different types of opt-LCLs we use shorthand such as obj-aggr-LCL (e.g., min--LCL for an opt-LCL where obj=min and aggr=), and we write -LCL for any opt-LCL where aggr=.

A solution s on a directed cycle G=(V,E) is a function s:VΓ. For the neighborhood Nr(vi) of a node vi, we define

s(Nr(vi))=(s(vi),s(v(i+1)modn),,s(v(i+r)modn)).

We let Sn be the solution space containing all solutions for a problem on a cycle with n nodes. A solution s to an opt-LCL Π=(Γ,r,c,aggr,obj) on a directed cycle G=(V,E) is said to be valid if and only if c(s(Nr(v))) holds for every vV.

The optimal cost of a problem Π=(Γ,r,c,aggr,obj) on a directed cycle G=(V,E) of length n is a function

OPT(n)=objsSnaggrvVc(s(Nr(v))).

Note that OPT depends only on n, since a directed cycle is determined by its length. A solution s is said to be optimal for some n if and only if

aggrvVc(s(Nr(v)))=OPT(n).

For an opt-LCL Π with obj=min, a solution of value v in an n-cycle is an α-approximation if αOPT(n)v. For obj=max, the condition is that 1αOPT(n)v. An algorithm is said to find an α-approximation for Π if it finds an α-approximation for any n.

2.2 The 𝗟𝗢𝗖𝗔𝗟 Model

We assume the usual distributed 𝖫𝖮𝖢𝖠𝖫 model of computation due to Linial [20]. The input graph G is perceived as defining a distributed network where nodes are computers and edges are communication links. The orientation of an edge is interpreted as input information and does not constrain communication in any way (i.e., it is allowed to send messages in the opposite direction of an edge). Each node is assigned a unique identifier from the set [nc], where c>1 is some constant of our choosing. Computation proceeds in rounds, where in each round a computer exchanges messages of unbounded size with its neighbors and executes a procedure where it decides whether to commit to an output or not and what to send on the next round of communication. We disregard time and space complexity of local computation and care only about the number of communication rounds required to solve the problem at hand. In this paper, our upper bounds give computable local procedures at every node whereas for the lower bounds we may as well assume an arbitrary procedure.

Whenever we address the randomized 𝖫𝖮𝖢𝖠𝖫 model, we say that an algorithm produces an α-approximation if, with high probability, the output is valid for the problem and additionally approximates the optimal solution within a factor of α. That is, with probability at most 1n, the result is a labeling that is either invalid or does not satisfy the allowed approximation ratio.

3 Defining Problem Parameters

In this section, we define the seven parameters βopt, βflex, βcoprime, δflex, βgap, δgap, and βconst that fully capture the complexity of producing approximations of an opt-LCL Π. Before we can do this in Definitions 3.5, 3.6, and 3.7, we first introduce some technical tools which we also require later in Section 4.

The main technical tool that we employ to analyze an opt-LCL Π is the de Bruijn graph associated with it. The values of the problem parameters for Π all come from characteristics of certain walks in the de Bruijn graph associated with it.

Definition 3.1 (De Bruijn Graph).

Given an opt-LCL Π=(Γ,r,c,aggr,obj), its de Bruijn graph is the digraph GdB=(VdB,AdB) where:

  • The nodes in VdB are all (r+1)-tuples (s1,,sr+1)Γr+1 fulfilling the condition c((s1,,sr+1)).

  • We have an edge from (s1,,sr+1) to (t1,,tr+1) if ti=si+1 holds for every i[r].

Moreover we equip this graph with the cost function c:VdBR(c), where R(c) is the range of c. The function c attributes to a particular neighborhood NΓr+1 its cost c(N). We denote by γ|Γ|r+1 the number of nodes in GdB.

Figures 1(a), 1(b), 1(c), 1(d), and 2 show the de Bruijn graphs associated with the problems from Section 1.3.

(a) The de Bruijn graph of Example 1.1.
(b) The de Bruijn graph of Example 1.2.
(c) The de Bruijn graph of Example 1.3. On the left-hand side is the general structure of the graph. On the right-hand side is the subgraph relevant for calculating the values in Table 1.
(d) The de Bruijn graph of Example 1.4.
Figure 1: De Bruijn graphs for Examples 1.1, 1.2, 1.3, and 1.4, with the cost of each neighborhood in blue. Neighborhoods with cost are not included in the De Bruijn graphs even if they occur in the respective problem.

Next we look at walks in the de Bruijn graph associated with Π as well as the total cost incurred by following such walks.

Definition 3.2 ((u,v)-walk).

A sequence of nodes W=w0,,wk with i,wiVdB is a (u,v)-walk if w0=u, v=wk and i,(wi,wi+1) is an arc of GdB. When u=v, we refer to W as a closed walk, or, in particular, as a closed v-walk.

Definition 3.3 (Cost of a walk).

Let W=w0wk be a walk with k1. If aggr=, then c(W)=1ki=1kc(wi). Otherwise, c(W)=aggri[1,,k]c(wi).

Note that we ignore the value of c(w0) systematically to avoid counting some element twice when W is a closed walk.

Finally, we define flexible nodes and components, which are relevant for flexible strategies.

Definition 3.4 (Flexible nodes and components).

A node v of GdB is flexible if there exists K such that kK, there exists a closed v-walk of length exactly k. We say in particular that v is K-flexible.

A strongly connected component SVdB of GdB is a flexible component if there is a node vS such that v is flexible.

We are now ready to define the seven parameters βopt, βflex, βcoprime, δflex, βgap, δgap, and βconst that capture the key properties of an opt-LCL Π. Recall OPT(n) denotes the value of an optimal solution for a cycle of n nodes. Basically, we relate each of the βa values to the optimal solution in a certain subclass of the solution space. These subclasses are implicitly defined by taking appropriate subgraphs of the de Bruijn graph associated with Π. Meanwhile, δflex and δgap are Boolean values that depend on certain specific properties of Π.

Definition 3.5 (Gopt, Gflex, Ggap, Gconst).

Let GdB be the de Bruijn graph associated with an opt-LCL. We define a set of (possibly empty) subgraphs 𝒢={Gopt,Gflex,Ggap,Gconst} for Π as follows:

  1. 1.

    Gopt=GdB.

  2. 2.

    GflexGopt containing only the flexible components of Gopt, excluding all edges between different components.

  3. 3.

    GgapGflex containing only components of Gflex that contain a self-loop.

  4. 4.

    GconstGgap containing only nodes of Ggap that have a self-loop and their self-loops.

Definition 3.6 (βopt, βflex, βgap, βconst, δflex, δgap).

Let Π=(Γ,r,c,aggr,obj) be an opt-LCL. For each Ga𝒢, we define 𝒞(Ga) as the set of cycles (including self-loops) in Ga and the real-valued parameter

βa=objC𝒞(Ga)c(C).

If 𝒞(Ga) is empty, then βa=. The parameters δflex and δgap are defined as follows:

  1. 1.

    δflex is false if there are two closed walks W1,W2 in Gflex of coprime lengths, sharing a node, so that c(W1)=c(W2)=βflex; otherwise, we set δflex to true.

  2. 2.

    δgap is false if βgap=βconst, and true otherwise.

As an example, for min--LCLs we can connect βa with concrete solutions to Π as follows:

  1. 1.

    βopt is the best possible relative value that can be achieved infinitely often (i.e., asymptotically) by the optimum over all cycle lengths;

  2. 2.

    βflex is the (again, asymptotically) worst such value;

  3. 3.

    βgap is the best possible value for solutions that admit constant neighborhoods, that is, where a single, fixed label appears;

  4. 4.

    βconst is the best possible (relative) value for constant solutions, that is, where the same fixed label is used everywhere.

Figure 2 shows how to obtain the relevant β and δ values for Example 1.5. Note that for -LCLs, if k is the length of the cycle providing βopt, then nmodk=0 implies OPT(k)/k=OPT(n)/n=βopt. In addition, for large n, we have βoptnOPT(n)βflexn+O(1). It is possible that δflex is false and βflex=βopt.

Definition 3.7 (βcoprime).

To define βcoprime, we let V(r)={vVdBc(v)r} and V(r)={vVdBc(v)r}. Let G(r) and G(r) be the subgraphs of the de Bruijn graph GdB induced by V(r) and V(r), respectively. We then define βcoprime for min-max-LCLs and max-min-LCLs as follows:

  • For min-max-LCLs, βcoprime is the smallest r0+ for which the graph G(r) contains two closed walks of coprime lengths that have a node in common.

  • For max-min-LCLs, βcoprime is the largest r0+ for which the graph G(r) contains two closed walks of coprime lengths that have a node in common.

Figure 2: An example for obtaining β and δ values from a de Bruijn graph GdB when considering Example 1.5, with the cost of each neighborhood in blue and the relevant subgraphs of GdB in green. After identifying these subgraphs one can obtain the β values by finding the cycles with the cheapest average cost in each subgraph. The edges of a selection of such cheapest cycles are shown in red. To illustrate why δflex is false, two cycles of the same cost are highlighted in GflexGgap; these cycles have coprime lengths. Since βgapβconst, δgap is true.

4 Problem Parameters Are Computable Efficiently

In this section, denote by GdB the de Bruijn graph of the given opt-LCL. Recall that γ is its number of nodes. Next we present the necessary result to prove that all seven problem parameters can be computed efficiently. Note that βgap, δflex and δgap are only relevant to -LCLs, so we do not need to compute them otherwise.

Let 𝒲opt be the set of closed walks that are optimal regarding the objective function. In the event that 𝒲opt=, it is impossible to form a closed walk in GdB, and hence it is impossible to solve Π. Hence we may assume 𝒲opt.

Lemma 4.1.

If 𝒲opt, there exists a closed walk in 𝒲opt of length at most γ.

Proof.

Let w𝒲opt be a closed v-walk of length at least γ+1. We prove that we can find a walk w𝒲opt such that |w|<|w|. We repeat this argument on w until we eventually find a walk of length at most γ.

As |w|γ+1, by the pigeonhole principle, either v is visited three times or there exists some vv visited twice. In the first case, we split w into two closed v-walks at v, L1 and L2. Observe that both L1 and L2 have lengths in [1,|w|1], and both have their cost equal to c(w), as otherwise w𝒲opt: If c(L1)<c(w) (resp., c(L1)>c(w) in the case of maximization), w is no longer optimal. Moreover, if c(L1)>c(w) (resp., c(L1)<c(w) for maximization), then for aggr= the other subwalk L2 must have better cost than w, while for aggr{min,max} this contradicts the definition of c(w) as the aggregate over w. Hence, both of them belong to 𝒲opt.

If there is some vv visited twice, the proof is quite similar. Split w into three walks at v: L1=v,,v, L2=v,,v and L3=v,,v. Now consider the two closed walks L1+L3 and L2. Their lengths are in [1,|w|1] and they both must have the same cost as w; hence they belong to 𝒲opt.

The definition of 𝒲opt suggests that the optimal elements of the de Bruijn graph are closed walks and not only simple cycles. In the next lemma we prove that we can actually think only about cycles:

Lemma 4.2.

If there is a closed walk W of length at most γ with cost β in the de Bruijn graph then there is a cycle of length at most γ with cost at least as good as β in terms of the optimization function.

Proof.

If W is a cycle, the proof is trivial. Otherwise, choose two occurrences of the same node in W that split W into two shorter closed walks W1 and W2. One of W1 and W2 has cost at least as good as W with respect to the objective function: for aggr= this follows from the fact that c(W) is the weighted average of c(W1) and c(W2), while for aggr{min,max} it follows directly from the aggregate over the concatenation. Continue with such a shorter closed walk until the remaining closed walk is a cycle.

In order to compute βflex and βcoprime it is crucial to identify flexible nodes of the graph. To figure out if a given node is flexible we prove the following lemma:

Lemma 4.3.

Node v is flexible if and only if v belongs to two closed walks of coprime lengths.

Proof.

Let A and B be two closed v-walks with a=|A| and b=|B| being coprime. By the Frobenius coin theorem, all integers at least (a1)(b1) can be obtained by a nonnegative integer combination of a and b. Concatenating the corresponding copies of A and B therefore gives closed v-walks of all sufficiently large lengths, hence v is flexible.

Now let v be a flexible node. Hence there exists K>0 such that for all kK, there exists a closed v-walk of length exactly k. Pick any aK and bK such that a and b are coprime. The corresponding closed v-walks have coprime lengths.

Figuring out the value of βflex is in principle easy: for each connected component, we can compute the cost of the optimum cycle in that component, and test which components are flexible. Testing the flexibility of a component can be achieved by considering walks of length O(γ) [9]. We then set βflex accordingly.

The hard part of our work is to decide whether δflex is false or not. Parameter δflex would be false only if there are two closed walks of coprime lengths, sharing a vertex and of cost βflex. We adapt the procedure for testing flexibility from [9] to take into account the cost function of the de Bruijn graph in the next two lemmas.

Definition 4.4.

If Π is an obj--LCL, then Gβflex=(Vβflex,Aβflex), where

Vβflex :={vVdBv belongs to a closed walk of cost βflex},
Aβflex :={(u,v)AdB(u,v) belongs to a closed walk of cost βflex}.

It would be reasonable to think that having very long walks might be good in order to optimize the cost of a solution. We prove in our next lemma that this is not the case: all the walks of Gβflex have cost exactly βflex, and hence the length of the walk does not matter. This lemma, coupled with the fact that testing flexibility can be done by considering walks of limited length will allow us to efficiently test whether δflex is false or not. As δflex is a parameter defined only for obj--LCLs, there is no analogous definition of Gβflex for min-max-LCLs and max-min-LCLs.

Lemma 4.5.

Let Π be an obj--LCL. Any closed walk in the graph Gβflex has cost exactly βflex.

Proof.

Define the shifted cost function cβflex:Vβflex by cβflex(v)=c(v)βflex. Call the cost of a walk W computed using cβflex instead of c the shifted cost of W. It is clear that a walk of cost βflex has shifted cost zero (as kβflex=i=1kc(wi) if and only if 0=i=1k(c(wi)βflex) for obj--LCLs). Moreover, by definition of βflex, all closed walks must have shifted cost at least 0 if obj=min (resp. at most 0 if obj=max). For simplicity, we can imagine that Gβflex is constructed by adding one by one a finite number of simple closed walks W1,W2, of shifted cost zero. (Note that we do not actually need to design an algorithm to build this sequence.)

We now prove the lemma by induction. Let W be a closed walk of Gβflex. Assume that W exists solely on the graph spanned by W1. It is clear that W1 being a simple closed walk, W consists of one or multiple instances of W1 and hence cβflex(W)=0.

Suppose now that for some k, any closed walk on the graph spanned by ikWi has shifted cost exactly zero. We want to prove that any closed walk W on ik+1Wi has shifted cost exactly 0. We use the following notation to represent W:

a11,,ak11=b11,,b11=a12,,ak22,,

where akii=b1i, bii=a1i+1 and bij are nodes of Wk+1, and aij are other nodes. The reader can refer to Figure 3 as a visual aid for the rest of the proof.

Consider a walk W defined in the following way: it is built following W but each time we reach a frontier node akii=b1i, we add a first segment going from akii to a1i+1 through nodes in j=1kWj then a second segment going from a1i+1=bii to b1i through nodes of Wk+1. Then, W resumes its visit of W. Note that it is possible to reorient the way in which W is travelled so that it can be seen as a visit of ik+1Wi, then when it reaches a frontier node b1i, it may go through Wk+1 entirely, then resume its visit of ikWi. Then it is clear that cβflex(W)=0 as it consists of multiple travels through Wk+1 and a walk on ikWi. Also, we have WW. If we had cβflex(W)>0 (resp., cβflex(W)<0), then in some way W would be compensating in some part of WW. However, WW consists of a collection of cycles of the form biib1i plus some segments going from akii to a1i+1 in ikWi. If cβflex(W)>0 (resp., cβflex(W)<0), then one of those cycles would be negative (resp., positive) shifted cost, which is a contradiction with the definition of βflex. Hence it follows that all closed walks on ik+1Wi have cost exactly βflex.

Figure 3: Illustration of the construction of W using W in Lemma 4.5. Each time W uses some part of Wk+1, we add at the end of it the unused part of Wk+1 (the green upper part) and any path going from akii to a1i+1 in ikWi. This way it is clear that W consists of possibly multiple instances of Wk+1 and some closed walk on ikWi, making the shifted cost of W necessarily 0.

For testing if a digraph has flexible nodes it is sufficient to look at closed walks of lengths at most 2γ+1, as shown in [9]. We give a version of this result tailored to our needs in the full version of this paper [6]:

Lemma 4.6.

Let H=(U,F) be a strongly connected digraph of n vertices. If H has flexible vertices, then there exists a node vU and two closed walks A,B containing v such that |A| and |B| are coprime and both have length 2n+1.

We can directly use this result to prove that in G(βcoprime) (resp. G(βcoprime)) there are two closed walks of coprime lengths that share a vertex, have length at most 2γ+1, and have cost at most (resp. at least) βcoprime as all walks of G(βcoprime) (resp. G(βcoprime)) must have cost at most (resp. at least) βcoprime. Thanks to Lemma 4.5, we know that all the walks of Gβflex have cost exactly βflex, so we can also apply the same result. Then we can test if there are two closed walks of coprime lengths sharing a vertex and having cost βflex by only looking at walks of length at most 2γ+1; hence deducing efficiently whether δflex is false or not.

Lemma 4.7.

Let GdB be the de Bruijn graph of a -LCL. If there are two closed walks of coprime lengths sharing a vertex and having cost βflex in GdB, then there are two such walks of length at most 2γ+1.

Proof.

Apply Lemma 4.6 to any strongly connected component of Gβflex containing such a shared vertex, where Gβflex is defined in Lemma 4.5, and note that those walks exist in Gβflex if and only if they exist in GdB.

Lemma 4.8.

For any opt-LCL, if they exist, we can compute the values of βopt, βflex, δflex, βcoprime, βgap, δgap, and βconst in time polynomial in |Γ|.

Proof.

First observe that for any opt-LCL we have r=O(1) and hence any algorithm running in time polynomial in |Γ| runs in time polynomial in γ=|Γ|r+1.

Let us start by addressing the more complex case of -LCLs. We will use the following subroutine: For sVdB, let A(s,K) be the following process: if it exists, denote by F(s,v,k) an optimal (s,v)-walk of length exactly k. Initialize as non-existing F(s,v,0) for vs and F(s,s,0)={s}. For all vVdB and k0 we can compute F(s,v,k+1) using the formula

argobjuVdBN(v)(c(F(s,u,k)+v)),

where N(v) represents the predecessors of v. Repeat this process until k=K. Then output the list of tuples (k,F(s,s,k)) for each value of k for which F(s,s,k) exists. Algorithm A(s,K) computes an optimal closed walk containing s of length k, for each kK for which such a walk exists. It can return multiple elements (at most K tuples).

When K is polynomial in γ, algorithm A(s,K) is clearly also polynomial in γ. We can then use it to obtain all the desired values:

  • βopt: Run A(s,γ) for all s and then figure for which s,k we have an optimal c(F(s,s,k)). Correctness is ensured by Lemma 4.1. Naturally, this procedure will also detect cycles among the closed walks. Lemma 4.2 ensures that the value we find will belong to at least one cycle, meeting the definition of βopt. The procedure takes polynomial time in γ as we call A(s,γ) only γ times.

  • βflex and δflex: First, determine the set of strongly connected components of GdB. For each component, flexibility can be determined by observing the output of A(s,2γ+1) for all s: if s is flexible, then it must belong to two closed walks of coprime lengths (Lemma 4.3), and there must be two such walks of lengths at most 2γ+1 (Lemma 4.6). Now, for nodes in flexible components, we can determine the value of βflex by running A(s,γ) for each s in a flexible component and picking the component that has the best cycle. As for δflex, it suffices to look at the output of A(s,2γ+1) in the good component and check if there are two closed walks containing s, of coprime lengths and cost βflex. In short, all we have to do is to run A(s,2γ+1) once for all nodes and all the information we need is in that output. Calling A(s,2γ+1) for all s is clearly polynomial in γ and all further analysis of the output is also polynomial.

  • βgap and δgap: Identify flexible components as explained above, but limit the considered components further by requiring them to contain one node with a self-loop. This procedure is clearly polynomial in γ and obtains βgap. Let Sloop be the set of all nodes with a self-loop. Then, δgap is false if and only if Sloop contains a node with a self-loop of cost βgap.

  • βconst: Run A(s,1) for every s. If A(s,1) is non-empty, take the walk of optimal value (depending on obj).

As for min-max-LCLs and max-min-LCLs, the procedure is much easier. Start by listing in L all costs c(s) for all nodes and sort L in ascending order if obj=min (resp. descending order if obj=max). Then, for all L, we can construct the graphs G() and G() in time polynomial in γ. Now, do the following: For all L in ascending order (resp. descending order), do the same procedure as for -LCLs to determine βopt. If the procedure fails (meaning we could not find a closed walk of length at most γ), continue through L. The first time the procedure succeeds, we have found βopt=. When βopt is found we can do a similar procedure to find βcoprime: build G() (resp. G()) and run A(s,2γ+1) for all s. As for obj--LCLs, this fully determines which components are flexible. The first time a flexible component is found, we have βcoprime=. We repeat the same idea to find βconst. Regarding the complexity, we can notice that all the procedures that were polynomial in γ for -LCLs are called at most |L|γ times, hence the whole algorithm runs in time polynomial in γ.

5 Problem Parameters Determine Distributed Complexity

In this final section, we show that the complexity of any given opt-LCL for any approximation ratio α can be determined using the parameters from Section 3. We do this by showing a series of lower and upper bounds in both deterministic and randomized 𝖫𝖮𝖢𝖠𝖫 based on these parameters. The result of all these observations is that every opt-LCL falls into one of five different classes based on its optimal complexity. An overview of these classes, which we label from A down to E, is provided in Table 5.

Throughout the previous sections we had to take into account that minimization and maximization problems behave in distinct ways, for example when referring to costs. However, the results we describe in this section can be readily ported from one case to the other by simple methods such as flipping signs or inverting fractions in the relevant inequalities. For example, an inequality αβopt>βflex in the minimization case becomes α/βopt>1/βflex in the maximization one (since we are interested in 1/α as a factor instead of α). Then we can simply recast the quantities 1/βopt and 1/βflex, respectively, as, say, βopt and βflex and proceed with the same argument. For this reason, in the text below we handle only the case of minimization problems and skip any mention of maximization problems altogether.

Table 5: Any α-approximation of an opt-LCL Π falls into one of these five complexity classes, depending on the structure of Π and the value of α.
Complexity
Class deterministic 𝖫𝖮𝖢𝖠𝖫 randomized 𝖫𝖮𝖢𝖠𝖫
A O(1) O(1)
B Θ(logn) O(1)
C Θ(logn) Θ(logn)
D Θ(n) Θ(n)
E unsolvable unsolvable

5.1 Lower Bounds

First we address lower bound results that show that Table 5 is complete, that is, that there is no opt-LCL whose optimal complexity lies in-between any of the classes A–E. For the full proofs of these results, please refer to the full version of this paper [6].

Lemma 5.1.

Assume there is a randomized 𝖫𝖮𝖢𝖠𝖫 algorithm that finds an α-approximation of some -LCL Π in o(n) rounds. Then αβoptβflex holds, where the inequality is strict if and only if δflex is true.

Since randomized 𝖫𝖮𝖢𝖠𝖫 is stronger than deterministic 𝖫𝖮𝖢𝖠𝖫, this directly implies a lower bound of Ω(n) for deterministic 𝖫𝖮𝖢𝖠𝖫 as well.

Lemma 5.2.

If there is a randomized 𝖫𝖮𝖢𝖠𝖫 algorithm that finds an α-approximation of some min-max-LCL Π in o(n) rounds, then αβoptβcoprime.

Lemma 5.3.

If an opt-LCL Π can be α-approximated in o(logn) rounds in the randomized 𝖫𝖮𝖢𝖠𝖫 model, then αβoptβgap. Additionally, if δgap is true, then this inequality is strict.

The proof follows more or less directly from well-known results on symmetry breaking and, in particular, the classical Ω(logn) lower bound for 3-coloring [21]. Note that this lower bound holds even when we are allowed to err, for instance, with probability 1/Ω(logn) (instead of 1/n as usual).

Lemma 5.4.

If an opt-LCL Π can be α-approximated in o(logn) rounds in the deterministic 𝖫𝖮𝖢𝖠𝖫 model, then αβoptβconst.

Our proof follows a similar idea to that of [11] based on an application of Ramsey’s theorem. For a set X and k+, let us write (Xk) for the subsets of X of size exactly k. For k,r,s+ where ks2, we let R(k;r,s) denote the size of the smallest n+ such that every r-coloring c:([n]s)[r] of subsets of size s of [n] admits a subset Y([n]k) for which c is monochromatic on (Ys). Such a number always exists by (the finite version of) Ramsey’s theorem (see, e.g., [15]) and can be asymptotically upper-bounded [12] by a power tower of height s. The core of the argument is that, in the o(logn) regime, we can use Ramsey’s theorem to pick identifiers that force the output of the deterministic 𝖫𝖮𝖢𝖠𝖫 algorithm to be constant almost everywhere (i.e., on all but o(n) nodes).

5.2 Upper Bounds

Finally we show upper bound results, that is, how each of the classes of Table 5 corresponds to obtaining an approximation to a desired factor α.

Lemma 5.5.

If αβoptβconst holds for some opt-LCL Π with obj=min, then the complexity class of an algorithm that finds an α-approximation of Π is A.

Proof.

αβoptβconst implies there exists a constant solution that approximates the best possible solution sufficiently well for any n, as there is a self-loop in the de Bruijn graph. An algorithm finding an α-approximation of Π can output a constant solution in O(1) rounds in deterministic 𝖫𝖮𝖢𝖠𝖫 and randomized 𝖫𝖮𝖢𝖠𝖫. Together with Lemma 5.4 this concludes the proof.

Lemma 5.6.

Suppose αβopt<βconst for some min--LCL Π and additionally one of the following two conditions holds:

  1. 1.

    αβoptβgap and δgap is false.

  2. 2.

    αβopt>βgap and δgap is true.

Then the complexity class of any algorithm finding an α-approximation of Π is B.

Proof.

As the definitions of βgap and δgap are simply more extensive than βflex and δflex, for deterministic 𝖫𝖮𝖢𝖠𝖫 we can proceed as in Lemma 5.8. For randomized 𝖫𝖮𝖢𝖠𝖫 it is known that we can find a ruling set in O(1) with high probability [5]. With a ruling set and K, as in Lemma 5.8, constant sequences can be labeled. Any superconstant-length sequences can be labeled with the self-loop included in the definition of βgap. With high probability the cycle will use the cost of the self-loop only rarely even if it is higher than βgap, so it is negligible on sufficiently large instances. Together with Lemma 5.3 this concludes the proof.

The connection between the definition of βconst and a solution of an opt-LCL in a directed cycle is comparatively easy to make. It is less obvious how the cheapest cycle in a flexible component is the deciding cost for finding labelings in complexity class C. Therefore we make a connection to Section 3 with Lemma 5.7 leading up to Lemma 5.8.

Lemma 5.7.

If βflex holds for a min--LCL Π, then there exists K0+ such that, for every KK0, the de Bruijn graph of Π admits a closed walk W of length K for which c(W)βflex+o(1).

Proof.

Since βflex, there exists at least one cycle of GdB through which βflex can be obtained. Let C be one of these cycles. Because C is in a flexible component S, there exists a K0 such that for all KK0 we can find a closed walk W of length K starting and ending in a node belonging to C. This is achieved by connecting two closed walks of coprime lengths in S to C. Now we compose a walk W of length K by finding a K[K0,K0+|C|1] so that KK is a multiple of |C|. The number of different walks we have to find to achieve this is clearly upper bounded by |C|. For large K, the cost contribution of W is negligible compared to the cost contribution of the copies of C, whose average cost is c(C)=βflex. Choosing K0 large enough handles the worst of the finitely many possible walks W, as their potential additional cost impact is o(1) compared to the length K.

Lemma 5.8.

Suppose one of the following two conditions holds for some min--LCL Π:

  1. 1.

    αβoptβflex and δflex is false.

  2. 2.

    αβopt>βflex and δflex is true.

Then the complexity class of any algorithm finding an α-approximation of Π is C, assuming Π does not belong to either class A or B.

Proof.

If the first condition holds, then δflex is false, so there are two closed walks of coprime lengths sharing a node and having cost exactly βflex. By the same Frobenius argument as in Lemma 4.3, concatenating these walks gives closed walks of cost exactly βflex for all sufficiently large lengths. If the second condition holds, Lemma 5.7 lets us choose K large enough so that the o(1) term is small enough. In either case, there exists a constant K for which any path of length K can be labeled so that the resulting solution approximates the best possible solution sufficiently well for any n. A 3-coloring can be found in Θ(logn) [10, 20]. Then this coloring can be used to find sequences of length at least K which can then be labeled to produce a sufficient approximation in each sequence of the graph and therefore also in the whole graph. Together with Lemma 5.1 this concludes the proof.

Lemma 5.9.

If αβoptβcoprime for some min-max-LCL Π, the complexity class of any algorithm finding an α-approximation of Π is C, assuming Π does not belong to class A.

Proof.

The proof is almost the same as the proof concerning the first condition in Lemma 5.8. The existence of the constant K for which any path of length K can be labeled according to requirements is here directly rooted in the definition of βcoprime. Together with Lemma 5.2 this concludes the proof.

Lemma 5.10.

If βflex= for a min--LCL Π or βcoprime= for a min-max-LCL Π, then the opt-LCL has no general solution and its complexity class is E. If the α-approximation algorithm is not in any of the classes A, B, C or E, then it is in the class D.

Proof.

If the relevant parameter is , then the problem does not have a valid solution on all instances. Otherwise, we know that αβoptn<βflexn+O(1) (αβopt<βcoprime) holds for infinitely many n, so on certain instances the whole graph must be known to all nodes in order to obtain an α-approximation.

6 Putting Things Together

To recap, given an opt-LCL Π, Section 4 allows us to calculate the values of all parameters βa and δa in polynomial time, and given these parameters, Section 5 lets us find out what is the round complexity of finding an α-approximation of Π in the deterministic 𝖫𝖮𝖢𝖠𝖫 and randomized 𝖫𝖮𝖢𝖠𝖫 models. We have:

Theorem 6.1.

Given an opt-LCL Π on directed cycles and a constant α, we can in polynomial time determine the asymptotic round complexity of finding an α-approximation of Π in the deterministic 𝖫𝖮𝖢𝖠𝖫 and randomized 𝖫𝖮𝖢𝖠𝖫 models.

References

  • [1] Amirreza Akbari, Xavier Coiteux-Roy, Francesco D’Amore, François Le Gall, Henrik Lievonen, Darya Melnyk, Augusto Modanese, Shreyas Pai, Marc-Olivier Renou, Václav Rozhon, and Jukka Suomela. Online locality meets distributed quantum computing. In Michal Koucký and Nikhil Bansal, editors, Proceedings of the 57th Annual ACM Symposium on Theory of Computing, STOC 2025, Prague, Czechia, June 23-27, 2025, pages 1295–1306. ACM, 2025. doi:10.1145/3717823.3718211.
  • [2] Alkida Balliu, Sebastian Brandt, Yi-Jun Chang, Dennis Olivetti, Mikaël Rabie, and Jukka Suomela. The distributed complexity of locally checkable problems on paths is decidable. In Peter Robinson and Faith Ellen, editors, Proceedings of the 2019 ACM Symposium on Principles of Distributed Computing, PODC 2019, Toronto, ON, Canada, July 29 - August 2, 2019, pages 262–271. ACM, 2019. doi:10.1145/3293611.3331606.
  • [3] Alkida Balliu, Sebastian Brandt, Yi-Jun Chang, Dennis Olivetti, Jan Studený, and Jukka Suomela. Efficient classification of locally checkable problems in regular trees. In Christian Scheideler, editor, 36th International Symposium on Distributed Computing, DISC 2022, October 25-27, 2022, Augusta, Georgia, USA, volume 246 of LIPIcs, pages 8:1–8:19. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2022. doi:10.4230/LIPIcs.DISC.2022.8.
  • [4] Alkida Balliu, Sebastian Brandt, Yi-Jun Chang, Dennis Olivetti, Jan Studený, Jukka Suomela, and Aleksandr Tereshchenko. Locally checkable problems in rooted trees. Distributed Computing, 36(3):277–311, 2023. doi:10.1007/S00446-022-00435-9.
  • [5] Alkida Balliu, Mohsen Ghaffari, Fabian Kuhn, and Dennis Olivetti. Node and edge averaged complexities of local graph problems. In Alessia Milani and Philipp Woelfel, editors, PODC ’22: ACM Symposium on Principles of Distributed Computing, Salerno, Italy, July 25 - 29, 2022, pages 4–14. ACM, 2022. doi:10.1145/3519270.3538419.
  • [6] Thomas Boudier, Fabian Kuhn, Augusto Modanese, Ronja Stimpert, and Jukka Suomela. Classification of local optimization problems in directed cycles. CoRR, abs/2602.13046, 2026. doi:10.48550/arXiv.2602.13046.
  • [7] Sebastian Brandt, Juho Hirvonen, Janne H. Korhonen, Tuomo Lempiäinen, Patric R. J. Östergård, Christopher Purcell, Joel Rybicki, Jukka Suomela, and Przemyslaw Uznanski. LCL problems on grids. In Elad Michael Schiller and Alexander A. Schwarzmann, editors, Proceedings of the ACM Symposium on Principles of Distributed Computing, PODC 2017, Washington, DC, USA, July 25-27, 2017, pages 101–110. ACM, 2017. doi:10.1145/3087801.3087833.
  • [8] Yi-Jun Chang, Tsvi Kopelowitz, and Seth Pettie. An exponential separation between randomized and deterministic complexity in the LOCAL model. SIAM Journal on Computing, 48(1):122–143, 2019. doi:10.1137/17M1117537.
  • [9] Yi-Jun Chang, Jan Studený, and Jukka Suomela. Distributed graph problems through an automata-theoretic lens. Theoretical Computer Science, 951:113710, 2023. doi:10.1016/J.TCS.2023.113710.
  • [10] Richard Cole and Uzi Vishkin. Deterministic coin tossing with applications to optimal parallel list ranking. Information and Control, 70(1):32–53, 1986. doi:10.1016/S0019-9958(86)80023-7.
  • [11] Andrzej Czygrinow, Michal Hanckowiak, and Wojciech Wawrzyniak. Fast distributed approximations in planar graphs. In Gadi Taubenfeld, editor, Distributed Computing, 22nd International Symposium, DISC 2008, Arcachon, France, September 22-24, 2008. Proceedings, volume 5218 of Lecture Notes in Computer Science, pages 78–92. Springer, 2008. doi:10.1007/978-3-540-87779-0_6.
  • [12] P. Erdős and R. Rado. Combinatorial theorems on classifications of subsets of a given set. Proceedings of the London Mathematical Society, s3-2(1):417–439, 1952. doi:10.1112/plms/s3-2.1.417.
  • [13] Mika Göös, Juho Hirvonen, and Jukka Suomela. Lower bounds for local approximation. Journal of the ACM, 60(5):39:1–39:23, 2013. doi:10.1145/2528405.
  • [14] Mika Göös and Jukka Suomela. No sublogarithmic-time approximation scheme for bipartite vertex cover. Distributed Computing, 27(6):435–443, 2014. doi:10.1007/S00446-013-0194-Z.
  • [15] Ronald L. Graham, Bruce L. Rothschild, and Joel H. Spencer. Ramsey Theory. Wiley Series in Discrete Mathematics and Optimization. Wiley, second edition, 1991.
  • [16] David G. Harris. Distributed local approximation algorithms for maximum matching in graphs and hypergraphs. SIAM Journal on Computing, 49(4):711–746, 2020. doi:10.1137/19M1279241.
  • [17] Ozan Heydt, Simeon Kublenz, Patrice Ossona de Mendez, Sebastian Siebertz, and Alexandre Vigny. Distributed domination on sparse graph classes. European Journal of Combinatorics, 123:103773, 2025. doi:10.1016/J.EJC.2023.103773.
  • [18] Fabian Kuhn, Thomas Moscibroda, and Roger Wattenhofer. Local computation: Lower and upper bounds. Journal of the ACM, 63(2):17:1–17:44, 2016. doi:10.1145/2742012.
  • [19] Christoph Lenzen and Roger Wattenhofer. Leveraging Linial’s locality limit. In Gadi Taubenfeld, editor, Distributed Computing, 22nd International Symposium, DISC 2008, Arcachon, France, September 22-24, 2008. Proceedings, volume 5218 of Lecture Notes in Computer Science, pages 394–407. Springer, 2008. doi:10.1007/978-3-540-87779-0_27.
  • [20] Nathan Linial. Locality in distributed graph algorithms. SIAM Journal on Computing, 21(1):193–201, 1992. doi:10.1137/0221015.
  • [21] Moni Naor. A lower bound on probabilistic algorithms for distributive ring coloring. SIAM Journal on Discrete Mathematics, 4(3):409–412, 1991. doi:10.1137/0404036.
  • [22] Moni Naor and Larry J. Stockmeyer. What can be computed locally? SIAM Journal on Computing, 24(6):1259–1277, 1995. doi:10.1137/S0097539793254571.
  • [23] David Peleg. Distributed Computing: A Locality-Sensitive Approach. Society for Industrial and Applied Mathematics, 2000. doi:10.1137/1.9780898719772.
  • [24] Aleksandr Tereshchenko and Jan Studený. LCL classifier, 2021. URL: https://lcl-classifier.cs.aalto.fi.