The Weighted Sparsity Problem: Complexity and Algorithms

Many optimization algorithms involve repeated processing of a fixed set of linear constraints. If the constraint matrix A is preprocessed to make it sparser, algebraic operations should become faster. In many applications there is a priori information about the likelihood that each column will appear in a basis, which can be expressed as weights on the columns. This leads to considering the weighted sparsity problem (WSP): Find a row-equivalent constraint matrix with as small a weight of nonzeros as possible. The WSP is shown to be NP-hard even with a nondegeneracy assumption, and even if restricted to instances with at most three nonzeros per both row and column. WSP is shown to have a polynomial algorithm when the number of nonzeros per either row or column is limited to at most two. This contrasts with previous results that, assuming only nondegeneracy, the unweighted version of WSP does have a polynomial algorithm (this has proven to be practically useful in tests on real data). The polynomial algorit...

[1]  A. Charnes Optimality and Degeneracy in Linear Programming , 1952 .

[2]  James G. Oxley,et al.  Matroid theory , 1992 .

[3]  R. Freund Projective Transformations for Interior Point Methods, Part II: Analysis of An Algorithm for Finding the Weighted Center of a Polyhedral System , 1988 .

[4]  Mauricio G. C. Resende,et al.  Data Structures and Programming Techniques for the Implementation of Karmarkar's Algorithm , 1989, INFORMS J. Comput..

[5]  S. Chang Increasing sparsity in matrices for large-scale optimization--theoretical properties and implementational aspects , 1989 .

[6]  A. Hoffman,et al.  A Fast Algorithm that makes Matrices Optimally Sparse , 1984 .

[7]  M. Iri,et al.  Structural solvability of systems of equations —A mathematical formulation for distinguishing accurate and inaccurate numbers in structural analysis of systems— , 1985 .

[8]  András Frank,et al.  Generalized polymatroids and submodular flows , 1988, Math. Program..

[9]  S. Making sparse matrices sparser : Computational results , 1988 .

[10]  David S. Johnson,et al.  The Rectilinear Steiner Tree Problem is NP Complete , 1977, SIAM Journal of Applied Mathematics.

[11]  Eugene Lawler,et al.  Combinatorial optimization , 1976 .

[12]  Nimrod Megiddo,et al.  Towards a Genuinely Polynomial Algorithm for Linear Programming , 1983, SIAM J. Comput..

[13]  S. T. McCormick,et al.  A combinatorial approach to some sparse matrix problems , 1983 .

[14]  J. Orlin On the simplex algorithm for networks and generalized networks , 1983 .

[15]  Martin Grötschel,et al.  The ellipsoid method and its consequences in combinatorial optimization , 1981, Comb..

[16]  Irvin Lustig,et al.  An analysis of an available set of linear programming test problems , 1989, Comput. Oper. Res..

[17]  S. Thomas McCormick,et al.  A hierarchical algorithm for making sparse matrices sparser , 1992, Math. Program..

[18]  S. Thomas McCormick,et al.  Making sparse matrices sparser: Computational results , 1990, Math. Program..

[19]  Andrew V. Goldberg,et al.  Combinatorial Algorithms for the Generalized Circulation Problem , 1991, Math. Oper. Res..

[20]  Thomas F. Coleman,et al.  Large Sparse Numerical Optimization , 1984, Lecture Notes in Computer Science.