Mapping constrained optimization problems to quantum annealing with application to fault diagnosis

Current quantum annealing (QA) hardware suffers from practical limitations such as finite temperature, sparse connectivity, small qubit numbers, and control error. We propose new algorithms for mapping boolean constraint satisfaction problems (CSPs) onto QA hardware mitigating these limitations. In particular we develop a new embedding algorithm for mapping a CSP onto a hardware Ising model with a fixed sparse set of interactions, and propose two new decomposition algorithms for solving problems too large to map directly into hardware. The mapping technique is locally-structured, as hardware compatible Ising models are generated for each problem constraint, and variables appearing in different constraints are chained together using ferromagnetic couplings. In contrast, global embedding techniques generate a hardware independent Ising model for all the constraints, and then use a minor-embedding algorithm to generate a hardware compatible Ising model. We give an example of a class of CSPs for which the scaling performance of D-Wave's QA hardware using the local mapping technique is significantly better than global embedding. We validate the approach by applying D-Wave's hardware to circuit-based fault-diagnosis. For circuits that embed directly, we find that the hardware is typically able to find all solutions from a min-fault diagnosis set of size N using 1000N samples, using an annealing rate that is 25 times faster than a leading SAT-based sampling method. Further, we apply decomposition algorithms to find min-cardinality faults for circuits that are up to 5 times larger than can be solved directly on current hardware.Comment: 22 pages, 4 figure

Kiertovaihtoalgoritmi ja muunnoksia yleistetylle ajoneuvoreititysongelmalle

Vehicle routing problems have numerous applications in fields such as transportation, supply logistics and network design. The optimal design of these routes fall in the category of NP-hard optimization problems, meaning that the computational complexity increases extremely fast with increasing problem size. The Generalized Vehicle Routing Problem (GVRP) is a general problem type that includes a broad variety of other problems as special cases. The main special feature of the GVRP is that the customers are grouped in clusters. For each cluster, only one customer is visited. In this thesis, we implement a heuristic algorithm to solve GVRP instances in reasonable time. Especially, we include a cyclic exchange method that considers a very large search neighborhood. In addition, we study the related Capacitated m-Ring-Star Problem (CmRSP). We present the Distance-Constrained Capacitated m-Ring-Star Problem (DCmRSP) and show that it contains the Multivehicle Covering Tour Problem (MCTP) as a special case. We show that DCmRSP instances can be transformed to (distance-constrained) GVRP with minor adaptations and solved with the same heuristic algorithm. Our algorithm is able to find best known solutions to all GVRP test instances; for two of them, our method shows strict improvement. The transformed CmRSP and MCTP instances are solved successfully by the same algorithm with adequate performance.Ajoneuvoreititysongelmilla on lukuisia sovelluksia muun muassa logistiikan ja verkostosuunnittelun aloilla. Tällaisten reittien optimaalinen ratkaiseminen kuuluu NP-vaikeiden optimointiongelmien kategoriaan, eli ratkaisuun vaadittava laskentateho kasvaa erittäin nopeasti ongelman koon suhteen. Yleistetty ajoneuvoreititysongelma (Generalized Vehicle Routing Problem, GVRP) on ongelmatyyppi, joka kattaa joukon muita reititysongelmia erikoistapauksina. GVRP:n selkein erityispiirre on asiakkaiden jako ryppäisiin: kussakin ryppäässä on käytävä tasan yhden asiakkaan luona. Tässä diplomityössä esitellään ja toteutetaan heuristinen algoritmi, joka etsii kohtalaisessa ajassa ratkaisuja GVRP-ongelmiin. Menetelmä sisältää kiertovaihtoalgoritmin, joka kykenee etsimään ratkaisuja hyvin laajasta ympäristöstä. Tutkimuksen kohteena on lisäksi m-rengastähtiongelma (Capacitated m-Ring-Star Problem, CmRSP). Esittelemme ongelman etäisyysrajoitetun version (DCmRSP), ja näytämme, että kyseiseen ongelmaan sisältyy usean ajoneuvon peittävän reitin ongelma (Multivehicle Covering Tour Problem). Näytämme, että DCmRSP-ongelman pystyy pienin muutoksin muuntamaan GVRP-ongelmaksi ja ratkaisemaan samalla heuristisella algoritmilla. Algoritmi löytää parhaat tunnetut ratkaisut kaikkiin GVRP-testitehtäviin. Kahdessa tapauksessa ratkaisu on parempi aiemmin löydettyihin nähden. Algoritmi kykenee ratkaisemaan muunnetut CmRSP- ja MCTP-testitehtävät kohtalaisella ratkaisulaadulla

Action planning for graph transition systems

Graphs are suitable modeling formalisms for software and hardware systems involving aspects such as communication, object orientation, concurrency, mobility and distribution. State spaces of such systems can be represented by graph transition systems, which are basically transition systems whose states and transitions represent graphs and graph morphisms. In this paper, we propose the modeling of graph transition systems in PDDL and the application of heuristic search planning for their analysis. We consider different heuristics and present experimental results

Fast Routing Table Construction Using Small Messages

We describe a distributed randomized algorithm computing approximate distances and routes that approximate shortest paths. Let n denote the number of nodes in the graph, and let HD denote the hop diameter of the graph, i.e., the diameter of the graph when all edges are considered to have unit weight. Given 0 < eps <= 1/2, our algorithm runs in weak-O(n^(1/2 + eps) + HD) communication rounds using messages of O(log n) bits and guarantees a stretch of O(eps^(-1) log eps^(-1)) with high probability. This is the first distributed algorithm approximating weighted shortest paths that uses small messages and runs in weak-o(n) time (in graphs where HD in weak-o(n)). The time complexity nearly matches the lower bounds of weak-Omega(sqrt(n) + HD) in the small-messages model that hold for stateless routing (where routing decisions do not depend on the traversed path) as well as approximation of the weigthed diameter. Our scheme replaces the original identifiers of the nodes by labels of size O(log eps^(-1) log n). We show that no algorithm that keeps the original identifiers and runs for weak-o(n) rounds can achieve a polylogarithmic approximation ratio. Variations of our techniques yield a number of fast distributed approximation algorithms solving related problems using small messages. Specifically, we present algorithms that run in weak-O(n^(1/2 + eps) + HD) rounds for a given 0 < eps <= 1/2, and solve, with high probability, the following problems: - O(eps^(-1))-approximation for the Generalized Steiner Forest (the running time in this case has an additive weak-O(t^(1 + 2eps)) term, where t is the number of terminals); - O(eps^(-2))-approximation of weighted distances, using node labels of size O(eps^(-1) log n) and weak-O(n^(eps)) bits of memory per node; - O(eps^(-1))-approximation of the weighted diameter; - O(eps^(-3))-approximate shortest paths using the labels 1,...,n.Comment: 40 pages, 2 figures, extended abstract submitted to STOC'1

A Distributed Newton Method for Network Utility Maximization

Most existing work uses dual decomposition and subgradient methods to solve Network Utility Maximization (NUM) problems in a distributed manner, which suffer from slow rate of convergence properties. This work develops an alternative distributed Newton-type fast converging algorithm for solving network utility maximization problems with self-concordant utility functions. By using novel matrix splitting techniques, both primal and dual updates for the Newton step can be computed using iterative schemes in a decentralized manner with limited information exchange. Similarly, the stepsize can be obtained via an iterative consensus-based averaging scheme. We show that even when the Newton direction and the stepsize in our method are computed within some error (due to finite truncation of the iterative schemes), the resulting objective function value still converges superlinearly to an explicitly characterized error neighborhood. Simulation results demonstrate significant convergence rate improvement of our algorithm relative to the existing subgradient methods based on dual decomposition.Comment: 27 pages, 4 figures, LIDS report, submitted to CDC 201

An asymmetric TSP with time windows and with time-dependent travel times and costs: An exact solution through a graph transformation

[EN] In this paper we deal with an extended version of the Asymmetric Traveling Salesman Problem with Time Windows (ATSPTW) that considers time-dependent travel times and costs, for a more accurate approximation of some routing problems inside large cities, in which the time or cost of traversing certain streets (e.g. main avenues) depends on the moment of the day (for example rush-hours). Unlike other existing papers about time-dependent routing problems, we focus on an exact method for solving this new problem. For this end we first transform the problem into an Asymmetric Generalized TSP and then into a Graphical Asymmetric TSP. In this way, we can apply a known exact algorithm for the Mixed General Routing Problem, which seems to run well with our resulting instances. Computational results are presented on a set of 270 adapted instances from benchmark ATSPTW instances.This work has been partially supported by the Ministerio de Ciencia y Tecnología of Spain (project TIC2003-05982-C05-01) and the Generalitat Valenciana (Ref: GRUPOS03/189). Thanks are due to Michel Gendreau, Alain Hertz, Gilbert Laporte and Mihnea Stan for providing us the set of benchmark ATSPTW instances, and to Matteo Fischetti and Norbert Ascheuer for their suggestions and help about the computational experiments. Last we are also indebted to the three anonymous referees for their valuable comments.Albiach, J.; Sanchís Llopis, JM.; Soler Fernández, D. (2008). An asymmetric TSP with time windows and with time-dependent travel times and costs: An exact solution through a graph transformation. European Journal of Operational Research. 189(3):789-802. https://doi.org/10.1016/j.ejor.2006.09.099S789802189

A new genetic algorithm for the asymmetric traveling salesman problem

The asymmetric traveling salesman problem (ATSP) is one of the most important combinatorial optimization problems. It allows us to solve, either directly or through a transformation, many real-world problems. We present in this paper a new competitive genetic algorithm to solve this problem. This algorithm has been checked on a set of 153 benchmark instances with known optimal solution and it outperforms the results obtained with previous ATSP heuristic methods. © 2012 Elsevier Ltd. All rights reserved.This work has been partially supported by the Ministerio de Educacion y Ciencia of Spain (Project No. TIN2008-06441-C02-01).Yuichi Nagata; Soler Fernández, D. (2012). A new genetic algorithm for the asymmetric traveling salesman problem. Expert Systems with Applications. 39(10):8947-8953. https://doi.org/10.1016/j.eswa.2012.02.029S89478953391