\(\mathtt{BAGEL}\): Projection-Free Algorithm for Adversarially
Constrained Online Convex Optimization
February 23, 2025
Projection-based algorithms for Constrained Online Convex Optimization (COCO) achieve optimal \(\mathcal{O}(T^{1/2})\) regret guarantees but face scalability challenges due to the computational complexity of projections. To circumvent this, projection-free methods utilizing Linear Optimization Oracles (LOO) have been proposed, albeit typically achieving slower \(\mathcal{O}(T^{3/4})\) regret rates. In this work, we examine whether the \(\mathcal{O}(T^{1/2})\) rate can be recovered in the projection-free setting by strengthening the oracle assumption. We introduce BAGEL, an algorithm utilizing a Separation Oracle (SO) that achieves \(\mathcal{O}(T^{1/2})\) regret and \(\tilde{\mathcal{O}}(T^{1/2})\) cumulative constraint violation (CCV) for convex cost functions. Our analysis shows that by leveraging an infeasible projection via SO, we can match the time-horizon dependence of projection-based methods with \(\tilde{\mathcal{O}}(T)\) oracle calls, provided dependence on the geometry of the action set. This establishes a specific regime where projection-free methods can attain the same convergence rates as projection-based counterparts.
Online Convex Optimization (OCO) has emerged as a foundational framework for sequential decision-making under uncertainty, with applications spanning resource allocation, real-time control systems, and adversarial robustness in machine learning [1]. A critical variant of this framework, Constrained Online Convex Optimization (COCO), requires agents to minimize cumulative costs while adhering to time-varying constraints. In this problem, the learner interacts with an environment over \(T\) rounds, where at round \(t\), the learner selects a action \(\mathbf{x}_t \in \mathcal{K}\), where \(\mathcal{K}\in\mathbb{R}^n\) is the action set. The adversary then reveals a convex cost function \(f_t: \mathcal{K}\to \mathbb{R}\) and \(k\) convex constraint functions \(g_{t,i}: \mathcal{K}\to \mathbb{R},\;i \in [k]\). Let \(\mathcal{K}^\star \subseteq \mathcal{X}\) denote the set of feasible actions satisfying all constraints \(g_{t,i}(\mathbf{x}) \leq 0\) for \(t \in [T],\;i \in [k]\). We assume \(\mathcal{K}^\star \neq \emptyset\) to ensure that at least one feasible solution exists. Then, we use the following performance metrics: \[\begin{align} \mathtt{Regret}(T) &= \sum\limits_{t=1}^Tf_t(\mathbf{x}_t) - \min_{\mathbf{x}^\star \in \mathcal{K}} \sum_{t=1}^T f_t(\mathbf{x}^\star), \\ \mathtt{CCV}(T) &= \max_{i \in [k]} \sum\limits_{t=1}^T\left(g_{t,i}(\mathbf{x}_t)\right)^+. \end{align}\] The goal is to design an algorithm that achieves sublinear \(\mathtt{Regret}(T)\) and minimize \(\mathtt{CCV}(T)\) simultaneously. COCO’s applications span diverse domains, including portfolio optimization under dynamic risk constraints [2], dynamic resource allocation for time-varying workloads [3], real-time dynamic pricing with inventory constraints [4], and collision-free trajectory planning under safety constraints [5].
A central challenge in online optimization is the design of computationally efficient, projection-free algorithms for the adversarial Constrained Online Convex Optimization (COCO) setting. We summarize the key related works in Table ¿tbl:tbl:comparison?1. We note that most works require solving a constrained convex optimization problem on each round (COO), or requires Euclidean projection operation on the convex set (PO). While projection-based algorithms, such as that of [6], establish the \(\tilde{O}(T^{1/2})\) bounds for both regret and CCV, this performance comes at the cost of a potentially intractable projection step in each round, limiting their applicability to domains like the spectahedron or flow polytopes. To circumvent this computational hurdle, current state-of-the-art projection-free methods relying on LOOs, such as [7] and [8], achieve regret rates of \(\mathcal{O}(T^{3/4})\) and CCV of \(\mathcal{O}(T^{7/8})\) or \(\tilde{\mathcal{O}}(T^{3/4})\). This gap with respect to horizon \(T\) raises a theoretical question: can projection-free methods match the horizon-dependent performance of projection-based methods in adversarial setting?
| Reference | Regret | CCV | \(\beta\) Range | Projection | |||
| type | |||||||
| to Oracle | Cost Function | ||||||
| [9] | \(\C{O}(T^{1/2})\) | \(\C{O}(T^{3/4})\) | - | Based | COO | - | Convex |
| \(\C{O}(\log(T))\) | \(\C{O}(T^{1/2} \sqrt{\log(T)})\) | - | \(\C{O}(T)\) | Strongly Convex | |||
| [10] | \(\C{O}(T^{\max(\beta,1-\beta)})\) | \(\C{O}(T^{1-\beta/2})\) | (0,1) | Based | COO | \(\C{O}(T)\) | Convex |
| \(\C{O}(\log(T))\) | \(\C{O}(T^{1/2} \sqrt{\log(T)})\) | - | \(\C{O}(T)\) | Strongly Convex | |||
| [6] | \(\C{O}(T^{1/2})\) | \(\C{O}(T^{1/2}\log(T))\) | - | Based | PO | \(\C{O}(T)\) | Convex |
| \(\C{O}(\log(T))\) | \(\C{O}(T^{1/2} \sqrt{\log(T)})\) | - | \(\C{O}(T)\) | Strongly Convex | |||
| [7] (Thm 4.1) | \(\C{O}(T^{3/4})\) | \(\C{O}(T^{7/8})\) | - | Free | LOO | \(\C{O}(T)\) | Convex |
| [8] (Thm 3.4) | \(\C{O}(T^{3/4})\) | \(\C{O}(T^{3/4})\) | - | Free | LOO | \(\C{O}(T)\) | Convex |
| This paper | \(\C{O}({T}^{1-\beta})\) | \(\C{O}({T}^{1-\beta}\log(T))\) | (0,1/2] | Free | SO | \(\tilde{\C{O}}(T^{2\beta})\) | Convex |
| \(\C{O}(T^{1-\beta}\log(T))\) | \(\C{O}(T^{1-\beta/2}\sqrt{\log(T)})\) | (0,1] | SO | \(\tilde{\C{O}}(T^{2\beta})\) | Strongly Convex |
In this paper, we demonstrate that this gap with respect to \(T\) can be closed when a Separation Oracle (SO) is available. While LOOs are computationally advantageous for some geometries like the nuclear norm ball, Separation Oracles can be more efficient for some others (e.g., spectral norm balls) and offer a natural alternative for accessing the action set. A detailed discussion on LOO and SO is provided in Section 3. We present \(\mathtt{BAGEL}\) (Blocked Adaptive online Gradient descEnt with infeasibLe projection), which achieves \(\mathcal{O}(T^{1/2})\) regret and \(\tilde{\mathcal{O}}(T^{1/2})\) CCV bounds in the projection-free setting for adversarial COCO. Our total calls to SO depends on geometric constants (specifically the ratio of enclosing to enclosed balls), that may scale with the ambient dimension.
In Section 4, we develop a projection-free algorithm for OCO based on Online Gradient Descent with adaptive step sizes and establish its regret guarantees. While adaptive step sizes do not improve the regret bounds for OCO, they play a crucial role in proving the guarantees for the COCO problem. In Section 5, we introduce an approach that constructs a surrogate function and applies base algorithm for OCO to it. Inspired by [6], we adopt a Lyapunov-like function to include cumulative constraint violation in weights when combining cost and constraint via a surrogate function. Furthermore, we extend our analysis to strongly convex cost functions in Appendix 10 and 12. The key contributions of this work are summarized as follows.
We propose \(\mathtt{BAGEL}\), a projection-free algorithm for Adversarial Constrained Online Convex Optimization (COCO), that achieves \(\tilde{\mathcal{O}}(T^{1/2})\) regret and CCV for convex optimization while maintaining in total \(\tilde{\mathcal{O}}(T)\) calls to the separation oracle (Theorem 2). Our results match those of the projection-based counterparts on time horizon \(T\), improving upon the \(T\)-dependence of recent LOO-based works (which typically scale as \(T^{3/4}\)), albeit dependent on the action set geometry and the availability of Separation Oracle.
We are the first to analyze projection-free OGD for COCO with a specific adaptive step-size of the form \(\eta_m \propto 1/\sqrt{\epsilon+\sum_{\tau=1}^m\|\bar{\nabla}_\tau \|^2}\). We show this seemingly minor modification is critical for computational feasibility. As detailed in Remark 2, this step-size reduces the total number of calls to the separation oracle from a prohibitive \(\mathcal{O}(T^2)\) to a near-linear \(\tilde{\mathcal{O}}(T)\), making \(\tilde{\mathcal{O}}(T^{1/2})\) bounds of regret and CCV practically achievable.
We adopt blocking technique and use a trade-off parameter \(\beta\) to allow adjustment of algorithm to computational limitations. The tradeoff is given in Table ¿tbl:tbl:comparison?, where we see that with decreased \(\beta\), we would need less calls to the separation oracle at the cost of higher regret and CCV.
We extend \(\mathtt{BAGEL}\) to handle strongly convex cost functions, achieving \(\mathcal{O}(\log T)\) regret and \(\tilde{\mathcal{O}}(T^{1/2})\) CCV (Theorem 4), matching the performance of projection-based algorithms when allowed more separation oracle calls. This extension demonstrates that \(\mathtt{BAGEL}\) is effective in a wider range of optimization settings.
In order to obtain the above results, we also study projection-free Online Convex Optimization (OCO) with adaptive step sizes (Theorem 1 for convex and Theorem 3 for strongly-convex). We note that previous work considered constant step-sizes [11], while the results for COCO needed such an analysis with adaptive step-sizes.
Adversarial COCO COCO problems have been widely studied, since they can handle complex action set [12]. Table ¿tbl:tbl:comparison? provides some of the key works for regret and CCV guarantees in the area. [9] proposed an algorithm that solves a constrained convex optimization problem in each round, and achieves \(\mathcal{O}(T^{1/2})\) regret with \(\mathcal{O}({T}^{3/4})\) CCV. This result is a special case of the result in [10]. Recently, [6] proposed a projection-based algorithm that achieves optimal regret and CCV of \(\tilde{\mathcal{O}}(T^{1/2})\). The results have been further extended to strongly convex cost functions, where all these works show a regret of \(\tilde{\mathcal{O}}(1)\) with a CCV of \(\tilde{\mathcal{O}}(T^{1/2})\). However, performing a constrained convex optimization or doing a orthogonal projection is computationally complex [11], [13], and thus we provide a projection-free algorithm in this paper.
Projection-free Algorithms for OCO Many projection-free algorithms have been proposed for classical OCO problems. However, regrets for Frank-Wolfe based approaches have been shown to achieve only \(\mathcal{O}(T^{3/4})\) for convex cost function [14], [15] and \(\mathcal{O}(T^{2/3})\) for strongly convex cost functions [16]. Even though [17] proves that on strongly convex set they can further improve the regrets to \(\mathcal{O}(T^{2/3})\) for convex cost function, these results are still not matching projection-based results. [18] proposes the concept of infeasible projection to replace the traditional orthogonal projection oracles, although it is only implemented for online linear optimization. In parallel with [19], [11] gives an approach to implement infeasible projection via separation oracle for OCO, which achieves an optimal regret of \(\mathcal{O}(T^{1/2})\) using online gradient descent based approach with constant step size. In their work, implementing infeasible projection via linear optimization oracle yields \(\mathcal{O}(T^{3/4})\) regret, which motivates the choice of separation oracle in our work. Our work builds upon the concept of infeasible projection, implemented via a separation oracle, pioneered by [11] for the unconstrained OCO setting. We specifically adopt the variant from [20] due to its favorable property of ensuring iterates remain within the constraint set \(\mathcal{K}\), which is crucial for our analysis. Our primary novelty lies in being the first to integrate this machinery with a carefully constructed adaptive step-size scheme and a surrogate framework inspired by [6] to solve the more complex constrained COCO problem, and in particular, to do so with near-optimal oracle complexity.
Projection-free Algorithms for Adversarial COCO Projection-free algorithms for online learning were proposed in response to the growing prevalence of high-dimensional data, which led to increased computational challenges for projection operators [11], [13], [15], [21], [22]. The primary challenge in projection-free COCO is to achieve strong guarantees using only computationally cheap oracles. As summarized in Table ¿tbl:tbl:comparison?, recent approaches based on Linear Optimization Oracles (LOOs) have made significant progress. [7] provides an algorithm with \(\mathcal{O}(T^{3/4})\) regret and \(\mathcal{O}(T^{7/8})\) CCV. Concurrently, [8] achieves \(\tilde{\mathcal{O}}(T^{3/4})\) for both metrics. While impressive, these results remain slower compared to projection-based methods. Our work differs by employing a Separation Oracle (SO), which is another natural projection-free setting in which the feasible action set \(\mathcal{K}\) can be given [11]. We demonstrate that this increase in oracle power is sufficient to close the performance gap entirely for convex functions.
We formally describe the problem as follows. Consider an agent is playing a consecutive online game with an adversary. The time horizon is \(T\) if the game stops after \(T\) steps. At each step \(t\), the agent chooses an action \(\mathbf{x}_t\), and the adversary chooses a cost function \(f_t\) and \(k\) constraint functions \(g_{t,i}\), for \(i \in [k]\), and reveals them to the agent. Note that the objective of the game is to minimize the cost \(f_t(\mathbf{x}_t)\) while satisfying the constraints \(g_{t,i}(\mathbf{x}_t)\leq 0, \forall i\in[k]\).
Preliminaries and Notations We say a function \(f\) over the action set \(\mathcal{K}\) is convex if for all \(\mathbf{x},\mathbf{y}\in \mathcal{K}\), \(f(\mathbf{y}) \geq f(\mathbf{x}) + \nabla f(\mathbf{x})^T(\mathbf{y}-\mathbf{x})\). We say a function \(f\) over the action set \(\mathcal{K}\) is \(\theta\)-strongly convex if for all \(\mathbf{x},\mathbf{y}\in \mathcal{K}\), \(f(\mathbf{y}) \geq f(\mathbf{x}) + \nabla f(\mathbf{x})(\mathbf{y}-\mathbf{x})+\frac{\theta}{2}\Vert \mathbf{y}-\mathbf{x}\Vert^2\). For a set \(\mathcal{K}\subseteq \mathbb{R}^d\), we define the set of \(\alpha\mathbf{x}+(1-\alpha)\mathbf{y}\) for all \(\mathbf{x},\mathbf{y}\in\mathcal{K}\) and \(\alpha\in\mathbb{R}\) as the affine hull of \(\mathcal{K}\), \(\operatorname{aff}(\mathcal{K})\). The diameter of the action set is defined as \(D \triangleq \max_{\mathbf{x},\mathbf{y}\in \mathcal{K}}\Vert \mathbf{x}-\mathbf{y}\Vert\). The relative interior of \(\mathcal{K}\), denoted by \(\operatorname{relint}(\mathcal{K})\), is defined as the set \(\{\mathbf{x}\in\mathcal{K}| \exists r > 0, \mathbb{B}_r(\mathbf{x}) \cap \operatorname{aff}(\mathcal{K}) \subseteq \mathcal{K}\}\). By definition, for any point \(\mathbf{c} \in \operatorname{relint}(\mathcal{K})\), there exists a real number \(r > 0\) such that \(\operatorname{aff}(\mathcal{K})\cap \mathbb{B}_r(\mathbf{c}) \subseteq \mathcal{K}\), where \(\mathbb{B}_r(\mathbf{c})\) represent the ball of radius \(r\) centering at \(\mathbf{c}\). We use \(r_{\mathcal{K}, \mathbf{c}}\) to denote the largest \(r\) satisfying this property. Given \(\mathbf{c}\) and \(r\), for any \(\delta \in [0,1)\), referred to as the shrinking parameter, we define the shrunk constrained set \(\mathcal{K}_\delta\) as the set \((1-\delta)\mathcal{K}+ \delta\mathbf{c}\).
Oracle Type In classical projection-based algorithms, Projection Oracle (PO/\(\mathcal{O}_{P}\)) and Convex Optimization Oracle (COO/\(\mathcal{O}_{CO}\)) are standard feasibility mechanisms for actions chosen by the algorithm. Given a closed convex set \(\mathcal{K} \subseteq \mathbb{R}^n\) and a query point \(y \in \mathbb{R}^n\), a PO returns the unique point \(\mathbf{x}\in \mathcal{K}\) that minimizes the Euclidean distance to \(\mathbf{y}\), i.e., \(\mathcal{O}_{P}(\mathbf{y}, \mathcal{K}) = \arg\min_{\mathbf{x}\in \mathcal{K}} \| \mathbf{y}- \mathbf{x}\|_2^2\). Given a general convex objective function \(h: \mathbb{R}^n \rightarrow \mathbb{R}\) over \(\mathcal{K}\), a COO returns a global minimizer, i.e., \(\mathcal{O}_{CO}(h, \mathcal{K}) \in \arg\min_{\mathbf{x}\in \mathcal{K}} \;h(\mathbf{x})\). In the context of COCO, COOs are often invoked to handle time-varying soft constraints \(g_t(\mathbf{x}) \le 0\) without performing a direct projection. If \(h(\mathbf{x})\) is quadratic, a COO acts as a PO. If \(h(\mathbf{x})\) includes complex non-linear constraints, a COO represents a more computationally intensive subroutine than a simple projection.
The Linear Optimization Oracle (LOO/\(\mathcal{O}_{LO}\)) and Separation Oracle (SO/\(\mathcal{O}_{S}\)) are established approaches to move from projection-based algorithm to projection-free paradigms. The LOO exploits the fact that for many complex geometries, minimizing a linear function is vastly cheaper than minimizing a quadratic function (PO) or a general convex function (COO). Given a convex set \(\mathcal{K}\subseteq \mathbb{R}^n\) and a linear objective vector \(\mathbf{c} \in \mathbb{R}^n\) (typically a gradient or a direction vector), a LOO returns a vertex or point in \(\mathcal{K}\) that minimizes the linear function: \(\mathcal{O}_{LO}(\mathbf{c}, \mathcal{K}) \in \arg\min_{\mathbf{x}\in \mathcal{K}} \;\langle \mathbf{c}, \mathbf{x}\rangle\). The SO is a more recent development, and another natural projection-free setting to access the action set. Given a convex set \(\mathcal{K} \subseteq \mathbb{R}^n\) and a query point \(\mathbf{y}\in \mathbb{R}^n\), it either asserts that \(\mathbf{y}\in \mathcal{K}\), or it returns a separating hyperplane vector \(\mathbf{g} \in \mathbb{R}^n\) such that for all \(\mathbf{x}\in \mathcal{K}\), \(\langle \mathbf{g}, \mathbf{y}- \mathbf{x}\rangle > 0\). Although LOOs are more prevalent in the literature, it is worth noting the complementing nature between the two oracles in efficiency. [11] gave the example of two unit balls: a nuclear norm ball \(\mathcal{B}_* = \{ \mathbf{x} \in \mathbb{R}^{m \times n} \mid \|\mathbf{x}\|_* \leq 1 \}\), and a spectral norm ball \(\mathcal{B}_2 = \{ \mathbf{x} \in \mathbb{R}^{m \times n} \mid \|\mathbf{x}\|_2 \leq 1 \}\). While the requirement of a full Singular Value Decomposition (SVD) makes Euclidean projection computationally prohibitive onto either set, the efficiency of projection-free oracles differs. The LOO is efficient for \(\mathcal{B}_*\) (requiring only a rank-one SVD) but expensive for \(\mathcal{B}_2\) where a full SVD is needed [23]. In contrast, the SO is expensive for \(\mathcal{B}_*\) due to the subgradient for the nuclear norm (full SVD), but efficient for \(\mathcal{B}_2\) since computing a subgradient for \(g_2(\mathbf{X}) := \|\mathbf{X}\|_2 - 1\) requires only the top singular vector pair (rank-one SVD).
Regrets and CCV For an online algorithm that solves the constrained optimization, there are two key metrics, regret and cumulative constraint violation (CCV). The regret represents the distance in the cumulative cost from the actions chosen by the agent to the optimal action in hindsight. Mathematically, we define the regret to be \(\mathtt{Regret}(T) \triangleq \sum_{t=1}^T f_t(\mathbf{x}_t) - \sum_{t=1}^T f_t(\mathbf{x}^*)\), where \(\mathbf{x}^* \triangleq \operatorname{argmax}_{\mathbf{x}\in\mathcal{K}}\sum_{t=1}^T f_t(\mathbf{x})\). The cumulative constraint violation (CCV) represents the cumulative violation to the constraint, and when there are more than one constraints we choose the one that is violated the most cumulatively. We define CCV to be \(\mathtt{CCV}(T) \triangleq \max_{i\in[k]} \sum_{t=1}^T (g_{t,i})^+\), where \((z)^+ \triangleq \max\{z, 0\}\). In each step, we can simplify the constraint by letting \(g_t(x)=\max_{i\in [k]}g_{t,i}(x)\). Note that if each constraint condition is satisfied, constraint \(g_t\) naturally follows. Without loss of generality, we can assume that in each step, the adversary only reveals one constraint function, namely \(g_t\).
Assumptions Under the given problem setting, we introduce a few assumptions that are fundamental to our proposed algorithm. The first assumption is about the action set \(\mathcal{K}\) that we can choose actions from.
Assumption 1. The action set \(\mathcal{K}\in \mathbb{R}^d\) is convex, closed and bounded with diameter \(D\), i.e., \(\max_{\mathbf{x},\mathbf{y}\in \mathcal{K}}\Vert \mathbf{x}-\mathbf{y}\Vert \leq D\). For simplicity, throughout this work, we will assume that \(\mathcal{K}\), \(\mathbf{c}\) and \(r\) are fixed, with \(\mathbf{c}\) and \(r\) defined in Preliminaries.
The second assumption concerns the cost and constraint functions at every time step.
Assumption 2. At each step \(t\), we assume the cost function \(f_t\) are either convex or strongly convex, and constraint function \(g_t\) are convex. Functions \(f_t\) and \(g_t\) are also assumed to be differentiable with their gradients bounded by \(M_1\), thus also \(M_1\)-Lipschitz continuous.
Besides assumption on action set and cost and constraint functions, we highlight a feasibility assumption that ensures the problem is solvable in Assumption 3.
Assumption 3. We say an action \(\mathbf{x}' \in \mathcal{K}\) is feasible if all constraint function values are non-positive, i.e., \(g_{t,i}(\mathbf{x}') \leq 0, \forall t,i\). We assume the feasible set \(\mathcal{K}'\), the set of all feasible actions, is not empty, i.e., \(\mathcal{K}'\neq \varnothing\).
Infeasible Projection via Separation Oracle To replace costly projection steps, we use the Infeasible Projection via Separation Oracle (IP-SO) to ensure feasibility of actions generated by the algorithm. Unlike a standard Euclidean projection which finds the closest point in \(\mathcal{K}\), the IP-SO finds a point \(\mathbf{y}'\) (not necessarily in \(\mathcal{K}\)) that is closer to all points in a slightly shrunk set \(\mathcal{K}_\delta\) than the original point \(\mathbf{y}\) was. The key is that this can be implemented efficiently: if an iterate is outside \(\mathcal{K}\), the SO provides a separating hyperplane, and we take a small step in a direction normal to that plane. Lemma 3 guarantees that this iterative process terminates quickly and yields a point with the desired properties, effectively mimicking a projection without the high computational cost.
Next, we formally introduce the concept of infeasible projection oracle. We say \(\mathbf{y}' \in \mathbb{R}^n\) is an infeasible projection of \(\mathbf{y}\) onto \(\mathcal{K}\) if for all \(\mathbf{x}\in \mathcal{K}, \Vert \mathbf{y}'-\mathbf{x}\Vert \leq \Vert \mathbf{y}-\mathbf{x}\Vert\). Note that, unlike projection, \(\mathbf{y}'\) is not assumed to be in \(\mathcal{K}\). We say an algorithm returns an infeasible projection for set \(\mathcal{K}\) if it takes input point \(\mathbf{y}\) and return \(\mathbf{y}'\) that is an infeasible projection of \(\mathbf{y}\). An infeasible projection could be implemented using a separation oracle [11]. For a convex set \(\mathcal{K}\subseteq \mathbb{R}^n\), a Separation Oracle (\(SO_\mathcal{K}\)) is an oracle that, given a point \(\mathbf{x}\in \mathbb{R}^n\), either confirms that \(\mathbf{x}\in \mathcal{K}\) or else gives a hyperplane separating \(\mathbf{x}\) from \(\mathcal{K}\). Such a separating hyperplane could be represented as a vector \(\mathbf{g}\) where if \(\mathbf{y}\notin \mathcal{K}\), then \(\langle \mathbf{y}- \mathbf{x}, \mathbf{g}\rangle > 0\) for all \(\mathbf{x}\in \mathcal{K}\).
Here we use a variant of the IP-SO algorithm described in [20], and included in Appendix 8 as Algorithm 1 for completeness. We use \(\mathcal{P}_{\mathcal{K}, \mathbf{c}, r, \delta}\) to denote this algorithm. This variant slightly differs from the original version described in [11] in that it does not require the assumption that \(\mathcal{K}\) contains a \(d\)-dimensional ball, neither does it require \(\boldsymbol{0} \in \mathcal{K}\). The key idea of the algorithm is to provide an infeasible projection onto \(\mathcal{K}_\delta\) that is contained within \(\mathcal{K}\). This guarantees that the infeasible projection is not too far from the projection, since any point in \(\mathcal{K}\) is not too far from \(\mathcal{K}_\delta\). Lemma 3 (given in Appendix 8) ensures that the algorithm output is an infeasible projection within certain number of iterations (or queries of separation oracle). When there is no ambiguity, we simply use \(\mathcal{P}\) and drop the subscripts.
In order to extend OCO methods to the COCO framework inspired by [6], we use a surrogate function that combines the cost and constraint, leveraging the cumulative constraint violation as weights. A key challenge is that the gradients can have vastly different magnitudes across iterations, as shown in Section 5. When constraint violations accumulate, the weight can become very large, leading to explosive gradients. An adaptive step-size dampens these large gradients by incorporating gradient magnitudes into the denominator, providing the stability required for the surrogate-based analysis to hold. Therefore, we first formally analyze a projection-free OGD with such an adaptive step-size (Algorithm [alg:general]), a necessary building block for the algorithm under adversarial constraints.
The proposed base algorithm is based on Online Gradient Descent with adaptive step-size. Although [11] proposed a projection free OGD using an infeasible projection oracle (Algorithm 1, Lemma 4), it only considered constant step size and did not consider feasibility constraints on the action set \(\mathcal{K}\). However, the guarantees of results for COCO algorithm require adaptive step-sizes and feasible solutions returned by infeasible projection oracle, which is why we propose and analyze the Base Algorithm described in Algorithm [alg:general]. Note that we add a constraint constant \(\epsilon\) so that the regret are guaranteed with only \(\tilde{\mathcal{O}}(T)\) total calls to the separation oracle.
Algorithm [alg:general] takes as input: action set \(\mathcal{K}\), time horizon \(T\), \(\mathbf{c} \in \operatorname{relint}(\mathcal{K})\), \(r = r_{\mathcal{K}, \mathbf{c}}\), block size \(K\), shrinking parameter \(\delta\in[0,1)\), diameter of the action set \(D\), an instance of IP-SO algorithm \(\mathcal{P}\) implemented by Algorithm 1, a constraint constant \(\epsilon>0\) for convex cost functions or the strongly-convex parameter \(\theta\) for strongly-convex cost functions. We initiate the program by choosing \(\mathbf{x}_1\in\mathcal{K}\), and pass \(\mathcal{K}, \mathbf{c}, r, \delta\) to \(\mathcal{P}\) to initiate algorithm instances. Within each block \(m\), there are \(K\) iterations, represented by \(\mathcal{T}_m\). At each step \(t\), we play action \(\mathbf{x}_t\), and observe the cost function \(f_t\) adversarially selected and revealed by the adversary, which we use to compute \(\nabla_t\), the gradient of \(f_t\) at \(\mathbf{x}_t\). After \(K\) iterations, we compute the average gradient in the block \(m\), \(\bar{\nabla}_m=\frac{1}{K} \sum_{t\in\mathcal{T}_m} \nabla_t\), and perform a gradient descent update, \(\mathbf{x}_m - \eta_m\bar{\nabla}_m\), with adaptive step-size \(\eta_m = \frac{D}{\sqrt{\epsilon + \sum^m_{\tau=1}\Vert \bar{\nabla}_\tau \Vert^2}}\) in the convex case, or \(\eta_m = \frac{1}{m \theta}\) in the strongly convex case. The gradient descent update \(\mathbf{x}_m-\eta_m\bar{\nabla}_m\) is passed to the IP-SO instance \(\mathcal{P}\) to obtain next action, \(\mathbf{x}_{m+1}\).
Action set \(\mathcal{K}\), time horizon \(T\), block size \(K\), \(\mathbf{c} \in \operatorname{relint}(\mathcal{K})\), \(r = r_{\mathcal{K}, \mathbf{c}}\), shrinking parameter \(\delta\in[0,1)\), diameter of the action set \(D\), an instance of IP-SO algorithm \(\mathcal{P}\) implemented by Algorithm 1, (in the convex case) constant \(\epsilon > 0\), (in the strongly-convex case) strongly-convex parameter \(\theta\) Initialization: Pick \(\mathbf{x}_1\in\mathcal{K}\); pass \(\mathcal{K},\mathbf{c},r,\delta\) to \(\mathcal{P}\) Play \(\mathbf{x}_m\) Adversary reveals \(f_t\) Compute \(\nabla_t \gets \nabla f_t(\mathbf{x}_m)\) \(\bar{\nabla}_m \gets \frac{1}{K}\sum_{t\in\mathcal{T}_m}\nabla_t\) Step size: \(\eta_m \gets \frac{D}{\sqrt{\epsilon + \sum_{\tau=1}^m \|\bar{\nabla}_\tau\|^2}}\) (Convex) \(\eta_m \gets \frac{1}{m\,\theta}\) (Strongly Convex) \(\mathbf{x}_{m+1}\gets \mathcal{P}(\mathbf{x}_m - \eta_m\bar{\nabla}_m)\)
Theorem 1. Assume functions \(f_t\) are convex and \(G_t\)-Lipschitz continuous. Then, Algorithm [alg:general] ensures that the regret is bounded as
\[\begin{align} \mathtt{Regret}(T) &= \sum_{t=1}^T f_t(\mathbf{x}_t) - \sum_{t=1}^T f_t(\mathbf{x}^*) \\ &\leq \frac{3}{2} D K \sqrt{\sum_{m=1}^{T/K} \Vert \bar{\nabla}_m \Vert ^2} + D \delta \sum\limits_{t=1}^TG_t, \end{align}\]
where \(\delta \in [0,1)\) is the shrinking parameter, \(r = r_{\mathcal{K}, \mathbf{c}}\), \(K\) is the block size, \(D\) is the diameter of \(\mathcal{K}\), and the definition of \(\bar{\nabla}_m\) is given in Algorithm [alg:general] (line 9). If \(\epsilon > 0\) is a constant independent of \(T\), then this is achieved with \(\mathcal{O}\left( (D/r)^2\delta^{-2} \log(T) + T K^{-1} \right)\) calls to a separation oracle. On the other hand, if we only assume \(\epsilon \geq 0\), then the same bounds for regret and CCV are achieved with \(\mathcal{O}\left( (D/r)^2T K^{-1} \delta^{-2}\right)\) calls to a separation oracle.
Remark 1. Note that, when block size \(K\) is equal to \(1\) and step-size is constant, this result reduces to a variant of Theorem 14 in [11]. Also note that the aforementioned lemma allows for the the actions to be outside the constraint set, which is why their bound does not include the term \(D \delta \sum_{t=1}^T G_t\).
Proof. Here is a high-level sketch proof, and more details are provided in Appendix 9.
We first demonstrate proof for bound of regret, and then give proof for bound of total calls to the separation oracle. We denote the actual optimal action \(\mathbf{x}^*\triangleq \operatorname{argmin}_{\mathbf{x}\in\mathcal{K}}\sum_{t=1}^T f_t(\mathbf{x})\) and the shrunk optimal action \(\tilde{\mathbf{x}}^*\triangleq (1-\delta)\mathbf{x}^*+\delta{\mathbf{c}} \in \mathcal{K}_\delta\). The high level idea is to use the shrunk optimal action as a bridge to break Regret into two parts: the distance between the action to the shrunk optimal action, and the distance between the shrunk optimal action and the actual optimal action, and we bound them separately. In specific, \(\mathtt{Regret}(T) = \sum\limits_{t=1}^T\left( f_t(\mathbf{x}_t) - f_t(\tilde{\mathbf{x}}^*) \right) + \sum\limits_{t=1}^T\left( f_t(\tilde{\mathbf{x}}^*) - f_t(\mathbf{x}^*) \right).\) For the first term, using Lemma 3, summing over telescoping term, and the convexity of the cost function \(f_t\), we have \[\sum\limits_{t=1}^Tf_t (\mathbf{x}_t) - \sum\limits_{t=1}^Tf_t (\tilde{\mathbf{x}}^*) \leq \frac{K D^2}{2\eta_{T/K}} + \frac{K}{2} \sum\limits_{m=1}^{T/K} \eta_m \Vert \bar{\nabla}_m \Vert^2.\] Although \(\eta_m = \frac{D}{\sqrt{\sum_{m=1}^{T/K} \Vert \bar{\nabla}_m \Vert^2}}\) would minimize the right hand side, we do not know apriori the future \(\bar{\nabla}_n\), \(n>m\) at any given \(m\). Instead, we use adaptive step size \(\eta_m = \frac{D}{\sqrt{\epsilon + \sum^m_{\tau=1}\Vert \bar{\nabla}_\tau \Vert^2}}\) for asymptotic approximation. Note that \(\epsilon\) is not present in previous works [6], [24], and we will discuss in Remark 2 why we need it for analysis of total calls to the separation oracles. By applying Lemma 1, we obtain that \(\sum_{m=1}^{T/K} \frac{\Vert \bar{\nabla}_m \Vert^2}{\sqrt{\epsilon + \sum^m_{\tau=1}\Vert \bar{\nabla}_\tau \Vert^2}} \leq 2 \sqrt{\sum_{m=1}^{T/K} \|\bar{\nabla}_m\|^2}.\) Thus, \[\sum\limits_{t=1}^Tf_t (\mathbf{x}_t) - \sum\limits_{t=1}^Tf_t (\tilde{\mathbf{x}}^*) \leq \frac{3 K D}{2}\sqrt{\sum\limits_{m=1}^{T/K}\Vert \bar{\nabla}_m \Vert^2}.\] For the second term, since \(f_t\) are \(G_t\)-Lipschitz continuous, and \(\Vert \tilde{\mathbf{x}}^*-\mathbf{x}^*\Vert \leq D\delta\), we have \[\sum\limits_{t=1}^T\left(f_t(\tilde{\mathbf{x}}^*)-f_t(\mathbf{x}^*)\right) \leq D \delta \sum\limits_{t=1}^TG_t.\] Combing the two, we have the regret bound in Theorem 1.
For total calls to the separation oracle, using triangle inequality, \(\forall \mathbf{y},\mathbf{y}_0 \in \mathcal{K}\), we have \(\operatorname{dist}(\mathbf{y}_0,\mathcal{K}_\delta)^2 - \operatorname{dist}(\mathbf{y},\mathcal{K}_\delta)^2 \leq \| \mathbf{y}- \mathbf{y}_0 \|^2 + 2D\delta \| \mathbf{y}- \mathbf{y}_0 \|,\) and for \(\mathbf{z}\in\mathcal{K}\), we have \(\operatorname{dist}(\mathbf{z},\mathcal{K}_\delta) \leq D\delta.\) Thus, using Lemma 3, we see the total calls to the separation oracle for Algorithm [alg:general] is bounded by \(\sum_{m = 1}^{T/K - 1}\left( \frac{\operatorname{dist}(\mathbf{x}_m - \eta_m \bar{\nabla}_m,\mathcal{K}_\delta)^2-\operatorname{dist}(\mathbf{x}_{m+1},\mathcal{K}_\delta)^2 }{\delta^2r^2}+1 \right) \leq \frac{1}{\delta^2} \sum_{m = 1}^{T/K} \left( \eta_m^2 \| \bar{\nabla}_m \|^2 + \frac{2 D \delta}{r} \eta_m \| \bar{\nabla}_m \| \right) + {\frac{D^2}{r^2}} + \frac{T}{K}.\)
If \(\epsilon>0\), then using Lemma 1, we have \(\sum_{m = 1}^{T/K} \frac{\| \bar{\nabla}_m \|^2}{\epsilon + \sum_{\tau=1}^{m}\| \bar{\nabla}_m \|^2} \leq \log\left( \frac{1}{\epsilon} \sum_{m = 1}^{T/K} \| \bar{\nabla}_m \|^2 \right),\) and therefore \(\sum_{m = 1}^{T/K} \eta_m^2 \| \bar{\nabla}_m \|^2 \leq 2D^2\log\left( \frac{M_1 T}{\epsilon K} \right).\) With Cauchy-Schwarz inequality, we have \(\left( \sum_{m = 1}^{T/K} \eta_m \| \bar{\nabla}_m \| \right)^2 \leq \frac{2D^2T}{K} \log\left( \frac{M_1 T}{\epsilon K} \right).\) Thus, the total number of calls is bounded by \[\begin{align} &\frac{2D^2}{\delta^2r^2} \log\left( \frac{M_1 T}{\epsilon K} \right) + \frac{2 \sqrt{2} D^2}{r^2 \delta} \sqrt{\frac{T}{K} \log\left( \frac{M_1 T}{\epsilon K} \right)} \\ & \quad + \frac{D^2}{r^2} + \frac{T}{K} = \mathcal{O}\left( (D/r)^2\delta^{-2} \log(T) + T K^{-1} \right). \end{align}\]
Remark 2. Note that if we if we only assume that \(\epsilon \geq 0\), instead of \(\epsilon>0\), then the above regret bound becomes infinite as it contains the term \(\log\left( \frac{M_1 T}{\epsilon K} \right)\). However, we may instead bound \(\sum_{m = 1}^{T/K} \frac{\| \bar{\nabla}_m \|^2}{\epsilon + \sum_{\tau=1}^{m}\| \bar{\nabla}_m \|^2}\) with \(T/K\), since each term is bounded by \(1\). Thus, with Cauchy-Schwartz inequality, we see that \(\left( \sum_{m = 1}^{T/K} \eta_m \| \bar{\nabla}_m \| \right)^2 \leq \frac{2D^2 T}{K} \cdot \frac{T}{K}\) and therefore the total number of calls to the separation oracle would be bounded by \[\begin{align} &\frac{2D^2}{\delta^2r^2} \frac{T}{K} + \frac{2 \sqrt{2} D^2}{r^2 \delta} \frac{T}{K} + \frac{D^2}{r^2} + \frac{T}{K} \\ & \quad= \mathcal{O}\left( (D/r)^2 T K^{-1} \delta^{-2} \right). \qedhere \end{align}\]
◻
Corollary 1. In Theorem 1, if we further assume that functions are \(M_1\)-Lipschitz continuous, then we see that \[\begin{align} \mathtt{Regret}(T) &\leq \frac{3}{2} D \sqrt{M_1 K T} + D \delta T M_1. \end{align}\] In particular, if \(K = \Theta(T^{1 - 2\beta})\) and \(\delta = \Theta(T^{-\beta})\) for some \(\beta \in [0, 1/2]\), then \(\mathtt{Regret}(T) = \mathcal{O}(DM_1^{1/2}T^{1 - \beta})\) with \(\tilde{\mathcal{O}}(T^{2\beta})\) calls to the separation oracle. On the other hand, if we assume that \(\epsilon = 0\), we obtain the same bounds for regret and CCV with \(\tilde{\mathcal{O}}(T^{4\beta})\) calls to the separation oracle.
We note that better regret bounds are achievable when the cost functions are chosen to be strongly convex. Thus, we give Theorem 3 to describe results for strongly convex functions. The proof details for Theorem 3 and Corollary 2 can be found in Appendix 10.
In this section, to solve Online Convex Optimization with Adversarial Constraint (Adversarial COCO), we propose Algorithm [alg:constrained], Blocked Adaptive online Gradient dEscent with infeasibLe projection (\(\mathtt{BAGEL}\)), where “infeasible projection” refers to the use of IP-SO described in Section 3. For computational convenience, we introduce a processing parameter \(\gamma\) and we let \(\tilde{f}_t = \gamma f_t\), \(\tilde{g}_{t,i}=\gamma(g_{t,i})^+\). To track the cumulative constraint violation, we let \(Q_t=Q_{t-1}+\tilde{g}_t(\mathbf{x}_t)\), and \(Q_0=0\). Inspired by [6], \(\mathtt{BAGEL}\) (Algorithm 2) combines costs and constraints into a single surrogate cost function, \(\hat{f}_t\), and applies Base Algorithm to it. This surrogate consists of the original cost function and the current constraint violation weighted by dynamic term. If violations have been low, the weight is small, and the algorithm focuses on minimizing cost. If violations have been high, the weight becomes large, forcing the algorithm to aggressively prioritize satisfying the constraint. Since our objective is to make the weighted cumulative constraint violation small, similar to [6], we introduce a potential function \(\Phi: \mathbb{R}^+\mapsto\mathbb{R}^+\) that is non-decreasing, differentiable, convex, and satisfies \(\Phi(0)=0\). Thus, we create the surrogate cost function to be \(\hat{f}= \tilde{f}_t + \Phi'(Q_t)\tilde{g}_t\) where \(\Phi'(Q_t)\) is the gradient of the pre-defined function \(\Phi\) at \(Q_t\).
Action set \(\mathcal{K}\), time horizon \(T\), \(\mathbf{c} \in \operatorname{relint}(\mathcal{K})\), \(r = r_{\mathcal{K}, \mathbf{c}}\), shrinking parameter \(\delta \in [0,1)\), block size \(K\), action set diameter \(D\), an instance of the Base Algorithm (Algorithm [alg:general]) \(\mathcal{A}\), processing parameter \(\gamma\), regularization parameter \(V\), Lyapunov function \(\Phi(\cdot)\), an instance of IP-SO (Algorithm 1) \(\mathcal{P}\), (in the convex case) constant \(\epsilon > 0\), (in the strongly-convex case) strongly-convex parameter \(\theta\). Initialization: Pick \(\mathbf{x}_1\in\mathcal{K}\); set \(Q_0\gets 0\); pass \(\mathcal{K},T,K,\mathbf{c},r,\delta,D,\mathcal{P}\), and \((\epsilon\text{ or }\theta)\) to \(\mathcal{A}\) Play \(\mathbf{x}_m\) Adversary reveals \(f_t\) and \(g_t\) Observe \(f_t(\mathbf{x}_m)\) and \((g_t(\mathbf{x}_m))^{+}\) \(\tilde{f}_t \gets \gamma\, f_t\); \(\tilde{g}_t \gets \gamma\, (g_t)^{+}\) \(Q_t \gets Q_{t-1} + (g_t(\mathbf{x}_m))^{+}\) \(\hat{f}_t \gets V\,\tilde{f}_t + \Phi'(Q_t)\,\tilde{g}_t\); \(\nabla_t \gets \nabla \hat{f}_t(\mathbf{x}_m)\) Pass \(\nabla_t\) to \(\mathcal{A}\) Receive \(\mathbf{x}_{m+1}\) from \(\mathcal{A}\)
The algorithm use the inputs: action set \(\mathcal{K}\), time horizon \(T\), \(\mathbf{c} \in \operatorname{relint}(\mathcal{K})\), \(r = r_{\mathcal{K}, \mathbf{c}}\), shrinking parameter \(\delta \in [0,1)\), block size \(K\), action set diameter \(D\), an instance of the Base Algorithm (Algorithm [alg:general]) \(\mathcal{A}\), processing parameter \(\gamma\), function \(\Phi(\cdot)\), an instance of IP-SO \(\mathcal{P}\), and constant \(\epsilon > 0\) for convex costs or strongly-convex parameter \(\theta\). We initiate the algorithm by picking the initial action \(\mathbf{x}_1\in\mathcal{K}\), and set \(Q_0=0\), and passing the required input to \(\mathcal{A}\) to initiate the algorithm instance. We break the horizon \(T\) into blocks with size \(K\). For each step \(t\) within block \(m\), represented by \(\mathcal{T}_m\), we play action \(\mathbf{x}_m\) determined by the algorithm, and then observe cost \(f_t(\mathbf{x}_t)\) and constraint violation \((g_t(\mathbf{x}_t))^+\), where the cost function and the maximum constraint function are chosen and revealed by an adversary. We transform them into \(\tilde{f}_t = \gamma f, \tilde{g}_t =\gamma (g_t(x_t))^+\), and update the cumulative constraint violation with \(Q_t = Q_{t-1}+\tilde{g}_t\). We compute the gradient \(\nabla_t\) of the surrogate function \(\hat{f}_t(x_t):=\tilde{f}_t + \Phi'(Q_t)\tilde{g}_t\), i.e., \(\nabla_t=\nabla \hat{f}_t\), and pass the gradient information to \(\mathcal{A}\). At the end of the each block, \(\mathcal{A}\) returns the action for the next block \(\mathbf{x}_{m+1}\). The following Theorem provides the regret and the CCV for the proposed algorithm in the convex case.
Theorem 2. Assume cost functions \(f_t\) and constraint functions \(g_t\) are convex and \(M_1\)-Lipschitz continuous over action set \(\mathcal{K}\), with diameter \(D\). Let \(\mathbf{c} \in \operatorname{relint}(\mathcal{K})\), and \(r = r_{\mathcal{K}, \mathbf{c}}\). With \(\beta\in[0,\frac{1}{2}]\), if we choose \(\gamma=(M_1R)^{-1}\), \(\delta=\Theta(T^{-\beta})\), \(K= \Theta(T^{1-2\beta})\), and \(\Phi(x)=e^{\lambda x}-1\) where \(\lambda=(2\delta T + 3\sqrt{2TK})^{-1}\), Algorithm [alg:constrained] ensures that: \[\begin{align} \mathtt{Regret}(T) &\leq D M_1 \delta T + \frac{3}{\sqrt{2}} DM_1\sqrt{TK} + DM_1 \\ &= \mathcal{O}(DM_1T^{1-\beta}), \end{align}\] and \[\begin{align} \mathtt{CCV}(T) &\leq 2 D M_1 (\delta T + \frac{3}{\sqrt{2}}\sqrt{TK}) \mathcal{O}(\log T) \\ &= \mathcal{O}(DM_1T^{1-\beta}\log T). \end{align}\] with \(\tilde{\mathcal{O}}((D/r)^2T^{2\beta})\) calls to the separation oracles.
Remark 3. Note that if we take \(\beta=\frac{1}{2}\), we match the \(\sqrt{T}\) for both regret and CCV of projection-based methods, with \(\tilde{\mathcal{O}}(T)\) total calls to the separation oracles. In this case, there will be no blocking effect because block size \(K=1\). However, considering the diverse application scenario, we introduce the trade-off parameter \(\beta\) to allow adjustment of the algorithm to real-life computational limitations. When \(\beta=\frac{1}{4}\), we have \(\mathcal{O}(T^{3/4})\) result for regret and CCV, which matches [7] on regret and [8] on regret and CCV, with at most \(\tilde{\mathcal{O}}(\sqrt{T})\) total calls to the separation oracles.
Remark 4. Note that our bounds of total calls to the SO contains the parameter \((D/r)^2\). The ratio \(D/r\) represents the condition number of the constraint set geometry (the ratio between the diameter and the radius of the largest enclosed ball). For “well-conditioned” sets (e.g., a ball or hypercube), this ratio is small or constant. However, for “ill-conditioned” sets in high dimensions, \(D/r\) can scale polynomially with the dimension \(d\) (e.g., \(\sqrt{d}\) or \(d\)). Thus, BAGEL offers a significant advantage in regimes where the time horizon \(T\) is the dominant factor (\(T \gg (D/r)^2\)), whereas LOO methods remain preferable for high-dimensional sets where \(D/r\) is prohibitively large.
Proof. A high-level sketch proof is given as follows, and more details can be found in Appendix 11.
We define \(\mathtt{R}(\hat{f}) \triangleq \sum_{t=1}^T \hat{f}_t(\mathbf{x}_t) - \sum_{t=1}^T \hat{f}_t(\mathbf{x}^*)\) and \(\mathtt{R}(\tilde{f}) \triangleq \sum_{t=1}^T \tilde{f}_t(\mathbf{x}_t) - \sum_{t=1}^T \tilde{f}_t(\mathbf{x}^*)\) where \(\mathbf{x}^* \triangleq \arg\min_{\mathbf{x}\in\mathcal{K}}\sum_{t=1}^T f_t(\mathbf{x})\). Since we pass gradient information of the surrogate functions \(\hat{f}_t\) to an instance of Algorithm [alg:general] at each iteration and then get updated actions back at each block, we can use Theorem 1 to derive \[\begin{align} \mathtt{R}(\hat{f}) \leq \frac{3}{2}DK\sqrt{\sum_{m=1}^{T/K} \Vert \bar{\nabla}_m \Vert ^2} + \delta D \sum\limits_{t=1}^TG_t, \end{align}\]
which is also formally described and proved in Lemma 4.
With \(f\) and \(g\) both \(M_1\)-Lipschitz continuous, we have \(\Vert \nabla_t \Vert\leq \gamma M_1(1+\Phi'(Q_t))\), and thus, surrogate functions \(\hat{f}_t\) are \(G_t\)-Lipschitz continuous where \(G_t=\gamma M_1(1+\Phi'(Q_t))\). Using Cauchy-Schwartz inequality, Jenson’s inequality, the convexity of \(\Phi(\cdot)\), and the non-decreasing nature of \(Q_t\), we have \(\sqrt{\sum_{m=1}^{T/K} \Vert \bar{\nabla}_m \Vert ^2} \leq \sqrt{ \frac{1}{K} \sum_{t=1}^{T} \Vert \nabla_t \Vert^2 } \leq \frac{\gamma M_1 \sqrt{2T}}{\sqrt{K}} (1 + \Phi'(Q_T))\). Substituting the \(\sqrt{\sum_{m=1}^{T/K} \Vert \bar{\nabla}_m \Vert ^2}\) term, we have \[\mathtt{R} (\hat{f}) \leq \gamma M_1 D (1+\Phi'(Q_T))(\delta T + \frac{3}{\sqrt{2}}\sqrt{TK}),\] which is formally described and proved in Equation 8 of Lemma 5.
Since \(\Phi(\cdot)\) is non-decreasing, differentiable and convex satisfying \(\Phi(0)=0\), we have \(\Phi(Q_t) - \Phi(Q_{t-1}) \leq \Phi'(Q_t)\tilde{g}_t(\mathbf{x}_t)\). Since \(\hat{f}_t(\mathbf{x}^*)=\tilde{f}_t(\mathbf{x}^*)+\Phi'(Q_t)\tilde{g}_t(\mathbf{x}^*)\) and \(\tilde{g}_t(\mathbf{x}^*)=\gamma (g_t(\mathbf{x}^*))^+ \leq 0\), plugging in and summing over telescoping term, we have \[\mathtt{R}(\hat{f}) \geq \mathtt{R}(\tilde{f}) + \Phi(Q_T),\] which is formally described and proved in Equation 8 of Lemma 5.
If we let \(\Phi(Q_T)=e^{\lambda Q_T}-1\), and denote \(S \triangleq \delta T + \frac{3}{\sqrt{2}}\sqrt{TK}\), we have \[\begin{align} \mathtt{R}(\tilde{f}) &\leq \mathtt{R} (\hat{f})-\Phi (Q_T) \\ &\leq \gamma M_1 D S + e^{\lambda Q_T} \left( \lambda \gamma M_1 D S-1 \right)+1. \end{align}\] Let \(\beta\in[0,\frac{1}{2}]\) be some real number. If we take \(\gamma=(M_1D)^{-1}\), \(\delta=\Theta(T^{-\beta})\), \(K=\Theta(T^{1-2\beta})\), \(\lambda=\frac{1}{2S}\), then we have \(\mathtt{R}(\tilde{f}) \leq S+1\). Since \(\tilde{f}_t=\gamma f_t\), we have \[\begin{align} \mathtt{Regret}(T) &= \gamma^{-1} \mathtt{R}(\tilde{f}) \\ &= {D M_1 \delta T + \frac{3}{\sqrt{2}} DM_1\sqrt{TK} + DM_1} \\ &= \mathcal{O}(D M_1 T^{1-\beta}). \end{align}\]
To investigate CCV, we first find the lower bound of \(\mathtt{R}(\tilde{f})\). Using the fact that \(\tilde{f}_t\) is \(\gamma M_1\)-Lipschitz continuous, i.e., \(\tilde{f}_t(\mathbf{x}_t)-\tilde{f}_t(\mathbf{x}^*) \geq -\gamma M_1 (\mathbf{x}_t-\mathbf{x}^*) \geq -1\), and summing over \(T\), we have \[\mathtt{R}(\tilde{f})=\sum\limits_{t=1}^T(\tilde{f}_t(\mathbf{x}_t)-\tilde{f}_t(\mathbf{x}^*)) \geq -T.\] Since \(\gamma=(M_1D)^{-1}\) and \(\lambda=\frac{1}{2S}\), we have \[\begin{align} \mathtt{R}(\tilde{f}) & \leq \gamma M_1 D S + e^{\lambda Q_T} \left( \lambda \gamma M_1 D S-1 \right)+1 \\ & \leq S - \frac{1}{2} e^{\frac{Q_T}{2S}}+1. \end{align}\] Thus, \(Q_T \leq 2S \log\left( 2(S + 1 + T )\right).\) As \(Q_T=\gamma \mathtt{CCV}(T)\) and \(S=\delta T + \frac{3}{\sqrt{2}}\sqrt{TK}\), we have \[\begin{align} \mathtt{CCV}(T) & = \gamma^{-1} Q_T \leq 2DM_1S * \log\left( 2(S + 1 + T )\right) \\ &\leq 2 D M_1 (\delta T + \frac{3}{\sqrt{2}}\sqrt{TK}) * \mathcal{O}(\log T) \\ & \leq \mathcal{O}(DM_1T^{1-\beta}\log T) \end{align}\]
Since calls to the separation oracles occur when \(\mathcal{A}\) is called, we can extend the results of Corollary 1 and derive that the oracle calls are bounded as \(\mathcal{O}((D/r)^2T^{2\beta}\log(T))=\tilde{\mathcal{O}}(T^{2\beta})\). ◻
In this paper, we investigated the gap regarding horizon \(T\) between projection-free and projection-based algorithms for for online convex optimization with adversarial constraints. We demonstrated that this gap can be closed for action sets with efficient Separation Oracle. Our algorithm, \(\mathtt{BAGEL}\), is the first projection-free algorithm to match the \(\tilde{\mathcal{O}}(T^{1/2})\) regret and CCV bounds of projection-based methods, using a number of separation oracle calls that is near-linear in the time horizon \(T\). We note that our algorithm assumes access to a separation oracle. The question remains open whether the bounds can be achieved with acess to LOO.
Here we introduce some technical lemmas that are used in our proofs.
Lemma 1 (Lemma 4.13 in [24]). Let \(a_0 \geq 0\) be a real number, \(N \geq 1\) be an integer, and \((a_t)_{t = 1}^N\) be a sequence of non-negative real numbers. For any non-increasing function \(f : [0, \infty) \to [0, \infty)\), we have \[\begin{align} \sum_{t = 1}^N a_t f\left( a_0 + \sum_{i=1}^t a_i \right) \leq \int_{a_0}^{\sum_{t = 1}^N a_t} f(x) dx. \end{align}\]
Lemma 2. Assume \(f\) is \(\theta\)-strongly convex and \(g\) is convex. Let \(a\) and \(b\) be some non-negative real number. Then \(a f + b g\) is \(a\theta\)-strongly convex.
Proof. Given \(f\) is \(\theta\)-strongly convex, its Hessian satisfies \(\nabla^2 f(x) \succeq \theta I\), where \(I\) stands for the identity matrix of appropriate dimensions. Since \(g\) is convex, its Hessian satisfies \(\nabla^2 g(x) \succeq 0\). Let \(h(x) = af(x) + bg(x)\). Then, its Hessian is given by \(\nabla^2 h(x) = a \nabla^2 f(x) + b \nabla^2 g(x) \succeq a\theta I\), where \(I\) stands for the identity matrix of appropriate dimensions. Thus, \(h(x)\) is strongly convex with parameter \(a\theta\). ◻
As mentioned in Section 2, our work builds upon the concept of infeasible projection, implemented via a separation oracle, pioneered by [11] for the unconstrained OCO setting. We specifically adopt the variant from [20] due to its favorable property of ensuring iterates remain within the constraint set \(\mathcal{K}\), which is crucial for our analysis. Here we formally introduce the details of IP-SO algorithm in Algorithm 1 and an important Lemma 3.
Constraint set \(\mathcal{K}\), \(\mathbf{c} \in \operatorname{relint}(\mathcal{K})\), \(r = r_{\mathcal{K}, \mathbf{c}}\), shrinking parameter \(\delta \in [0,1)\), initial point \(\mathbf{y}_0\) \(\mathbf{y}_1 \gets \mathbf{P}_{\operatorname{aff}(\mathcal{K})}(\mathbf{y}_0)\) \(\mathbf{y}_2 \gets \mathbf{c} + \frac{\mathbf{y}_1 - \mathbf{c}}{\max\{1, \|\mathbf{y}_1\|/D \}}\) Call \(\operatorname{SO}_{\mathcal{K}}\) with input \(\mathbf{y}_i\) Set \(\mathbf{g}_i\) to be the hyperplane returned by \(\operatorname{SO}_{\mathcal{K}}\) \(\mathbf{g}'_i \gets \mathbf{P}_{\operatorname{aff}(\mathcal{K})-\mathbf{c}}(\mathbf{g}_i)\) \(\mathbf{y}_{i+1} \gets \mathbf{y}_i - \delta r\, \dfrac{\mathbf{g}'_i}{\|\mathbf{g}'_i\|}\) \(\mathbf{y}\gets \mathbf{y}_i\)
Lemma 3 (Lemma 5 in [20]). Given a constrained set \(\mathcal{K}\), a shrinking parameter \(\delta \in [0, 1)\) (where \(r\) is as previously described), and an initial point \(\mathbf{y}_0\), Algorithm 1 stops after at most \(\frac{\operatorname{dist}(\mathbf{y}_0,\mathcal{K}_\delta)^2-\operatorname{dist}(\mathbf{y},\mathcal{K}_\delta)^2}{\delta^2r^2}+1\) iterations and returns \(\mathbf{y}\in \mathcal{K}\) such that for all \(\mathbf{x}\in \mathcal{K}_\delta\), we have \(\Vert \mathbf{y}-\mathbf{x}\Vert \leq \Vert \mathbf{y}_0 - \mathbf{x}\Vert\).
Proof. Let \(\mathbf{x}^*\in \operatorname{argmin}_{\mathbf{x}\in\mathcal{K}}\sum_{t=1}^T f_t(\mathbf{x})\) and \(\tilde{\mathbf{x}}^*=(1-\delta)\mathbf{x}^*+\delta{\mathbf{c}} \in \mathcal{K}_\delta\). Thus, the regret of the proposed algorithm is given as \[\begin{align} \label{eq:decomposeRegret} \mathtt{Regret}(T) &= \sum\limits_{t=1}^T\left( f_t(\mathbf{x}_t) - f_t(\mathbf{x}^*) \right) \nonumber \\ & = \sum\limits_{t=1}^T\left( f_t(\mathbf{x}_t) - f_t(\tilde{\mathbf{x}}^*) \right) \nonumber \\ & + \sum\limits_{t=1}^T\left( f_t(\tilde{\mathbf{x}}^*) - f_t(\mathbf{x}^*) \right). \end{align}\tag{1}\] We bound the two terms separately.
First we consider the term \(\sum_{t=1}^T \left( f_t(\mathbf{x}_t) - f_t(\tilde{\mathbf{x}}^*) \right)\). From Lemma 3, we have for \(\tilde{\mathbf{x}}^*\in \mathcal{K}_\delta\), \[\begin{align} \Vert \mathbf{x}_{m+1}-\tilde{\mathbf{x}}^*\Vert ^2 & \leq \Vert \mathbf{x}_m - \eta_m \bar{\nabla}_m -\tilde{\mathbf{x}}^*\Vert ^2 \\ &\leq \Vert \mathbf{x}_m - \tilde{\mathbf{x}}^*\Vert ^2 + \eta_m^2 \Vert \bar{\nabla}_m \Vert^2 \\ & - 2\eta_m \langle \bar{\nabla}_m, \mathbf{x}_m - \tilde{\mathbf{x}}^*\rangle. \end{align}\] Rearranging, we have \[\begin{align} \langle \bar{\nabla}_m, \mathbf{x}_m - \tilde{\mathbf{x}}^*\rangle &\leq \frac{1}{2\eta_m} (\Vert \mathbf{x}_m - \tilde{\mathbf{x}}^*\Vert ^2- \Vert \mathbf{x}_{m+1}-\tilde{\mathbf{x}}^*\Vert ^2) \nonumber \\ &+ \frac{\eta_m}{2} \Vert \bar{\nabla}_m \Vert^2 \end{align}\] Summing this over \(m\) from \(1\) to \(T/K\), we have \[\begin{align} &\sum\limits_{m=1}^{T/K} \langle \bar{\nabla}_m , \mathbf{x}_m - \tilde{\mathbf{x}}^*\rangle \leq \frac{\Vert \mathbf{x}_1 - \tilde{\mathbf{x}}^*\Vert ^2}{2\eta_1} \\ & \quad + \sum\limits_{m=2}^{T/K} (\frac{1}{2\eta_m}-\frac{1}{2\eta_{m-1}}) \Vert \mathbf{x}_m - \tilde{\mathbf{x}}^*\Vert ^2 + \sum\limits_{m=1}^{T/K} \frac{\eta_m}{2} \Vert \bar{\nabla}_m \Vert^2 \\ &\quad \leq \frac{D^2}{2\eta_1} + D^2 \sum\limits_{m=2}^{T/K} \left| \frac{1}{2\eta_m}-\frac{1}{2\eta_{m-1}} \right| + \sum\limits_{m=1}^{T/K} \frac{\eta_m}{2} \Vert \bar{\nabla}_m \Vert^2 \\ &\quad = \frac{D^2}{2\eta_1} + D^2 \sum\limits_{m=2}^{T/K} \left( \frac{1}{2\eta_m}-\frac{1}{2\eta_{m-1}} \right) + \sum\limits_{m=1}^{T/K} \frac{\eta_m}{2} \Vert \bar{\nabla}_m \Vert^2 \\ &\quad = \frac{D^2}{2\eta_{T/K}} + \sum\limits_{m=1}^{T/K} \frac{\eta_m}{2} \Vert \bar{\nabla}_m \Vert^2, \end{align}\] where we used the facts that \(\Vert \mathbf{x}_m - \tilde{\mathbf{x}}^*\Vert \leq D\) and that \(\eta_m\) is non-increasing. By the convexity of \(f_t\), we have \(f_t (\mathbf{x}_t) - f_t(\tilde{\mathbf{x}}^*) \leq \langle \nabla_t, \mathbf{x}_t - \mathbf{z}\rangle\). Summing this over \(t\), we have \[\begin{align} \label{eq:xt2txs} &\sum\limits_{t=1}^Tf_t (\mathbf{x}_t) - \sum\limits_{t=1}^Tf_t (\tilde{\mathbf{x}}^*) \leq \sum\limits_{t=1}^T\langle \nabla_t, \mathbf{x}_t - \tilde{\mathbf{x}}^*\rangle \nonumber \\ &\quad \quad = \sum\limits_{m=1}^{T/K} \sum\limits_{t\in\mathcal{T}_m} \langle \nabla_t, \mathbf{x}_t - \tilde{\mathbf{x}}^*\rangle \nonumber\\ &\quad \quad= \sum\limits_{m=1}^{T/K} K \langle \bar{\nabla}_m, \mathbf{x}_m - \tilde{\mathbf{x}}^*\rangle \nonumber \\ &\quad \quad \leq \frac{K D^2}{2\eta_{T/K}} + \frac{K}{2} \sum\limits_{m=1}^{T/K} \eta_m \Vert \bar{\nabla}_m \Vert^2. \end{align}\tag{2}\] Note that if \(\eta_t\) is chosen to be a constant, \(\eta_m = \frac{D}{\sqrt{\sum_{m=1}^{T/K} \Vert \bar{\nabla}_m \Vert^2}}\) minimizes the right hand side. However, this step size could not be chosen since at a given \(m\), the future \(\bar{\nabla}_m\)’s are not known. Thus, we use adaptive step size in this work, where \(\eta_m = \frac{D}{\sqrt{\epsilon + \sum^m_{\tau=1}\Vert \bar{\nabla}_\tau \Vert^2}}\). Using this, Equation 2 is further bounded as \[\begin{align} & \sum\limits_{t=1}^Tf_t (\mathbf{x}_t) - \sum\limits_{t=1}^Tf_t (\tilde{\mathbf{x}}^*) \leq \frac{K D^2}{2 \eta_{T/K}} + \frac{K}{2} \sum\limits_{m=1}^{T/K} \eta_m \Vert \bar{\nabla}_m \Vert^2 \nonumber \\ & \leq \frac{K D}{2}\sqrt{\sum\limits_{m=1}^{T/K}\Vert \bar{\nabla}_m \Vert^2} + \frac{K}{2} \sum\limits_{m=1}^{T/K} \frac{D \Vert \bar{\nabla}_m \Vert^2}{\sqrt{\epsilon + \sum^m_{\tau=1}\Vert \bar{\nabla}_\tau \Vert^2}}. \end{align}\] By applying Lemma 1 to \(f(x) = 1/\sqrt{x}\), \(a_0 = \epsilon\), \(N = T/K\), and \(a_m = \Vert \bar{\nabla}_m \Vert^2\) for \(m\ge 1\), we see that \[\begin{align} &\sum_{m=1}^{T/K} \frac{\Vert \bar{\nabla}_m \Vert^2}{\sqrt{\epsilon + \sum^m_{\tau=1}\Vert \bar{\nabla}_\tau \Vert^2}} = \sum_{m=1}^{T/K} a_m f\left( a_0 + \sum^m_{\tau=1} a_\tau \right) \\ &\qquad \leq \int_{a_0}^{\sum_{t = \tau}^{T/K} a_\tau} f(x) dx = 2 \sqrt{\sum_{m=1}^{T/K} \|\bar{\nabla}_m\|^2} - 2 \sqrt{\epsilon} \\ &\qquad \leq 2 \sqrt{\sum_{m=1}^{T/K} \|\bar{\nabla}_m\|^2}. \end{align}\] Thus \[\begin{align} \label{eq:xt2txs1} \sum\limits_{t=1}^Tf_t (\mathbf{x}_t) - \sum\limits_{t=1}^Tf_t (\tilde{\mathbf{x}}^*) &\leq \frac{3 K D}{2}\sqrt{\sum\limits_{m=1}^{T/K}\Vert \bar{\nabla}_m \Vert^2}. \end{align}\tag{3}\]
On the other hand, the expression \(\sum\limits_{t=1}^T\left(f_t(\tilde{\mathbf{x}}^*)-f_t(\mathbf{x}^*)\right)\) may be bounded as \[\begin{align} \label{eq:txs2xs} \sum\limits_{t=1}^T\left(f_t(\tilde{\mathbf{x}}^*)-f_t(\mathbf{x}^*)\right) &\overset{(a)}{\leq} \sum\limits_{t=1}^TG_t \Vert \tilde{\mathbf{x}}^*-\mathbf{x}^*\Vert \nonumber\\ & \overset{(b)}{\leq} \delta D \sum\limits_{t=1}^TG_t, \end{align}\tag{4}\] where (a) is because \(f_t\) are \(G_t\)-Lipschitz continuous, and (b) is due to the fact that \[\begin{align} \Vert \tilde{\mathbf{x}}^*-\mathbf{x}^*\Vert= \Vert (1-\delta)\mathbf{x}^*+\delta\mathbf{c}-\mathbf{x}^*\Vert = \delta\Vert \mathbf{x}^*- \mathbf{c} \Vert \leq \delta D. \end{align}\] The claim now follows by using the bounds in 3 and 4 together with 1 .
Next we examine the number of calls to the separation oracles through the infeasible projection algorithm. Using triangle inequality, \(\forall \mathbf{y}, \mathbf{y}_0 \in \mathcal{K}\), we have \[\begin{align} \operatorname{dist}(\mathbf{y}_0,\mathcal{K}_\delta) &= \min_{\mathbf{z}\in \mathcal{K}_\delta} \| \mathbf{y}_0 - \mathbf{z}\| \\ &\leq \min_{\mathbf{z}\in \mathcal{K}_\delta} \left( \| \mathbf{y}- \mathbf{z}\| + \| \mathbf{y}- \mathbf{y}_0 \| \right) \\ &= \operatorname{dist}(\mathbf{y},\mathcal{K}_\delta) + \| \mathbf{y}- \mathbf{y}_0 \|. \end{align}\] Therefore \[\begin{align} &\operatorname{dist}(\mathbf{y}_0,\mathcal{K}_\delta)^2 - \operatorname{dist}(\mathbf{y},\mathcal{K}_\delta)^2 \\ &\leq \left( \operatorname{dist}(\mathbf{y},\mathcal{K}_\delta) + \| \mathbf{y}- \mathbf{y}_0 \| \right)^2 - \operatorname{dist}(\mathbf{y},\mathcal{K}_\delta)^2 \\ &= \| \mathbf{y}- \mathbf{y}_0 \|^2 + 2 \| \mathbf{y}- \mathbf{y}_0 \| \operatorname{dist}(\mathbf{y},\mathcal{K}_\delta) \\ &\leq \| \mathbf{y}- \mathbf{y}_0 \|^2 + 2 D \delta \| \mathbf{y}- \mathbf{y}_0 \|, \end{align}\] where the last inequality follows from that fact that, for all \(\mathbf{z}\in \mathcal{K}\), we have \(\tilde{\mathbf{z}} \triangleq (1-\delta)\mathbf{z}+\delta{\mathbf{c}} \in \mathcal{K}_\delta\) and therefore \[\begin{align} \operatorname{dist}(\mathbf{z},\mathcal{K}_\delta) &\leq \| \mathbf{z}- \tilde{\mathbf{z}} \| = \| \mathbf{z}- (1-\delta)\mathbf{z}-\delta{\mathbf{c}} \| \\ &= \delta \| \mathbf{z}- \mathbf{c} \| \leq D \delta. \end{align}\] Thus, using Lemma 3, we see that the number of calls in Algorithm [alg:general] is bounded by \[\begin{align} \label{eq:oracle-calls} &\sum_{m = 1}^{T/K - 1}\left( \frac{\operatorname{dist}(\mathbf{x}_m - \eta_m \bar{\nabla}_m,\mathcal{K}_\delta)^2-\operatorname{dist}(\mathbf{x}_{m+1},\mathcal{K}_\delta)^2 }{\delta^2r^2}+1 \right) \nonumber \\ &= \frac{1}{\delta^2r^2} \sum_{m = 1}^{T/K - 1}\left( \operatorname{dist}(\mathbf{x}_m - \eta_m \bar{\nabla}_m,\mathcal{K}_\delta)^2-\operatorname{dist}(\mathbf{x}_{m+1},\mathcal{K}_\delta)^2 \right) \nonumber \\ &\qquad+ \frac{T}{K} - 1 \nonumber \\ &\leq \frac{1}{\delta^2r^2} \sum_{m = 1}^{T/K - 1} ( \eta_m^2 \| \bar{\nabla}_m \|^2 + 2D\delta \eta_m \| \bar{\nabla}_m \| \nonumber \\ &\qquad+ \operatorname{dist}(\mathbf{x}_m, \mathcal{K}_\delta)^2-\operatorname{dist}(\mathbf{x}_{m+1},\mathcal{K}_\delta)^2 ) + \frac{T}{K} \nonumber \\ &= \frac{1}{\delta^2r^2} \sum_{m = 1}^{T/K - 1} \left( \eta_m^2 \| \bar{\nabla}_m \|^2 + 2D\delta \eta_m \| \bar{\nabla}_m \| \right) \nonumber \\ &\qquad + \frac{1}{\delta^2r^2} \left( \operatorname{dist}(\mathbf{x}_1,\mathcal{K}_\delta)^2-\operatorname{dist}(\mathbf{x}_{T/K},\mathcal{K}_\delta)^2 \right) + \frac{T}{K} \nonumber \\ &\leq \frac{1}{\delta^2r^2} \sum_{m = 1}^{T/K} \left( \eta_m^2 \| \bar{\nabla}_m \|^2 + 2D\delta \eta_m \| \bar{\nabla}_m \| \right) \nonumber \\ & \qquad + {\color{blue}}\frac{D^2}{r^2} + \frac{T}{K}, \end{align}\tag{5}\] where we used the fact that \(\operatorname{dist}(\mathbf{x}_1,\mathcal{K}_\delta) \leq D\delta\) and that \(\left( \eta_m^2 \| \nabla_m \|^2 + 2D\delta \eta_m \| \nabla_m \| \right)\) is positive when \(m=\frac{T}{K}\) in the last inequality.
If \(\epsilon > 0\), then using Lemma 1 for \(f(x) = 1/x\), \(N = T/K\), \(a_0 = \epsilon\), and \(a_i = \| \bar{\nabla}_i \|^2\), we see that \[\begin{align} \sum_{m = 1}^{T/K} \frac{\| \bar{\nabla}_m \|^2}{\epsilon + \sum_{\tau=1}^{m}\| \bar{\nabla}_m \|^2} &= \sum_{m = 1}^{T/K} a_m f\left( a_0 + \sum_{\tau=1}^{m} a_\tau \right) \\ &\leq \int_{a_0}^{\sum_{m = 1}^{T/K} a_m} f(x) dx \\ &= \log\left( \sum_{m = 1}^{T/K} \| \bar{\nabla}_m \|^2 \right) - \log(\epsilon) \\ &= \log\left( \frac{1}{\epsilon} \sum_{m = 1}^{T/K} \| \bar{\nabla}_m \|^2 \right). \end{align}\] Therefore \[\begin{align} \sum_{m = 1}^{T/K} \eta_m^2 \| \bar{\nabla}_m \|^2 &= D^2 \sum_{m = 1}^{T/K} \frac{\| \bar{\nabla}_m \|^2}{\epsilon + \sum_{\tau=1}^{m}\| \bar{\nabla}_m \|^2} \\ &\leq D^2\log\left( \frac{1}{\epsilon} \sum_{m = 1}^{T/K} \| \bar{\nabla}_m \|^2 \right) \\ &\leq 2D^2\log\left( \frac{M_1 T}{\epsilon K} \right), \end{align}\] where we used the fact that \(\bar{\nabla}_m=\frac{1}{K}\sum_{t\in\mathcal{T}_m}\nabla_t\) is bounded by \(M_1\). Thus, using Cauchy–Schwarz inequality, we see that \[\begin{align} \left( \sum_{m = 1}^{T/K} \eta_m \| \bar{\nabla}_m \| \right)^2 &\leq \left( \sum_{m = 1}^{T/K} 1 \right) \left( \sum_{m = 1}^{T/K} \eta_m^2 \| \bar{\nabla}_m \|^2 \right) \\ &\leq \frac{2D^2T}{K} \log\left( \frac{M_1 T}{\epsilon K} \right). \end{align}\] Therefore, the number of calls is bounded by \[\begin{align} &\frac{1}{\delta^2r^2} \sum_{m = 1}^{T/K - 1} \left( \eta_m^2 \| \nabla_m \|^2 + 2D\delta \eta_m \| \nabla_m \| \right) + \frac{D^2}{r^2} + \frac{T}{K} \\ &\leq \frac{2D^2}{\delta^2r^2} \log\left( \frac{M_1 T}{\epsilon K} \right) + \frac{2 \sqrt{2} D^2}{r^2 \delta} \sqrt{\frac{T}{K} \log\left( \frac{M_1 T}{\epsilon K} \right)} \\ & \qquad + \frac{D^2}{r^2} + \frac{T}{K} \\ &= \mathcal{O}( (D/r)^2\delta^{-2} \log(T) \\ & \qquad + (D/r)\delta^{-1} (\log(T))^{1/2} (T K^{-1})^{1/2} + T K^{-1} ) \\ &= \mathcal{O}\left( (D/r)^2\delta^{-2} \log(T) + T K^{-1} \right), \end{align}\] where we used the fact that \(\epsilon\) is a constant independent of \(T\) and \(\sqrt{a b} = \mathcal{O}(a + b)\) in the last equality.
Note that, if we only assume that \(\epsilon \geq 0\), then the above regret bound becomes infinite as it contains the term \(\log\left( \frac{M_1 T}{\epsilon K} \right)\). However, we may instead bound \(\sum_{m = 1}^{T/K} \frac{\| \bar{\nabla}_m \|^2}{\epsilon + \sum_{\tau=1}^{m}\| \bar{\nabla}_m \|^2}\) with \(T/K\), since each term is bounded by \(1\). Thus, using Cauchy–Schwarz inequality, we see that \[\begin{align} \left( \sum_{m = 1}^{T/K} \eta_m \| \bar{\nabla}_m \| \right)^2 &\leq \left( \sum_{m = 1}^{T/K} 1 \right) \left( \sum_{m = 1}^{T/K} \eta_m^2 \| \bar{\nabla}_m \|^2 \right) \\ &\leq \frac{2D^2 T}{K} \cdot \frac{T}{K}. \end{align}\] Therefore, if we allow \(\epsilon=0\), i.e., if we only assume \(\epsilon \geq 0\), the number of calls is bounded by \[\begin{align} &\frac{1}{\delta^2r^2} \sum_{m = 1}^{T/K - 1} \left( \eta_m^2 \| \nabla_m \|^2 + 2D \delta \eta_m \| \nabla_m \| \right) + \frac{D^2}{r^2} + \frac{T}{K} \\ &\quad\leq \frac{2D^2}{\delta^2r^2} \frac{T}{K} + \frac{2 \sqrt{2} D^2}{r^2 \delta} \frac{T}{K} + \frac{D^2}{r^2} + \frac{T}{K} \\ &\quad= \mathcal{O}\left( (D/r)^2\delta^{-2} T K^{-1} + (D/r)\delta^{-1} T K^{-1} + T K^{-1} \right) \\ &\quad= \mathcal{O}\left((D/r)^2 T K^{-1} \delta^{-2} \right). \qedhere \end{align}\] ◻
Theorem 3. Assume functions \(f_t\) are \(G_t\)-Lipschitz continuous and \(\theta\)-strongly convex for some constant \(\theta > 0\). Then Algorithm [alg:general] ensures that the regret is bounded as \[\begin{align} \mathtt{Regret}(T) \leq \frac{1}{2 \theta} \sum\limits_{m=1}^{T/K} \sum\limits_{t\in\mathcal{T}_m} \frac{\Vert \nabla_t \Vert^2}{m} + \frac{\delta }{r} D \sum\limits_{t=1}^TG_t, \end{align}\] where \(\delta\), \(K\), and \(D\) are defined as in Theorem 1. Moreover, if we assume all \(G_t \leq M_1\) for some constant \(M_1 \geq 0\), then this requires \(\mathcal{O}\left( \delta^{-2} + \delta^{-1}\log(T) + T K^{-1} \right)\) calls to the separation oracle.
Proof. Let \(\mathbf{x}^*\in \operatorname{argmin}_{\mathbf{x}\in\mathcal{K}}\sum_{t=1}^T f_t(\mathbf{x})\) and \(\tilde{\mathbf{x}}^*=(1-\frac{\delta}{r})\mathbf{x}^*+\frac{\delta}{r}{\mathbf{c}} \in \mathcal{K}_\delta\). We may decompose the regret in the same manner as Equation 1 to see that \[\begin{align} \mathtt{Regret}(T) &= \sum\limits_{t=1}^T\left( f_t(\mathbf{x}_t) - f_t(\mathbf{x}^*) \right) \\ &= \sum\limits_{t=1}^T\left( f_t(\mathbf{x}_t) - f_t(\tilde{\mathbf{x}}^*) \right) \\ &+ \sum\limits_{t=1}^T\left( f_t(\tilde{\mathbf{x}}^*) - f_t(\mathbf{x}^*) \right). \end{align}\] Since we assume \(f_t\) to be \(G_t\)-Lipschitz continuous, similar to Equation 4 , we have \[\begin{align} \label{eq:strconv-lip-cons} \sum\limits_{t=1}^T\left( f_t(\tilde{\mathbf{x}}^*) - f_t(\mathbf{x}^*) \right) \leq \frac{\delta D}{r} \sum\limits_{t=1}^TG_t. \end{align}\tag{6}\] The infeasible projection operator ensures that \(\Vert \mathbf{x}_{t+1}-\tilde{\mathbf{x}}^*\Vert \leq \Vert \mathbf{y}_t - \tilde{\mathbf{x}}^*\Vert\), extending Theorem 4.1 by [25], we have \[\begin{align} \sum\limits_{t=1}^T\left( f_t(\mathbf{x}_t) - f_t(\tilde{\mathbf{x}}^*) \right) \leq \frac{1}{2} \sum\limits_{m=1}^{T/K} \sum\limits_{t\in\mathcal{T}_m} \frac{G_t^2}{m \theta} \end{align}\] Thus, we conclude that \[\begin{align} \mathtt{Regret}(T) \leq \frac{1}{2 \theta} \sum\limits_{m=1}^{T/K} \sum\limits_{t\in\mathcal{T}_m} \frac{G_t^2}{m} + \frac{\delta D}{r} \sum\limits_{t=1}^TG_t \end{align}\] Note that if we assume all \(G_t\leq M_1\), thus all functions are \(M_1\)-Lipschitz, then \[\begin{align} \mathtt{Regret}(T) &\leq \frac{M_1^2 K}{2 \theta} \sum\limits_{m=1}^{T/K} \frac{1}{m} + \frac{\delta D M_1 T}{r} \\ &\leq \frac{M_1^2 K}{2 \theta} (1+\log(\frac{T}{K})) + \frac{\delta D M_1 T}{r} \end{align}\] where the last is due to the sum of harmonic series. If \(\delta = \Theta ( T^{-\beta}\operatorname{log}(T))\), \(K=\Theta(T^{1-\beta})\), then \(\mathtt{Regret}(T)=\mathcal{O}(T^{1-\beta}\log(T))\).
Next we examine the number of calls to the separation oracles. Similar to the convex case, the number of calls to the separation oracle is bounded by Equation 5 : \[\begin{align} &\sum_{m = 1}^{T/K - 1}\left( \frac{\operatorname{dist}(\mathbf{x}_m - \eta_m \bar{\nabla}_m,\mathcal{K}_\delta)^2-\operatorname{dist}(\mathbf{x}_{m+1},\mathcal{K}_\delta)^2 }{\delta^2}+1 \right) \nonumber \\ &\quad\leq \frac{1}{\delta^2} \sum_{m = 1}^{T/K} \left( \eta_m^2 \| \bar{\nabla}_m \|^2 + \frac{2 D \delta}{r} \eta_m \| \bar{\nabla}_m \| \right) + \frac{D}{r \delta} + \frac{T}{K}, \end{align}\] but for the strongly convex functions, \(\eta_m=\frac{1}{m\theta}\). Thus, if we assume all Lipschitz parameters are bounded by \(M_1\), then we have \[\begin{align} \sum_{m = 1}^{T/K} \eta_m^2 \| \bar{\nabla}_m \|^2 \leq \frac{M_1^2}{\theta^2} \sum_{m = 1}^{T/K} \frac{1}{m^2} \leq \frac{2 M_1^2 }{\theta^2} \end{align}\] where the last inequality follows from the fact that \(\sum_{i = 1}^\infty i^{-2} < 2\). On the other hand, \[\begin{align} \sum_{m = 1}^{T/K} \eta_m \| \bar{\nabla}_m \| \leq \frac{M_1}{\theta} \sum_{m = 1}^{T/K} \frac{1}{m} \leq \frac{M_1}{\theta}\left( \log\left( \frac{T}{K} \right) + 1 \right) \end{align}\] where the last inequality follows from the fact that \(\sum_{i = 1}^N i^{-1} < \log(N) + 1\).
Combining the results, we see that if all functions are \(M_1\)-Lipchitz, we have the number of calls to the separation oracles may be bounded as \[\begin{align} &\sum_{m = 1}^{T/K - 1}\left( \frac{\operatorname{dist}(\mathbf{x}_m - \eta_m \bar{\nabla}_m,\mathcal{K}_\delta)^2-\operatorname{dist}(\mathbf{x}_{m+1},\mathcal{K}_\delta)^2 }{\delta^2}+1 \right) \nonumber \\ &\quad\leq \frac{1}{\delta^2} \left( \frac{2 M_1^2}{\theta^2} + \frac{2 D \delta}{r} \frac{M_1}{\theta}\left( \log\left( \frac{T}{K} \right) + 1 \right) \right) + \frac{D}{r \delta} + \frac{T}{K} \\ &\quad = \mathcal{O}\left( \delta^{-2} + \delta^{-1}\log(T) + \delta^{-1} + T K^{-1} \right) \\ &\quad = \mathcal{O}\left( \delta^{-2} + \delta^{-1}\log(T) + T K^{-1} \right). \qedhere \end{align}\] ◻
Corollary 2. In Theorem 3, if we assume all \(G_t\leq M_1\) for some constant \(M_1\), then we have \[\begin{align} \mathtt{Regret}(T) \leq \frac{M_1^2 K}{2 \theta} (1+\log(\frac{T}{K})) + \frac{\delta D M_1 T}{r}. \end{align}\] In particular, if we choose \(\delta=\Theta(T^{-\beta}\operatorname{log}(T))\) and \(K=\Theta(T^{1-\beta})\) for some \(\beta \in (0, 1]\), then \(\mathtt{Regret}(T)= \mathcal{O}(T^{1-\beta}\log(T))\) with \(\mathcal{O}(T^{2\beta}(\log(T))^{-2})\) calls to the separation oracle.
The proof proceeds by connecting the performance of \(\mathtt{BAGEL}\) on the COCO problem to the performance of its base OCO algorithm (Theorem 1) on the sequence of surrogate costs. The first step is bounding surrogate regret. We apply Theorem 1 directly to the sequence of surrogate functions \(\{\hat{f}_t\}\) to obtain a bound on the surrogate regret, \(\mathtt{R}(\hat{f})\). This bound depends on the norms of the surrogate gradients. The second step is bounding surrogate gradients. We then bound the surrogate gradient norms in terms of the Lyapunov function’s derivative, \(\Phi'(Q_T)\), which effectively captures the worst-case cumulative violation. The third step is relating surrogate and true regret. We leverage the definition of \(\hat{f}_t\) and the convexity of \(\Phi\) to establish the key inequality: \(\mathtt{R}(\hat{f})\geq \mathtt{R}(\tilde{f})+\Phi(Q_T)\). This inequality forms a bridge, linking the (boundable) surrogate regret to the two quantities we ultimately wish to control: the true regret \(\mathtt{R}(\tilde{f})\) and the total violation (captured by \(\Phi(Q_T)\)). Step 4: Parameter Tuning. By combining these inequalities and carefully choosing our hyperparameters (\(\lambda,\gamma, \mathcal{K}, \delta\)) to balance the resulting terms, we can solve for the final bounds on \(\mathtt{Regret}(T)\) and \(\mathtt{CCV}(T)\).
We start with bounding \(\mathtt{R}(\hat{f})\):
Lemma 4. Assume surrogate functions \(\hat{f}_t\) are convex and \(G_t\)-Lipschitz, then Algorithm [alg:constrained] ensures that the term \(\mathtt{R}(\hat{f})\) defined as follows is bounded: \[\begin{align} \mathtt{R}(\hat{f}) &\triangleq \sum_{t=1}^T \hat{f}_t(\mathbf{x}_t) - \sum_{t=1}^T \hat{f}_t(\mathbf{x}^*) \\ &\leq \frac{3}{2}DK\sqrt{\sum_{m=1}^{T/K} \Vert \bar{\nabla}_m \Vert ^2} + D \delta \sum\limits_{t=1}^TG_t, \end{align}\] with \(\tilde{\mathcal{O}}(T^{2\beta})\) total calls to the SO, where \(\mathbf{x}^*=\arg\min_{\mathbf{x}\in\mathcal{K}}\sum_{t=1}^T f_t(\mathbf{x})\), \(\bar{\nabla}_m\) is defined in Algorithm [alg:general](line 9), \(D\) is the diameter of the action set \(\mathcal{K}\), \(K\) is the block size, \(\delta \in [0,1)\) is the shrinking parameter.
Proof. Let \(\mathbf{x}^*=\arg\min_{\mathbf{x}\in\mathcal{K}}\sum_{t=1}^T f_t(\mathbf{x})\) and \(\mathbf{x}'=\arg\min_{\mathbf{x}\in\mathcal{K}}\sum_{t=1}^T \hat{f}_t(\mathbf{x})\). Proof of Lemma 4 follows from that of Theorem 1, by replacing \(\mathbf{x}'\) with \(\mathbf{x}^*\). Note that \(\hat{f}_t\) and \(\tilde{f}_t\) share the same action set \(\mathcal{K}\). ◻
Next, we aim to bound \(\mathtt{R}(\tilde{f})\), using the relationship between \(\hat{f}_t\) and \(\tilde{f}_t\), and replacing the \(\Vert \bar{\nabla}_m \Vert ^2\) term in Lemma 4. The results are provided below in Lemma 5 with proof.
Lemma 5. Assume \(\Phi(\cdot):\mathbb{R}^+\mapsto\mathbb{R}^+\) is a non-decreasing, differentiable, and convex potential function satisfies \(\Phi(0)=0\). Assume cost functions \(f_t\) and constraint function \(g_t\) are \(M_1\)-Lipschitz. Define \(\mathtt{R}(\tilde{f}) \triangleq \sum_{t=1}^T \tilde{f}_t(\mathbf{x}_t) - \sum_{t=1}^T \tilde{f}_t(\mathbf{x}^*)\) and Algorithm [alg:constrained] ensures that \[\label{eq:sur2real:1} \Phi (Q_T) + \mathtt{R}(\tilde{f}) \leq \mathtt{R} (\hat{f})\tag{7}\] and \[\label{eq:sur2real:2} \mathtt{R} (\hat{f}) \leq \gamma M_1 D (1+\Phi'(Q_T)) (\delta T + \frac{3}{\sqrt{2}}\sqrt{TK})\tag{8}\] where \(\delta \in [0,1)\) is the shrinking parameter, \(K\) is the block size, \(D\) is the diameter of \(\mathcal{K}\), and the definition of \(\bar{\nabla}_m\) is given in Algorithm [alg:general] (line 9), \(\gamma\) is the processing parameter.
Proof. Note that \(\hat{f}_t(\mathbf{x}_t)=\tilde{f}_t(\mathbf{x}_t)+\Phi'(Q_t)\tilde{g}_t(\mathbf{x}_t)\), where \(\tilde{f}_t=\gamma f_t\), \(\tilde{g}_t=\gamma (g_t)^+\), and \(\Phi(\cdot)\) is a non-decreasing, differentiable, convex function satisfying \(\Phi(0)=0\). By convexity of \(\Phi(\cdot)\), we have \[\begin{align} \Phi(Q_t) &\leq \Phi(Q_{t-1}) + \Phi'(Q_t)(Q_t-Q_{t-1}) \\ &= \Phi(Q_{t-1}) + \Phi'(Q_t)\tilde{g}_t(\mathbf{x}_t). \end{align}\] Rearranging, we have \[\begin{align} \Phi(Q_t) - \Phi(Q_{t-1}) \leq \Phi'(Q_t)\tilde{g}_t(\mathbf{x}_t). \end{align}\] Since \(\hat{f}_t(\mathbf{x}^*)=\tilde{f}_t(\mathbf{x}^*)+\Phi'(Q_t)\tilde{g}_t(\mathbf{x}^*)\) and \(\tilde{g}_t(\mathbf{x}^*)=\gamma (g_t(\mathbf{x}^*))^+ \leq 0\), plugging in the above inequality, we have \[\begin{align} &\hat{f}_t(\mathbf{x}_t) - \hat{f}_t(\mathbf{x}^*) = \left(\tilde{f}_t(\mathbf{x}_t)-\tilde{f}_t(\mathbf{x}^*) \right) \\ & \qquad+ \Phi'(Q_t)(\tilde{g}_t(\mathbf{x}_t)-\tilde{g}_t(\mathbf{x}^*)) \\ & \geq \left(\tilde{f}_t(\mathbf{x}_t) - \tilde{f}_t(\mathbf{x}^*) \right) + \Phi'(Q_t)\tilde{g}_t(\mathbf{x}_t) \\ & \geq \left(\tilde{f}_t(\mathbf{x}_t) - \tilde{f}_t(\mathbf{x}^*) \right) + \Phi(Q_t) - \Phi(Q_{t-1}). \end{align}\] Summing the above inequality over \(t\in[T]\), with \(Q_0=0\) and \(\Phi(0)=0\), we have \[\begin{align} \mathtt{R}(\hat{f}) \geq \mathtt{R}(\tilde{f}) + \Phi(Q_T), \end{align}\] which gives Equation 7 in Lemma 5.
Since both \(f\) and \(g\) are assumed to be \(M_1\)-Lipschitz continuous, we have \[\begin{align} \label{eq:Gt} \Vert \nabla_t \Vert&= \Vert \nabla \hat{f}_t\Vert \leq \Vert \nabla \tilde{f}_t\Vert + \Phi'(Q_t) \Vert \nabla \tilde{g}_t\Vert \nonumber \\ &\leq \gamma M_1(1+\Phi'(Q_t)). \end{align}\tag{9}\] Thus, surrogate functions \(\hat{f}_t\) are \(G_t\)-Lipschitz continuous, where \(G_t = \gamma M_1(1+\Phi'(Q_t))\). Since \(\bar{\nabla}_m=\frac{1}{K} \sum_{t\in\mathcal{T}_m} \nabla_t\), using Cauchy-Schwartz inequality, we have \(\Vert \bar{\nabla}_m \Vert ^2 = \frac{1}{K^2} \Vert \sum_{t\in\mathcal{T}_m} \nabla_t \Vert ^2 \leq \frac{1}{K} \sum_{t\in\mathcal{T}_m}\Vert \nabla_t \Vert^2\). Thus, \[\begin{align} \label{eq:sqrt-sum-gm} \sqrt{\sum_{m=1}^{T/K} \Vert \bar{\nabla}_m \Vert ^2} &\leq \sqrt{\frac{1}{K} \sum_{m=1}^{T/K} \sum_{t\in\mathcal{T}_m}\Vert \nabla_t \Vert^2} \nonumber\\ &= \sqrt{ \frac{1}{K} \sum_{t=1}^{T} \Vert \nabla_t \Vert^2 } \nonumber \\ &\overset{(a)}\leq \frac{\gamma M_1}{\sqrt{K}}\sqrt{ \sum_{t=1}^{T} (1+\Phi'(Q_t))^2 } \nonumber \\ &\overset{(b)}\leq \frac{\gamma M_1}{\sqrt{K}}\sqrt{ \sum_{t=1}^{T} 2(1+\Phi'(Q_t)^2) } \nonumber \\ &\overset{(c)}\leq \frac{\gamma M_1 \sqrt{2}}{\sqrt{K}} (\sqrt{T}+\sqrt{\sum_{t=1}^{T}\Phi'(Q_t)^2})\nonumber \\ &\overset{(d)}\leq \frac{\gamma M_1 \sqrt{2T}}{\sqrt{K}} (1 + \Phi'(Q_T)) \end{align}\tag{10}\] where (a) follows from Equation 9 and the fact that \(G_t=\gamma M_1(1+\Phi'(Q_t))\), (b) follows from Cauchy-Schwartz inequality, (c) follows from Jensen’s inequality, and (d) is because \(\Phi'(\cdot)\) is non-decreasing as \(\Phi(\cdot)\) is convex and \(Q_t\) is non-decreasing, thus \(\Phi'(Q_t) \leq \Phi'(Q_T), \forall t\in[T]\). By plugging Equation 9 and Equation 10 into Lemma 4 and rearranging terms, we have Equation 8 in Lemma 5. ◻
Having necessary lemmas at hand, we now move towards parameter tuning. If we set \(\Phi(Q_T)=e^{\lambda Q_T}-1\) with some real number \(\lambda\), we have \(\Phi'(Q_T)=\lambda e^{\lambda Q_T}\). Letting \(S \triangleq \delta T + \frac{3}{\sqrt{2}}\sqrt{TK}\), we have \[\begin{align} \label{eq:regretTft} &\mathtt{R}(\tilde{f}) \overset{(a)}{\leq} \mathtt{R} (\hat{f})-\Phi (Q_T) \nonumber \\ &\quad\overset{(b)}{\leq} \gamma M_1 D (V+\lambda e^{\lambda Q_T})(\delta T + \frac{3}{\sqrt{2}}\sqrt{TK})-e^{\lambda Q_T}+1\nonumber\\ &\quad\overset{(c)} = \gamma M_1 D S + e^{\lambda Q_T} \left( \lambda \gamma M_1 D S-1 \right)+1 \end{align}\tag{11}\] where (a) uses Equation 7 , (b) uses Equation 8 , and (c) rearranges terms and adopts \(S\).
We make a few observation of Equation 11 . Note that if we take \(\gamma=(M_1D)^{-1}\), we can offset \(\gamma M_1 D\) term entirely. If we choose \(\lambda=\frac{1}{2S}\), then \(\lambda \gamma M_1 D S -1 = -\frac{1}{2}<0\). Thus, based on Equation 11 , we have \[\begin{align} &\mathtt{R}(\tilde{f}) \leq S+1= \delta T + \frac{3}{\sqrt{2}}\sqrt{TK}+1. \end{align}\] Since \(\tilde{f}_t=\gamma f_t\), we have \(\mathtt{R}(\tilde{f}) = \gamma\mathtt{Regret}(T)\). Let \(\beta\in[0,\frac{1}{2}]\) be some real number. Thus, if we choose \(\delta=\Theta(T^{-\beta})\) and \(K=\Theta(T^{1-2\beta})\), \[\begin{align} \mathtt{Regret}(T) &= \gamma^{-1} \mathtt{R}(\tilde{f}) \\ &= D M_1 \delta T + \frac{3}{\sqrt{2}}DM_1\sqrt{TK}+DM_1 \\ &= \mathcal{O}(DM_1 T^{1-\beta}). \end{align}\]
We now take a closer look at the Cumulative Constraint Violation (CCV). With \(\tilde{f}_t=\gamma f_t\) and \(f_t\) being \(M_1\)-Lipschitz, we have \(\tilde{f}_t\) is \(\gamma M_1\)-Lipscitz continuous, and thus \(\tilde{f}_t(\mathbf{x}_t)-\tilde{f}_t(\mathbf{x}^*) \geq -\gamma M_1 (\mathbf{x}_t-\mathbf{x}^*) \geq -1\). Taking sum over \(t\in[T]\), we have \[\begin{align} \label{eq:lowbon4regret} \mathtt{R}(\tilde{f})=\sum\limits_{t=1}^T(\tilde{f}_t(\mathbf{x}_t)-\tilde{f}_t(\mathbf{x}^*)) \geq -T \end{align}\tag{12}\] Replacing \(\mathtt{R}(\tilde{f})\) in Equation 11 with the lower bound in Equation 12 and rearranging terms, we have \[\begin{align} e^{\lambda Q_T} \leq \frac{VS + 1 + VT }{1-\lambda S}. \end{align}\] Since \(\lambda=\frac{1}{2S}\), we have \[\begin{align} Q_T &\leq 2S \log\left( 2(S + \frac{1}{V} + T)\right) \end{align}\] Considering that \(Q_T=\sum\limits_{t=1}^T(\tilde{g}_t)^+=\sum\limits_{t=1}^T\gamma(g_t)^+=\gamma \mathtt{CCV}(T)\), with \(S = \delta T + \frac{3}{\sqrt{2}}\sqrt{TK}\) and \(\gamma=(DM_1)^{-1}\), if we choose \(\delta=\Theta(T^{-\beta})\) and \(K=\Theta(T^{1-2\beta})\) for some \(\beta\in(0,\frac{1}{2}]\), we have \[\begin{align} \mathtt{CCV}(T) &= \gamma^{-1} Q_T \leq 2DM_1S * \mathcal{O}(\log T) \\ &\leq 2 D M_1 (\delta T + \frac{3}{\sqrt{2}}\sqrt{TK}) * \mathcal{O}(\log T) \\ &\leq \mathcal{O}(DM_1T^{1-\beta}\log T) \end{align}\]
Since calls to the infeasible projection and separation oracles are done in \(\mathcal{A}\), we can extend results from Corollary 1 and derive that the oracle calls are bounded as \(\tilde{\mathcal{O}}(T^{2\beta})\).
We note that the regret bounds can improve when the cost functions are strongly convex. In this section, we describe how the proposed algorithm, \(\mathtt{BAGEL}\), can be implemented when the cost functions are strongly convex, and analyze the corresponding \(\mathtt{Regret}\), \(\mathtt{CCV}\), and total number of calls to the separation oracle guarantees.
Let the cost functions \(f_t\) be \(\theta\)-strongly convex while the constraint functions are considered to be convex, not necessarily strongly convex. Let them all be \(M_1\)-Lipschitz continuous. In this case, we still apply \(\mathtt{BAGEL}\) (Algorithm [alg:constrained]) with different parameters. In line 10, we add a regularization parameter \(V\) to control the influence of costs on the surrogate function, i.e., we let \(\hat{f}_t \gets V\,\tilde{f}_t + \Phi'(Q_t)\,\tilde{g}_t\). Besides, we choose \(\gamma = 1\), \(V = \frac{8M_1^2 K \log(Te/K)}{\theta}\), and \(\Phi(\cdot) = (\cdot)^2\). Further, the Base Algorithm instance \(\mathcal{A}\) uses step size \(\eta_m = \frac{1}{m\theta}\). For clarity, we include the updated BAGEL algorithm for strongly-convex costs in Algorithm [alg:strongly-convex].
Action set \(\mathcal{K}\), time horizon \(T\), \(\mathbf{c} \in \operatorname{relint}(\mathcal{K})\), \(r = r_{\mathcal{K}, \mathbf{c}}\), shrinking parameter \(\delta \in [0,1)\), block size \(K\), action set diameter \(D\), an instance of the Base Algorithm (Algorithm [alg:general]) \(\mathcal{A}\), processing parameter \(\gamma\), regularization parameter \(V\), Lyapunov function \(\Phi(\cdot)\), an instance of IP-SO (Algorithm 1) \(\mathcal{P}\), strongly-convex parameter \(\theta\). Initialization: Pick \(\mathbf{x}_1\in\mathcal{K}\); set \(Q_0\gets 0\); pass \(\mathcal{K},T,K,\mathbf{c},r,\delta,D,\mathcal{P}\), and \((\epsilon\text{ or }\theta)\) to \(\mathcal{A}\) Play \(\mathbf{x}_m\) Adversary reveals \(f_t\) and \(g_t\) Observe \(f_t(\mathbf{x}_m)\) and \((g_t(\mathbf{x}_m))^{+}\) \(\tilde{f}_t \gets \gamma\, f_t\); \(\tilde{g}_t \gets \gamma\, (g_t)^{+}\) \(Q_t \gets Q_{t-1} + (g_t(\mathbf{x}_m))^{+}\) \(\hat{f}_t \gets V\,\tilde{f}_t + \Phi'(Q_t)\,\tilde{g}_t\); \(\nabla_t \gets \nabla \hat{f}_t(\mathbf{x}_m)\) Pass \(\nabla_t\) to \(\mathcal{A}\) Receive \(\mathbf{x}_{m+1}\) from \(\mathcal{A}\)
Having studied the regret for OCO for strongly convex functions in Section 4, we now use the result to study the regret and CCV for the adversarial COCO problem for strongly convex costs, with similar approach to the proofs in Section 5. The results are given in Theorem 4, and the detailed proof can be found in Appendix 12.
Theorem 4. Assume cost functions \(f_t\) are \(\theta\)-strongly convex and constraint functions \(g_t\) are convex, and they are all \(M_1\)-Lipschitz continuous. With \(\beta\in (0,1]\), if we choose \(\gamma=1\), \(\delta= \Theta( T^{-\beta} \log(T) )\), \(K= \Theta(T^{1-\beta})\), \(V=\frac{8 M_1^2 K \log(Te/K)}{\theta}\), and \(\Phi(\cdot)=(\cdot)^2\), Algorithm [alg:constrained] ensures that \[\begin{align} \mathtt{Regret}(T) &= \mathcal{O}(T^{1-\beta}\log(T)), \\ \mathtt{CCV}(T) &= \mathcal{O}(T^{1-\frac{\beta}{2}}\sqrt{\log(T)}), \end{align}\] with at most \(\mathcal{O}(T^{2\beta}(\log(T))^{-2})\) calls to the separation oracle.
Remark 5. If we take \(\beta=1\), we obtain optimal results \(\log(T)\) for \(\mathtt{Regret}\) and \(\sqrt{T\log(T)}\) for \(\mathtt{CCV}\) if permitted \(\mathcal{O}(T^2)\) calls to the separation oracle, which matches the optimal result for projection-based algorithm in [6]. If only \(\mathcal{O}(T)\) calls permitted, we take \(\beta=\frac{1}{2}\), and still have \(\mathcal{O}(\sqrt{T})\) regret and \(\mathcal{O}(T^{\frac{3}{4}})\) CCV, which shows significant improvement compared to \(\mathcal{O}(T^\frac{3}{4})\) regret and \(\mathcal{O}(T^{\frac{7}{8}})\) CCV in [7].
Proof. We observe that Equation 7 still holds when \(f_t\) are strongly convex. Similar to proof of Theorem 2, we first derive bound on \(\mathtt{R}(\hat{f})\), and move to \(\mathtt{R}(\tilde{f})\), before obtaining results on \(\mathtt{Regret}(T)\). Then we bound regret from below and choose parameters appropriately to obtain results for \(\mathtt{CCV}(T)\).
If \(\gamma=1\), we have \(\tilde{f}_t=f_t\), \(\tilde{g}_t=(g_t)^+\). Thus, \(\hat{f}_t(\mathbf{x}_t)=V\tilde{f}_t(\mathbf{x}_t)+\Phi'(Q_t)\tilde{g}_t=Vf_t+\Phi'(Q_t)(g_t)^+\), where \(V\) is a positive real number, and \(\Phi(\cdot)\) is a Lyapunov function as described in Section 5. Since \(g_t\) are convex, \((g_t)^+\) are also convex. Since \(\Phi(\cdot)\) is non-decreasing, \(\Phi'(Q_t)\) is non-negative. As \(f_t\) are \(\theta\)-strongly convex, we can apply Lemma 2 and derive that \(\hat{f}_t(\mathbf{x}_t)\) are \(V\theta\)-strongly convex. Since both \(f\) and \(g\) are assumed to be \(M_1\)-Lipschitz continuous, Equation 9 still holds: \(\Vert \nabla_t \Vert\leq \gamma M_1(V+\Phi'(Q_t))\). Thus, surrogate functions \(\hat{f}_t\) are \(G_t\)-Lipschitz continuous, where \(G_t = \gamma M_1(V+\Phi'(Q_t))\). Similar to how we extend Theorem 1 to get Lemma 4, we can extend Theorem 3 and get: \[\begin{align} \mathtt{R}(\hat{f}) & \overset{(a)}\leq \frac{1}{2V\theta}\sum_{m=1}^{T/K}\sum_{t\in\mathcal{T}_m} \frac{\Vert \nabla_t \Vert^2}{m} + \frac{\delta D}{r} \sum\limits_{t=1}^TG_t\\ & \overset{(b)} \leq \frac{1}{2V\theta} \sum_{m=1}^{T/K}\sum_{t\in\mathcal{T}_m} \frac{M_1^2(V+\Phi'(Q_t))^2}{m} \\ & \qquad+ \frac{\delta D}{r} \sum\limits_{t=1}^TM_1(V+\Phi'(Q_t)) \\ & \overset{(c)} \leq \frac{VM_1^2K}{\theta}(1+\log(\frac{T}{K})) \\ & \qquad+\frac{M_1^2}{V\theta } \sum_{m=1}^{T/K}\sum_{t\in\mathcal{T}_m} \frac{(\Phi'(Q_t))^2}{m} \\ & \qquad+ \frac{\delta D}{r} M_1 \sum\limits_{t=1}^T(V+\Phi'(Q_t)) \end{align}\] where (a) extends Theorem 3 and uses the fact that \(\hat{f}_t(\mathbf{x}_t)\) are \(V\theta\)-strongly convex, (b) follows from Equation 9 and the fact that \(G_t = \gamma M_1(V+\Phi'(Q_t))\), and (c) uses the \(1+\log(n)\) bound on sum of harmonic series with \(n\) terms. Since Equation 7 still holds, taking \(\Phi(\cdot)=(\cdot)^2\), we have \[\begin{align} &Q_T^2+V\mathtt{Regret}(T) \leq \mathtt{R}(\hat{f}) \leq \frac{VM_1^2K}{\theta}(1+\log(\frac{T}{K})) \\ & \qquad + \frac{4M_1^2}{\theta V } \sum_{m=1}^{T/K}\sum_{t\in\mathcal{T}_m} \frac{Q_t^2}{m} + \frac{\delta D}{r} M_1 \sum\limits_{t=1}^T(V+2Q_t). \end{align}\] Since \(Q_t\) is non-decreasing, we have \(Q_T \geq Q_t \geq 0, \forall t\in[T]\). With the sum of harmonic series \(\sum\limits_{t=1}^T\frac{1}{t} = \log(Te)\), and \(V\) being a positive real number, we have \[\begin{align} \label{eq:stronvex-reg} &\mathtt{Regret}(T) \leq \frac{M_1^2K}{\theta}(1+\log(\frac{T}{K})) \nonumber \\ & \quad +\frac{4M_1^2KQ_T^2 \log(Te/K)}{\theta V^2} + \frac{\delta D}{r} M_1 T (1+\frac{2Q_T}{V}) - \frac{Q_T^2}{V} \nonumber\\ & \leq \frac{M_1^2K}{\theta}(1+\log(\frac{T}{K})) +\frac{4M_1^2KQ_T^2 \log(Te/K)}{\theta V^2} \nonumber \\ & \quad + \frac{\delta D}{r} M_1 T (1+\frac{2Q_T}{V}). \end{align}\tag{13}\] If we choose \(\delta=\Theta(T^{-\beta} \log(T))\) and \(K=\Theta(T^{1-\beta})\), where \(\beta \in (0,1]\) is a trade-off parameter,we have \[\begin{align} \mathtt{Regret}(T) = \mathcal{O}(T^{1-\beta}\log(T)), \end{align}\] and hence we prove the result for \(\mathtt{Regret}\) in Theorem 4.
Next we take a look at \(\mathtt{CCV}(T)\). Given that the cost functions are \(M_1\)-Lipschitz continuous, we have \[\begin{align} \mathtt{Regret}(T)&=\sum\limits_{t=1}^T(f_t(\mathbf{x}_t)-f_t(\mathbf{x}^*)) \\ &\geq \sum\limits_{t=1}^T-M_1\Vert \mathbf{x}_t - \mathbf{x}^*\Vert \\ &\geq -M_1DT. \end{align}\] Replacing \(\mathtt{Regret}(T)\) in Equation 13 with the above lower bound, and choose \(V=\frac{8M_1^2K\log(Te/K)}{\theta}\), we have, \[\begin{align} \frac{1}{2}Q_T^2-\frac{2\delta D}{r}M_1 T Q_T &\leq V(\frac{\delta D}{r} M_1 T + D M_1 T \\ &\quad+ \frac{M_1K}{\theta}(1+\log(\frac{T}{K}))). \end{align}\] Let \(B_T\triangleq V(\frac{\delta D}{r} M_1 T + D M_1 T + \frac{M_1K}{\theta}(1+\log(\frac{T}{K})))\), then we have \[\begin{align} Q_T^2- \frac{4\delta D}{r} M_1 T Q_T - 2B_T \leq 0. \end{align}\] Note that the left hand side is a quadratic function of \(Q_T\). Solving the inequality, we have \[\begin{align} Q_T \leq \frac{\frac{4\delta D}{r} M_1 T + \sqrt{(\frac{4\delta D}{r} M_1 T)^2 + 8B_T}}{2}. \end{align}\] Since \(B_T = \mathcal{O}(T^{2-\beta}\log(T))\), we have \(Q_T = \mathcal{O}(T^{1-\frac{\beta}{2}}\sqrt{\log(T)})\), and \(\mathtt{CCV}(T)= \sum\limits_{t=1}^T(g_t)^+ = Q_T = \mathcal{O}(T^{1-\frac{\beta}{2}}\sqrt{\log(T)})\) for \(\theta\)-strongly convex cost.
Since calls to the infeasible projection and separation oracles are done in \(\mathcal{A}\), we can directly extend results from Corollary 2 and bound oracle calls as \(\mathcal{O}(T^{2\beta}(\log(T))^{-2}))\). ◻
For simplicity, in Table ¿tbl:tbl:comparison?, we only compare results with respect to time horizon \(T\); for more detailed discussion, we refer to our theorem statements and remarks.↩︎