2,237 research outputs found
A lagrangian relaxation-based heuristic to solve large extended graph partitioning problems
© Springer International Publishing Switzerland 2016. The paper is concerned with the planning of training sessions in large organisations requiring periodic retraining of their staff. The allocation of students must take into account student preferences as well as the desired composition of study groups. The paper presents a bicriteria Quadratic Multiple Knapsack formulation of the considered practical problem, and a novel solution procedure based on Lagrangian relaxation. The paper presents the results of computational experiments aimed at testing the optimisation procedure on real world data originating from Australia’s largest electricity distributor. Results are compared and validated against a Genetic Algorithm based matheuristic
Improving the Asymmetric TSP by Considering Graph Structure
Recent works on cost based relaxations have improved Constraint Programming
(CP) models for the Traveling Salesman Problem (TSP). We provide a short survey
over solving asymmetric TSP with CP. Then, we suggest new implied propagators
based on general graph properties. We experimentally show that such implied
propagators bring robustness to pathological instances and highlight the fact
that graph structure can significantly improve search heuristics behavior.
Finally, we show that our approach outperforms current state of the art
results.Comment: Technical repor
Using a conic bundle method to accelerate both phases of a quadratic convex reformulation
We present algorithm MIQCR-CB that is an advancement of method
MIQCR~(Billionnet, Elloumi and Lambert, 2012). MIQCR is a method for solving
mixed-integer quadratic programs and works in two phases: the first phase
determines an equivalent quadratic formulation with a convex objective function
by solving a semidefinite problem , and, in the second phase, the
equivalent formulation is solved by a standard solver. As the reformulation
relies on the solution of a large-scale semidefinite program, it is not
tractable by existing semidefinite solvers, already for medium sized problems.
To surmount this difficulty, we present in MIQCR-CB a subgradient algorithm
within a Lagrangian duality framework for solving that substantially
speeds up the first phase. Moreover, this algorithm leads to a reformulated
problem of smaller size than the one obtained by the original MIQCR method
which results in a shorter time for solving the second phase.
We present extensive computational results to show the efficiency of our
algorithm
A Lagrangian relaxation approach to the edge-weighted clique problem
The -clique polytope is the convex hull of the node and edge incidence vectors of all subcliques of size at most of a complete graph on nodes. Including the Boolean quadric polytope as a special case and being closely related to the quadratic knapsack polytope, it has received considerable attention in the literature. In particular, the max-cut problem is equivalent with optimizing a linear function over . The problem of optimizing linear functions over has so far been approached via heuristic combinatorial algorithms and cutting-plane methods. We study the structure of in further detail and present a new computational approach to the linear optimization problem based on Lucena's suggestion of integrating cutting planes into a Lagrangian relaxation of an integer programming problem. In particular, we show that the separation problem for tree inequalities becomes polynomial in our Lagrangian framework. Finally, computational results are presented. \u
An Exact Algorithm for Side-Chain Placement in Protein Design
Computational protein design aims at constructing novel or improved functions
on the structure of a given protein backbone and has important applications in
the pharmaceutical and biotechnical industry. The underlying combinatorial
side-chain placement problem consists of choosing a side-chain placement for
each residue position such that the resulting overall energy is minimum. The
choice of the side-chain then also determines the amino acid for this position.
Many algorithms for this NP-hard problem have been proposed in the context of
homology modeling, which, however, reach their limits when faced with large
protein design instances.
In this paper, we propose a new exact method for the side-chain placement
problem that works well even for large instance sizes as they appear in protein
design. Our main contribution is a dedicated branch-and-bound algorithm that
combines tight upper and lower bounds resulting from a novel Lagrangian
relaxation approach for side-chain placement. Our experimental results show
that our method outperforms alternative state-of-the art exact approaches and
makes it possible to optimally solve large protein design instances routinely
Efficient Semidefinite Branch-and-Cut for MAP-MRF Inference
We propose a Branch-and-Cut (B&C) method for solving general MAP-MRF
inference problems. The core of our method is a very efficient bounding
procedure, which combines scalable semidefinite programming (SDP) and a
cutting-plane method for seeking violated constraints. In order to further
speed up the computation, several strategies have been exploited, including
model reduction, warm start and removal of inactive constraints.
We analyze the performance of the proposed method under different settings,
and demonstrate that our method either outperforms or performs on par with
state-of-the-art approaches. Especially when the connectivities are dense or
when the relative magnitudes of the unary costs are low, we achieve the best
reported results. Experiments show that the proposed algorithm achieves better
approximation than the state-of-the-art methods within a variety of time
budgets on challenging non-submodular MAP-MRF inference problems.Comment: 21 page
Eigenvector Synchronization, Graph Rigidity and the Molecule Problem
The graph realization problem has received a great deal of attention in
recent years, due to its importance in applications such as wireless sensor
networks and structural biology. In this paper, we extend on previous work and
propose the 3D-ASAP algorithm, for the graph realization problem in
, given a sparse and noisy set of distance measurements. 3D-ASAP
is a divide and conquer, non-incremental and non-iterative algorithm, which
integrates local distance information into a global structure determination.
Our approach starts with identifying, for every node, a subgraph of its 1-hop
neighborhood graph, which can be accurately embedded in its own coordinate
system. In the noise-free case, the computed coordinates of the sensors in each
patch must agree with their global positioning up to some unknown rigid motion,
that is, up to translation, rotation and possibly reflection. In other words,
to every patch there corresponds an element of the Euclidean group Euc(3) of
rigid transformations in , and the goal is to estimate the group
elements that will properly align all the patches in a globally consistent way.
Furthermore, 3D-ASAP successfully incorporates information specific to the
molecule problem in structural biology, in particular information on known
substructures and their orientation. In addition, we also propose 3D-SP-ASAP, a
faster version of 3D-ASAP, which uses a spectral partitioning algorithm as a
preprocessing step for dividing the initial graph into smaller subgraphs. Our
extensive numerical simulations show that 3D-ASAP and 3D-SP-ASAP are very
robust to high levels of noise in the measured distances and to sparse
connectivity in the measurement graph, and compare favorably to similar
state-of-the art localization algorithms.Comment: 49 pages, 8 figure
Railway Crew Rescheduling with Retiming
Railway operations are disrupted frequently, e.g. the Dutch railway network experiences about three large disruptions per day on average. In such a disrupted situation railway operators need to quickly adjust their resource schedules. Nowadays, the timetable, the rolling stock and the crew schedule are recovered in a sequential way. In this paper, we model and solve the crew rescheduling problem with retiming. This problem extends the crew rescheduling problem by the possibility to delay the departure of some trains. In this way we partly integrate timetable adjustment and crew rescheduling. The algorithm is based on column generation techniques combined with Lagrangian heuristics. In order to prevent a large increase in computational time, retiming is allowed only for a limited number of trains where it seems very promising. Computational experiments with real-life disruption data show that, compared to the classical approach, it is possible to find better solutions by using crew rescheduling with retiming.
- …