skip to main content


Title: Sum of squares generalizations for conic sets
Abstract

Polynomial nonnegativity constraints can often be handled using thesum of squarescondition. This can be efficiently enforced using semidefinite programming formulations, or as more recently proposed by Papp and Yildiz (Papp D in SIAM J O 29: 822–851, 2019), using the sum of squares cone directly in an interior point algorithm. Beyond nonnegativity, more complicated polynomial constraints (in particular, generalizations of the positive semidefinite, second order and$$\ell _1$$1-norm cones) can also be modeled through structured sum of squares programs. We take a different approach and propose using more specialized cones instead. This can result in lower dimensional formulations, more efficient oracles for interior point methods, or self-concordant barriers with smaller parameters.

 
more » « less
Award ID(s):
1835443
PAR ID:
10408377
Author(s) / Creator(s):
; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Mathematical Programming
Volume:
199
Issue:
1-2
ISSN:
0025-5610
Page Range / eLocation ID:
p. 1417-1429
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    A well-known open problem of Meir and Moser asks if the squares of sidelength 1/nfor$$n\ge 2$$n2can be packed perfectly into a rectangle of area$$\sum _{n=2}^\infty n^{-2}=\pi ^2/6-1$$n=2n-2=π2/6-1. In this paper we show that for any$$1/21/2<t<1, and any$$n_0$$n0that is sufficiently large depending on t, the squares of sidelength$$n^{-t}$$n-tfor$$n\ge n_0$$nn0can be packed perfectly into a square of area$$\sum _{n=n_0}^\infty n^{-2t}$$n=n0n-2t. This was previously known (if one packs a rectangle instead of a square) for$$1/21/2<t2/3(in which case one can take$$n_0=1$$n0=1).

     
    more » « less
  2. Abstract

    A graphGisH-freeif it has no induced subgraph isomorphic toH. We prove that a$$P_5$$P5-free graph with clique number$$\omega \ge 3$$ω3has chromatic number at most$$\omega ^{\log _2(\omega )}$$ωlog2(ω). The best previous result was an exponential upper bound$$(5/27)3^{\omega }$$(5/27)3ω, due to Esperet, Lemoine, Maffray, and Morel. A polynomial bound would imply that the celebrated Erdős-Hajnal conjecture holds for$$P_5$$P5, which is the smallest open case. Thus, there is great interest in whether there is a polynomial bound for$$P_5$$P5-free graphs, and our result is an attempt to approach that.

     
    more » « less
  3. Abstract

    We construct an example of a group$$G = \mathbb {Z}^2 \times G_0$$G=Z2×G0for a finite abelian group $$G_0$$G0, a subsetEof $$G_0$$G0, and two finite subsets$$F_1,F_2$$F1,F2of G, such that it is undecidable in ZFC whether$$\mathbb {Z}^2\times E$$Z2×Ecan be tiled by translations of$$F_1,F_2$$F1,F2. In particular, this implies that this tiling problem isaperiodic, in the sense that (in the standard universe of ZFC) there exist translational tilings ofEby the tiles$$F_1,F_2$$F1,F2, but no periodic tilings. Previously, such aperiodic or undecidable translational tilings were only constructed for sets of eleven or more tiles (mostly in $$\mathbb {Z}^2$$Z2). A similar construction also applies for$$G=\mathbb {Z}^d$$G=Zdfor sufficiently large d. If one allows the group$$G_0$$G0to be non-abelian, a variant of the construction produces an undecidable translational tiling with only one tile F. The argument proceeds by first observing that a single tiling equation is able to encode an arbitrary system of tiling equations, which in turn can encode an arbitrary system of certain functional equations once one has two or more tiles. In particular, one can use two tiles to encode tiling problems for an arbitrary number of tiles.

     
    more » « less
  4. Abstract

    LetXbe ann-element point set in thek-dimensional unit cube$$[0,1]^k$$[0,1]kwhere$$k \ge 2$$k2. According to an old result of Bollobás and Meir (Oper Res Lett 11:19–21, 1992) , there exists a cycle (tour)$$x_1, x_2, \ldots , x_n$$x1,x2,,xnthrough thenpoints, such that$$\left( \sum _{i=1}^n |x_i - x_{i+1}|^k \right) ^{1/k} \le c_k$$i=1n|xi-xi+1|k1/kck, where$$|x-y|$$|x-y|is the Euclidean distance betweenxandy, and$$c_k$$ckis an absolute constant that depends only onk, where$$x_{n+1} \equiv x_1$$xn+1x1. From the other direction, for every$$k \ge 2$$k2and$$n \ge 2$$n2, there existnpoints in$$[0,1]^k$$[0,1]k, such that their shortest tour satisfies$$\left( \sum _{i=1}^n |x_i - x_{i+1}|^k \right) ^{1/k} = 2^{1/k} \cdot \sqrt{k}$$i=1n|xi-xi+1|k1/k=21/k·k. For the plane, the best constant is$$c_2=2$$c2=2and this is the only exact value known. Bollobás and Meir showed that one can take$$c_k = 9 \left( \frac{2}{3} \right) ^{1/k} \cdot \sqrt{k}$$ck=9231/k·kfor every$$k \ge 3$$k3and conjectured that the best constant is$$c_k = 2^{1/k} \cdot \sqrt{k}$$ck=21/k·k, for every$$k \ge 2$$k2. Here we significantly improve the upper bound and show that one can take$$c_k = 3 \sqrt{5} \left( \frac{2}{3} \right) ^{1/k} \cdot \sqrt{k}$$ck=35231/k·kor$$c_k = 2.91 \sqrt{k} \ (1+o_k(1))$$ck=2.91k(1+ok(1)). Our bounds are constructive. We also show that$$c_3 \ge 2^{7/6}$$c327/6, which disproves the conjecture for$$k=3$$k=3. Connections to matching problems, power assignment problems, related problems, including algorithms, are discussed in this context. A slightly revised version of the Bollobás–Meir conjecture is proposed.

     
    more » « less
  5. Abstract

    This paper studies several solution paths of sparse quadratic minimization problems as a function of the weighing parameter of the bi-objective of estimation loss versus solution sparsity. Three such paths are considered: the “$$\ell _0$$0-path” where the discontinuous$$\ell _0$$0-function provides the exact sparsity count; the “$$\ell _1$$1-path” where the$$\ell _1$$1-function provides a convex surrogate of sparsity count; and the “capped$$\ell _1$$1-path” where the nonconvex nondifferentiable capped$$\ell _1$$1-function aims to enhance the$$\ell _1$$1-approximation. Serving different purposes, each of these three formulations is different from each other, both analytically and computationally. Our results deepen the understanding of (old and new) properties of the associated paths, highlight the pros, cons, and tradeoffs of these sparse optimization models, and provide numerical evidence to support the practical superiority of the capped$$\ell _1$$1-path. Our study of the capped$$\ell _1$$1-path is interesting in its own right as the path pertains to computable directionally stationary (= strongly locally minimizing in this context, as opposed to globally optimal) solutions of a parametric nonconvex nondifferentiable optimization problem. Motivated by classical parametric quadratic programming theory and reinforced by modern statistical learning studies, both casting an exponential perspective in fully describing such solution paths, we also aim to address the question of whether some of them can be fully traced in strongly polynomial time in the problem dimensions. A major conclusion of this paper is that a path of directional stationary solutions of the capped$$\ell _1$$1-regularized problem offers interesting theoretical properties and practical compromise between the$$\ell _0$$0-path and the$$\ell _1$$1-path. Indeed, while the$$\ell _0$$0-path is computationally prohibitive and greatly handicapped by the repeated solution of mixed-integer nonlinear programs, the quality of$$\ell _1$$1-path, in terms of the two criteria—loss and sparsity—in the estimation objective, is inferior to the capped$$\ell _1$$1-path; the latter can be obtained efficiently by a combination of a parametric pivoting-like scheme supplemented by an algorithm that takes advantage of the Z-matrix structure of the loss function.

     
    more » « less