150 research outputs found

    Avoiding unnecessary demerging and remerging of multi‐commodity integer flows

    No full text
    Resource flows may merge and demerge at a network node. Sometimes several demerged flows may be immediately merged again, but in different combinations compared to before they were demerged. However, the demerging is unnecessary in the first place if the total resources at each of the network nodes involved remains unchanged. We describe this situation as “unnecessary demerging and remerging (UDR)” of flows, which would incur unnecessary operations and costs in practice. Multi‐commodity integer flows in particular will be considered in this paper. This deficiency could be theoretically overcome by means of fixed‐charge variables, but the practicality of this approach is restricted by the difficulty in solving the corresponding integer linear program (ILP). Moreover, in a problem where the objective function has many cost elements, it would be helpful if such operational costs are optimized implicitly. This paper presents a heuristic branching method within an ILP solver for removing UDR without the use of fixed‐charge variables. We use the concept of “flow potentials” (different from “flow residues” for max‐flows) guided by which underutilized arcs are heuristically banned, thus reducing occurrences of UDR. Flow connection bigraphs and flow connection groups (FCGs) are introduced. We prove that if certain conditions are met, fully utilizing an arc will guarantee an improvement within an FCG. Moreover, a location sub‐model is given when the former cannot guarantee an improvement. More importantly, the heuristic approach can significantly enhance the full fixed‐charge model by warm‐starting. Computational experiments based on real‐world instances have shown the usefulness of the proposed methods

    Decomposition Methods and Network Design Problems

    Get PDF
    Decomposition based approaches are recalled from primal and dual point of view. The possibility of building partially disaggregated reduced master problems is investigated. This extends the idea of aggregated-versus-disaggregated formulation to a gradual choice of alternative level of aggregation. Partial aggregation is applied to the linear multicommodity minimum cost flow problem. The possibility of having only partially aggregated bundles opens a wide range of alternatives with different trade-offs between the number of iterations and the required computation for solving it. This trade-off is explored for several sets of instances and the results are compared with the ones obtained by directly solving the natural node-arc formulation. An iterative solution process to the route assignment problem is proposed, based on the well-known Frank Wolfe algorithm. In order to provide a first feasible solution to the Frank Wolfe algorithm, a linear multicommodity min-cost flow problem is solved to optimality by using the decomposition techniques mentioned above. Solutions of this problem are useful for network orientation and design, especially in relation with public transportation systems as the Personal Rapid Transit. A single-commodity robust network design problem is addressed. In this, an undirected graph with edge costs is given together with a discrete set of balance matrices, representing different supply/demand scenarios. The goal is to determine the minimum cost installation of capacities on the edges such that the flow exchange is feasible for every scenario. A set of new instances that are computationally hard for the natural flow formulation are solved by means of a new heuristic algorithm. Finally, an efficient decomposition-based heuristic approach for a large scale stochastic unit commitment problem is presented. The addressed real-world stochastic problem employs at its core a deterministic unit commitment planning model developed by the California Independent System Operator (ISO)

    Learning-Based Matheuristic Solution Methods for Stochastic Network Design

    Full text link
    Cette dissertation consiste en trois Ă©tudes, chacune constituant un article de recherche. Dans tous les trois articles, nous considĂ©rons le problĂšme de conception de rĂ©seaux multiproduits, avec coĂ»t fixe, capacitĂ© et des demandes stochastiques en tant que programmes stochastiques en deux Ă©tapes. Dans un tel contexte, les dĂ©cisions de conception sont prises dans la premiĂšre Ă©tape avant que la demande rĂ©elle ne soit rĂ©alisĂ©e, tandis que les dĂ©cisions de flux de la deuxiĂšme Ă©tape ajustent la solution de la premiĂšre Ă©tape Ă  la rĂ©alisation de la demande observĂ©e. Nous considĂ©rons l’incertitude de la demande comme un nombre fini de scĂ©narios discrets, ce qui est une approche courante dans la littĂ©rature. En utilisant l’ensemble de scĂ©narios, le problĂšme mixte en nombre entier (MIP) rĂ©sultant, appelĂ© formulation Ă©tendue (FE), est extrĂȘmement difficile Ă  rĂ©soudre, sauf dans des cas triviaux. Cette thĂšse vise Ă  faire progresser le corpus de connaissances en dĂ©veloppant des algorithmes efficaces intĂ©grant des mĂ©canismes d’apprentissage en matheuristique, capables de traiter efficacement des problĂšmes stochastiques de conception pour des rĂ©seaux de grande taille. Le premier article, s’intitulĂ© "A Learning-Based Matheuristc for Stochastic Multicommodity Network Design". Nous introduisons et dĂ©crivons formellement un nouveau mĂ©canisme d’apprentissage basĂ© sur l’optimisation pour extraire des informations concernant la structure de la solution du problĂšme stochastique Ă  partir de solutions obtenues avec des combinaisons particuliĂšres de scĂ©narios. Nous proposons ensuite une matheuristique "Learn&Optimize", qui utilise les mĂ©thodes d’apprentissage pour dĂ©duire un ensemble de variables de conception prometteuses, en conjonction avec un solveur MIP de pointe pour rĂ©soudre un problĂšme rĂ©duit. Le deuxiĂšme article, s’intitulĂ© "A Reduced-Cost-Based Restriction and Refinement Matheuristic for Stochastic Network Design". Nous Ă©tudions comment concevoir efficacement des mĂ©canismes d’apprentissage basĂ©s sur l’information duale afin de guider la dĂ©termination des variables dans le contexte de la conception de rĂ©seaux stochastiques. Ce travail examine les coĂ»ts rĂ©duits associĂ©s aux variables hors base dans les solutions dĂ©terministes pour guider la sĂ©lection des variables dans la formulation stochastique. Nous proposons plusieurs stratĂ©gies pour extraire des informations sur les coĂ»ts rĂ©duits afin de fixer un ensemble appropriĂ© de variables dans le modĂšle restreint. Nous proposons ensuite une approche matheuristique utilisant des techniques itĂ©ratives de rĂ©duction des problĂšmes. Le troisiĂšme article, s’intitulĂ© "An Integrated Learning and Progressive Hedging Method to Solve Stochastic Network Design". Ici, notre objectif principal est de concevoir une mĂ©thode de rĂ©solution capable de gĂ©rer un grand nombre de scĂ©narios. Nous nous appuyons sur l’algorithme Progressive Hedging (PHA), ou les scĂ©narios sont regroupĂ©s en sous-problĂšmes. Nous intĂ©grons des methodes d’apprentissage au sein de PHA pour traiter une grand nombre de scĂ©narios. Dans notre approche, les mĂ©canismes d’apprentissage developpĂ©s dans le premier article de cette thĂšse sont adaptĂ©s pour rĂ©soudre les sous-problĂšmes multi-scĂ©narios. Nous introduisons une nouvelle solution de rĂ©fĂ©rence Ă  chaque Ă©tape d’agrĂ©gation de notre ILPH en exploitant les informations collectĂ©es Ă  partir des sous problĂšmes et nous utilisons ces informations pour mettre Ă  jour les pĂ©nalitĂ©s dans PHA. Par consĂ©quent, PHA est guidĂ© par les informations locales fournies par la procĂ©dure d’apprentissage, rĂ©sultant en une approche intĂ©grĂ©e capable de traiter des instances complexes et de grande taille. Dans les trois articles, nous montrons, au moyen de campagnes expĂ©rimentales approfondies, l’intĂ©rĂȘt des approches proposĂ©es en termes de temps de calcul et de qualitĂ© des solutions produites, en particulier pour traiter des cas trĂšs difficiles avec un grand nombre de scĂ©narios.This dissertation consists of three studies, each of which constitutes a self-contained research article. In all of the three articles, we consider the multi-commodity capacitated fixed-charge network design problem with uncertain demands as a two-stage stochastic program. In such setting, design decisions are made in the first stage before the actual demand is realized, while second-stage flow-routing decisions adjust the first-stage solution to the observed demand realization. We consider the demand uncertainty as a finite number of discrete scenarios, which is a common approach in the literature. By using the scenario set, the resulting large-scale mixed integer program (MIP) problem, referred to as the extensive form (EF), is extremely hard to solve exactly in all but trivial cases. This dissertation is aimed at advancing the body of knowledge by developing efficient algorithms incorporating learning mechanisms in matheuristics, which are able to handle large scale instances of stochastic network design problems efficiently. In the first article, we propose a novel Learning-Based Matheuristic for Stochastic Network Design Problems. We introduce and formally describe a new optimizationbased learning mechanism to extract information regarding the solution structure of a stochastic problem out of the solutions of particular combinations of scenarios. We subsequently propose the Learn&Optimize matheuristic, which makes use of the learning methods in inferring a set of promising design variables, in conjunction with a state-ofthe- art MIP solver to address a reduced problem. In the second article, we introduce a Reduced-Cost-Based Restriction and Refinement Matheuristic. We study on how to efficiently design learning mechanisms based on dual information as a means of guiding variable fixing in the context of stochastic network design. The present work investigates how the reduced cost associated with non-basic variables in deterministic solutions can be leveraged to guide variable selection within stochastic formulations. We specifically propose several strategies to extract reduced cost information so as to effectively identify an appropriate set of fixed variables within a restricted model. We then propose a matheuristic approach using problem reduction techniques iteratively (i.e., defining and exploring restricted region of global solutions, as guided by applicable dual information). Finally, in the third article, our main goal is to design a solution method that is able to manage a large number of scenarios. We rely on the progressive hedging algorithm (PHA) where the scenarios are grouped in subproblems. We propose a two phase integrated learning and progressive hedging (ILPH) approach to deal with a large number of scenarios. Within our proposed approach, the learning mechanisms from the first study of this dissertation have been adapted as an efficient heuristic method to address the multi-scenario subproblems within each iteration of PHA.We introduce a new reference point within each aggregation step of our proposed ILPH by exploiting the information garnered from subproblems, and using this information to update the penalties. Consequently, the ILPH is governed and guided by the local information provided by the learning procedure, resulting in an integrated approach capable of handling very large and complex instances. In all of the three mentioned articles, we show, by means of extensive experimental campaigns, the interest of the proposed approaches in terms of computation time and solution quality, especially in dealing with very difficult instances with a large number of scenarios

    Scheduled service network design for integrated planning of rail freight transportation

    Get PDF
    Cette thĂšse Ă©tudie une approche intĂ©grant la gestion de l’horaire et la conception de rĂ©seaux de services pour le transport ferroviaire de marchandises. Le transport par rail s’articule autour d’une structure Ă  deux niveaux de consolidation oĂč l’affectation des wagons aux blocs ainsi que des blocs aux services reprĂ©sentent des dĂ©cisions qui complexifient grandement la gestion des opĂ©rations. Dans cette thĂšse, les deux processus de consolidation ainsi que l’horaire d’exploitation sont Ă©tudiĂ©s simultanĂ©ment. La rĂ©solution de ce problĂšme permet d’identifier un plan d’exploitation rentable comprenant les politiques de blocage, le routage et l’horaire des trains, de mĂȘme que l’habillage ainsi que l’affectation du traffic. Afin de dĂ©crire les diffĂ©rentes activitĂ©s ferroviaires au niveau tactique, nous Ă©tendons le rĂ©seau physique et construisons une structure de rĂ©seau espace-temps comprenant trois couches dans lequel la dimension liĂ©e au temps prend en considĂ©ration les impacts temporels sur les opĂ©rations. De plus, les opĂ©rations relatives aux trains, blocs et wagons sont dĂ©crites par diffĂ©rentes couches. Sur la base de cette structure de rĂ©seau, nous modĂ©lisons ce problĂšme de planification ferroviaire comme un problĂšme de conception de rĂ©seaux de services. Le modĂšle proposĂ© se formule comme un programme mathĂ©matique en variables mixtes. Ce dernie r s’avĂšre trĂšs difficile Ă  rĂ©soudre en raison de la grande taille des instances traitĂ©es et de sa complexitĂ© intrinsĂšque. Trois versions sont Ă©tudiĂ©es : le modĂšle simplifiĂ© (comprenant des services directs uniquement), le modĂšle complet (comprenant des services directs et multi-arrĂȘts), ainsi qu’un modĂšle complet Ă  trĂšs grande Ă©chelle. Plusieurs heuristiques sont dĂ©veloppĂ©es afin d’obtenir de bonnes solutions en des temps de calcul raisonnables. PremiĂšrement, un cas particulier avec services directs est analysĂ©. En considĂ©rant une cara ctĂ©ristique spĂ©cifique du problĂšme de conception de rĂ©seaux de services directs nous dĂ©veloppons un nouvel algorithme de recherche avec tabous. Un voisinage par cycles est privilĂ©giĂ© Ă  cet effet. Celui-ci est basĂ© sur la distribution du flot circulant sur les blocs selon les cycles issus du rĂ©seau rĂ©siduel. Un algorithme basĂ© sur l’ajustement de pente est dĂ©veloppĂ© pour le modĂšle complet, et nous proposons une nouvelle mĂ©thode, appelĂ©e recherche ellipsoidale, permettant d’amĂ©liorer davantage la qualitĂ© de la solution. La recherche ellipsoidale combine les bonnes solutions admissibles gĂ©nĂ©rĂ©es par l’algorithme d’ajustement de pente, et regroupe les caractĂ©ristiques des bonnes solutions afin de crĂ©er un problĂšme Ă©lite qui est rĂ©solu de facon exacte Ă  l’aide d’un logiciel commercial. L’heuristique tire donc avantage de la vitesse de convergence de l’algorithme d’ajustement de pente et de la qualitĂ© de solution de la recherche ellipsoidale. Les tests numĂ©riques illustrent l’efficacitĂ© de l’heuristique proposĂ©e. En outre, l’algorithme reprĂ©sente une alternative intĂ©ressante afin de rĂ©soudre le problĂšme simplifiĂ©. Enfin, nous Ă©tudions le modĂšle complet Ă  trĂšs grande Ă©chelle. Une heuristique hybride est dĂ©veloppĂ©e en intĂ©grant les idĂ©es de l’algorithme prĂ©cĂ©demment dĂ©crit et la gĂ©nĂ©ration de colonnes. Nous proposons une nouvelle procĂ©dure d’ajustement de pente oĂč, par rapport Ă  l’ancienne, seule l’approximation des couts liĂ©s aux services est considĂ©rĂ©e. La nouvelle approche d’ajustement de pente sĂ©pare ainsi les dĂ©cisions associĂ©es aux blocs et aux services afin de fournir une dĂ©composition naturelle du problĂšme. Les rĂ©sultats numĂ©riques obtenus montrent que l’algorithme est en mesure d’identifier des solutions de qualitĂ© dans un contexte visant la rĂ©solution d’instances rĂ©elles.This thesis studies a scheduled service network design problem for rail freight transportation planning. Rails follow a special two level consolidation organization, and the car-to-block, block-to-service handling procedure complicates daily operations. In this research, the two consolidation processes as well as the operation schedule are considered simultaneously, and by solving this problem, we provide an overall cost-effective operating plan, including blocking policy, train routing, scheduling, make-up policy and traffic distribution. In order to describe various rail operations at the tactical level, we extend the physical network and construct a 3-layer time-space structure, in which the time dimension takes into consideration the temporal impacts on operations. Furthermore, operations on trains, blocks, and cars are described in different layers. Based on this network structure, we model the rail planning problem to a service network design formulation. The proposed model relies on a complex mixed-integer programming formulation. The problem is very hard to solve due to the computational difficulty as well as the tremendous size of the application instances. Three versions of the problem are studied, which are the simplified model (with only non-stop services), complete model (with both non-stop and multi-stop services) and very-large-scale complete model. Heuristic algorithms are developed to provide good feasible solutions in reasonable computing efforts. A special case with non-stop services is first studied. According to a specific characteristic of the direct service network design problem, we develop a tabu search algorithm. The tabu search moves in a cycle-based neighborhood, where flows on blocks are re-distributed according to the cycles in a conceptual residual network. A slope scaling based algorithm is developed for the complete model, and we propose a new method, called ellipsoidal search, to further improve the solution quality. Ellipsoidal search combines the good feasible solutions generated from the slope scaling, and collects the features of good solutions into an elite problem, and solves it with exact solvers. The algorithm thus takes advantage of the convergence speed of slope scaling and solution quality of ellipsoidal search, and is proven effective. The algorithm also presents an alternative for solving the simplified problem. Finally, we work on the very-large-size complete model. A hybrid heuristic is developed by integrating the ideas of previous research with column generation. We propose a new slope scaling scheme where, compared with the previous scheme, only approximate service costs instead of both service and block costs are considered. The new slope scaling scheme thus separates the block decisions and service decisions, and provide a natural decomposition of the problem. Experiments show the algorithm is good to solve real-life size instances

    Models and Methods for Merge-In-Transit Operations

    Get PDF
    We develop integer programming formulations and solution methods for addressing operational issues in merge-in-transit distribution systems. The models account for various complex problem features including the integration of inventory and transportation decisions, the dynamic and multimodal components of the application, and the non-convex piecewise linear structure of the cost functions. To accurately model the cost functions, we introduce disaggregation techniques that allow us to derive a hierarchy of linear programming relaxations. To solve these relaxations, we propose a cutting-plane procedure that combines constraint and variable generation with rounding and branch-and-bound heuristics. We demonstrate the effectiveness of this approach on a large set of test problems with instances with up to almost 500,000 integer variables derived from actual data from the computer industry. Key words : Merge-in-transit distribution systems, logistics, transportation, integer programming, disaggregation, cutting-plane method

    Learning to compare nodes in branch and bound with graph neural networks

    Full text link
    En informatique, la rĂ©solution de problĂšmes NP-difficiles en un temps raisonnable est d’une grande importance : optimisation de la chaĂźne d’approvisionnement, planification, routage, alignement de sĂ©quences biologiques multiples, inference dans les modĂšles graphiques pro- babilistes, et mĂȘme certains problĂšmes de cryptographie sont tous des examples de la classe NP-complet. En pratique, nous modĂ©lisons beaucoup d’entre eux comme un problĂšme d’op- timisation en nombre entier, que nous rĂ©solvons Ă  l’aide de la mĂ©thodologie sĂ©paration et Ă©valuation. Un algorithme de ce style divise un espace de recherche pour l’explorer rĂ©cursi- vement (sĂ©paration), et obtient des bornes d’optimalitĂ© en rĂ©solvant des relaxations linĂ©aires sur les sous-espaces (Ă©valuation). Pour spĂ©cifier un algorithme, il faut dĂ©finir plusieurs pa- ramĂštres, tel que la maniĂšre d’explorer les espaces de recherche, de diviser une recherche l’espace une fois explorĂ©, ou de renforcer les relaxations linĂ©aires. Ces politiques peuvent influencer considĂ©rablement la performance de rĂ©solution. Ce travail se concentre sur une nouvelle maniĂšre de dĂ©river politique de recherche, c’est Ă  dire le choix du prochain sous-espace Ă  sĂ©parer Ă©tant donnĂ© une partition en cours, en nous servant de l’apprentissage automatique profond. PremiĂšrement, nous collectons des donnĂ©es rĂ©sumant, sur une collection de problĂšmes donnĂ©s, quels sous-espaces contiennent l’optimum et quels ne le contiennent pas. En reprĂ©sentant ces sous-espaces sous forme de graphes bipartis qui capturent leurs caractĂ©ristiques, nous entraĂźnons un rĂ©seau de neurones graphiques Ă  dĂ©terminer la probabilitĂ© qu’un sous-espace contienne la solution optimale par apprentissage supervisĂ©. Le choix d’un tel modĂšle est particuliĂšrement utile car il peut s’adapter Ă  des problĂšmes de diffĂ©rente taille sans modifications. Nous montrons que notre approche bat celle de nos concurrents, consistant Ă  des modĂšles d’apprentissage automatique plus simples entraĂźnĂ©s Ă  partir des statistiques du solveur, ainsi que la politique par dĂ©faut de SCIP, un solveur open-source compĂ©titif, sur trois familles NP-dures: des problĂšmes de recherche de stables de taille maximum, de flots de rĂ©seau multicommoditĂ© Ă  charge fixe, et de satisfiabilitĂ© maximum.In computer science, solving NP-hard problems in a reasonable time is of great importance, such as in supply chain optimization, scheduling, routing, multiple biological sequence align- ment, inference in probabilistic graphical models, and even some problems in cryptography. In practice, we model many of them as a mixed integer linear optimization problem, which we solve using the branch and bound framework. An algorithm of this style divides a search space to explore it recursively (branch) and obtains optimality bounds by solving linear relaxations in such sub-spaces (bound). To specify an algorithm, one must set several pa- rameters, such as how to explore search spaces, how to divide a search space once it has been explored, or how to tighten these linear relaxations. These policies can significantly influence resolution performance. This work focuses on a novel method for deriving a search policy, that is, a rule for select- ing the next sub-space to explore given a current partitioning, using deep machine learning. First, we collect data summarizing which subspaces contain the optimum, and which do not. By representing these sub-spaces as bipartite graphs encoding their characteristics, we train a graph neural network to determine the probability that a subspace contains the optimal so- lution by supervised learning. The choice of such design is particularly useful as the machine learning model can automatically adapt to problems of different sizes without modifications. We show that our approach beats the one of our competitors, consisting of simpler machine learning models trained from solver statistics, as well as the default policy of SCIP, a state- of-the-art open-source solver, on three NP-hard benchmarks: generalized independent set, fixed-charge multicommodity network flow, and maximum satisfiability problems

    Iterative restricted space search : a solving approach based on hybridization

    Get PDF
    Face Ă  la complexitĂ© qui caractĂ©rise les problĂšmes d'optimisation de grande taille l'exploration complĂšte de l'espace des solutions devient rapidement un objectif inaccessible. En effet, Ă  mesure que la taille des problĂšmes augmente, des mĂ©thodes de solution de plus en plus sophistiquĂ©es sont exigĂ©es afin d'assurer un certain niveau d 'efficacitĂ©. Ceci a amenĂ© une grande partie de la communautĂ© scientifique vers le dĂ©veloppement d'outils spĂ©cifiques pour la rĂ©solution de problĂšmes de grande taille tels que les mĂ©thodes hybrides. Cependant, malgrĂ© les efforts consentis dans le dĂ©veloppement d'approches hybrides, la majoritĂ© des travaux se sont concentrĂ©s sur l'adaptation de deux ou plusieurs mĂ©thodes spĂ©cifiques, en compensant les points faibles des unes par les points forts des autres ou bien en les adaptant afin de collaborer ensemble. Au meilleur de notre connaissance, aucun travail Ă  date n'Ă  Ă©tĂ© effectuĂ© pour dĂ©velopper un cadre conceptuel pour la rĂ©solution efficace de problĂšmes d'optimisation de grande taille, qui soit Ă  la fois flexible, basĂ© sur l'Ă©change d'information et indĂ©pendant des mĂ©thodes qui le composent. L'objectif de cette thĂšse est d'explorer cette avenue de recherche en proposant un cadre conceptuel pour les mĂ©thodes hybrides, intitulĂ© la recherche itĂ©rative de l'espace restreint, ±Iterative Restricted Space Search (IRSS)>>, dont, la principale idĂ©e est la dĂ©finition et l'exploration successives de rĂ©gions restreintes de l'espace de solutions. Ces rĂ©gions, qui contiennent de bonnes solutions et qui sont assez petites pour ĂȘtre complĂštement explorĂ©es, sont appelĂ©es espaces restreints "Restricted Spaces (RS)". Ainsi, l'IRSS est une approche de solution gĂ©nĂ©rique, basĂ©e sur l'interaction de deux phases algorithmiques ayant des objectifs complĂ©mentaires. La premiĂšre phase consiste Ă  identifier une rĂ©gion restreinte intĂ©ressante et la deuxiĂšme phase consiste Ă  l'explorer. Le schĂ©ma hybride de l'approche de solution permet d'alterner entre les deux phases pour un nombre fixe d'itĂ©rations ou jusqu'Ă  l'atteinte d'une certaine limite de temps. Les concepts clĂ©s associĂ©es au dĂ©veloppement de ce cadre conceptuel et leur validation seront introduits et validĂ©s graduellement dans cette thĂšse. Ils sont prĂ©sentĂ©s de maniĂšre Ă  permettre au lecteur de comprendre les problĂšmes que nous avons rencontrĂ©s en cours de dĂ©veloppement et comment les solutions ont Ă©tĂ© conçues et implĂ©mentĂ©es. À cette fin, la thĂšse a Ă©tĂ© divisĂ©e en quatre parties. La premiĂšre est consacrĂ©e Ă  la synthĂšse de l'Ă©tat de l'art dans le domaine de recherche sur les mĂ©thodes hybrides. Elle prĂ©sente les principales approches hybrides dĂ©veloppĂ©es et leurs applications. Une brĂšve description des approches utilisant le concept de restriction d'espace est aussi prĂ©sentĂ©e dans cette partie. La deuxiĂšme partie prĂ©sente les concepts clĂ©s de ce cadre conceptuel. Il s'agit du processus d'identification des rĂ©gions restreintes et des deux phases de recherche. Ces concepts sont mis en oeuvre dans un schĂ©ma hybride heuristique et mĂ©thode exacte. L'approche a Ă©tĂ© appliquĂ©e Ă  un problĂšme d'ordonnancement avec deux niveaux de dĂ©cision, reliĂ© au contexte des pĂątes et papier: "Pulp Production Scheduling Problem". La troisiĂšme partie a permit d'approfondir les concepts dĂ©veloppĂ©s et ajuster les limitations identifiĂ©es dans la deuxiĂšme partie, en proposant une recherche itĂ©rative appliquĂ©e pour l'exploration de RS de grande taille et une structure en arbre binaire pour l'exploration de plusieurs RS. Cette structure a l'avantage d'Ă©viter l'exploration d 'un espace dĂ©jĂ  explorĂ© prĂ©cĂ©demment tout en assurant une diversification naturelle Ă  la mĂ©thode. Cette extension de la mĂ©thode a Ă©tĂ© testĂ©e sur un problĂšme de localisation et d'allocation en utilisant un schĂ©ma d'hybridation heuristique-exact de maniĂšre itĂ©rative. La quatriĂšme partie gĂ©nĂ©ralise les concepts prĂ©alablement dĂ©veloppĂ©s et conçoit un cadre gĂ©nĂ©ral qui est flexible, indĂ©pendant des mĂ©thodes utilisĂ©es et basĂ© sur un Ă©change d'informations entre les phases. Ce cadre a l'avantage d'ĂȘtre gĂ©nĂ©ral et pourrait ĂȘtre appliquĂ© Ă  une large gamme de problĂšmes

    Matheuristics: using mathematics for heuristic design

    Get PDF
    Matheuristics are heuristic algorithms based on mathematical tools such as the ones provided by mathematical programming, that are structurally general enough to be applied to different problems with little adaptations to their abstract structure. The result can be metaheuristic hybrids having components derived from the mathematical model of the problems of interest, but the mathematical techniques themselves can define general heuristic solution frameworks. In this paper, we focus our attention on mathematical programming and its contributions to developing effective heuristics. We briefly describe the mathematical tools available and then some matheuristic approaches, reporting some representative examples from the literature. We also take the opportunity to provide some ideas for possible future development

    On High-Performance Benders-Decomposition-Based Exact Methods with Application to Mixed-Integer and Stochastic Problems

    Get PDF
    RÉSUMÉ : La programmation stochastique en nombres entiers (SIP) combine la difficultĂ© de l’incertitude et de la non-convexitĂ© et constitue une catĂ©gorie de problĂšmes extrĂȘmement difficiles Ă  rĂ©soudre. La rĂ©solution efficace des problĂšmes SIP est d’une grande importance en raison de leur vaste applicabilitĂ©. Par consĂ©quent, l’intĂ©rĂȘt principal de cette dissertation porte sur les mĂ©thodes de rĂ©solution pour les SIP. Nous considĂ©rons les SIP en deux Ă©tapes et prĂ©sentons plusieurs algorithmes de dĂ©composition amĂ©liorĂ©s pour les rĂ©soudre. Notre objectif principal est de dĂ©velopper de nouveaux schĂ©mas de dĂ©composition et plusieurs techniques pour amĂ©liorer les mĂ©thodes de dĂ©composition classiques, pouvant conduire Ă  rĂ©soudre optimalement divers problĂšmes SIP. Dans le premier essai de cette thĂšse, nous prĂ©sentons une revue de littĂ©rature actualisĂ©e sur l’algorithme de dĂ©composition de Benders. Nous fournissons une taxonomie des amĂ©liorations algorithmiques et des stratĂ©gies d’accĂ©lĂ©ration de cet algorithme pour synthĂ©tiser la littĂ©rature et pour identifier les lacunes, les tendances et les directions de recherche potentielles. En outre, nous discutons de l’utilisation de la dĂ©composition de Benders pour dĂ©velopper une (mĂ©ta- )heuristique efficace, dĂ©crire les limites de l’algorithme classique et prĂ©senter des extensions permettant son application Ă  un plus large Ă©ventail de problĂšmes. Ensuite, nous dĂ©veloppons diverses techniques pour surmonter plusieurs des principaux inconvĂ©nients de l’algorithme de dĂ©composition de Benders. Nous proposons l’utilisation de plans de coupe, de dĂ©composition partielle, d’heuristiques, de coupes plus fortes, de rĂ©ductions et de stratĂ©gies de dĂ©marrage Ă  chaud pour pallier les difficultĂ©s numĂ©riques dues aux instabilitĂ©s, aux inefficacitĂ©s primales, aux faibles coupes d’optimalitĂ© ou de rĂ©alisabilitĂ©, et Ă  la faible relaxation linĂ©aire. Nous testons les stratĂ©gies proposĂ©es sur des instances de rĂ©fĂ©rence de problĂšmes de conception de rĂ©seau stochastique. Des expĂ©riences numĂ©riques illustrent l’efficacitĂ© des techniques proposĂ©es. Dans le troisiĂšme essai de cette thĂšse, nous proposons une nouvelle approche de dĂ©composition appelĂ©e mĂ©thode de dĂ©composition primale-duale. Le dĂ©veloppement de cette mĂ©thode est fondĂ© sur une reformulation spĂ©cifique des sous-problĂšmes de Benders, oĂč des copies locales des variables maĂźtresses sont introduites, puis relĂąchĂ©es dans la fonction objective. Nous montrons que la mĂ©thode proposĂ©e attĂ©nue significativement les inefficacitĂ©s primales et duales de la mĂ©thode de dĂ©composition de Benders et qu’elle est Ă©troitement liĂ©e Ă  la mĂ©thode de dĂ©composition duale lagrangienne. Les rĂ©sultats de calcul sur divers problĂšmes SIP montrent la supĂ©rioritĂ© de cette mĂ©thode par rapport aux mĂ©thodes classiques de dĂ©composition. Enfin, nous Ă©tudions la parallĂ©lisation de la mĂ©thode de dĂ©composition de Benders pour Ă©tendre ses performances numĂ©riques Ă  des instances plus larges des problĂšmes SIP. Les variantes parallĂšles disponibles de cette mĂ©thode appliquent une synchronisation rigide entre les processeurs maĂźtre et esclave. De ce fait, elles souffrent d’un important dĂ©sĂ©quilibre de charge lorsqu’elles sont appliquĂ©es aux problĂšmes SIP. Cela est dĂ» Ă  un problĂšme maĂźtre difficile qui provoque un important dĂ©sĂ©quilibre entre processeur et charge de travail. Nous proposons une mĂ©thode Benders parallĂšle asynchrone dans un cadre de type branche-et-coupe. L’assouplissement des exigences de synchronisation entraine des problĂšmes de convergence et d’efficacitĂ© divers auxquels nous rĂ©pondons en introduisant plusieurs techniques d’accĂ©lĂ©ration et de recherche. Les rĂ©sultats indiquent que notre algorithme atteint des taux d’accĂ©lĂ©ration plus Ă©levĂ©s que les mĂ©thodes synchronisĂ©es conventionnelles et qu’il est plus rapide de plusieurs ordres de grandeur que CPLEX 12.7.----------ABSTRACT : Stochastic integer programming (SIP) combines the difficulty of uncertainty and non-convexity, and constitutes a class of extremely challenging problems to solve. Efficiently solving SIP problems is of high importance due to their vast applicability. Therefore, the primary focus of this dissertation is on solution methods for SIPs. We consider two-stage SIPs and present several enhanced decomposition algorithms for solving them. Our main goal is to develop new decomposition schemes and several acceleration techniques to enhance the classical decomposition methods, which can lead to efficiently solving various SIP problems to optimality. In the first essay of this dissertation, we present a state-of-the-art survey of the Benders decomposition algorithm. We provide a taxonomy of the algorithmic enhancements and the acceleration strategies of this algorithm to synthesize the literature, and to identify shortcomings, trends and potential research directions. In addition, we discuss the use of Benders decomposition to develop efficient (meta-)heuristics, describe the limitations of the classical algorithm, and present extensions enabling its application to a broader range of problems. Next, we develop various techniques to overcome some of the main shortfalls of the Benders decomposition algorithm. We propose the use of cutting planes, partial decomposition, heuristics, stronger cuts, and warm-start strategies to alleviate the numerical challenges arising from instabilities, primal inefficiencies, weak optimality/feasibility cuts, and weak linear relaxation. We test the proposed strategies with benchmark instances from stochastic network design problems. Numerical experiments illustrate the computational efficiency of the proposed techniques. In the third essay of this dissertation, we propose a new and high-performance decomposition approach, called Benders dual decomposition method. The development of this method is based on a specific reformulation of the Benders subproblems, where local copies of the master variables are introduced and then priced out into the objective function. We show that the proposed method significantly alleviates the primal and dual shortfalls of the Benders decomposition method and it is closely related to the Lagrangian dual decomposition method. Computational results on various SIP problems show the superiority of this method compared to the classical decomposition methods as well as CPLEX 12.7. Finally, we study parallelization of the Benders decomposition method. The available parallel variants of this method implement a rigid synchronization among the master and slave processors. Thus, it suffers from significant load imbalance when applied to the SIP problems. This is mainly due to having a hard mixed-integer master problem that can take hours to be optimized. We thus propose an asynchronous parallel Benders method in a branchand- cut framework. However, relaxing the synchronization requirements entails convergence and various efficiency problems which we address them by introducing several acceleration techniques and search strategies. In particular, we propose the use of artificial subproblems, cut generation, cut aggregation, cut management, and cut propagation. The results indicate that our algorithm reaches higher speedup rates compared to the conventional synchronized methods and it is several orders of magnitude faster than CPLEX 12.7
    • 

    corecore