Lower Bounds and Optimal Algorithms for Non-Smooth Convex Decentralized Optimization over Time-Varying Networks

Authors: dmitry  kovalev
ekaterina  borodich
alexander  gasnikov
dmitrii  feoktistov
Venue: NeurIPS_cc_2024_Conference
Published: 2024-05-14 16:33:30

Abstract: We consider the task of minimizing the sum of convex functions stored in a decentralized manner across the nodes of a communication network. This problem is relatively well-studied in the scenario when the objective functions are smooth, or the links of the network are fixed in time, or both. In particular, lower bounds on the number of decentralized communications and (sub)gradient computations required to solve the problem have been established, along with matching optimal algorithms. However, the remaining and most challenging setting of non-smooth decentralized optimization over time-varying networks is largely underexplored, as neither lower bounds nor optimal algorithms are known in the literature. We resolve this fundamental gap with the following contributions: (i) we establish the first lower bounds on the communication and subgradient computation complexities of solving non-smooth convex decentralized optimization problems over time-varying networks; (ii) we develop the first optimal algorithm that matches these lower bounds and offers substantially improved theoretical performance compared to the existing state of the art.

Main paper: 
 introduction: 
 In this paper, we study the decentralized optimization problem. Specifically, given a set of n compute nodes connected through a communication network, our goal is to solve the following finite-sum optimization problem with quadratic regularization:EQUATIONwhere r ≥ 0 is a regularization parameter, and each function f i (x) : R d → R is stored on the corresponding node i ∈ {1, . . . , n}. Each node i can perform computations based on its local state and data, and can directly communicate with other nodes through the links in the communication network.Decentralized optimization problems find applications in a wide variety of fields. These include network resource allocation (Beck et al., 2014) , distributed model predictive control (Giselsson et al., 2013) , power system control (Gan et al., 2012) , distributed spectrum sensing (Bazerque and Giannakis, 2009) , and optimization in sensor networks (Rabbat and Nowak, 2004) . In addition, such problems cover the supervised training of machine learning models through empirical risk minimization, thus attracting significant interest from the machine learning community (Lian et al., 2017; Ryabinin et al., 2021; Ryabinin and Gusev, 2020) .In our paper, we focus on the setting in which the links in the communication network are allowed to change over time. Such time-varying networks (Zadeh, 1961; Kolar et al., 2010) hold significant relevance to many practical applications. For instance, in sensor networks, changes in the link structure can be caused by the motion of sensors and disturbances in the wireless signal connecting pairs of sensors. Similarly, in distributed machine learning, connections between compute nodes can intermittently appear and disappear due to network unreliability (Ryabinin and Gusev, 2020) . Lastly, we anticipate that the time-varying setting will be supported by future-generation federated learning systems (Konecnỳ et al., 2016; McMahan et al., 2017) , where communication between pairs of mobile devices or between mobile devices and servers will be affected by their physical proximity, which naturally changes over time.In this work, we consider the decentralized optimization problem in the case when the objective function is convex (or strongly convex). At first glance, it may seem that the convexity assumption is restrictive and should not be considered. However, as we will see further, even in this fundamental setting, the existing algorithmic developments are limited and have significant gaps that need to be closed. Moreover, considering the convex optimization setting offers important benefits compared to general non-convex functions. One such benefit is that convex optimization often serves as a source of inspiration for the development of algorithms that turn out to be highly effective in solving practical problems, even non-convex ones.For example, state-of-the-art optimization algorithms such as Adam (Kingma and Ba, 2014) and RMSProp (Hinton et al., 2012) employ the momentum trick, which is observed to be efficient for numerous tasks, including the training of deep neural networks. However, from the perspective of non-convex optimization theory, momentum is useless because, for non-convex problems, the iteration complexity of the standard gradient method cannot be improved (Carmon et al., 2020) . On the other hand, it was theoretically proven that momentum substantially boosts the convergence speed of the gradient method when applied to convex functions (Nesterov, 1983) . In other words, convex optimization theory suggests that the momentum trick should be used, while non-convex theory suggests that it should not, and the former aligns much more closely with practical observations. A similar situation can be seen with other state-of-the-art optimization methods, including distributed local gradient methods (Mishchenko et al., 2022; Sadiev et al., 2022; Karimireddy et al., 2020) , adaptive gradient methods (Duchi et al., 2011) , etc. Such inconsistency between non-convex theoretical convergence guarantees for optimization algorithms and their actual performance in practice can be attributed to the fact that the class of non-convex functions is far too broad. This is why many optimization research papers try to narrow down this class by considering additional assumptions such as Polyak-Łojasiewicz condition (Karimi et al., 2016) , bounded non-convexity (Carmon et al., 2018; Allen-Zhu, 2018) , quasi-strong convexity (Necoara et al., 2019) , etc. However, these assumptions can be seen as relaxations of the standard convexity property. Therefore, we naturally opt to focus on the convex decentralized optimization problem, leaving potential generalizations for future work.Decentralized optimization has been attracting a lot of attention for more than a decade. Plenty of algorithms have been developed, including EXTRA (Shi et al., 2015) , DIGing (Nedic et al., 2017) , SONATA (Scutari and Sun, 2019) , NIDS (Li et al., 2019) , APM-C (Li et al., 2018; Rogozin et al., 2021) , and many others. In recent years, the focus of the research community has shifted towards the more complex task of finding, in some sense, the best possible algorithms for solving decentralized Table 1 : Summary of the existing state-of-the-art results in decentralized convex optimization. Multiple paper references are provided for each problem setting: papers marked with * provide lower complexity bounds, and papers marked with † provide optimal algorithms that match the corresponding lower bounds.Non-Smooth Setting Fixed Networks Algorithm 1 (this paper) † Theorems 1 and 2 (this paper) * optimization problems (Scaman et al., 2017 (Scaman et al., , 2018;; Lan et al., 2020; Kovalev et al., 2020 Kovalev et al., , 2021b Kovalev et al., ,a, 2022;; Hendrikx et al., 2021; Li et al., 2022; Li and Lin, 2021; Metelev et al., 2024) . This task consists of finding a lower bound on the complexity5 of solving a given subclass of decentralized problems and finding an algorithm whose complexity matches this lower bound. Such algorithms are called optimal because their complexity cannot be improved for a given problem class due to the established lower bounds.We discuss the four main classes of decentralized optimization problems that cover smooth6 and nonsmooth objective functions, and fixed and time-varying communication networks. We reference the existing state-of-the-art research papers that collectively solve the task of finding optimal algorithms for these classes. These papers are summarized in Table 1 . In the case of smooth and strongly convex objective functions and fixed communication networks, Scaman et al. (2017) established the lower bounds on the number of communication rounds and the number of local gradient computations required to find the solution. These lower bounds were matched by OPAPC algorithm of Kovalev et al. (2020) . In the case of smooth and strongly convex problems over time-varying networks, lower complexity bounds were provided by Kovalev et al. (2021a) , and two optimal algorithms were developed: ADOM+ (Kovalev et al., 2021a) and AccGT (Li and Lin, 2021) . In the case of non-smooth convex problems over fixed networks, lower bounds were established by Scaman et al. (2018) , and two optimal algorithms were proposed: DCS (Lan et al., 2020) and MSPD (Scaman et al., 2018) .Our paper primarily focuses on the remaining and most challenging setting of non-smooth convex decentralized optimization problems over time-varying networks. Only a few algorithms have been developed for this setting, including the distributed subgradient method (D-SubGD) by Nedic and Ozdaglar (2009) , the subgradient-push method (SubGD-Push) by Nedić and Olshevsky (2014) , and ZOSADOM by Lobanov et al. (2023) . Moreover, to the best of our knowledge, neither lower complexity bounds nor optimal algorithms have been proposed in this setting. Consequently, in this work, we close this significant gap with the following key contributions:(i) We establish the first lower bounds on the number of decentralized communications and local subgradient computations required to solve problem (1) in the non-smooth convex setting over time-varying networks, (ii) We show that our lower bounds are tight by developing the first optimal algorithm that matches these lower bounds. The proposed algorithm has state-of-the-art theoretical communication complexity, which outclasses the existing methods described in the literature. 
notation and assumptions: 
 In this paper, we are going to use the following notations: ⊗ denotes the Kronecker matrix product, I p denotes a p × p identity matrix, 1 p = (1, . . . , 1) ⊤ ∈ R p , e p j ∈ R p for j ∈ {1, . . . , p} denotes the j-th unit basis vector, where p ∈ {1, 2, . . .}. In addition, ∥•∥ denotes the standard Euclidean norm of a vector, and ⟨•, •⟩ denotes the standard scalar product of two vectors.Further, we describe the assumptions that we impose on problem 1. As discussed in Section 1.2, we assume the convexity of the objective function in problem (1). In particular, we assume that functions f 1 (x), . . . , f n (x) are convex, which is formally described in Assumption 1.Assumption 1. Each function f i (x) is convex. That is, for all x ′ , x ∈ R d and τ ∈ [0, 1], the following inequality holds:EQUATIONIn addition, we assume that the objective functions f 1 (x), . . . , f n (x) are Lipschitz continuous, which is formalized in Assumption 2. This property is widely used in the theoretical analysis of non-smooth optimization algorithms, such as the subgradient method (Nesterov, 2013) , dual extrapolation method (Nesterov, 2009) , etc. Assumption 2. Each functionf i (x) is M -Lipschitz continuous for M ≥ 0. That is, for all x ′ , x ∈ R d ,the following inequality holds:EQUATIONWe also need the following Assumption 3, which ensures the existence of a solution to problem (1). Note that in the strongly convex case (r > 0), the solution always exists and is unique. However, in the convex case (r = 0), we need to explicitly assume the existence of a solution. Assumption 3. There exists a solution x * ∈ R d to problem (1) and a distance R > 0 such that∥x * ∥ ≤ R. (4)Next, we formally describe the decentralized communication setting. The communication network is typically represented by a graph G(V, E), where V = {1, . . . , n} is the set of compute nodes and E ⊂ V × V is the set of links in the network. As mentioned earlier, we allow the communication links to change over time. Thus, we introduce the continuous time parameter τ ≥ 0 and a set-valued function E(τ ) : R + → 2 V×V , which represents the time-varying set of edges.7 Our time-varying network is then denoted as G(τ ) = (V, E(τ )).Decentralized communication is typically represented via a matrix-vector multiplication with the so-called gossip matrix associated with the communication network (Scaman et al., 2017; Kovalev et al., 2021a) . In the time-varying setting, we represent the gossip matrix by a matrix-valued function W(τ ) : R + → R n×n , which satisfies the following Assumption 4. Assumption 4. For all τ ≥ 0, the gossip matrix W(τ ) ∈ R n×n associated with the time-varying communication network G(V, E(τ )) satisfies the following properties:(i) W(τ ) ij = 0 if i ̸ = j and (j, i) / ∈ E(τ ), (ii) W(τ )1 n = 0 and W(τ ) ⊤ 1 n = 0.We also define the so-called condition number of the network χ ≥ 1, which indicates how well the network G(τ ) is connected (Scaman et al., 2017; Kovalev et al., 2021a) . In particular, the communication complexity of most decentralized optimization algorithms depends on χ. Assumption 5 provides the formal definition of this quantity. Assumption 5. There exists a constant χ ≥ 1 such that the following inequality holds for all τ ≥ 0:EQUATION)3 Lower Complexity BoundsQuestion: Does the paper discuss the limitations of the work performed by the authors? Answer: [Yes] Justification: The theoretical results provided in the paper require certain assumptions described in Section 2. Guidelines:• The answer NA means that the paper has no limitation while the answer No means that the paper has limitations, but those are not discussed in the paper. • The authors are encouraged to create a separate "Limitations" section in their paper.• The paper should point out any strong assumptions and how robust the results are to violations of these assumptions (e.g., independence assumptions, noiseless settings, model well-specification, asymptotic approximations only holding locally). The authors should reflect on how these assumptions might be violated in practice and what the implications would be. • The authors should reflect on the scope of the claims made, e.g., if the approach was only tested on a few datasets or with a few runs. In general, empirical results often depend on implicit assumptions, which should be articulated. • The authors should reflect on the factors that influence the performance of the approach.For example, a facial recognition algorithm may perform poorly when image resolution is low or images are taken in low lighting. Or a speech-to-text system might not be used reliably to provide closed captions for online lectures because it fails to handle technical jargon. • The authors should discuss the computational efficiency of the proposed algorithms and how they scale with dataset size. • If applicable, the authors should discuss possible limitations of their approach to address problems of privacy and fairness. • While the authors might fear that complete honesty about limitations might be used by reviewers as grounds for rejection, a worse outcome might be that reviewers discover limitations that aren't acknowledged in the paper. The authors should use their best judgment and recognize that individual actions in favor of transparency play an important role in developing norms that preserve the integrity of the community. Reviewers will be specifically instructed to not penalize honesty concerning limitations. 
decentralized subgradient optimization algorithms: 
 In this section, we present the lower bounds on the number of decentralized communications and the number of local subgradient computations required to solve problem (1). These lower bounds apply to a particular class of algorithms, which we refer to as the class of decentralized subgradient optimization algorithms. This class can be seen as an adaptation of black-box optimization procedures (Scaman et al., 2018) to the time-varying network setting, or an adaptation of first-order decentralized optimization algorithms (Kovalev et al., 2021a) to the non-smooth optimization setting.Non-smooth optimization algorithms typically perform incremental updates by computing the subgradient of a given objective function. The set of all subgradients of a convex function, called the subdifferential, can be multivalued in general. Thus, it is necessary to select the specific subgradient that the algorithm will use. This is done by the subgradient oracle, which is described by Definition 1.Definition 1. For each i ∈ V, a function ∇f i (x) : R d → R d is called a subgradient oracle associated with the function f i (x) if, for all x ∈ R d , it satisfies ∇f i (x) ∈ ∂f i (x). That is, for each i ∈ V and for all x, x ′ ∈ R d , the following inequality holds:EQUATIONFurther, we provide the formal description of the class of decentralized subgradient optimization algorithms in the following Definition 2. Definition 2. An algorithm is called a decentralized subgradient optimization algorithm with the subgradient computation time τ sub > 0 and decentralized communication time τ com > 0 if it satisfies the following constraints:(i) Internal memory. At any time τ ≥ 0, each node i ∈ V maintains an internal memory, which is represented by a set-valued function M i (τ ) : R + → 2 R d . The internal memory can be updated by subgradient computation or decentralized communication, which is formally represented by the following inclusion:EQUATIONwhere set-valued functions M sub i (τ ), M com i (τ ) : R + → 2 R d are defined below. (ii) Subgradient computation. At any time τ ≥ 0, each node i ∈ V can update its internal memory M i (τ ) by computing the subgradient ∇f i (x) of the function f i (x), which takes time τ sub . That is, for all τ ≥ 0, the set M sub i (τ ) is defined as follows:EQUATION)(iii) Decentralized communication. At any time τ ≥ 0, each node i ∈ V can update its internal memory M i (τ ) by performing decentralized communication across the communication network, which takes time τ com . That is, for all τ ≥ 0, the set M com i (τ ) is defined as follows:EQUATION(iv) Initialization and output. At time τ = 0, each node i ∈ V must initialize its internal memory with the zero vector, that is, M i (0) = {0}. At any time τ ≥ 0, each node i ∈ V must specify a single output vector from its internal memory, x o,i (τ ) ∈ M i (τ ).Now, we are ready to present the lower bounds on the execution time τ ≥ 0 required to find an ϵ-approximate solution8 to problem (1) by any algorithm satisfying Definition 2. Theorem 1 provides the lower bound in the strongly convex case (r > 0), and Theorem 2 provides the lower bound in the convex case (r = 0). These lower bounds naturally depend on the precision ϵ > 0, the parameters of the problem, including the Lipschitz constant M > 0, the regularization parameter r ≥ 0, the distance R > 0, and the parameters of the network, including the condition number χ ≥ 1, communication time τ com > 0, and subgradient computation time τ sub > 0.Table 2 : Lower bounds on the communication complexity of solving problem (1) in the centralized (Arjevani and Shamir, 2015) , decentralized fixed network (Scaman et al., 2018) , and decentralized time-varying network (Theorems 1 and 2) settings.SettingStrongly convex Ω (M/ √ rϵ) Ω √ χM/ √ rϵ Ω (χM/ √ rϵ) Convex Ω (M R/ϵ) Ω √ χM R/ϵ Ω (χM R/ϵ)Centralized Fixed networks9 Time-varying networks Theorem 1. For arbitrary parameters M, r, ϵ, τ com , τ sub > 0 and χ ≥ 1, there exists an optimization problem of the form (1) satisfying Assumptions 1, 2 and 3, corresponding subgradient oracles given by Definition 1, a time varying network G(τ ) = (V, E(τ )), and a corresponding time-varying gossip matrix W(τ ) satisfying Assumptions 4 and 5, such that at least the following time τ is required to reach precision p(x o,i (τ )) -p(x * ) ≤ ϵ by any decentralized subgradient optimization algorithm satisfying Definition 2:EQUATION)Theorem 2. For arbitrary parameters M, R, ϵ, τ com , τ sub > 0 and χ ≥ 1, there exists an optimization problem of the form (1) with zero regularization (r = 0) satisfying Assumptions 1, 2 and 3, corresponding subgradient oracles given by Definition 1, a time varying network G(τ ) = (V, E(τ )), and a corresponding time-varying gossip matrix W(τ ) satisfying Assumptions 4 and 5, such that at least the following time τ is required to reach precision p(x o,i (τ )) -p(x * ) ≤ ϵ by any decentralized subgradient optimization algorithm satisfying Definition 2:EQUATION)The proofs of Theorems 1 and 2 can be found in Appendix B. Further, we provide a brief and informal description of the main theoretical ideas that underlie these proofs:(i) We select a specific "hard" instance of problem (1). In particular, we choose the objective function of the form p(x) = a d-1 j=1 |⟨e d j+1 -e d j , x⟩| -a⟨e d 1 , x⟩ + r 2 ∥x∥ 2 , which was used by Arjevani and Shamir (2015) ; Scaman et al. (2018) in the proof of lower bounds on the communication complexity in centralized and fixed-network settings. One can show that the gap p(x) -p(x * ) is lower-bounded by a positive constant as long as the last component of the vector x is zero, and it takes Ω(τ sub • d) time to break this bound due to the constraint on the subgradient updates (8).(ii) We split the objective function between two nodes of a star-topology network with a timevarying central node, which was previously utilized by Kovalev et al. (2021a) in the proof of lower bounds for optimizing smooth functions. One can show that it takes Ω(n) = Ω(χ) communications to exchange information between the two selected nodes due to the timevarying center. This contrasts with the fixed path-topology network used by Scaman et al. (2017 Scaman et al. ( , 2018)) , where such an exchange would take Ω(n) = Ω( √ χ) communications.Moreover, using the constraint (8), we can show that it takes Ω(τ com • nd) time to make the last component of the vector x nonzero and break the lower bound on the gap p(x) -p(x * ), thanks to the way we split the objective function. (iii) Based on the above considerations, we show that the total execution time required to solve the problem is lower-bounded by Ω (τ com • nd + τ sub • d). Thus, we obtain the desired results by making a specific choice of the dimension d, network size n, and other parameters of problem (1).We compare the lower complexity bounds for solving non-smooth convex optimization problems in the three main distributed optimization settings: centralized, decentralized fixed network, andAlgorithm 11: input: x 0 = x -1 = x0 ∈ (R d ) n , y 0 = y 0 ∈ (R d ) n , z 0 = z 0 ∈ L ⊥ , m 0 ∈ (R d ) n 2: parameters: K, T ∈ {1, 2, . . .}, {(α k , β k , γ k , σ k , λ k , τ k x , η k x , η k y , η k z , θ k z )} K-1 k=0 ⊂ R 10 + 3: for k = 0, 1, . . . , K -1 do 4: y k = α k y k + (1 -α k )y k , z k = α k z k + (1 -α k )z k 5: g k y = ∇ y G(y k , z k ), g k z = ∇ z G(y k , z k ), where function G(y, z) is defined in eq. ( 12)6:gk z = (W k ⊗ I d )g k z , ĝk z = (W k ⊗ I d )(g k z + m k ), where W k denotes the gossip matrix W(τ ) at the current time τ 7:y k+1 = y k -η k y (g k y + xk+1 ), z k+1 = z k -η k z ĝk z , xk+1 = x k + γ k (x k -x k-1 )8:y k+1 = y k + α k (y k+1 -y k ), z k+1 = z k -θ k z gk z , m k+1 = (η k z /η k+1 z )(m k + g k z -ĝk z ) 9:x k,0 = x k 10:for t = 0, 1, . . . , T -1 do 11:g k,t x = ( ∇f 1 (x k,t 1 ), . . . , ∇f n (x k,t n ))12:x k,t+1 = x k,t -η k x g k,t x + β k x k,t+1 -y k+1 + τ k x (x k,t+1 -x k )13:x k+1 = σ k x k,T + (1 -σ k )x k+1 , xk+1 = 1 T T t=1 x k,t , x k+1 = α k xk+1 + (1 -α k )x k 14: (x K a , y K a , z K a ) = ( K k=1 λ k ) -1 K k=1 λ k (x k , y k , z k ) 15: output: x K o = 1 n n i=1 x K a,i ∈ R d , where (x K a,1 , . . . , x K a,n ) = x K a ∈ (R d ) ndecentralized time-varying network. The lower subgradient computation complexity bounds coincide in these cases (Nesterov (2013) , Scaman et al. (2018) ,Theorems 1 and 2). However, the situation with the communication complexity is different. See Table 2 for a summary.Theorems 1 and 2 imply that the communication complexity in the decentralized time-varying network setting is proportional to the network condition number χ. In contrast, the communication complexity in the fixed network setting is proportional to √ χ, which reflects the fact that time-varying networks are more difficult to deal with compared to fixed networks. In particular, there was a long-standing conjecture that the "upgrade" from the factor χ to the factor √ χ in communication complexity is impossible in the time-varying network setting. Only recently, this conjecture was proved for smooth functions by Kovalev et al. (2021a) , and now we resolve this open question in the non-smooth case as well.Question: For each theoretical result, does the paper provide the full set of assumptions and a complete (and correct) • The experimental setting should be presented in the core of the paper to a level of detail that is necessary to appreciate the results and make sense of them. • The full details can be provided either with the code, in appendix, or as supplemental material. 
optimal algorithm: 
 In this section, we develop an optimal algorithm for solving the non-smooth convex decentralized optimization problem (1) over time-varying networks. The design of our algorithm relies on a specific saddle-point reformulation of the problem, which we describe in the following section.Let functionsF (x) : (R d ) n → R and G(y, z) : (R d ) n × (R d) n → R be defined as follows:EQUATION)where x = (x 1 , . . . , x n ) ∈ (R d ) n , and r x , r yz > 0 are some constants that satisfy r x + 1/r yz = r.(13) Consider the following saddle-point problem:EQUATIONwhere L ⊥ ⊂ (R d ) n is the orthogonal complement to the so-called consensus space L ⊂ (R d ) n , defined as follows:L = {(x 1 , . . . , x n ) : x 1 = . . . = x n }, L ⊥ = {(x 1 , . . . , x n ) : n i=1 x i = 0}.(15) One can show that the saddle-point problem ( 14) is equivalent to the minimization problem (1). This is justified by the following Lemma 1. The proof of the lemma can be found in the Appendix A.Lemma 1. Problem ( 14) is equivalent to problem (1) in the following sense:min x∈(R d ) n max y∈(R d ) n max z∈L ⊥ Q(x, y, z) = n • min x∈R d p(x).(16)The saddle-point reformulation of the form ( 14) was first introduced by Kovalev et al. (2020 Kovalev et al. ( , 2021a) ) to develop optimal decentralized algorithms for optimizing smooth functions. However, these are not applicable to the non-smooth case. To the best of our knowledge, the only attempt to adapt the reformulation ( 14) to the non-smooth setting was made by Lobanov et al. (2023) . However, their results have significant downsides, which we discuss in Section 4.3.Now, we present Algorithm 1 for solving problem (1). We provide upper bounds on the number of decentralized communications K and the number of subgradient computations K × T required to find an ϵ-approximate solution to the problem. Theorems 3 and 4 provide the upper bounds in the strongly convex (r > 0) and convex (r = 0) cases, respectively. The proofs can be found in Appendix D. The total execution time of Algorithm 1 is upper-bounded asτ = O (τ com • K + τ sub • K × T ), where the communication time τ com > 0 and the subgradient computation time τ sub > 0 are described in Definition 2. This upper-bound on the execution time cannot be improved because of the lower bounds established in the previous Section 3. Therefore, Algorithm 1 is an optimal algorithm for solving problem (1). Theorem 3. Under Assumptions 1, 2, 3, 4 and 5, let r > 0 (strongly convex case). Then Algorithm 1 requires K = O χM √ rϵ decentralized communications (line 6 of Algorithm 1) and Assumptions 1, 2, 3, 4 and 5 K ×T = O M 2 rϵ subgradient computations (line 11 of Algorithm 1) to reach precision p(x K o ) -p(x * ) ≤ ϵ. Theorem 4. Under, let r = 0 (convex case). Then Algorithm 1 requires K = O χM R ϵ decentralized communications (line 6 of Algorithm 1) and K × T = O M 2 R 2 ϵ 2 subgradient computations (line 11 of Algorithm 1) to reach precision p(x K o ) -p(x * ) ≤ ϵ.The design of Algorithm 1 is based on the fundamental Forward-Backward algorithm (Bauschke and Combettes, 2011) . Let E = (R d ) n × (R d ) n × L ⊥ be a Euclidean space, and consider a monotone operator A(u) : E → E and a maximally-monotone multivalued operator B(u) : E → 2 E defined as follows:EQUATIONwhere u = (x, y, z) ∈ E, and P = ( 14) is equivalent to the following monotone inclusion problem:I n -(1/n)1 n 1 ⊤ n ) ⊗ I d ∈ R nd×nd is the orthogonal projection matrix onto L ⊥ . Then problem (EQUATIONThe basic Forward-Backward algorithm iteratesu k+1 = (id + B) -1 (u k -A(u k )), where id is the identity operator and (id + B) -1 denotes the inverse of the operator id(u) + B(u), which is called resolvent. Algorithm 1 can be obtained by making the following major modifications to these iterations:(i) We accelerate the convergence of the Forward-Backward algorithm using Nesterov acceleration (Nesterov, 1983) . Although this mechanism cannot be applied to the general monotone inclusion problem (18), Kovalev et al. (2020) showed that it can be used when the operator A(u) is equal to the gradient of a smooth convex function, which is true in our case.(ii) Computation of the operator A(u) requires multiplication with the matrix P. This, in turn, requires an exact averaging of a vector, which is difficult to do over the time-varying network. Kovalev et al. (2021b) showed that this obstacle can be tackled with the Error-Feedback mechanism for decentralized communication, which we also utilize. (iii) At each iteration of the algorithm, we have to compute the resolvent, which requires solving an auxiliary subproblem min x max y τx2 ∥x -x k ∥ 2 + F (x) -⟨y, x⟩ - τy 2 ∥y -y k ∥ 2. This problem cannot be solved exactly, so we have to find an approximate solution using an Table 3 : The execution time τ required to find an ϵ-approximate solution to the decentralized optimization problem (1) by the following algorithms: D-SubGD (Nedic and Ozdaglar, 2009) , SubGD-Push (Nedić and Olshevsky, 2014) , ZO-SADOM (Lobanov et al., 2023) , and Algorithm 1 (this paper). Decentralized communication and subgradient computation complexities are marked with green and yellow colors, respectively. For D-SubGD, the complexity is not provided because the algorithm converges only to a neighborhood of the solution. For SubGD-Push, poly(M, R, d) denotes a certain polynomial in M, R, d. For ZO-SADOM, the differences from the optimal complexities are highlighted in red color.Strongly-convex case complexity Convex case complexityD-SubGD SubGD-Push ZO-SADOM τ com • χM d 1/4 log 1 ϵ √ rϵ + τ sub • M 2 d log 1 ϵ rϵ τ com • χM Rd 1/4 log 1 ϵ ϵ + τ sub • M 2 R 2 d log 1 ϵ ϵ 2 Algorithm 1 τ com • χM √ rϵ + τ sub • M 2 rϵ τ com • χM R ϵ + τ sub • M 2 R 2 ϵ 2 Lower Bounds τ com • χM √ rϵ + τ sub • M 2 rϵ τ com • χM R ϵ + τ sub • M 2 R 2 ϵ 2 N/A τ com • poly(M, R, d) • n 2n log 2 1 ϵ ϵ 2 + τ sub • poly(M, R, d) • n 2n log 2 1 ϵ ϵ 2additional "inner" algorithm based on the subgradient method (Nesterov, 2013) and the Chambolle-Pock operator splitting (Chambolle and Pock, 2011) . We also have to conduct a careful analysis to find an efficient way to combine the inner and the "outer" Forward-Backward algorithms and avoid unnecessary waste of subgradient calls.The design of Algorithm 1 shares some similarities with the algorithm of Kovalev et al. (2021a) such as (i) and (ii) above. However, Kovalev et al. (2021a) simply add the gradient ∇F (x) to the operator A(u) and use the accelerated version of the Forward-Backward algorithm, which we obviously cannot do as the function F (x) is not smooth. Instead, we have to put the subdifferential ∂F (x) into the operator B(u) and follow (iii) above. Part (iii), in turn, shares some similarities with the algorithm of Lan et al. (2020) . However, Lan et al. (2020) simply have a zero operator A(u) = 0, which makes (i) and (ii) above unnecessary in their case. In contrast, we cannot make such simplifications because we work in the much more complicated setting of time-varying networks.One could naturally expect that the existing optimal algorithms, originally developed for fixed networks, such as DCS (Lan et al., 2020) and MSPD (Scaman et al., 2018) , could be applied to solve problem (1) over time-varying networks. However, this is not the case, which is justified by the lack of corresponding theoretical guarantees and was shown empirically by Kovalev et al. (2021b) . Therefore, we have to consider only those algorithms that were specifically developed for the time-varying network setting.We provide a comparison of our Algorithm 1 with the existing state-of-the-art decentralized methods for solving convex non-smooth optimization problems over time-varying networks in Table 3 . 10These include D-SubGD (Nedic and Ozdaglar, 2009) , SubGD-Push (Nedić and Olshevsky, 2014) , and ZO-SADOM (Lobanov et al., 2023) . The first two algorithms have poor performance: D-SubGD converges only to limited precision, and SubGD-Push converges at a slow rate of O(log 2 (1/ϵ)/ϵ 2 ), which does not match even the iteration complexity of the standard centralized subgradient method, let alone the improved complexity of Algorithm 1. The complexity of ZO-SADOM is also worse than the lower bounds. Moreover, the theoretical results of Lobanov et al. (2023) have substantial drawbacks compared to ours:(i) Lobanov et al. (2023) do not provide any theoretical insights or innovations in the analysis of their algorithm. In particular, they use the randomized smoothing technique (Duchi et al., 2012) to obtain a smooth approximation of the objective p(x), and apply the existing algorithm of Kovalev et al. (2021a) to minimize this approximation. In contrast, we develop a new algorithm that directly works with the original non-smooth objective p(x).(ii) ZO-SADOM has extra factors d 1/4 log(1/ϵ) and d log(1/ϵ) in the decentralized communication and subgradient computation complexities, respectively, compared to the optimal complexity of our Algorithm 1. Thus, the performance of ZO-SADOM can be poor when applied, for instance, to large-scale machine learning problems in which the dimension d can be huge.The orthogonal complement L ⊥ to the consensus space L is given as follows:EQUATION)Let us perform the maximization of Q(x, y, z) in the variable y ∈ (R d ) n :max y∈(R d ) n Q(x, y, z) (a) = max y∈(R d ) n F (x) + ⟨y, x⟩ -G(y, z) (b) = F (x) + max y∈(R d ) n ⟨y, x⟩ - r yz 2 ∥y + z∥ 2 = F (x) + 1 2r yz ∥x∥ 2 -⟨x, z⟩,where (a) uses the definition of Q(x, y, z) eq. ( 14); (b) uses the definition of G(y, z) in eq. ( 12). Next, we perform maximization in the variable z ∈ L ⊥ :max z∈L ⊥ max y∈(R d ) n Q(x, y, z) = max z∈L ⊥ F (x) + 1 2r yz ∥x∥ 2 -⟨x, z⟩ = F (x) + 1 2r yz ∥x∥ 2 + max z∈L ⊥ [-⟨x, z⟩] = F (x) + 1 2r yz ∥x∥ 2 + I L (x),where I L (x) : (R d ) n → R is the indicator function, which is defined as follows:EQUATION)Now, we can rewrite the saddle-point problem (14) as followsmin x∈(R d ) n max y∈(R d ) n max z∈L ⊥ Q(x, y, z) (a) = min x∈(R d ) n max z∈L ⊥ max y∈(R d ) n Q(x, y, z) (b) = min x∈(R d ) n F (x) + 1 2r yz ∥x∥ 2 + I L (x) (c) = min x∈(R d ) n n i=1 f i (x i ) + r x + 1/r yz 2 ∥x i ∥ 2 + I L (x) (d) = min x∈(R d ) n n i=1 f i (x i ) + r 2 ∥x i ∥ 2 + I L (x) (e) = n • min x∈R d p(x).where (a) uses the fact that we can exhange the order of the two consecutive maximizations; (b) uses the previous equation; (c) uses the definition of F (x) in eq. ( 12); (d) uses eq. ( 13); (e) uses the definition of p(x) in eq. ( 1) and the definition of I L (x).B.1 The Hard Instance of Problem (1)Compute nodes. In this proof, we consider the case when χ ≥ 3. The case χ < 3 can be proven using the fixed-network argument of Scaman et al. (2018) . We choose n = 3⌊χ/3⌋, which implies that n ≥ 3 and n mod 3 = 0. We also divide the set of nodes V = {1, . . . , n} into the following three disjoint subsets: V 1 = {1, . . . , n/3}, V 2 = {n/3 + 1, . . . , 2n/3} and V 3 = {2n/3 + 1, . . . , n}.Objective functions. We fix an arbitrary odd integer d ∈ {3, 5, 7, . . .} and define functions f 1 (x), . . . , f n (x) : R d → R as follows:EQUATION)where a > 0 is an arbitrary constant and functions h 1 (x), . . . , h d-1 (x) : R d → R are defined as follows:EQUATION)Consequently, the objective function p(x) in problem ( 1) is given as follows:EQUATION)We also define the subgradient oracles ∇f 1 (x), . . . , ∇f n (x) : R d → R d as follows:EQUATIONwhere ∇h 1 (x), . . . , ∇h d-1 (x) : R d → R d are the subgradient oracles associated with functions h 1 (x), . . . , h d-1 (x), defined as follows: ∇h j (x) =    e d j+1 -Time-varying network. We choose the time-varying network G(τ ) = (V, E(τ )) to be a startopology undirected graph with the time-varying center node i c (τ ) ∈ V. Formally, we define the edges of the time-varying network E(τ ) ⊂ V × V as follows:EQUATION)We also specify the center node i c (τ ) at a given time τ ≥ 0 as follows:EQUATION)We choose the time-varying gossip matrix W(τ ) ∈ R n×n to be the Laplacian matrix of the graph G(τ ). Formally, W(τ ) is defined as follows:EQUATION)where deg i (τ ) denotes the degree of the node i ∈ V in the graph G(τ ), i.e.,EQUATION)One can observe, that the time-varying gossip matrix W(τ ) satisfies Assumption 4, in particular, ker W(τ ) = ker W(τ ) ⊤ = span({1 n }). Moreover, one can show that W(τ ) is a symmetric matrix, and λ max (W(τ )) = 1 and λ + min (W(τ )) = 1/n ≥ 1/χ. Hence, W(τ ) satisfies Assumption 5.Further, we define linear spaces K 0 , . . . , K d ⊂ R d as follows:EQUATIONIn order to prove Theorems 1 and 2, we will use the following auxiliary lemmas. The proofs of these lemmas can be found in Appendix C. Furthermore, the proof of Theorem 1 is contained in Appendix B.3, and the proof of Theorem 2 is contained in Appendix B.4. Lemma 2. For all τ ≥ 0, the following statements hold:(i) Let i ∈ V 1 .Then, for all j ∈ {1, . . . , (d -1)/2},EQUATION) (ii) Let i ∈ V 2 .Then, for all j ∈ {0, . . . , (d -1)/2},EQUATION) (ii) Let i ∈ V 3 .Then, for all j ∈ {0, . . . , d},EQUATION)The proof of Lemma 2 is contained in Appendix C.1. Lemma 3. Let k ∈ {0, . . . , n(d -1)/6 -1}. Then, for all τ < (k + 1)τ com , the following inclusion holds:EQUATIONwhere p = ⌊3k/n⌋ and q = k mod (n/3).The proof of Lemma 3 is contained in Appendix C.2. Lemma 4. Let functions f 1 , . . . , f n (x) be defined by eq. ( 21). Then problem eq.(1) has a unique solution x * ∈ R d , which is given as follows:EQUATION)Moreover, for all x ∈ K d-1 , the following inequality holds:EQUATION)The proof of Lemma 4 is contained in Appendix C.3.Decentralized communication. Lemma 3 implies that M i (τ ) ⊂ K d-1 as long as τ < τ com • n(d -1)/6. Hence, Lemma 4 implies eq. ( 36) for all x ∈ M i (τ ) as long as τ < τ com • n(d -1)/6. Let the constant a > 0 be chosen as follows:EQUATIONThen, each function f i (x) defined by eq. ( 21) is M -Lipschitz. Indeed, the case i ∈ V 3 is trivial. In the case when i ∈ V 1 , we can prove the M -Lipschitz continuity of f i (x) as follows:f i (x) -f i (x ′ ) = a (d-1)/2 j=1 ⟨x, e d 2j -e d 2j-1 ⟩ -⟨x ′ , e d 2j -e d 2j-1 ⟩ -a⟨x -x ′ , e d 1 ⟩ ≤ a (d-1)/2 j=1 ⟨x -x ′ , e d 2j -e d 2j-1 ⟩ + a ⟨x -x ′ , e d 1 ⟩ ≤ a (d-1)/2 j=1 ⟨x -x ′ , e d 2j ⟩ + ⟨x -x ′ , e d 2j-1 ⟩ + a ⟨x -x ′ , e d 1 ⟩ = a d-1 j=1 ⟨x -x ′ , e d j ⟩ + a ⟨x -x ′ , e d 1 ⟩ ≤ 2a d j=1 ⟨x -x ′ , e d j ⟩ ≤ 2a √ d∥x -x ′ ∥ ≤ M ∥x -x ′ ∥.In the case when i ∈ V 2 , we can prove the M -Lipschitz continuity of f i (x) similarly.Without loss of generality, we assume ϵ ≤ M 2 /(576r) and define d ∈ {3, 5, . . .} as follows:EQUATION)Using eqs. ( 37) and ( 38), for all τ < τ com • n(d -1)/6 and x ∈ M i (τ ), we obtainp(x) -p(x * ) ≥ M 2 36rd 2 > ϵ.Hence, to reach precision p(x) -p(x * ) ≤ ϵ for some x ∈ M i (τ ), it is necessary that τ satisfiesEQUATIONSubgradient computation. We also need to prove that to reach precision p(x) -p(x * ) ≤ ϵ for some x ∈ M i (τ ), it is necessary that τ satisfiesEQUATIONWe can do this by providing an extended version of our hard problem instance, described in Appendix B.1. In particular, we consider the following instance of problem (1):EQUATION)where functions f 1 (x), . . . , f n (x) : R d → R are defined in Appendix B.1 by eq. ( 21), and functionsf ′ 1 (x ′ ), . . . , f ′ n (x ′ ) : R d ′ →R are defined as follows:EQUATIONwhere b > 0 is some constant. Then, by choosing an appropriate subgradient oracle ∇f ′ i (x ′ ) associated with each function f ′ i (x ′ ) (see Section 3.2.1 of Nesterov ( 2013)) we can obtain both lower bounds (39) and (40), which concludes the proof.Our proof of Theorem 2 is very similar to the proof of Theorem 1 with the following differences. Let function h δ (x) : R d → R be the Huber function, which is defined as follows:EQUATION)Note that function h δ (x) is continuously differentiable and ( √ dδ)-Lipschitz continuous.In the proof of Theorem 3 we used functions f 1 (x), . . . , f n (x) defined in eq. ( 21) of Appendix B.1.Here we use a slightly different choice, that is, functions f 1 (x), . . . , f n (x) are defined as follows:EQUATION)Consequently, our hard instance of problem (1), which is described in Appendix B.1, turns into the following:EQUATIONwhere c > 0 is some constant, and functions h 1 (x), . . . , h d-1 (x) are defined in eq. ( 22).One can show that Lemmas 2 and 3 still hold true. We can also replace Lemma 4 with the following Lemma 5. The proof of this lemma is a trivial extension of the proof of Lemma 4, which uses the fact that ∇( 1 2 ∥•∥ 2 )(x * ) = ∇h δ (x * ) as long as δ and x * are defined by eq. ( 46) and eq. ( 47), respectively. Lemma 5. Let δ be defined as follows:EQUATIONProblem eq. ( 45) has a solution x * ∈ R d , which is given as follows:EQUATION)Moreover, for all x ∈ K d-1 , the following inequality holds:EQUATION)One can also show that each function f i (x) defined in eq. ( 44) is M f -Lipschitz continuous, where M f is defined as follows:M f = 2a √ d + cδ √ d = 2a √ d + a/(3 √ d) ≤ 3a √ d.(49) Let us choose a and c as follows:EQUATION)This choice of a and c implies M f ≤ M and ∥x * ∥ ≤ R. Moreover, eq. ( 48) impliesEQUATIONas long as x * ∈ K d-1 . Next, without loss of generality we can assume ϵ ≤ (M R)/72 and choose d ∈ {3, 5, . . .} as follows:EQUATIONwhich, for all x ∈ M i (τ ), implies p(x) -p(x * ) > ϵ as long as τ satisfiesEQUATIONwhich concludes the proof.C Proofs of Lemmas from Section B.2C.1 Proof of Lemma 2Statement (i). Let i ∈ V 1 and x ⊂ K 2j for j ∈ {1, . . . , (d -1)/2}. Then for l ≥ 2j + 1 we obtain ⟨e d l+1 -e d l , x⟩ = 0, which implies ∇h l (x) = 0 due to eq. ( 25). Hence, we obtain the following:1 a ∇f i (x) (a) = ∇h 1 (x) + ∇h 3 (x) + • • • + ∇h d-2 (x) -e d 1 (b) = ∇h 1 (x) + ∇h 3 (x) + • • • + ∇h 2j-1 (x) -e d 1 (c) ⊂ span {e d 1 , e d 2 } ∪ • • • ∪ {e d 2j-1 , e d 2j } (d) ⊂ K 2j ,where (a) uses eq. ( 24); (b) uses the fact that ∇h l (x) = 0 for l ≥ 2j + 1; (c) uses eq. ( 25); (d) uses the definition of K 2j in eq. ( 30). Hence,M i (τ ) ⊂ K 2j implies M sub i (τ + τ sub ) ⊂ K 2j by the definition of M sub i (•) in eq. (8).Statement (ii). Let i ∈ V 2 and x ⊂ K 2j+1 for j ∈ {0, . . . , (d -1)/2}. Then for l ≥ 2j + 2 we obtain ⟨e d l+1 -e d l , x⟩ = 0, which implies ∇h l (x) = 0 due to eq. ( 25). Hence, we obtain the following:1 a ∇f i (x) (a) = ∇h 2 (x) + ∇h 4 (x) + • • • + ∇h d-1 (x) (b) = ∇h 2 (x) + ∇h 4 (x) + • • • + ∇h 2j (x) (c) ⊂ span {e d 2 , e d 3 } ∪ • • • ∪ {e d 2j , e d 2j+1 } (d) ⊂ K 2j ,where (a) uses eq. ( 24); (b) uses the fact that ∇h l (x) = 0 for l ≥ 2j + 2; (c) uses eq. ( 25); (d) uses the definition of K 2j+1 in eq. (30). Hence, M i (τ ) ⊂ K 2j+1 implies M sub i (τ + τ sub ) ⊂ K 2j+1 by the definition of M sub i (•) in eq. ( 8).Statement (iii). This statement is trivially implied by the definition of ∇f i (x) in eq. ( 24) and the definition of M sub i (•) in eq. ( 8).We prove the lemma using the induction on k.Base case: k = 0. In this case, we assume τ < (k + 1)τ com = τ com . Hence, for all i ∈ V, we obtainM com i (τ ) = ∅ and M i (τ ) ⊂ M sub i (τ ).Using Lemma 2 and the fact that M com i (τ ) = ∅, we can easily obtainM i (τ ) ⊂ M sub i (τ ) ⊂    K 2 i ∈ V 1 K 1 i ∈ V 2 K 0 i ∈ V 3, which implies the desired eq. ( 34) for k = p = q = 0.Induction hypothesis. Let k ′ ∈ {0, 1, 2, . . .}. We assume that eq. ( 34) holds for all τ < (k ′ + 1)τ com , that is,EQUATIONwhere p ′ = ⌊3k ′ /n⌋ and q ′ = k ′ mod (n/3).Induction step. We assume that the induction hypothesis (54) is true. Our goal is to prove that eq. ( 34) holds for k = k ′ + 1. When 0 ≤ τ < kτ com , the desired eq. ( 34) is implied by the induction hypothesis (54). Thus, we can assume kτ com ≤ τ < (k + 1)τ com . Further, we consider two cases: q ̸ = 0 and q = 0.Induction step, case q ̸ = 0. In this case, p = p ′ and q = q ′ + 1.Part (i). First, we consider the caseEQUATIONEquation ( 55) implies τ -τ sub < (k ′ + 1)τ com and τ -τ com < (k ′ + 1)τ com . Using the induction hypothesis (54) and the fact that p ′ = p and q ′ = q -1, we getEQUATION)Hence, using Lemma 2, we obtainEQUATION)Equations ( 27) and ( 55) imply i c (τ ) = 2n/3 + q + 1. Hence, using eq. ( 56), we getM ic(τ ) (τ -τ com ) ⊂ K 2p+1 .For i ̸ = i c (τ ), using eqs. ( 9) and ( 56), we getEQUATIONFor i = i c (τ ) = 2n/3 + q + 1, using eqs. ( 9) and ( 56), we getEQUATIONHence, using eqs. ( 58) and ( 59), for all i ∈ V, we obtainEQUATIONNow, we combine eqs. ( 57) and (60), and obtainEQUATIONThus, we were able to prove eq. ( 34) for τ satisfying (55).Part (ii). We can prove the general casekτ com ≤ τ < min{(k + 1)τ com , kτ com + lτ sub }for arbitrary l ∈ {1, 2, . . .} using the induction on l. The only difference compared to the proof in the previous part is in eq. ( 56), which will change toM i (τ -τ sub ) ⊂ K 2p+2 i ∈ V 1 or (i ∈ V 3 and i ≤ 2n/3 + q + 1) K 2p+1 i ∈ V 2 or (i ∈ V 3 and i > 2n/3 + q + 1) ,and eq. ( 57) will change as follows due to Lemma 2:M sub i (τ ) ⊂ K 2p+2 i ∈ V 1 or (i ∈ V 3 and i ≤ 2n/3 + q + 1) K 2p+1 i ∈ V 2 or (i ∈ V 3 and i > 2n/3 + q + 1) .However, the rest of the proof, including eq. ( 61) will remain unchanged.Induction step, case q = 0. In this case p = p ′ + 1 and q ′ = n/3 -1.Part (i). First, we consider the caseEQUATIONEquation ( 62) implies τ -τ sub < (k ′ + 1)τ com and τ -τ com < (k ′ + 1)τ com . Using the induction hypothesis (54) and the fact that p ′ = p -1 and q ′ = n/3 -1, we getEQUATION)Equations ( 27) and ( 62) imply i c (τ ) = 2n/3 + 1. Using eq. ( 63), we getM ic(τ ) (τ -τ com ) ⊂ K 2p .For i ̸ = i c (τ ), using eqs. ( 9) and ( 63), we getEQUATION)For i = i c (τ ) = 2n/3 + 1, using eqs. ( 9) and ( 63), we getEQUATIONHence, using eqs. ( 64) and ( 65), for all i ∈ V, we obtainEQUATIONUsing Lemma 2, from eq. ( 63) we obtainEQUATIONHence, using eqs. ( 66) and ( 67), for all i ∈ V, we obtainEQUATIONwhich implies eq. ( 34) for τ satisfying (62).Part (ii). Next, we consider the caseEQUATIONEquation (66) still holds for all i ∈ V and τ satisfying eq. ( 69). From eqs. ( 68) and ( 69), for all i ∈ V, we obtainM i (τ -τ sub ) ⊂ K 2p ,which, due to Lemma 2, implies the following:EQUATION)Hence, using eqs. ( 66) and (70), we obtainEQUATIONwhich implies eq. ( 34) for τ satisfying (69).Part(iii). We can prove the general caseEQUATIONfor l ∈ {2, 3, . . .} using the induction on l. There will be no differences compared to the proof in the previous part. Indeed, eqs. ( 66) and (71) will still hold for all i ∈ V and τ satisfying eq. ( 72).One can show, that x * defined in eq. ( 35) is indeed the unique minimizer of the function p(x) defined in eq. ( 23). Moreover, we can obtain the following:p(x * ) = - a 2 18rd.We can lower-bound function p(x) as follows:p(x) = r 2 ∥x∥ 2 - a 3 ⟨e d 1 , x⟩ + a 3 d-1 j=1 ⟨x, e d j+1 -e d j ⟩ ≥ - a 3 ⟨e d 1 , x⟩ + a 3 d-1 j=1 ⟨x, e d j ⟩ -⟨x, e d j+1 ⟩ = - a 3 ⟨e d d , x⟩ = 0as long as x ∈ K d-1 . Hence, for all x ∈ K j , we obtainp(x) -p(x * ) ≥ a 2 18rd ,which concludes the proof.In order to prove Theorems 3 and 4, we will use the following auxiliary lemmas. The proofs of these lemmas can be found in Appendix E. Furthermore, the proof of Theorem 3 is contained in Appendix D.2, and the proof of Theorem 4 is contained in Appendix D.3. Lemma 6. Under Assumptions 1, 2 and 3, let r > 0 (strongly convex case). Then there exists a solution (w * , y * , z * ) ∈ L × (R d ) n × L ⊥ to problem (1), which satisfies the following conditions0 ∈ ∂ x Q(w * , y * , z * ), 0 = ∇ y Q(w * , y * , z * ), L ∋ ∇ z Q(w * , y * , z * ).(73) Moreover, the following inequalities hold:EQUATION)The proof of Lemma 6 is contained in Appendix E.1. Lemma 7. Under Assumptions 1 and 2, let η 0 x , . . . , η K-1x and β 0 , . . . , β K-1 be chosen as follows:η k x = 1/(τ k x T ), β k = r x , σ k = τ k x /(2τ k x + β k ) for k ∈ {0, . . . , K -1}.(75) Then, for all x ∈ (R d ) n and k ∈ {0, . . . , K -1}, the following inequality holds:EQUATION)The proof of Lemma 7 is contained in Appendix E.2. Lemma 8. Under Assumption 4, for all k ∈ {0, . . . , K -1}, the iterates of Algorithm 1 satisfyEQUATIONwhere P ∈ R nd×nd is the orthogonal projection matrix onto L ⊥ , which is given as follows:EQUATIONThe proof of Lemma 8 is contained in Appendix E.3. Lemma 9. Under Assumptions 4 and 5, for all k ∈ {0, . . . , K -1} the following inequality holds:EQUATION)The proof of Lemma 9 is contained in Appendix E.4. Lemma 10. Under Assumptions 4 and 5, let parameters θ 0 z , . . . , θ K-1 z be chosen as follows: θ k z = 1/(2r yz ) for k = 0, . . . , K -1.(80) Then, for all k ∈ {0, . . . , K -1}, the following inequality holds:EQUATIONThe proof of Lemma 10 is contained in Appendix E.5. Lemma 11. Under Assumptions 1 and 2 and under conditions of Lemmas 7 and 10, let parameters α 0 , . . . , α K-1 and γ 0 , . . . , γ K-1 be chosen as follows:α k = 3/(k + 3), γ k = (k + 2)/(k + 3) for k = 0, . . . , K -1. (82) Let parameters τ 0 x , . . . , τ K-1 x , η 0 y , . . . , η K-1 y, and η 0 z , . . . , η K-1 z be chosen as follows:EQUATION) where τ x , η y and η z are defined as follows:τ x = 1 2 r x , η y = (4r yz ) -1 , η z = (10r yz χ 2 ) -1 , r x = 2 3 r, r yz = 3/r.(84) Let parameters λ 1 , . . . , λ K be chosen as follows:λ K = α -2 K-1 and λ k = α -2 k-1 + α -1 k -α -2 k for k = 1, . . . , K -1.(85) Let the input of Algorithm 1 be chosen as follows:x 0 = 0, y 0 = 0, z 0 = 0, m 0 = 0. (86) Then, for all x, y ∈ (R d ) n and z ∈ L ⊥ , the following inequality holds:EQUATION)The proof of Lemma 11 is contained in Appendix E.6.We can upper-bound rx 2 ∥x K a -w * ∥ 2 , where w * is defined in Lemma 6, as follows:r x 2 ∥x K a -w * ∥ 2 (a) ≤ Q(x K a , y * , z * ) -Q(w * , y * , z * ) (b) ≤ Q(x K a , y * , z * ) -Q(w * , y K a , z K a ) (c) ≤ 2 K 2 r∥w * ∥ 2 + 18 r ∥y * ∥ 2 + 45χ 2 r ∥z * ∥ 2 + 72nM 2 rKT (d) ≤ 2 K 2 nM 2 r + 18(1 + r x /r) 2 nM 2 r + 180nχ 2 M 2 r + 72nM 2 rKTwhere (a) uses Lemma 6 and the strong convexity of Q(x, y, z) in x; (b) and (d) use Lemma 6; (c) uses Lemma 11. Using the definition of r x in eq. ( 84)r∥x K a -w * ∥ 2 ≤ 6 K 2 51nM 2 r + 180nχ 2 M 2 r + 72nM 2 rKT ≤ 1386nχ 2 M 2 rK 2 + 72nM 2 rKT .Next, we can upper-bound n(p(x K o ) -p(x * )) as follows:n(p(x K o ) -p(x * )) (a) = n i=1 f i (x K o ) -f i (x * ) + r 2 ∥x K o ∥ 2 - r 2 ∥x * ∥ 2 (b) = n i=1 f i (x K o ) -f i (x * ) + r 2 ∥ 1 n n j=1 x K a,j ∥ 2 - r 2 ∥x * ∥ 2 (c) ≤ n i=1 f i (x K o ) -f i (x * ) + r 2 ∥x K a,i ∥ 2 - r 2 ∥x * ∥ 2 (d) = n i=1 f i (x K o ) -f i (x * ) + r 2 ∥x K a ∥ 2 - r 2 ∥w * ∥ 2 (e) ≤ n i=1 f i (x K a,i ) -f i (x * ) + M ∥x K a,i -x K o ∥ + r 2 ∥x K a ∥ 2 - r 2 ∥w * ∥ 2 (f ) = F (x K a ) -F (w * ) + 1 2r yz ∥x K a ∥ 2 - 1 2r yz ∥w * ∥ 2 + n i=1 M ∥x K a,i -x K o ∥ (g) ≤ F (x K a ) -F (w * ) + 1 2r yz ∥x K a ∥ 2 - 1 2r yz ∥w * ∥ 2 + n i=1 M 2 n i=1 ∥x K a,i -x K o ∥ 2 (h) = F (x K a ) -F (w * ) + 1 2r yz ∥x K a ∥ 2 - 1 2r yz ∥w * ∥ 2 + √ nM ∥x K a ∥ Pwhere (a) uses the definition of p(x) in eq. ( 1); (b) uses the definition of x K o on line 15 of Algorithm 1; (c) uses the convexity of ∥•∥ 2 ; (d) uses the definition of w * in eq. ( 95); (e) uses Assumption 2; (f) uses the definition of function F (x) in eq. ( 12) and eq. ( 13); (g) uses the Cauchy-Schwarz inequality; (h) uses the definition of P in eq. ( 78).Next, for arbitrary z ∈ L ⊥ we define y = -r -1 yz x K a -z. Then, we get∇ y Q(x K a , y, z) = 0 and Q(x K a , y, z) = F (x K a ) + 1 2ryz ∥x K a ∥ 2 + ⟨x K a , z⟩.Plugging this into the previous upper-bound gives the following:n(p(x K o ) -p(x * )) ≤ Q(x K a , y, z) -F (w * ) - 1 2r yz ∥w * ∥ 2 -⟨x K a , z⟩ + √ nM ∥x K a ∥ P (a) = Q(x K a , y, z) -Q(w * , y * , z * ) -⟨x K a , z⟩ + √ nM ∥x K a ∥ P (b) ≤ Q(x K a , y, z) -Q(w * , y K a , z K a ) -⟨x K a , z⟩ + √ nM ∥x Ka ∥ P where (a) uses the definition of y * in eq. ( 98) and the definition of z * in eq. ( 99); (b) uses Lemma 6.Next, we choose z ∈ L ⊥ as follows:EQUATION)Then, ⟨x K a , z⟩ = + √ nM ∥x K a ∥ P and we obtain the following:n(p(x K o ) -p(x * )) ≤ Q(x K a , y, z) -Q(w * , y K a , z K a ) (a) ≤ 2 K 2 r∥w * ∥ 2 + 18 r ∥y∥ 2 + 45χ 2 r ∥z∥ 2 + 72nM 2 rKT (b) = 2 K 2 r∥w * ∥ 2 + 18 r ∥r -1 yz x K a + z∥ 2 + 45χ 2 r ∥z∥ 2 + 72nM 2 rKT = 2 K 2 r∥w * ∥ 2 + 18 r ∥r -1 yz (x K a -w * + w * ) + z∥ 2 + 45χ 2 r ∥z∥ 2 + 72nM 2 rKT (c) ≤ 2 K 2 r∥w * ∥ 2 + 54 rr 2 yz ∥x K a -w * ∥ 2 + 54 rr 2 yz ∥w * ∥ 2 + 54 r ∥z∥ 2 + 45χ 2 r ∥z∥ 2 + 72nM 2 rKT ≤ 2 K 2 r∥w * ∥ 2 + 54 rr 2 yz ∥x K a -w * ∥ 2 + 54 rr 2 yz ∥w * ∥ 2 + 99χ 2 r ∥z∥ 2 + 72nM 2 rKT (d) ≤ 2 K 2 r∥w * ∥ 2 + 54 rr 2 yz ∥x K a -w * ∥ 2 + 54 rr 2 yz ∥w * ∥ 2 + 99nχ 2 M 2 r + 72nM 2 rKT (e) ≤ 2 K 2 nM 2 r + 54nM 2 r 3 r 2 yz + 54 rr 2 yz ∥x K a -w * ∥ 2 + 99nχ 2 M 2 r + 72nM 2 rKT (f ) = 2 K 2 7nM 2 r + 6r∥x K a -w * ∥ 2 + 99nχ 2 M 2 r + 72nM 2 rKT ≤ 212nχ 2 M 2 rK 2 + 72nM 2 rKT + 12r K 2 ∥x K a -w * ∥ 2 (g) ≤ 212nχ 2 M 2 rK 2 + 72nM 2 rKT + 12 K 2 1386nχ 2 M 2 rK 2 + 72nM 2 rKT ,where 84); (g) uses the previously obtained upper-bound on r∥x K a -w * ∥ 2 . Dividing both sides of the inequality by n gives the following:p(x K o ) -p(x * ) ≤ 212χ 2 M 2 rK 2 + 72M 2 rKT + 12 K 2 1386χ 2 M 2 rK 2 + 72M 2 rKT .Hence, choosing the parameters K and T such thatK ≥ O χM √ rϵ and K × T ≥ O M 2 rϵ implies p(x K o ) -p(x *) ≤ ϵ, which concludes the proof.With r = 0, the original problem (1) turns into the following problem:EQUATION)Let x * ∈ R d be the solution to problem (89), such that ∥x * ∥ ≤ R, which always exists due to Assumption 3. Let r > 0 be an arbitrary regularization parameter. We can upper-bound function f (x) using the regularized objective function p(x) defined in eq. ( 1) as follows:f (x) ≤ f (x) + r 2 ∥x∥ 2 = p(x).On the other hand, we can lower-bound f (x * ) as follows:f (x * ) = p(x * ) - r 2 ∥x * ∥ 2 ≥ min x ′ ∈R d p(x ′ ) - r 2 ∥x * ∥ 2 ≥ min x ′ ∈R d p(x ′ ) - rR 2 2 .Hence, we can upper-bound the function suboptimality gap in problem (89) as follows:f (x) -f (x * ) ≤ p(x) -min x ′ ∈R d p(x ′ ) + rR 2 2 .Let the regularization parameter r > 0 be chosen as follows:EQUATION)Then, we obtain the following:EQUATION)We can apply Algorithm 1 to solving the regularized problem (1) with the regularization parameter r defined in eq. ( 90). Theorem 3 implies that, to reach precisionp(x K o ) -min x ′ ∈R d p(x ′ ) ≤ ϵ 2 (92)it is sufficient to perform the following number of decentralized communications:EQUATIONand the following number of subgradient computations:EQUATION)where (a) and (b) use the definition of r in eq. ( 90). Using eqs. ( 91) and ( 92), we also obtain the desired precision f (x K o ) -f (x * ) ≤ ϵ, which concludes the proof.First, we pick the solution x * ∈ R d to problem (1), which is unique due to Assumption 3 and the fact that r > 0. Next, we define w * ∈ L as follows:EQUATION)From Assumptions 1 and 2 it follows that dom p(x) = R d and dom f i (x) = R d for all i ∈ {1, . . . , n}, which implies the following:EQUATION)Hence, there exists a vector∆ * = (∆ * 1 , . . . , ∆ * n ) ∈ (R d ) n such that ∆ * i ∈ ∂f i (x *) for all i ∈ {1, . . . , n}, and the following relation holds:EQUATION)Next, we define y * ∈ (R d ) n as follows:EQUATION)From Assumptions 1 and 2 it follows that domF (x) = (R d ) n , which implies y * ∈ ∂F (w * ) and 0 ∈ ∂(F (•) -⟨y * , •⟩)(w * ) = ∂ x Q(w * , y * , z * ).Next, we define z * ∈ L ⊥ as follows:EQUATION)Note that the inclusion z * ∈ L ⊥ is implied by eq. ( 97). Further, we get∇ z Q(w * , y * , z * ) (a) = -r yz (y * + z * ) (b) = -r yz (r x -r)w * ∈ L,where (a) uses the definition of Q(x, y, z) in eq. ( 14); (b) uses the definition of y * and z * , and the last inclusion follows from the definition of w * . Moreover, we obtain the following∇ y Q(w * , y * , z * ) (a) = -w * -r yz (y * + z * ) (b) = -r yz (r -1 yz + r x -r)w * (c) = 0,where (a) uses the definition of Q(x, y, z) in eq. ( 14); (b) uses the definition of y * and z * ; (c) uses eq. ( 13).From Assumption 2 it follows that ∥∆ * i ∥ ≤ M for all i ∈ {1, . . . , n}. Hence, using eq. ( 97), we get r∥x * ∥ ≤ M , which implies ∥w * ∥ 2 ≤ nM 2 /r 2 . Moreover, we get∥y * ∥ ≤ ∥∆ * ∥ + r x ∥w * ∥ ≤ √ n(M + r x M/r) = (1 + r x /r) √ nM, which implies ∥y * ∥ 2 ≤ (1 + r x /r) 2 nM 2 . Finally, we obtain ∥z * ∥ ≤ r∥w * ∥ + ∥∆ * ∥ ≤ 2 √ nM,which implies ∥z * ∥ 2 ≤ 4nM 2 and concludes the proof.We start with the following upper-bound on1 2η k x ∥x k,t+1 -x∥ 2 : 1 2η k x ∥x k,t+1 -x∥ 2 (a) = 1 2η k x ∥x k,t -x∥ 2 - 1 2η k x ∥x k,t+1 -x k,t ∥ 2 + 1 η k x ⟨x k,t+1 -x k,t , x k,t+1 -x⟩ (b) = 1 2η k x ∥x k,t -x∥ 2 - 1 2η k x ∥x k,t+1 -x k,t ∥ 2 -⟨g k,t x + β k x k,t+1 -y k+1 + τ k x (x k,t+1 -x k ), x k,t+1 -x⟩ = 1 2η k x ∥x k,t -x∥ 2 - 1 2η k x ∥x k,t+1 -x k,t ∥ 2 + ⟨y k+1 , x k,t+1 -x⟩ -⟨β k x k,t+1 + τ k x (x k,t+1 -x k ), x k,t+1 -x⟩ -⟨g k,t x , x k,t+1 -x⟩ (c) ≤ 1 2η k x ∥x k,t -x∥ 2 - 1 2η k x ∥x k,t+1 -x k,t ∥ 2 + ⟨y k+1 , x k,t+1 -x⟩ - τ k x 2 ∥x k,t+1 -x k ∥ 2 - τ k x 2 ∥x k,t+1 -x∥ 2 + τ k x 2 ∥x k -x∥ 2 - β k 2 ∥x k,t+1 ∥ 2 - β k 2 ∥x k,t+1 -x∥ 2 + β k 2 ∥x∥ 2 -⟨g k,t x , x k,t -x⟩ -⟨g k,t x , x k,t+1 -x k,t ⟩ (d) ≤ 1 2η k x ∥x k,t -x∥ 2 - 1 2η k x ∥x k,t+1 -x k,t ∥ 2 + ⟨y k+1 , x k,t+1 -x⟩ - τ k x 2 ∥x k,t+1 -x k ∥ 2 - τ k x 2 ∥x k,t+1 -x∥ 2 + τ k x 2 ∥x k -x∥ 2 - β k 2 ∥x k,t+1 ∥ 2 - β k 2 ∥x k,t+1 -x∥ 2 + β k 2 ∥x∥ 2 + n i=1 (f i (x i ) -f i (x k,t i ) -⟨g k,t x,i , x k,t+1 i -x k,t i ⟩),where (g k,tx,1 , . . . , g k,t x,n ) = ( ∇f 1 (x k,t 1 ), . . . , ∇f n (x k,t n )) = g k,t x ∈ (R d ) n , (1 2η k x ∥x k,t+1 -x∥ 2 (a) ≤ 1 2η k x ∥x k,t -x∥ 2 - 1 2η k x ∥x k,t+1 -x k,t ∥ 2 + ⟨y k+1 , x k,t+1 -x⟩ - τ k x 2 ∥x k,t+1 -x k ∥ 2 - τ k x 2 ∥x k,t+1 -x∥ 2 + τ k x 2 ∥x k -x∥ 2 - β k 2 ∥x k,t+1 ∥ 2 - β k 2 ∥x k,t+1 -x∥ 2 + β k 2 ∥x∥ 2 + n i=1 (f i (x i ) -f i (x k,t+1 i ) + M ∥x k,t+1 i -x k,t i ∥ + ∥g k,t x,i ∥∥x k,t+1 i -x k,t i ∥) (b) ≤ 1 2η k x ∥x k,t -x∥ 2 - 1 2η k x ∥x k,t+1 -x k,t ∥ 2 + ⟨y k+1 , x k,t+1 -x⟩ - τ k x 2 ∥x k,t+1 -x k ∥ 2 - τ k x 2 ∥x k,t+1 -x∥ 2 + τ k x 2 ∥x k -x∥ 2 - β k 2 ∥x k,t+1 ∥ 2 - β k 2 ∥x k,t+1 -x∥ 2 + β k 2 ∥x∥ 2 + n i=1 (f i (x i ) -f i (x k,t+1 i ) + 2M ∥x k,t+1 i -x k,t i ∥) (c) ≤ 1 2η k x ∥x k,t -x∥ 2 - 1 2η k x ∥x k,t+1 -x k,t ∥ 2 + ⟨y k+1 , x k,t+1 -x⟩ - τ k x 2 ∥x k,t+1 -x k ∥ 2 - τ k x 2 ∥x k,t+1 -x∥ 2 + τ k x 2 ∥x k -x∥ 2 + F (x) -F (x k,t+1 ) - r x 2 ∥x k,t+1 -x∥ 2 + n i=1 1 2η k x ∥x k,t+1 i -x k,t i ∥ 2 + 2η k x M 2 = 1 2η k x ∥x k,t -x∥ 2 + τ k x 2 ∥x k -x∥ 2 + ⟨y k+1 , x k,t+1 -x⟩ + 2nη k x M 2 - τ k x 2 ∥x k,t+1 -x k ∥ 2 - τ k x + r x 2 ∥x k,t+1 -x∥ 2 -F (x k,t+1 ) + F (x),where (a) uses Assumption 2 and the Cauchy-Schwarz inequality; (b) uses the inequality ∥g k,t x,i ∥ ≤ M , which follows from Assumption 2; (c) uses the definition of β k in eq. ( 75), the definition of F (x) in eq. ( 12) and Young's inequality. After rearranging, we obtain1 2η k x ∥x k,t+1 -x∥ 2 ≤ 1 2η k x ∥x k,t -x∥ 2 + τ k x 2 ∥x k -x∥ 2 + 2nη k x M 2 -∆ k,t+1 ,where ∆ k,t+1 is defined as∆ k,t+1 = F (x k,t+1 ) -F (x) -⟨y k+1 , x k,t+1 -x⟩ + τ k x + r x 2 ∥x k,t+1 -x∥ 2 + τ k x 2 ∥x k,t+1 -x k ∥ 2Now, we sum these inequalities for t = 0, . . . , T -1 and obtain1 2η k x ∥x k,T -x∥ 2 ≤ 1 2η k x ∥x k,0 -x∥ 2 + τ k x T 2 ∥x k -x∥ 2 + 2nη k x M 2 T - T t=1 ∆ k,t .Dividing both sides of the inequality by T gives1 2η k x T ∥x k,T -x∥ 2 ≤ 1 2η k x T ∥x k,0 -x∥ 2 + τ k x 2 ∥x k -x∥ 2 + 2nη k x M 2 - 1 T T t=1 ∆ k,t .Using the definition of ∆ k,t , the definition of xk on line 13 of Algorithm 1 and Assumption 1, we obtain1 2η k x T ∥x k,T -x∥ 2 ≤ 1 2η k x T ∥x k,0 -x∥ 2 + τ k x 2 ∥x k -x∥ 2 + 2nη k x M 2 -F (x k+1 ) -F (x) -⟨y k+1 , xk+1 -x⟩ - τ k x + r x 2 ∥x k+1 -x∥ 2 + τ k x 2 ∥x k+1 -x k ∥ 2 .Using the definition of η k x and β k in eq. ( 75), we obtainτ k x 2 ∥x k,T -x∥ 2 + τ k x + β k 2 ∥x k+1 -x∥ 2 ≤ τ k x 2 ∥x k,0 -x∥ 2 + τ k x 2 ∥x k -x∥ 2 + 2nη k x M 2 -F (x k+1 ) -F (x) -⟨y k+1 , xk+1 -x⟩ + τ k x 2 ∥x k+1 -x k ∥ 2 .Using the definition of x k,0 on line 9 of Algorithm 1, the definition of x k+1 on line 13 of Algorithm 1, the definition of η k x , β k and σ k in eq. ( 75) and the convexity of ∥•∥, we obtain(τ k x + 1 2 r x )∥x k+1 -x∥ 2 ≤ τ k x ∥x k -x∥ 2 + 2nM 2 τ k x T -F (x k+1 ) -F (x) -⟨y k+1 , xk+1 -x⟩ + τ k x 2 ∥x k+1 -x k ∥ 2 ,which concludes the proof.Using Assumption 4, and the definition of P in eq. ( 78), we obtainEQUATIONThen, the desired relations can be trivially obtained by analyzing the lines of Algorithm 1.We can upper-bound ∥η k+1 z m k+1 ∥ 2 P as follows:∥η k+1 z m k+1 ∥ 2 P (a) = ∥η k z (m k + g k z -ĝk z )∥ 2 P (b) = ∥η k z (m k + g k z -(W k ⊗ I d )(m k + g k z ))∥ 2 P (c) = ∥η k z (P(m k + g k z ) -(W k ⊗ I d )P(m k + g k z ))∥ 2 (d) ≤ (1 -1/χ)∥η k z P(m k + g k z )∥ 2 (e) ≤ (1 -1/χ) (1 + 1/(2χ))∥η k z m k ∥ 2 P + (1 + 2χ)∥η k z g k z ∥ 2 P ≤ (1 -1/(2χ))∥η k z m k ∥ 2 P + 2χ∥η k z g k z ∥ 2 Pwhere (a) uses line 8 of Algorithm 1; (b) uses line 6 of Algorithm 1; (c) uses eq. ( 100); (d) uses Assumption 5; (e) uses the parallelogram rule and Young's inequality. Using this, we obtain∥η k z m k ∥ 2 P ≤ 2χ∥η k z m k ∥ 2 P -2χ∥η k+1 z m k+1 ∥ 2 P + 4χ 2 ∥η k z g k z ∥ 2 P , which concludes the proof.We can upper bound ∥g k z -Pg k z ∥ 2 as follows:∥g k z -Pg k z ∥ 2 (a) = ∥(W k ⊗ I d )g k z -Pg k z ∥ 2 (b) = ∥(W k ⊗ I d )Pg k z -Pg k z ∥ 2 (c) ≤ (1 -1/χ)∥g k z ∥ 2 Pwhere (a) uses line 6 of Algorithm 1; (b) uses eq. ( 100); (c) uses Assumption 5. On the other hand, ∥g k z -Pg k z ∥ 2 is equal to the following:∥g k z -Pg k z ∥ 2 (a) = ∥g k z ∥ 2 + ∥g k z ∥ 2 P -2⟨g k z , Pg k z ⟩ (b) = 1 (θ k z ) 2 ∥z k+1 -z k ∥ 2 + ∥g k z ∥ 2 P + 2 θ k z ⟨z k+1 -z k , Pg k z ⟩.where (a) uses the parallelogram rule; (b) uses line 8 of Algorithm 1. Hence, we obtain the following1 (θ k z ) 2 ∥z k+1 -z k ∥ 2 + 2 θ k z ⟨z k+1 -z k , Pg k z ⟩ + 1 χ ∥g k z ∥ 2 P ≤ 0.After rearranging and multiplying both sides of the inequality byθ k z 2α k , we obtain 0 ≥ α -1 k ⟨z k+1 -z k , Pg k z ⟩ + 1 2θ k z ∥z k+1 -z k ∥ 2 + θ k z 2α k χ ∥g k z ∥ 2 P (a) = α -1 k ⟨P(z k+1 -z k ), g k z ⟩ + r yz ∥z k+1 -z k ∥ 2 + 1 4α k χr yz ∥g k z ∥ 2 P (b) = α -1 k ⟨z k+1 -z k , g k z ⟩ + r yz ∥z k+1 -z k ∥ 2 + 1 4α k χr yz ∥g k z ∥ 2 Pwhere (a) uses eq. ( 80); (b) uses Lemma 8, which concludes the proof.Using eq. ( 101) and lines 7 and 8 of Algorithm 1, we obtainẑk+1 = ẑk + z k+1 -z k -P(η k+1 z m k+1 -η k z m k ) = ẑk -P(η k z ĝk z + η k+1 z m k+1 -η k z m k ) = ẑk -P(η k z ĝk z + η k+1 z (η k z /η k+1 z )(m k + g k z -ĝk z ) -η k z m k ) = ẑk -η k z Pg k z .Hence, we can upper-bound 1 2η k z ∥ẑ k+1 -z∥ 2 as follows:1 2η k z ∥ẑ k+1 -z∥ 2 (a) = 1 2η k z ∥ẑ k -z∥ 2 + 1 2η k z ∥ẑ k+1 -ẑk ∥ 2 + 1 η k z ⟨ẑ k+1 -ẑk , ẑk -z⟩ (b) = 1 2η k z ∥ẑ k -z∥ 2 + η k z 2 ∥g k z ∥ 2 P -⟨Pg k z , ẑk -z⟩ = 1 2η k z ∥ẑ k -z∥ 2 + η k z 2 ∥g k z ∥ 2 P -⟨Pg k z , z k -z⟩ + ⟨Pg k z , z k -ẑk ⟩ (c) = 1 2η k z ∥ẑ k -z∥ 2 + η k z 2 ∥g k z ∥ 2 P -⟨Pg k z , z k -z⟩ + η k z ⟨Pg k z , Pm k ⟩,where (a) uses the parallelogram rule; (b) uses the update rule for ẑk which we previously obtained; (c) uses eq. ( 101). Further, we can upper-bound the term η k z ⟨Pg k z , Pm k ⟩ as follows Plugging this into the previous inequality gives η k z ⟨Pg k z , Pm k ⟩ (a) ≤ 1 η k z ∥η k z g k z ∥ P ∥η k z m k ∥ P (b) ≤ 1 2η k z 2χ∥η k z g k z ∥ 2 P + 1 2χ ∥η k z m k ∥ 2 P (c) ≤ 1 2η k z 4χ∥η k z g k z ∥ 2 P + ∥η k z m k ∥ 2 P -1 2η k z ∥ẑ k+1 -z∥ 2 ≤ 1 2η k z ∥ẑ k -z∥ 2 + η k z 2 ∥g k z ∥ 2 P -⟨Pg k z , z k -z⟩ + 1 2η k z 4χ∥η k z g k z ∥ 2 P + ∥η k z m k ∥ 2 P -∥η k+1 z m k+1 ∥ 2 P (a) = 1 2η k z ∥ẑ k -z∥ 2 + η k z (1 + 4χ) 2 ∥g k z ∥ 2 P -⟨g k z , z k -z k + z k -z⟩ + 1 2η k z ∥η k z m k ∥ 2 P -∥η k+1 z m k+1 ∥ 2 P (b) = 1 2η k z ∥ẑ k -z∥ 2 + η k z (1 + 4χ) 2 ∥g k z ∥ 2 P + 1 2η k z ∥η k z m k ∥ 2 P -∥η k+1 z m k+1 ∥ 2 P + ⟨g k z , z -z k ⟩ + (1 -α k )α -1 k ⟨g k z , z k -z k ⟩ (c) ≤ 1 2η k z ∥ẑ k -z∥ 2 + η k z (1 + 4χ) 2 ∥g k z ∥ 2 P + 1 2η k z ∥η k z m k ∥ 2 P -∥η k+1 z m k+1 ∥ 2 P + ⟨g k z , z -z k ⟩ + (1 -α k )α -1 k ⟨g k z , z k -z k ⟩ -α -1 k ⟨z k+1 -z k , g k z ⟩ + r yz ∥z k+1 -z k ∥ 2 - 1 4α k χr yz ∥g k z ∥ 2 P (d) ≤ 1 2η k z ∥ẑ k -z∥ 2 + 1 2η k z ∥η k z m k ∥ 2 P -∥η k+1 z m k+1 ∥ 2 P + ⟨g k z , z -z k ⟩ + (1 -α k )α -1 k ⟨g k z , z k -z k ⟩ -α -1 k ⟨z k+1 -z k , g k z ⟩ + r yz ∥z k+1 -z k ∥ 21 2η k y ∥y k+1 -y∥ 2 + 1 2η k z ∥ẑ k+1 -z∥ 2 ≤ 1 2η k y ∥y k -y∥ 2 + 1 2η k z ∥ẑ k -z∥ 2 + 2η k y γ 2 k ∥x k-1 -xk ∥ 2 -⟨x k+1 , y k+1 -y⟩ + γ k ⟨x k-1 -xk , y k -y⟩ -⟨x k -xk+1 , y k+1 -y⟩ + 1 2η k z ∥η k z m k ∥ 2 P -∥η k+1 z m k+1 ∥ 2 P + ⟨g k y , y -y k ⟩ + ⟨g k z , z -z k ⟩ + (1 -α k )α -1 k ⟨g k y , y k -y k ⟩ + ⟨g k z , z k -z k ⟩ -α -1 k ⟨g k y , y k+1 -y k ⟩ + ⟨z k+1 -z k , g k z ⟩ + r yz ∥y k+1 -y k ∥ 2 + r yz ∥z k+1 -z k ∥ 2 (a) ≤ 1 2η k y ∥y k -y∥ 2 + 1 2η k z ∥ẑ k -z∥ 2 + 2η k y γ 2 k ∥x k-1 -xk ∥ 2 -⟨x k+1 , y k+1 -y⟩ + γ k ⟨x k-1 -xk , y k -y⟩ -⟨x k -xk+1 , y k+1 -y⟩ + 1 2η k z ∥η k z m k ∥ 2 P -∥η k+1 z m k+1 ∥ 2 P + G(y, z) -G(y k , z k ) + (1 -α k )α -1 k G(y k , z k ) -G(y k , z k ) -α -1 k G(y k+1 , z k+1 ) -G(y k , z k ) = 1 2η k y ∥y k -y∥ 2 + 1 2η k z ∥ẑ k -z∥ 2 + 2η k y γ 2 k ∥x k-1 -xk ∥ 2 -⟨x k+1 , y k+1 -y⟩ + γ k ⟨x k-1 -xk , y k -y⟩ -⟨x k -xk+1 , y k+1 -y⟩ + 1 2η k z ∥η k z m k ∥ 2 P -∥η k+1 z m k+1 ∥ 2 P + (1 -α k )α -1 k G(y k , z k ) -G(y, z) -α -1 k G(y k+1 , z k+1 ) -G(y, z) ,where (a) uses the definition of g k y and g k z on line 5 of Algorithm 1 and the convexity and (2r yz )smoothness of the function G(y, z). Further, we divide both sides of the inequality by α k and, using eq. ( 83), obtain the following:1 2η y ∥y k+1 -y∥ 2 + 1 2η z ∥ẑ k+1 -z∥ 2 + 1 2η z ∥η k+1 z m k+1 ∥ 2 P ≤ 1 2η y ∥y k -y∥ 2 + 1 2η z ∥ẑ k -z∥ 2 + 1 2η z ∥η k z m k ∥ 2 P + 2η y α -2 k γ 2 k ∥x k-1 -xk ∥ 2 + γ k α -1 k ⟨x k-1 -xk , y k -y⟩ -α -1 k ⟨x k -xk+1 , y k+1 -y⟩ -α -1 k ⟨x k+1 , y k+1 -y⟩ + (α -2 k -α -1 k ) G(y k , z k ) -G(y, z) -α -2 k G(y k+1 , z k+1 ) -G(y, z)Next, we divide the inequality in Lemma 7 by α k and, using the definition of τ k x and τ x in eqs. ( 83) and ( 84), obtain the following:τ x (α -2 k + α -1 k )∥x k+1 -x∥ 2 ≤ τ x α -2 k ∥x k -x∥ 2 - τ x α -2 k 2 ∥x k+1 -x k ∥ 2 + 2nM 2 τ x T -α -1 k F (x k+1 ) -F (x) -⟨y k+1 , xk+1 -x⟩ .Combining this inequality with the previous upper-bound gives the following:τ x (α -2 k + α -1 k )∥x k+1 -x∥ 2 + 1 2η y ∥y k+1 -y∥ 2 + 1 2η z ∥ẑ k+1 -z∥ 2 + 1 2η z ∥η k+1 z m k+1 ∥ 2 P ≤ τ x α -2 k ∥x k -x∥ 2 + 1 2η y ∥y k -y∥ 2 + 1 2η z ∥ẑ k -z∥ 2 + 1 2η z ∥η k z m k ∥ 2 P + 2nM 2 τ x T - τ x α -2 k 2 ∥x k+1 -x k ∥ 2 + 2η y α -2 k γ 2 k ∥x k-1 -xk ∥ 2 -α -1 k ⟨x k -xk+1 , y k+1 -y⟩ + γ k α -1 k ⟨x k-1 -xk , y k -y⟩ -α -1 k F (x k+1 ) -F (x) + ⟨y k+1 , x⟩ -⟨x k+1 , y⟩ + (α -2 k -α -1 k ) G(y k , z k ) -G(y, z) -α -2 k G(y k+1 , z k+1 ) -G(y, z) (a) = τ x α -2 k ∥x k -x∥ 2 + 1 2η y ∥y k -y∥ 2 + 1 2η z ∥ẑ k -z∥ 2 + 1 2η z ∥η k z m k ∥ 2 P + 2nM 2 τ x T - τ x α -2 k 2 ∥x k+1 -x k ∥ 2 + 2η y α -2 k γ 2 k ∥x k-1 -xk ∥ 2 -α -1 k ⟨x k -xk+1 , y k+1 -y⟩ + γ k α -1 k ⟨x k-1 -xk , y k -y⟩ -α -1 k F (x k+1 ) -⟨x k+1 , y⟩ -G(y, z) + (α -2 k -α -1 k ) G(y k , z k ) + ⟨y k , x⟩ -F (x) -α -2 k G(y k+1 , z k+1 ) + ⟨y k+1 , x⟩ -F (x) (b) = τ x α -2 k ∥x k -x∥ 2 + 1 2η y ∥y k -y∥ 2 + 1 2η z ∥ẑ k -z∥ 2 + 1 2η z ∥η k z m k ∥ 2 P + 2nM 2 τ x T - τ x α -2 k 2 ∥x k+1 -x k ∥ 2 + 2η y α -2 k γ 2 k ∥x k-1 -xk ∥ 2 -α -1 k ⟨x k -xk+1 , y k+1 -y⟩ + γ k α -1 k ⟨x k-1 -xk , y k -y⟩ -(α -2 k -α -1 k )Q(x, y k , z k ) + α -2 k Q(x, y k+1 , z k+1 ) -α -1 k Q(x k+1 , y, z) (c) ≤ τ x α -2 k ∥x k -x∥ 2 + 1 2η y ∥y k -y∥ 2 + 1 2η z ∥ẑ k -z∥ 2 + 1 2η z ∥η k z m k ∥ 2 P + 2nM 2 τ x T - τ x α -2 k 2 ∥x k+1 -x k ∥ 2 + 2η y α -2 k γ 2 k ∥x k-1 -xk ∥ 2 -α -1 k ⟨x k -xk+1 , y k+1 -y⟩ + γ k α -1 k ⟨x k-1 -xk , y k -y⟩ -(α -2 k -α -1 k )Q(x, y k , z k ) + α -2 k Q(x, y k+1 , z k+1 ) -α -2 k Q(x k+1 , y, z) + (α -2 k -α -1 k )Q(x k , y, z) = τ x α -2 k ∥x k -x∥ 2 + 1 2η y ∥y k -y∥ 2 + 1 2η z ∥ẑ k -z∥ 2 + 1 2η z ∥η k z m k ∥ 2 P + 2nM 2 τ x T - τ x α -2 k 2 ∥x k+1 -x k ∥ 2 + 2η y α -2 k γ 2 k ∥x k-1 -xk ∥ 2 -α -1 k ⟨x k -xk+1 , y k+1 -y⟩ + γ k α -1 k ⟨x k-1 -xk , y k -y⟩ + (α -2 k -α -1 k ) Q(x k , y, z) -Q(x, y k , z k ) -α -2 k Q(x k+1 , y, z) -Q(x, y k+1 , z k+1 )where (a) uses the fact that y k+1 = α -1 k y k+1 -(1 -α k )α -1 k y k , which follows from lines 4 and 8 of Algorithm 1; (b) uses the definition of Q(x, y, z) in eq. ( 14); (c) uses line 13 of Algorithm 1 and Assumption 1.Further, let α K = 3/(K + 3). Then from eq. ( 82) it follows that α -2 k + α -1 k ≥ α -2 k+1 , γ k α -1 k = (k + 2)/3 and α -1 k = (k + 3)/3. Hence, we obtain the following:τ x α -2 k+1 ∥x k+1 -x∥ 2 + 1 2η y ∥y k+1 -y∥ 2 + 1 2η z ∥ẑ k+1 -z∥ 2 + 1 2η z ∥η k+1 z m k+1 ∥ 2 P ≤ τ x α -2 k ∥x k -x∥ 2 + 1 2η y ∥y k -y∥ 2 + 1 2η z ∥ẑ k -z∥ 2 + 1 2η z ∥η k z m k ∥ 2 P + 2nM 2 τ x T - τ x (k + 3) 2 18 ∥x k+1 -x k ∥ 2 + 4η y (k + 2) 2 18 ∥x k-1 -xk ∥ 2 - k + 3 3 ⟨x k -xk+1 , y k+1 -y⟩ + k + 2 3 ⟨x k-1 -xk , y k -y⟩ + (α -2 k -α -1 k ) Q(x k , y, z) -Q(x, y k , z k ) -α -2 k Q(x k+1 , y, z) -Q(x, y k+1 , z k+1 ) (a) = τ x α -2 k ∥x k -x∥ 2 + 1 2η y ∥y k -y∥ 2 + 1 2η z ∥ẑ k -z∥ 2 + 1 2η z ∥η k z m k ∥ 2 P + 2nM 2 τ x T - 2η y (k + 3) 2 9 ∥x k+1 -x k ∥ 2 + 2η y (k + 2) 2 9 ∥x k-1 -xk ∥ 2 - k + 3 3 ⟨x k -xk+1 , y k+1 -y⟩ + k + 2 3 ⟨x k-1 -xk , y k -y⟩ + (α -2 k -α -1 k ) Q(x k , y, z) -Q(x, y k , z k ) -α -2 k Q(x k+1 , y, z) -Q(x, y k+1 , z k+1 ) ,where (a) uses eq. ( 84).Next, we sum these inequalities for k = 0, . . . , K -1 and obtain the following:τ x α -2 K ∥x K -x∥ 2 + 1 2η y ∥y K -y∥ 2 + 1 2η z ∥ẑ K -z∥ 2 + 1 2η z ∥η K z m K ∥ 2 P ≤ τ x α -2 0 ∥x 0 -x∥ 2 + 1 2η y ∥y 0 -y∥ 2 + 1 2η z ∥ẑ 0 -z∥ 2 + 1 2η z ∥η 0 z m 0 ∥ 2 P + 2nM 2 K τ x T + 8η y 9 ∥x 0 -x -1 ∥ 2 + 2 3 ⟨x -1 -x0 , y 0 -y⟩ - 2η y (K + 2) 2 9 ∥x K -x K-1 ∥ 2 -1 3 (K + 2)⟨x K-1 -xK , y K -y⟩ + K-1 k=0 (α -2 k -α -1 k ) Q(x k , y, z) -Q(x, y k , z k ) - K-1 k=0 α -2 k Q(x k+1 , y, z) -Q(x, y k+1 , z k+1 ) (a) ≤ τ x α -2 0 ∥x 0 -x∥ 2 + 1 2η y ∥y 0 -y∥ 2 + 1 2η z ∥ẑ 0 -z∥ 2 + 1 2η z ∥η 0 z m 0 ∥ 2 P + 2nM 2 K τ x T - 2η y (K + 2) 2 9 ∥x K -x K-1 ∥ 2 + η y (K + 2) 2 9 ∥x K-1 -xK ∥ 2 + 1 4η y ∥y K -y∥ 2 + K-1 k=0 (α -2 k -α -1 k ) Q(x k , y, z) -Q(x, y k , z k ) - K-1 k=0 α -2 k Q(x k+1 , y, z) -Q(x, y k+1 , z k+1 ) (b) = τ x ∥x 0 -x∥ 2 + 1 2η y ∥y 0 -y∥ 2 + 1 2η z ∥ẑ 0 -z∥ 2 + 1 2η z ∥η 0 z m 0 ∥ 2 P + 2nM 2 K τ x T - η y (K + 2) 2 9 ∥x K -x K-1 ∥ 2 + 1 4η y ∥y K -y∥ 2 -α -2 K-1 Q(x K , y, z) -Q(x, y K , z K ) + K-1 k=1 (α -2 k -α -1 k -α -2 k-1 ) Q(x k , y, z) -Q(x, y k , z k ) (c) = τ x ∥x 0 -x∥ 2 + 1 2η y ∥y 0 -y∥ 2 + 1 2η z ∥ẑ 0 -z∥ 2 + 1 2η z ∥η 0 z m 0 ∥ 2 P + 2nM 2 K τ x T - η y (K + 2) 2 9 ∥x K -x K-1 ∥ 2 + 1 4η y ∥y K -y∥ 2 - K k=1 λ k Q(x k , y, z) -Q(x, y k , z k ) ,where (a) uses the Cauchy-Schwarz inequality, Young's inequality, and the initialization x0 = x -1 on line 1 of Algorithm 1; (b) uses the fact that α 0 = 1, which follows from eq. ( 82); (c) uses the definition of λ k in eq. ( 85). Further, we obtain the following:τ x α -2 K ∥x K -x∥ 2 + 1 2η y ∥y K -y∥ 2 + 1 2η z ∥ẑ K -z∥ 2 + 1 2η z ∥η K z m K ∥ 2 P (a) ≤ τ x ∥x 0 -x∥ 2 + 1 2η y ∥y 0 -y∥ 2 + 1 2η z ∥ẑ 0 -z∥ 2 + 1 2η z ∥η 0 z m 0 ∥ 2 P + 2nM 2 K τ x T - η y (K + 2) 2 9 ∥x K -x K-1 ∥ 2 + 1 4η y ∥y K -y∥ 2 - K k=1 λ k Q(x K a , y, z) -Q(x, y K a , z K a ) (b) = τ x ∥x 0 -x∥ 2 + 1 2η y ∥y 0 -y∥ 2 + 1 2η z ∥ẑ 0 -z∥ 2 + 1 2η z ∥η 0 z m 0 ∥ 2 P + 2nM 2 K τ x T - η y (K + 2) 2 9 ∥x K -x K-1 ∥ 2 + 1 4η y ∥y K -y∥ 2 - K-1 k=0 α -1 k Q(x K a , y, z) -Q(x, y K a , z K a ) ,where (a) uses the convexity of Q(x, y, z) in x (follows from Assumption 1) and the concavity of Q(x, y, z) in (y, z), line 14 of Algorithm 1, and the fact that λ k ≥ 0, which follows from eqs. ( 82) and ( 85); (b) use the definition of λ k in eq. ( 85) and the fact that α 0 = 1, which follows from eq. ( 85). Next, we do rearranging and use eqs. ( 84) and ( 86), which gives the following:K-1 k=0 α -1 k Q(x K a , y, z) -Q(x, y K a , z K a ) ≤ r 3 ∥x∥ 2 + 6 r ∥y∥ 2 + 15χ 2 r ∥z∥ 2 + 12nM 2 K rT .Next, we divide both sides of the inequality by K-1 k=0 α -1 k , which gives the following:Q(x K a , y, z) -Q(x, y K a , z K a ) ≤ K-1 k=0 α -1 k -1 r 3 ∥x∥ 2 + 6 r ∥y∥ 2 + 15χ 2 r ∥z∥ 2 + 12nM 2 K rT .Further, we can estimateK-1 k=0 α -1 k as follows: K-1 k=0 α -1 k (a) = K-1 k=0 k + 3 3 = K + 1 3 K-1 k=0 k = K + K(K -1) 6 = K(K + 5) 6 ≥ K 2 6 ,where (a) uses eq. ( 82). Plugging this into the previous inequality givesQ(x K a , y, z) -Q(x, y K a , z K a ) ≤ 6 K 2 r 3 ∥x∥ 2 + 6 r ∥y∥ 2 + 15χ 2 r ∥z∥ 2 + 12nM 2 K rT = 1 K 2 2r∥x∥ 2 + 36 r ∥y∥ 2 + 90χ 2 r ∥z∥ 2 + 72nM 2 rKTwhich concludes the proof.1. Claims Question: Do the main claims made in the abstract and introduction accurately reflect the paper's contributions and scope? Answer: [Yes] Justification: See the abstract and the introduction (Section 1). Guidelines:• The answer NA means that the abstract and introduction do not include the claims made in the paper. • The abstract and/or introduction should clearly state the claims made, including the contributions made in the paper and important assumptions and limitations. A No or NA answer to this question will not be perceived well by the reviewers. • The claims made should match theoretical and experimental results, and reflect how much the results can be expected to generalize to other settings. • It is fine to include aspirational goals as motivation as long as it is clear that these goals are not attained by the paper. 
experiment statistical significance: 
 Question: Does the paper report error bars suitably and correctly defined or other appropriate information about the statistical significance of the experiments? Answer: [NA] Justification: The paper does not contain experimental results. Guidelines:• The answer NA means that the paper does not include experiments.• The authors should answer "Yes" if the results are accompanied by error bars, confidence intervals, or statistical significance tests, at least for the experiments that support the main claims of the paper. • The factors of variability that the error bars are capturing should be clearly stated (for example, train/test split, initialization, random drawing of some parameter, or overall run with given experimental conditions). • The method for calculating the error bars should be explained (closed form formula, call to a library function, bootstrap, etc.)• The assumptions made should be given (e.g., Normally distributed errors).• It should be clear whether the error bar is the standard deviation or the standard error of the mean. • It is OK to report 1-sigma error bars, but one should state it. The authors should preferably report a 2-sigma error bar than state that they have a 96% CI, if the hypothesis of Normality of errors is not verified. • For asymmetric distributions, the authors should be careful not to show in tables or figures symmetric error bars that would yield results that are out of range (e.g. negative error rates). Guidelines:• The answer NA means that the authors have not reviewed the NeurIPS Code of Ethics.• If the authors answer No, they should explain the special circumstances that require a deviation from the Code of Ethics. • The authors should make sure to preserve anonymity (e.g., if there is a special consideration due to laws or regulations in their jurisdiction). 10. Broader Impacts Question: Does the paper discuss both potential positive societal impacts and negative societal impacts of the work performed? Answer: [NA] Justification: The paper provides theoretical research and there is no societal impact from the work performed. Guidelines:• The answer NA means that there is no societal impact of the work performed.• If the authors answer NA or No, they should explain why their work has no societal impact or why the paper does not address societal impact. • Examples of negative societal impacts include potential malicious or unintended uses (e.g., disinformation, generating fake profiles, surveillance), fairness considerations (e.g., deployment of technologies that could make decisions that unfairly impact specific groups), privacy considerations, and security considerations.• The conference expects that many papers will be foundational research and not tied to particular applications, let alone deployments. However, if there is a direct path to any negative applications, the authors should point it out. For example, it is legitimate to point out that an improvement in the quality of generative models could be used to generate deepfakes for disinformation. On the other hand, it is not needed to point out that a generic algorithm for optimizing neural networks could enable people to train models that generate Deepfakes faster. • The authors should consider possible harms that could arise when the technology is being used as intended and functioning correctly, harms that could arise when the technology is being used as intended but gives incorrect results, and harms following from (intentional or unintentional) misuse of the technology. • If there are negative societal impacts, the authors could also discuss possible mitigation strategies (e.g., gated release of models, providing defenses in addition to attacks, mechanisms for monitoring misuse, mechanisms to monitor how a system learns from feedback over time, improving the efficiency and accessibility of ML). 
safeguards: 
 Question: Does the paper describe safeguards that have been put in place for responsible release of data or models that have a high risk for misuse (e.g., pretrained language models, image generators, or scraped datasets)?Answer: [NA] Justification: The paper does not pose the risks described in the question.Guidelines:• The answer NA means that the paper poses no such risks.• Released models that have a high risk for misuse or dual-use should be released with necessary safeguards to allow for controlled use of the model, for example by requiring that users adhere to usage guidelines or restrictions to access the model or implementing safety filters. • Datasets that have been scraped from the Internet could pose safety risks. The authors should describe how they avoided releasing unsafe images. • We recognize that providing effective safeguards is challenging, and many papers do not require this, but we encourage authors to take this into account and make a best faith effort.12. Licenses for existing assets Question: Are the creators or original owners of assets (e.g., code, data, models), used in the paper, properly credited and are the license and terms of use explicitly mentioned and properly respected?Answer: [NA] Justification: The paper does not use existing assets.Guidelines:• The answer NA means that the paper does not use existing assets.• The authors should cite the original paper that produced the code package or dataset.• The authors should state which version of the asset is used and, if possible, include a URL. • The name of the license (e.g., CC-BY 4.0) should be included for each asset.• For scraped data from a particular source (e.g., website), the copyright and terms of service of that source should be provided. • If assets are released, the license, copyright information, and terms of use in the package should be provided. For popular datasets, paperswithcode.com/datasets has curated licenses for some datasets. Their licensing guide can help determine the license of a dataset. • For existing datasets that are re-packaged, both the original license and the license of the derived asset (if it has changed) should be provided.By complexity, we mean, depending on the context, the number of subgradient computations or decentralized communications required to solve the problem.A function is called smooth if it is continuously differentiable and has a Lipschitz-continuous gradient.By 2 V×V = {E : E ⊂ V × V} we denote the set of all subsets of V × V.A vector x ∈ R d is called an ϵ-approximate solution to problem (1) if p(x) -p(x * ) ≤ ϵ.Scaman et al. (2018) do not provide any lower complexity bounds in the strongly convex setting. However, the desired lower bound on the communication complexity can be obtained by extending their analysis.We ignore universal constants in Table3like in the O(•) and Ω(•) notation. 
 
 Appendix:  


Decision: This paper derives lower bounds on the communication and computation complexities of solving non-smooth convex decentralized optimization problems over time-varying networks and designs an algorithm that achieves these lower bounds. All reviewers appreciate the contributions of this work, recognizing its novelty and nontrivial nature. The only minor weakness lies in the somewhat limited simulation setups.

In the final version I recommend the authors to discuss the time-varying assumption. In general, time-varying graphs are assumed to be contractive over a finite length, but this work seem to assume contraction for every graph.

Reviews: review: paper_summary: paper_summary: This paper derives lower bounds on the communication and computation complexities of solving non-smooth convex decentralized optimization problems over time-varying networks and designs and algorithm that achieves these lower bounds.
strength_weakness: strength_weakness: None.
questions: questions: None.
limitations: limitations: Yes.
paper_summary: paper_summary: The paper studies non-smooth decentralized optimization with time-varying communication networks. The paper presents execution time lower bounds of subgradient algorithms for strongly-convex and convex cases. Then, the paper develops algorithms achieving matching execution time.
strength_weakness: strength_weakness: There is no major flaw I see in this work, but again I am not very familiar with the field. One minor weakness might be a lack of numerical experiments validating the convergence rate of the proposed algorithm. Nonetheless, I think it is not a big issue for a theory paper.
questions: questions: Why does the objective in eq (1) have a quadratic regularization?

The theorem 1 and 2 state lower bounds of execution times but the authors interpret these results as communication and computation bounds. I wonder can you make this argument formal in the classical definition of communication complexity (Yao, 1979) and computation complexity (Papadimitriou, 2003).

References:

Yao, A. C. C. (1979). Some complexity questions related to distributive computing (preliminary report). In Proceedings of the eleventh annual ACM symposium on Theory of computing (pp. 209-213).

Papadimitriou, C. H. (2003). Computational complexity. In Encyclopedia of computer science (pp. 260-265).
limitations: limitations: Yes.
paper_summary: paper_summary: The authors introduce an algorithm which optimally bounds the complexity of algorithms for non-smooth convex decentralised optimisation over time varying networks.
strength_weakness: strength_weakness: The paper can be quite dense and while the authors provide slight intuition about the proof sketch in the main paper, most of the actual paper lies in the appendix which can make it quite inaccessible.
questions: questions: Based on the above comments, I would ask the authors the following question,

 - Would it be possible to make the main paper less dense and focus more on guiding the user through the intuition of the proofs?
limitations: limitations: The authors have sufficiently addressed any limitations that might arise from the submitted work.
paper_summary: paper_summary: This paper studies the problem of decentralized optimization for non-smooth convex opjectives and time-varying networks. The paper introduces an algorithm to solve this problem together with matching lower bounds on the required communication and subgradient computations, thereby proving that the proposed algorithm is optimal in this sense.
strength_weakness: strength_weakness: In my experience, optimization algorithms can sometimes yield unsatisfactory performance despite having good theoretical guarantees. Therefore, I believe that adding (even a small) simulation example showcasing the proposed algorithm's performance on a relevant problem would greatly improve the paper. 

Furthermore, while the paper is generally very well written, the algorithm could be explained a little bit better. Even just mentioning which lines of the algorithm correspond to which step in lines 282 - 296 would greatly help a reader who is not necessarily familiar with each of the references. In exchange, I believe that the justification for considering convex cost functions (Section 1.2) could be shortened. 

Finally, there are a few small typos that should be corrected:
I think that $\mathbf{1}_p$ should be defined as $(1,\dots,1)^\top$; Algorithm 1 requires an initialization for $\bar{z}^0$ and $\bar{y}^0$, and in the proofs (e.g., after line 765) you need to define w.l.o.g. that $\tilde{x}^0 = x^{-1}$; in line 657, I think you should have $+\langle x_a^K,z \rangle$ and you need to flip the sign in (87); in lines 692-693, the function $F(x)$ maps from $(\mathbb{R}^n)^d$, i.e., the argument should be $w^*$ instead of $x^*$; in Appendix E.3, E.4, and E:5, you previously used the notation $W(k)$ instead of $W_k$; in the first equations of Appendix E.6 (after (b)), I believe it should be $y^k$ instead of $y^{k+1}$, which is corrected in the next equation, and I do not see why you need line 7 of Algorithm 1 for (d); after line 773, there should be $k$ instead of $K$ in the third sum.
questions: questions: Why does (57) not change in Part(ii) (line 591)? I believe the result is correct, but I am not sure whether this statement is true.
limitations: limitations: I do not fully agree with the authors' answer to Question 2 in the NeurIPS Checklist, as limitations of the presented results and avenues for future research are not clearly discussed in the paper. More specifically, the authors only provide the required assumptions, but do not discuss whether these assumptions are restrictive in certain scenarios (except for Assumption 1) and whether these assumptions could be relaxed in future work.

Comments: title: Read your rebuttal
comment: I thank the authors for their clarifications. Based on all the information presented I will be keeping my score as-is. Best of luck to the authors for the final evaluation.
Thank you for your clarifications. I read them carefully and decided to keep my original score, which is already relatively high.
Thank you for the clarification, I decide to keep my positive rating.