带有限反馈下层的二层规划问题的部分合作模型

对下层最优反馈为离散有限多个的二层规划问题的部分合作模型进行探讨. 当下层的合作程度依赖于上层的决策变量时, 给出一个确定合作系数函数的一般方法, 进而得到一个新的部分合作模型. 在适当地假设下, 可保证所给的部分合作模型一定可以找到比悲观解要好的解, 并结合新的部分合作模型对原不适定问题进行分析, 得到了一些有益的结论. 最后以实际算例说明了所给部分合作模型的可行性.     

Linear bilevel multi-follower programming with independent followers

This paper considers a particular case of linear bilevel programming problems with one leader and multiple followers. In this model, the followers are independent, meaning that the objective function and the set of constraints of each follower only include the leader’s variables and his own variables. We prove that this problem can be reformulated into a linear bilevel problem with one leader and one follower by defining an adequate second level objective function and constraint region. In the second part of the paper we show that the results on the optimality of the linear bilevel problem with multiple independent followers presented in Shi et al. [The kth-best approach for linear bilevel multi-follower programming, J. Global Optim. 33, 563–578 (2005)] are based on a misconstruction of the inducible region.     

An efficient evolutionary algorithm for the ring star problem

This paper addresses the ring star problem (RSP). The goal is to locate a cycle through a subset of nodes of a network aiming to minimize the sum of the cost of installing facilities on the nodes on the cycle, the cost of connecting them and the cost of assigning the nodes not on the cycle to their closest node on the cycle. A fast and efficient evolutionary algorithm is developed which is based on a new formulation of the RSP as a bilevel programming problem with one leader and two independent followers. The leader decides which nodes to include in the ring, one follower decides about the connections of the cycle and the other follower decides about the assignment of the nodes not on the cycle. The bilevel approach leads to a new form of chromosome encoding in which genes are associated to values of the upper level variables. The quality of each chromosome is evaluated by its fitness, by means of the objective function of the RSP. Hence, in order to compute the value of the lower level variables, two optimization problems are solved for each chromosome. The computational results show the efficiency of the algorithm in terms of the quality of the solutions yielded and the computing time. A study to select the best configuration of the algorithm is presented. The algorithm is tested on a set of benchmark problems providing very accurate solutions within short computing times. Moreover, for one of the problems a new best solution is found.     

Nash bargaining and leader–follower models in water allocation: Application to the Zarrinehrud River basin,Iran

Mathematical models for conflict resolution are very important in integrated water resources and environmental management. This study proposes a new methodology to resolve conflicts among different water users and water suppliers while considering environmental requirements and the system’s constraints. A two-level leader–follower model is applied to maximize the net benefit with the Iran Water Resources Management Company as the leader and agricultural, domestic, and industrial users as followers subject to the system’s constraints. As a comparison, the Nash bargaining solution is also used to find a solution when simultaneous moves are assumed by the participants. The suggested method is then applied to the real case of the Zarrinehrud River basin that is one of the areas facing water shortages in Iran. For the actual optimization, Genetic Algorithm is used in order to avoid local optimum. As the contribution of this study, the results show that benefits for the leader in the leader–follower model increased in comparison with the Nash bargaining solutions.     

A decomposition approach for global optimum search in QP,NLP and MINLP problems

Nonconvex programming problems are frequently encountered in engineering and operations research. A large variety of global optimization algorithms have been proposed for the various classes of programming problems. A new approach for global optimum search is presented in this paper which involves a decomposition of the variable set into two sets —complicating and noncomplicating variables. This results in a decomposition of the constraint set leading to two subproblems. The decomposition of the original problem induces special structure in the resulting subproblems and a series of these subproblems are then solved, using the Generalized Benders' Decomposition technique, to determine the optimal solution. The key idea is to combine a judicious selection of the complicating variables with suitable transformations leading to subproblems which can attain their respective global solutions at each iteration. Mathematical properties of the proposed approach are presented. Even though the proposed approach cannot guarantee the determination of the global optimum, computational experience on a number of nonconvex QP, NLP and MINLP example problems indicates that a global optimum solution can be obtained from various starting points.     

On the use of intersection cuts for bilevel optimization

We address a generic mixed-integer bilevel linear program (MIBLP), i.e., a bilevel optimization problem where all objective functions and constraints are linear, and some/all variables are required to take integer values. We first propose necessary modifications needed to turn a standard branch-and-bound MILP solver into an exact and finitely-convergent MIBLP solver, also addressing MIBLP unboundedness and infeasibility. As in other approaches from the literature, our scheme is finitely-convergent in case both the leader and the follower problems are pure integer. In addition, it is capable of dealing with continuous variables both in the leader and in follower problems—provided that the leader variables influencing follower’s decisions are integer and bounded. We then introduce new classes of linear inequalities to be embedded in this branch-and-bound framework, some of which are intersection cuts based on feasible-free convex sets. We present a computational study on various classes of benchmark instances available from the literature, in which we demonstrate that our approach outperforms alternative state-of-the-art MIBLP methods.     

Branch-and-Price Algorithms for the One-Dimensional Cutting Stock Problem

We compare two branch-and-price approaches for the cutting stock problem. Each algorithm is based on a different integer programming formulation of the column generation master problem. One formulation results in a master problem with 0–1 integer variables while the other has general integer variables. Both algorithms employ column generation for solving LP relaxations at each node of a branch-and-bound tree to obtain optimal integer solutions. These different formulations yield the same column generation subproblem, but require different branch-and-bound approaches. Computational results for both real and randomly generated test problems are presented.     

Improve-and-Branch Algorithm for the Global Optimization of Nonconvex NLP Problems

A new algorithm to solve nonconvex NLP problems is presented. It is based on the solution of two problems. The reformulated problem RP is a suitable reformulation of the original problem and involves convex terms and concave univariate terms. The main problem MP is a nonconvex NLP that outer-approximates the feasible region and underestimate the objective function. MP involves convex terms and terms which are the products of concave univariate functions and new variables. Fixing the variables in the concave terms, a convex NLP that overestimates the feasible region and underestimates the objective function is obtained from the MP. Like most of the deterministic global optimization algorithms, bounds on all the variables in the nonconvex terms must be provided. MP forces the objective value to improve and minimizes the difference of upper and lower bound of all the variables either to zero or to a positive value. In the first case, a feasible solution of the original problem is reached and the objective function is improved. In general terms, the second case corresponds to an infeasible solution of the original problem due to the existence of gaps in some variables. A branching procedure is performed in order to either prove that there is no better solution or reduce the domain, eliminating the local solution of MP that was found. The MP solution indicates a key point to do the branching. A bound reduction technique is implemented to accelerate the convergence speed. Computational results demonstrate that the algorithm compares very favorably to other approaches when applied to test problems and process design problems. It is typically faster and it produces very accurate results.     

A decomposition approach to solve a bilevel capacitated facility location problem with equity constraints

In this paper, we study a capacitated facility location problem with two decision makers. One (say, the leader) decides on which subset of facilities to open and the capacity to be installed in each facility with the goal of minimizing the overall costs; the second decision maker (say, the follower), once the facilities have been designed, aims at maximizing the profit deriving from satisfying the demands of a given set of clients beyond a certain threshold imposed by the leader. The leader can foresee but cannot control the follower’s behavior. The resulting mathematical formulation is a discrete–continuous bilevel optimization problem. We propose a decomposition approach to cope with the bilevel structure of the problem and the integrality of a subset of variables under the control of the leader. Such a proposal has been tested on a set of benchmark instances available in the literature.     

Post-optimality sensitivity analysis in abstract spaces with applications to continuous-time programming problems

The programming problem under consideration consists in maximizing a concave objective functional, subject to convex operator inequality contraints. The assumptions include the existence of an optimum solution, Fréchet differentiability of all operators involved, and the existence of the topological complement of the null space of the Fréchet derivative of the constraint operator. It is shown that the rate of change of the optimum value of the objective functional due to the perturbation is measured by the dual. The optimum values of the primal variables are locally approximated as linear functions of the perturbation; the theory of generalized inverse operators is used in the approximation. We give an approximation to the primal variables if the problem is perturbed. The results are specialized for some continuous-time and finite-dimensional cases. Two examples for finite-dimensional problems are given. We apply the theory to the continuous-time linear programming problem and prove some continuity results for the optimal primal and dual objective functionals.The authors are indebted to the Natural Sciences and Engineering Research Council of Canada for financial support through Grants A4109 and A7329, respectively. They would also like to thank the referee for his comments.     

A discrete bilevel brain storm algorithm for solving a sales territory design problem: a case study

A sales territory design problem faced by a manufacturing company that supplies products to a group of customers located in a service region is addressed in this paper. The planning process of designing the territories has the objective to minimizing the total dispersion of the customers without exceeding a limited budget assigned to each territory. Once territories have been determined, a salesperson has to define the day-by-day routes to satisfy the demand of customers. Currently, the company has established a service level policy that aims to minimize total waiting times during the distribution process. Also, each territory is served by a single salesperson. A novel discrete bilevel optimization model for the sales territory design problem is proposed. This problem can be seen as a bilevel problem with a single leader and multiple independent followers, in which the leader’s problem corresponds to the design of territories (manager of the company), and the routing decision for each territory corresponds to each follower. The hierarchical nature of the current company’s decision-making process triggers some particular characteristics of the bilevel model. A brain storm algorithm that exploits these characteristics is proposed to solve the discrete bilevel problem. The main features of the proposed algorithm are that the workload is used to verify the feasibility and to cluster the leader’s solutions. In addition, four discrete mechanisms are used to generate new solutions, and an elite set of solutions is considered to reduce computational cost. This algorithm is used to solve a real case study, and the results are compared against the current solution given by the company. Results show a reduction of more than 20% in the current costs with the solution obtained by the proposed algorithm. Furthermore, a sensitivity analysis is performed, providing interesting managerial insights to improve the current operations of the company.     

The bilevel linear/linear fractional programming problem

Bilevel programming involves two optimization problems where the constraint region of the first level problem is implicitly determined by another optimization problem. In this paper we consider the bilevel linear/linear fractional programming problem in which the objective function of the first level is linear, the objective function of the second level is linear fractional and the feasible region is a polyhedron. For this problem we prove that an optimal solution can be found which is an extreme point of the polyhedron. Moreover, taking into account the relationship between feasible solutions to the problem and bases of the technological coefficient submatrix associated to variables of the second level, an enumerative algorithm is proposed that finds a global optimum to the problem.     

Controlled predatory pricing in a multiperiod Stackelberg game: an MPEC approach

We analyze a multiperiod oligopolistic market where each period is a Stackelberg game between a leader firm and multiple follower firms. The leader chooses his production level first, taking into account the reaction of the followers. Then, the follower firms decide their production levels after observing the leader’s decision. The difference between the proposed model and other models discussed in literature is that the leader firm has the power to force the follower firms out of business by preventing them from achieving a target sales level in a given time period. The leader firm has an incentive to lower the market prices possibly lower than the Stackelberg equilibrium in order to push the followers to sell less and eventually go out of business. Intentionally lowering the market prices to force competitors to fail is known as predatory pricing, and is illegal under antitrust laws since it negatively affects consumer welfare. In this work, we show that there exists a predatory pricing strategy where the market price is above the average cost and consumer welfare is preserved. We develop a mixed integer nonlinear problem (MINLP) that models the multiperiod Stackelberg game. The MINLP problem is transformed to a mixed integer linear problem (MILP) by using binary variables and piecewise linearization. A cutting plane algorithm is used to solve the resulting MILP. The results show that firms can engage in predatory pricing even if the average market price is forced to remain higher than the average cost. Furthermore, we show that in order to protect the consumers, antitrust laws can control predatory pricing by setting rules on consumer welfare.     

Evolutionary hybrid approaches for generation scheduling in power systems

Generation scheduling (GS) in power systems is a tough optimisation problem which continues to present a challenge for efficient solution techniques. The solution is to define on/off decisions and generation levels for each electricity generator of a power system for each scheduling interval. The solution procedure requires simultaneous consideration of binary decision and continuous variables. In recent years researchers have focused much attention on developing new hybrid approaches using evolutionary and traditional exact methods for this type of mixed-integer problems. This paper investigates how the optimum or near optimum solution for the GS problem may be quickly identified. A design is proposed which uses a variety of metaheuristic, heuristics and mathematical programming techniques within a hybrid framework. The results obtained for two case studies are promising and show that the hybrid approach offers an effective alternative for solving the GS problems within a realistic timeframe.     

A mathematical programming model and solution for scheduling production orders in Shanghai Baoshan Iron and Steel Complex

In this paper, we investigate the production order scheduling problem derived from the production of steel sheets in Shanghai Baoshan Iron and Steel Complex (Baosteel). A deterministic mixed integer programming (MIP) model for scheduling production orders on some critical and bottleneck operations in Baosteel is presented in which practical technological constraints have been considered. The objective is to determine the starting and ending times of production orders on corresponding operations under capacity constraints for minimizing the sum of weighted completion times of all orders. Due to large numbers of variables and constraints in the model, a decomposition solution methodology based on a synergistic combination of Lagrangian relaxation, linear programming and heuristics is developed. Unlike the commonly used method of relaxing capacity constraints, this methodology alternatively relaxes constraints coupling integer variables with continuous variables which are introduced to the objective function by Lagrangian multipliers. The Lagrangian relaxed problem can be decomposed into two sub-problems by separating continuous variables from integer ones. The sub-problem that relates to continuous variables is a linear programming problem which can be solved using standard software package OSL, while the other sub-problem is an integer programming problem which can be solved optimally by further decomposition. The subgradient optimization method is used to update Lagrangian multipliers. A production order scheduling simulation system for Baosteel is developed by embedding the above Lagrangian heuristics. Computational results for problems with up to 100 orders show that the proposed Lagrangian relaxation method is stable and can find good solutions within a reasonable time.     

On the integer-valued variables in the linear vertex packing problem

Given a graph with weights on vertices, the vertex packing problem consists of finding a vertex packing (i.e. a set of vertices, no two of them being adjacent) of maximum weight. A linear relaxation of one binary programming formulation of this problem has these two well-known properties: (i) every basic solution is (0, 1/2, 1)-valued, (ii) in an optimum linear solution, an integer-valued variable keeps the same value in an optimum binary solution.As an answer to an open problem from Nemhauser and Trotter, it is shown that there is a unique maximal set of variables which are integral in optimal (VLP) solutions.This research was supported by National Research Council of Canada GRANT A8528 and RD 804.     

A Stackelberg equilibrium for a missile procurement problem

This paper deals with a procurement problem of missiles involving the efficient assignment of the missiles to some targets. Within a fixed amount of budget, a leader purchases several types of missiles, by which he aims to damage as much value as possible a follower hides into some facilities later. The effectiveness of the missile depends on the type of missile and facility. A payoff of the game is the expected amount of destroyed value. The problem is generalized as a two-person zero-sum game of distributing discrete resources with a leader and a follower. Our problem is to derive a Stackelberg equilibrium for the game. This type of game has an abundance of applications. The problem is first formulated into an integer programming problem with a non-separable objective function of variables and it is further equivalently transformed into a maximin integer knapsack problem. We propose three exacts methods and an approximation method for an optimal solution.     

Multiple Neutral Data Fitting

A method is proposed for estimating the relationship between a number of variables; this differs from regression where the emphasis is on predicting one of the variables. Regression assumes that only one of the variables has error or natural variability, whereas our technique does not make this assumption; instead, it treats all variables in the same way and produces models which are units invariant – this is important for ensuring physically meaningful relationships. It is thus superior to orthogonal regression in that it does not suffer from being scale-dependent. We show that the solution to the estimation problem is a unique and global optimum. For two variables the method has appeared under different names in various disciplines, with two Nobel laureates having published work on it.     

A linear differential game on a plane with a minimum-type functional

A differential game on a plane with a functional in the form of the minimum, with respect to time, of a certain prescribed phase vector function (quality function) is considered. It is proved that the game value is constant outside a certain bounded region, consisting of two parts. In the first subregion, the value is equal to the quality function, and in the second it satisfies Bellman's equation. For the constant-value region, where the players’ optimum strategies are not unique, single-valued guaranteeing players’ strategies are proposed. The results of a numerical investigation of the problem are presented.     

运筹学中若干线性目标规划和线性规划的人工智能-代数解法

运筹学中的线性目标规划和线性规划问题一直分别采用修正单纯形法和单纯形法求解.当变量稍多时计算还是有些繁琐、费时,最近作者通过研究发现,可应用人工智能-代数方法求得这两类问题的解,而且具有相当广泛的适用性.若干例题说明,本法的结果和传统方法的结果由于传统算法在计算中发生的错误,除少数例外大都是一致的.本文的一个 重要目的是希望和广大读者一起研究该方法是否具有通用性.