A “joint + marginal” heuristic for 0/1 programs

Abstract

We present a new algorithm for solving a polynomial program P based on the recent joint + marginal approach of the first author for, parametric optimization. The idea is to first consider the variable x1 as a parameter and solve the associated (n-1)-variable (x2,...,xn) problem P(x1) where the parameter x1 is fixed and takes values in some interval Y1 with some probability uniformly distributed on Y1. Then one considers the hierarchy of what we call joint+marginal semidefinite relaxations, whose duals provide a sequence of univariate polynomial approximations that converges to the optimal value function J(x1) of problem P(x1), as k increases. Then with k fixed a priori, one computes a minimizer of the univariate polynomial pk(x1) on the interval Y1, which reduces to solving a single semidefinite program. One iterates the procedure with now an (n-2)-variable problem P(x2) with parameter x2 in some new interval Y2, etc. The quality of the approximation depends on how large k can be chosen (in general for significant size problems, k=1 is the only choice). Preliminary numerical results are provided