The quadratic assignment problem is easy for Robinsonian matrices with Toeplitz structure

We present a new polynomially solvable case of the Quadratic Assignment Problem in Koopmans-Beckman form $QAP(A,B)$, by showing that the identity permutation is optimal when $A$ and $B$ are respectively a Robinson similarity and dissimilarity matrix and one of $A$ or $B$ is a Toeplitz matrix. A Robinson (dis)similarity matrix is a symmetric matrix whose entries (increase) decrease monotonically along rows and columns when moving away from the diagonal, and such matrices arise in the classical seriation problem.Comment: 15 pages, 2 figure

The linearization problem of a binary quadratic problem and its applications

We provide several applications of the linearization problem of a binary quadratic problem. We propose a new lower bounding strategy, called the linearization-based scheme, that is based on a simple certificate for a quadratic function to be non-negative on the feasible set. Each linearization-based bound requires a set of linearizable matrices as an input. We prove that the Generalized Gilmore-Lawler bounding scheme for binary quadratic problems provides linearization-based bounds. Moreover, we show that the bound obtained from the first level reformulation linearization technique is also a type of linearization-based bound, which enables us to provide a comparison among mentioned bounds. However, the strongest linearization-based bound is the one that uses the full characterization of the set of linearizable matrices. Finally, we present a polynomial-time algorithm for the linearization problem of the quadratic shortest path problem on directed acyclic graphs. Our algorithm gives a complete characterization of the set of linearizable matrices for the quadratic shortest path problem

New special cases of the quadratic assignment problem with diagonally structured coefficient matrices

We consider new polynomially solvable cases of the well-known Quadratic Assignment Problem involving coefficient matrices with a special diagonal structure. By combining the new special cases with polynomially solvable special cases known in the literature we obtain a new and larger class of polynomially solvable special cases of the QAP where one of the two coefficient matrices involved is a Robinson matrix with an additional structural property: this matrix can be represented as a conic combination of cut matrices in a certain normal form. The other matrix is a conic combination of a monotone anti-Monge matrix and a down-benevolent Toeplitz matrix. We consider the recognition problem for the special class of Robinson matrices mentioned above and show that it can be solved in polynomial time

Linear Programming Methods for Identifying Solvable Cases of the Quadratic Assignment Problem

This research effort is concerned with identifying and characterizing families of polynomially solvable instances of the celebrated NP-hard quadratic assignment problem (qap). The approach is novel in that it uses polyhedral methods based on an equivalent mixed 0-1 linear reformulation of the problem. The continuous relaxation of this mixed 0-1 form yields a feasible region having extreme points that are both binary and fractional. The solvable instances of concern essentially possess objective function structures that ensure a binary extreme point must be optimal, so that the linear program solves the qap. The ultimate contribution of this work is the unification and subsumption of a variety of known solvable instances of the qap, and the development of a theoretical framework for identifying richer families of solvable instances. The qap was introduced over 50 years ago in the context of facility layout and location. The underlying mathematical structure, from which the problem draws its name, consists of the minimization of a quadratic function of binary variables over an assignment polytope. Since its inception, this structure has received considerable attention from various researchers, both practitioners and theoreticians alike, due to the diversity of practical applications and the resistance to exact solution procedures. Unfortunately, the combinatorial explosion of feasible solutions to the qap, in terms of the number of binary variables, creates a significant gap between the sizes of the motivating applications and the instances that can be solved by state-of-the-art solution algorithms. The most successful algorithms rely on linear forms of the qap to compute bounds within enumerative schemes. The inability to solve large qap instances has motivated researchers to seek special objective function structures that permit polynomial solvability. Various, seemingly unrelated, structures are found in the literature. This research shows that many such structures can be explained in terms of the linear reformulation which results from applying the level-1 reformulation-linearization technique (RLT) to the qap. In fact, the research shows that the level-1 RLT not only serves to explain many of these instances, but also allows for simplifications and/or generalizations. One important structure centers around instances deemed to be linearizable, where a qap instance is defined to be linearizazble if it can be equivalently rewritten as a linear assignment problem that preserves the objective function value at all feasible points. A contribution of this effort is that the constraint structure of a relaxed version of the continuous relaxation of the level-1 RLT form gives rise to a necessary and sufficient condition for an instance of the qap to be linearizable. Specifically, an instance of the qap is linearizable if and only if the given relaxed level-1 RLT form has a finite optimal solution. For all such cases, an optimal solution must occur at a binary extreme point. As a consequence, all linearizable qap instances are solvable via the level-1 RLT. The converse, however is not true, as the continuous relaxation of the level-1 RLT form can have a binary optimal solution when the qap is not linearizable. Thus, the linear program available from the level-1 RLT theoretically identifies a richer family of solvable instances. Notably, and as a consequence of this study, the level-1 RLT serves as a unifying entity in that it integrates the computation of linear programming-based bounds with the identification of polynomially solvable special cases, a relationship that was previously unnoticed

Learning to Rank and Quadratic Assignment

In this paper we show that the optimization of several ranking-based performance measures, such as precision-at-k and average-precision, is intimately related to the solution of quadratic assignment problems. Both the task of test-time prediction of the best ranking and the task of constraint generation in estimators based on structured support vector machines can all be seen as special cases of quadratic assignment problems. Although such problems are in general NP-hard, we identify a polynomially-solvable subclass (for both inference and learning) that still enables the modeling of a substantial number of pairwise rank interactions. We show preliminary results on a public benchmark image annotation data set, which indicates that this model can deliver higher performance over ranking models without pairwise rank dependencies

Combinatorial algorithms for the seriation problem

In this thesis we study the seriation problem, a combinatorial problem arising in data analysis, which asks to sequence a set of objects in such a way that similar objects are ordered close to each other. We focus on the combinatorial structure and properties of Robinsonian matrices, a special class of structured matrices which best achieve the seriation goal. Our contribution is both theoretical and practical, with a particular emphasis on algorithms. In Chapter 2 we introduce basic concepts about graphs, permutations and proximity matrices used throughout the thesis. In Chapter 3 we present Robinsonian matrices, discussing their characterizations and recognition algorithms existing in the literature. In Chapter 4 we discuss Lexicographic Breadth-First search (Lex-BFS), a special graph traversal algorithm used in multisweep algorithms for the recognition of several classes of graphs. In Chapter 5 we introduce a new Lex-BFS based algorithm to recognize Robinsonian matrices, which is derived from a new characterization of Robinsonian matrices in terms of straight enumerations of unit interval graphs. In Chapter 6 we introduce the novel Similarity-First Search algorithm (SFS), a weighted version of Lex-BFS which we use in a multisweep algorithm for the recognition of Robinsonian matrices. In Chapter 7 we model the seriation problem as an instance of Quadratic Assignment Problem (QAP) and we show that if the data has a Robinsonian structure, then one can find an optimal solution for QAP using a Robinsonian recognition algorithm. In Chapter 8 we discuss how to solve the seriation problem when the data does not have a Robinsonian structure, by finding a Robinsonian approximation of the original data. Finally, in Chapter 9 we discuss some experiments which we have carried out in order to compare the performance of the algorithms introduced in the thesis

Graph Similarity and Approximate Isomorphism

The graph similarity problem, also known as approximate graph isomorphism or graph matching problem, has been extensively studied in the machine learning community, but has not received much attention in the algorithms community: Given two graphs G,H of the same order n with adjacency matrices A_G,A_H, a well-studied measure of similarity is the Frobenius distance dist(G,H):=min_{pi}|A_G^{pi}-A_H|_F, where pi ranges over all permutations of the vertex set of G, where A_G^pi denotes the matrix obtained from A_G by permuting rows and columns according to pi, and where |M |_F is the Frobenius norm of a matrix M. The (weighted) graph similarity problem, denoted by GSim (WSim), is the problem of computing this distance for two graphs of same order. This problem is closely related to the notoriously hard quadratic assignment problem (QAP), which is known to be NP-hard even for severely restricted cases. It is known that GSim (WSim) is NP-hard; we strengthen this hardness result by showing that the problem remains NP-hard even for the class of trees. Identifying the boundary of tractability for WSim is best done in the framework of linear algebra. We show that WSim is NP-hard as long as one of the matrices has unbounded rank or negative eigenvalues: hence, the realm of tractability is restricted to positive semi-definite matrices of bounded rank. Our main result is a polynomial time algorithm for the special case where the associated (weighted) adjacency graph for one of the matrices has a bounded number of twin equivalence classes. The key parameter underlying our algorithm is the clustering number of a graph; this parameter arises in context of the spectral graph drawing machinery