Complementary column generation and bounding approaches for set partitioning formulations

Nowadays, large-scale failures in elastic optical inter-data center networks (EO-DCNs) are critical issues for network survivability. Once the original connection or the primary DC fails, an alternative data center (DC) can serve the cloud services with replicated content. In this paper, we study disaster-resilient cloud services provisioning with distance-adaptive modulation formats, considering content placement, routing with modulation format adaption, protection of path and content, and spectrum allocation simultaneously. Both dedicated end-to-content backup path protection (DEBPP) and shared end-to-content backup path protection (SEBPP) are studied. First, integer linear program (ILP) models are formulated with the objective of minimum spectrum usage. For the sake of scalability, a column generation (CG) approach is developed to decompose the ILP model into a master problem and a pricing problem. The CG approach is shown to achieve comparable results to ILP models. Numerical simulations are conducted to investigate the spectrum usage of DEBPP and SEBPP with different traffic parameters and different network topologies.

A column generation (CG) approach for the solution of timetabling problems is presented. This methodology could be used for various instances of the timetabling problem, although in this paper the solution of the high-school situation in Greece is presented. The results obtained show clearly that the CG approach that has been extremely successful in recent years in the solution of airline crew scheduling problems could also be very efficient and robust for the solution of timetabling problems. Several large timetabling problems corresponding to real problems have been successfully solved, with the solutions obtained feasible and of very high quality in accordance with the problem definition. In addition, none of the solutions contained any idle hour for any of the teachers, which was one of the main goals of this optimization effort.

A column generation approach to capacitated p-median problems

The problem of scheduling and tree routing in WiMAX/802.16 based mesh networks were not defined in the standard and are thus subject to extensive research. In this paper, we consider the problem of joint routing and scheduling in 802.16-based wireless mesh network, with the objective of determining a minimum length schedule that satisfies a given (uplink/downlink) end-to-end traffic demand. Minimizing the schedule length amounts to maximizing the spectrum spatial reuse by concurrently transmitting on as many links as possible, which we refer to as a transmission configuration (a group of links that can simultaneously transmit without violating the signal-to-interferenceplus- noise ratio (SINR) requirement). Our model is referred to as maximum spatial reuse (MSR). Since there is an overwhelming number of possible transmission configurations to be assigned to time slots, we adopt the column generation technique to construct our MSR model. We present two formulations for modeling MSR, namely the link-based column generation (CGLink) formulation and the path-based column generation (CGPath) formulation. These two formulations differ mainly in the number of routing decision variables. Our experimental results indicate that the path-based formulation needs much less computational (CPU) time than the link-based formulation in order to determine the (same) optimized solution with the same spatial reuse gain.

We consider two approaches for solving the classical minimum vertex coloring problem—that is, the problem of coloring the vertices of a graph so that adjacent vertices have different colors and minimizing the number of used colors—namely, constraint programming and column generation. Constraint programming is able to solve very efficiently many of the benchmarks but suffers from a lack of effective bounding methods. On the contrary, column generation provides tight lower bounds by solving the fractional vertex coloring problem exploited in a branch-and-price algorithm, as already proposed in the literature. The column generation approach is here enhanced by using constraint programming to solve the pricing subproblem and to compute heuristic solutions. Moreover, new techniques are introduced to improve the performance of the column generation approach in solving both the linear relaxation and the integer problem. We report extensive computational results applied to the benchmark instances: we are able to prove optimality of 11 new instances and to improve the best-known lower bounds on 17 other instances. Moreover, we extend the solution approaches to a generalization of the problem known as the minimum vertex graph multicoloring problem, where a given number of colors has to be assigned to each vertex.

When a column generation approach is applied to decomposable mixed integer programming problems, it is standard to formulate and solve the master problem as a linear program. Seen in the dual space, this results in the algorithm known in the nonlinear programming community as the cutting-plane algorithm of Kelley and Cheney-Goldstein. However, more stable methods with better theoretical convergence rates are known and have been used as alternatives to this standard. One of them is the bundle method; our aim is to illustrate its differences with Kelley’s method. In the process we review alternative stabilization techniques used in column generation, comparing them from both primal and dual points of view. Numerical comparisons are presented for five applications: cutting stock (which includes bin packing), vertex coloring, capacitated vehicle routing, multi-item lot sizing, and traveling salesman. We also give a sketchy comparison with the volume algorithm.

In intensity-modulated radiation therapy (IMRT), a network flow is adopted to solve the pricing problem of the generic column generation approach in order to obtain a deliverable aperture. However, excessive computation results from the direct use of a network flow. In addition, a decline in plan quality may result from the direct determination of the leaf position using the gradient information. To overcome these problems, a column generation approach based on region growth is proposed. The proposed method is designed to reduce the computational cost of solving the pricing problem and improve the IMRT plan quality. First, the gradients of the beamlets are obtained by an objective function constructed under the constraint conditions of the organs. Second, the gradients are transformed nonlinearly. Third, the positions of the continuous negative gradient regions in each row of the aperture are determined and stored. Fourth, these gradients are taken as a whole and added to the aperture network flow, which is solved as a shortest-path problem. Finally, the deliverable aperture is obtained and added to the treatment plan. To verify the effectiveness of the proposed method, experiments involving five five-field prostate cancer cases and five nine-field head and neck cancer cases were conducted. Compared with the generic column generation method, the dose distribution of the target is ensured by the proposed method, which also effectively protects organs at risk and reduces the running time. Specifically, in ten groups of comparative experiments, the normal tissue complication probability of the proposed method is reduced by up to 3.37%, and the maximum acceleration rate is 20.44%. According to the experimental results, the proposed method is more consistent with clinical requirements compared with the generic column generation method.

The noisy intermediate-scale quantum (NISQ) devices that have just emerged in recent years provide an opportunity for the physical realization of quantum circuits. As the first step of mapping quantum circuits to these devices, qubit allocation imposes a great impact on the number of additional quantum operations required throughout the mapping. The global optimization of qubit allocation is NP-hard, but since many current NISQ devices have very few qubits, it is possible to solve this problem exactly. In this paper, we propose an exact approach for qubit allocation, which takes advantage of the branch and bound algorithm as the basic framework. In order to further prune the considerably huge state space, three techniques have been proposed and integrated into the exact approach, including an optimized lower bound for quadratic assignment problem, prioritization of logical qubits, and branch pruning by structural symmetry of physical qubits. Moreover, based on the multiple optimal solutions obtained by the exact approach, we give an error-aware qubit allocation method. For problems that are too large to be solved by the exact approach, we also give a heuristic qubit allocation approach with polynomial time complexity. Experimental evaluations show that the exact approach proposed in this paper is feasible for the qubit allocation of current NISQ architectures. For all benchmarks considered, this exact approach can find multiple optimal solutions to the qubit allocation problem on IBM Melbourne, a 16-qubit NISQ architecture, in no more than 20 min. It can serve as a evaluation baseline of any heuristic approach. Experimental evaluations also show that the proposed heuristic approach is near-optimal in terms of routing cost. Both the exact approach and the heuristic approach can be integrated into existing quantum circuit mapping approaches.

One-dimensional bin packing with pattern-dependent processing time

A column generation approach for the unconstrained binary quadratic programming problem

To develop a novel aperture-based algorithm for volumetric modulated are therapy (VMAT) treatment plan optimization with high quality and high efficiency. The VMAT optimization problem is formulated as a large-scale convex programming problem solved by a column generation approach. The authors consider a cost function consisting two terms, the first enforcing a desired dose distribution and the second guaranteeing a smooth dose rate variation between successive gantry angles. A gantry rotation is discretized into 180 beam angles and for each beam angle, only one MLC aperture is allowed. The apertures are generated one by one in a sequential way. At each iteration of the column generation method, a deliverable MLC aperture is generated for one of the unoccupied beam angles by solving a subproblem with the consideration of MLC mechanic constraints. A subsequent master problem is then solved to determine the dose rate at all currently generated apertures by minimizing the cost function. When all 180 beam angles are occupied, the optimization completes, yielding a set of deliverable apertures and associated dose rates that produce a high quality plan. The algorithm was preliminarily tested on five prostate and five head-and-neck clinical cases, each with one full gantry rotation without any couch/collimator rotations. High quality VMAT plans have been generated for all ten cases with extremely high efficiency. It takes only 5-8 min on CPU (MATLAB code on an Intel Xeon 2.27 GHz CPU) and 18-31 s on GPU (CUDA code on an NVIDIA Tesla C1060 GPU card) to generate such plans. The authors have developed an aperture-based VMAT optimization algorithm which can generate clinically deliverable high quality treatment plans at very high efficiency.

A hybrid exact approach for maximizing lifetime in sensor networks with complete and partial coverage constraints

We propose a column generation algorithm for solving sequence dependent flowshop scheduling problems (SDST flowshop) to minimize the total weighted tardiness. The continuous relaxation of the set partitioning formulation of the Dantzig-Wolfe decomposition for SDST flowshop is solved by the column generation. The pricing subproblem for the column generation is decomposed into each job-level subproblem. In order to strengthen the lower bound, the additional constraints are imposed to the pricing problems in the proposed algorithm. The cut generation algorithm is developed for the column generation. Computational experiments show that the proposed method can derive the solutions with a smaller duality gap compared with those of the conventional algorithm of ordinary column generation and Lagrangian relaxation with cuts.

Electric bus fleet composition and scheduling

The problem of scheduling and routing tree construction in WiMAX/802.16 based mesh networks is not defined in the standard and has thus been the subject to extensive research. We consider the problem of joint routing and scheduling in WiMAX-based mesh networks, with the objective of determining a minimum schedule period that satisfies a given (uplink/downlink) traffic demand. Minimizing the length of a schedule amounts to maximizing the spectrum spatial reuse by activating concurrently as many links. This group of transmission links active concurrently is referred to as the transmission group and refers to the set of wireless links that can simultaneously transmit without violating the signal-to-interference-plus-noise ratio (SINR) requirement. Our model is referred to as maximum spatial reuse (MSR). We assume centralized scheduling at the base station and attempt to maximize the system throughput through appropriate routing tree selection and achieving efficient spectrum reuse through opportunistic link scheduling. We present an ILP optimization model for the joint problem, which relies on the enumeration of all possible link schedules. Given its complexity, we decompose the problem using a column generation (CG) approach. We present two formulations for modeling MSR, namely the link-based (CGLink) and the path-based (CGPath) formulation. These two formulations differ mainly in the number of routing decision variables. Our experimental results indicate that the path-based formulation needs much less computational (CPU) time than the link-based in order to determine the (same) optimal solution with the same spatial reuse gain.