Source author record

Kyriakos Axiotis

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

6works
5topics
4close collaborators

Published work

6 publications

preprint2022arXiv

Iterative Hard Thresholding with Adaptive Regularization: Sparser Solutions Without Sacrificing Runtime

We propose a simple modification to the iterative hard thresholding (IHT) algorithm, which recovers asymptotically sparser solutions as a function of the condition number. When aiming to minimize a convex function f(x)f(x) with condition number κκ subject to xx being an ss-sparse vector, the standard IHT guarantee is a solution with relaxed sparsity O(sκ2)O(sκ^2), while our proposed algorithm, regularized IHT, returns a solution with sparsity O(sκ)O(sκ). Our algorithm significantly improves over ARHT which also finds a solution of sparsity O(sκ)O(sκ), as it does not require re-optimization in each iteration (and so is much faster), is deterministic, and does not require knowledge of the optimal solution value f(x∗)f(x^*) or the optimal sparsity level ss. Our main technical tool is an adaptive regularization framework, in which the algorithm progressively learns the weights of an ℓ2\ell_2 regularization term that will allow convergence to sparser solutions. We also apply this framework to low rank optimization, where we achieve a similar improvement of the best known condition number dependence from κ2κ^2 to κκ.

preprint2021arXiv

Decomposable Submodular Function Minimization via Maximum Flow

This paper bridges discrete and continuous optimization approaches for decomposable submodular function minimization, in both the standard and parametric settings. We provide improved running times for this problem by reducing it to a number of calls to a maximum flow oracle. When each function in the decomposition acts on O(1)O(1) elements of the ground set VV and is polynomially bounded, our running time is up to polylogarithmic factors equal to that of solving maximum flow in a sparse graph with O(∣V∣)O(\vert V \vert) vertices and polynomial integral capacities. We achieve this by providing a simple iterative method which can optimize to high precision any convex function defined on the submodular base polytope, provided we can efficiently minimize it on the base polytope corresponding to the cut function of a certain graph that we construct. We solve this minimization problem by lifting the solutions of a parametric cut problem, which we obtain via a new efficient combinatorial reduction to maximum flow. This reduction is of independent interest and implies some previously unknown bounds for the parametric minimum s,ts,t-cut problem in multiple settings.

preprint2021arXiv

Local Search Algorithms for Rank-Constrained Convex Optimization

We propose greedy and local search algorithms for rank-constrained convex optimization, namely solving min⁡rank(A)≤r∗ R(A)\underset{\mathrm{rank}(A)\leq r^*}{\min}\, R(A) given a convex function R:Rm×n→RR:\mathbb{R}^{m\times n}\rightarrow \mathbb{R} and a parameter r∗r^*. These algorithms consist of repeating two steps: (a) adding a new rank-1 matrix to AA and (b) enforcing the rank constraint on AA. We refine and improve the theoretical analysis of Shalev-Shwartz et al. (2011), and show that if the rank-restricted condition number of RR is κκ, a solution AA with rank O(r∗⋅min⁡{κlog⁡R(0)−R(A∗)ε,κ2})O(r^*\cdot \min\{κ\log \frac{R(\mathbf{0})-R(A^*)}ε, κ^2\}) and R(A)≤R(A∗)+εR(A) \leq R(A^*) + ε can be recovered, where A∗A^* is the optimal solution. This significantly generalizes associated results on sparse convex optimization, as well as rank-constrained convex optimization for smooth functions. We then introduce new practical variants of these algorithms that have superior runtime and recover better solutions in practice. We demonstrate the versatility of these methods on a wide range of applications involving matrix completion and robust principal component analysis.

preprint2020arXiv

Circulation Control for Faster Minimum Cost Flow in Unit-Capacity Graphs

We present an m4/3+o(1)log⁡Wm^{4/3+o(1)}\log W-time algorithm for solving the minimum cost flow problem in graphs with unit capacity, where WW is the maximum absolute value of any edge weight. For sparse graphs, this improves over the best known running time for this problem and, by well-known reductions, also implies improved running times for the shortest path problem with negative weights, minimum cost bipartite b\boldsymbol{\mathit{b}}-matching when ∥b∥1=O(m)\|\boldsymbol{\mathit{b}}\|_1 = O(m), and recovers the running time of the currently fastest algorithm for maximum flow in graphs with unit capacities (Liu-Sidford, 2020). Our algorithm relies on developing an interior point method-based framework which acts on the space of circulations in the underlying graph. From the combinatorial point of view, this framework can be viewed as iteratively improving the cost of a suboptimal solution by pushing flow around circulations. These circulations are derived by computing a regularized version of the standard Newton step, which is partially inspired by previous work on the unit-capacity maximum flow problem (Liu-Sidford, 2019), and subsequently refined based on the very recent progress on this problem (Liu-Sidford, 2020). The resulting step problem can then be computed efficiently using the recent work on ℓp\ell_p-norm minimizing flows (Kyng-Peng-Sachdeva-Wang, 2019). We obtain our faster algorithm by combining this new step primitive with a customized preconditioning method, which aims to ensure that the graph on which these circulations are computed has sufficiently large conductance.

preprint2020arXiv

Sparse Convex Optimization via Adaptively Regularized Hard Thresholding

The goal of Sparse Convex Optimization is to optimize a convex function ff under a sparsity constraint s≤s∗γs\leq s^*γ, where s∗s^* is the target number of non-zero entries in a feasible solution (sparsity) and γ≥1γ\geq 1 is an approximation factor. There has been a lot of work to analyze the sparsity guarantees of various algorithms (LASSO, Orthogonal Matching Pursuit (OMP), Iterative Hard Thresholding (IHT)) in terms of the Restricted Condition Number κκ. The best known algorithms guarantee to find an approximate solution of value f(x∗)+εf(x^*)+ε with the sparsity bound of γ=O(κmin⁡{log⁡f(x0)−f(x∗)ε,κ})γ= O\left(κ\min\left\{\log \frac{f(x^0)-f(x^*)}ε, κ\right\}\right), where x∗x^* is the target solution. We present a new Adaptively Regularized Hard Thresholding (ARHT) algorithm that makes significant progress on this problem by bringing the bound down to γ=O(κ)γ=O(κ), which has been shown to be tight for a general class of algorithms including LASSO, OMP, and IHT. This is achieved without significant sacrifice in the runtime efficiency compared to the fastest known algorithms. We also provide a new analysis of OMP with Replacement (OMPR) for general ff, under the condition s>s∗κ24s > s^* \frac{κ^2}{4}, which yields Compressed Sensing bounds under the Restricted Isometry Property (RIP). When compared to other Compressed Sensing approaches, it has the advantage of providing a strong tradeoff between the RIP condition and the solution sparsity, while working for any general function ff that meets the RIP condition.

preprint2016arXiv

On the Size and the Approximability of Minimum Temporally Connected Subgraphs

We consider temporal graphs with discrete time labels and investigate the size and the approximability of minimum temporally connected spanning subgraphs. We present a family of minimally connected temporal graphs with nn vertices and Ω(n2)Ω(n^2) edges, thus resolving an open question of (Kempe, Kleinberg, Kumar, JCSS 64, 2002) about the existence of sparse temporal connectivity certificates. Next, we consider the problem of computing a minimum weight subset of temporal edges that preserve connectivity of a given temporal graph either from a given vertex r (r-MTC problem) or among all vertex pairs (MTC problem). We show that the approximability of r-MTC is closely related to the approximability of Directed Steiner Tree and that r-MTC can be solved in polynomial time if the underlying graph has bounded treewidth. We also show that the best approximation ratio for MTC is at least O(2log⁡1−εn)O(2^{\log^{1-ε} n}) and at most O(min⁡{n1+ε,(ΔM)2/3+ε})O(\min\{n^{1+ε}, (ΔM)^{2/3+ε}\}), for any constant ε>0ε> 0, where MM is the number of temporal edges and ΔΔ is the maximum degree of the underlying graph. Furthermore, we prove that the unweighted version of MTC is APX-hard and that MTC is efficiently solvable in trees and 22-approximable in cycles.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.