Aninteger programming problem is amathematical optimization orfeasibility program in which some or all of the variables are restricted to beintegers. In many settings the term refers tointegerlinear programming (ILP), in which theobjective function and the constraints (other than the integer constraints) arelinear.

Integer programming isNP-complete. In particular, the special case of 0–1 integer linear programming, in which unknowns are binary, and only the restrictions must be satisfied, is one ofKarp's 21 NP-complete problems.^[1]

If some decision variables are not discrete, the problem is known as amixed-integer programming problem.^[2]

In integer linear programming, thecanonical form is distinct from thestandard form. An integer linear program in canonical form is expressed thus (note that it is the${\displaystyle \mathbf{x}}$ vector which is to be decided):^[3]

and an ILP in standard form is expressed as

where${\displaystyle \mathbf{c}\in {\mathbb{R}}^n,\mathbf{b}\in {\mathbb{R}}^m}$ are vectors and${\displaystyle A\in {\mathbb{R}}^{m\times n}}$ is a matrix. As with linear programs, ILPs not in standard form can beconverted to standard form by eliminating inequalities, introducing slack variables (${\displaystyle \mathbf{s}}$) and replacing variables that are not sign-constrained with the difference of two sign-constrained variables.

The plot on the right shows the following problem.

The feasible integer points are shown in red, and the red dashed lines indicate their convex hull, which is the smallest convex polyhedron that contains all of these points. The blue lines together with the coordinate axes define the polyhedron of the LP relaxation, which is given by the inequalities without the integrality constraint. The goal of the optimization is to move the black dashed line as far upward while still touching the polyhedron. The optimal solutions of the integer problem are the points${\displaystyle (1,2)}$ and${\displaystyle (2,2)}$ that both have an objective value of 2. The unique optimum of the relaxation is${\displaystyle (1.8,2.8)}$ with objective value of 2.8. If the solution of the relaxation is rounded to the nearest integers, it is not feasible for the ILP.

The following is a reduction from minimumvertex cover to integer programming that will serve as the proof of NP-hardness.

Let${\displaystyle G=(V,E)}$ be an undirected graph. Define a linear program as follows:

Given that the constraints limit${\displaystyle y_v}$ to either 0 or 1, any feasible solution to the integer program is a subset of vertices. The first constraint implies that at least one end point of every edge is included in this subset. Therefore, the solution describes a vertex cover. Additionally given some vertex cover C,${\displaystyle y_v}$ can be set to 1 for any${\displaystyle v\in C}$ and to 0 for any${\displaystyle v\notin C}$ thus giving us a feasible solution to the integer program. Thus we can conclude that if we minimize the sum of${\displaystyle y_v}$ we have also found the minimum vertex cover.^[4]

Mixed-integer linear programming (MILP) involves problems in which only some of the variables,${\displaystyle x_i}$, are constrained to be integers, while other variables are allowed to be non-integers.

Zero–one linear programming (orbinary integer programming) involves problems in which the variables are restricted to be either 0 or 1. Any bounded integer variable can be expressed as a combination ofbinary variables.^[5] For example, given an integer variable,${\displaystyle 0\le x\le U}$, the variable can be expressed using${\displaystyle \lfloor {\log }_{2}U\rfloor +1}$ binary variables:

$${\displaystyle x=x_1+2x_2+4x_3+\cdots +2^{\lfloor {\log }_{2}U\rfloor }{x}_{\lfloor {\log }_{2}U\rfloor +1}.}$$

There are two main reasons for using integer variables when modeling problems as a linear program:

1. The integer variables represent quantities that can only be integer. For example, it is not possible to build 3.7 cars.
2. The integer variables represent decisions (e.g. whether to include an edge in agraph) and so should only take on the value 0 or 1.

These considerations occur frequently in practice and so integer linear programming can be used in many applications areas, some of which are briefly described below.

Mixed-integer programming has many applications in industrial productions, including job-shop modelling. One important example happens in agriculturalproduction planning and involves determining production yield for several crops that can share resources (e.g. land, labor, capital, seeds, fertilizer, etc.). A possible objective is to maximize the total production, without exceeding the available resources. In some cases, this can be expressed in terms of a linear program, but the variables must be constrained to be integer.

These problems involve service and vehicle scheduling in transportation networks. For example, a problem may involve assigning buses or subways to individual routes so that a timetable can be met, and also to equip them with drivers. Here binary decision variables indicate whether a bus or subway is assigned to a route and whether a driver is assigned to a particular train or subway. The zero–one programming technique has been successfully applied to solve a project selection problem in which projects are mutually exclusive and/or technologically interdependent.

Territorial partitioning or districting problems consist of partitioning a geographical region into districts in order to plan some operations while considering different criteria or constraints. Some requirements for this problem are: contiguity, compactness, balance or equity, respect of natural boundaries, and socio-economic homogeneity. Some applications for this type of problem include: political districting, school districting, health services districting and waste management districting.

The goal of these problems is to design a network of lines to install so that a predefined set of communication requirements are met and the total cost of the network is minimal.^[6] This requires optimizing both the topology of the network along with setting the capacities of the various lines. In many cases, the capacities are constrained to be integer quantities. Usually there are, depending on the technology used, additional restrictions that can be modeled as linear inequalities with integer or binary variables.

The task of frequency planning inGSM mobile networks involves distributing available frequencies across the antennas so that users can be served and interference is minimized between the antennas.^[7] This problem can be formulated as an integer linear program in which binary variables indicate whether a frequency is assigned to an antenna.

• Cash flow matching
• Energy system optimization^[8][9]
• UAVguidance^[10][11]
• Transit maplayouting^[12]

The naive way to solve an ILP is to simply remove the constraint thatx is integer, solve the corresponding LP (called theLP relaxation of the ILP), and then round the entries of the solution to the LP relaxation. But, not only may this solution not be optimal, it may not even be feasible; that is, it may violate some constraint.

While in general the solution to LP relaxation will not be guaranteed to be integral, if the ILP has the form${\displaystyle max{\mathbf{c}}^{\mathrm{T}}\mathbf{x}}$ such that${\displaystyle A\mathbf{x}=\mathbf{b}}$ where${\displaystyle A}$ and${\displaystyle \mathbf{b}}$ have all integer entries and${\displaystyle A}$ istotally unimodular, then every basic feasible solution is integral. Consequently, the solution returned by thesimplex algorithm is guaranteed to be integral. To show that every basic feasible solution is integral, let${\displaystyle \mathbf{x}}$ be an arbitrary basic feasible solution . Since${\displaystyle \mathbf{x}}$ is feasible,we know that${\displaystyle A\mathbf{x}=\mathbf{b}}$. Let${\displaystyle {\mathbf{x}}_0=[x_{n_1},x_{n_2},\cdots ,x_{n_j}]}$ be the elements corresponding to the basis columns for the basic solution${\displaystyle \mathbf{x}}$. By definition of a basis, there is some square submatrix${\displaystyle B}$ of${\displaystyle A}$ with linearly independent columns such that${\displaystyle B{\mathbf{x}}_0=\mathbf{b}}$.

Since the columns of${\displaystyle B}$ are linearly independent and${\displaystyle B}$ is square,${\displaystyle B}$ is nonsingular,and therefore by assumption,${\displaystyle B}$ isunimodular and so${\displaystyle det(B)=\pm 1}$. Also, since${\displaystyle B}$ is nonsingular, it is invertible and therefore${\displaystyle {\mathbf{x}}_0=B^{-1}\mathbf{b}}$. By definition,${\displaystyle B^{-1}=\frac{B^{\mathrm{a}\mathrm{d}\mathrm{j}}}{det(B)}=\pm B^{\mathrm{a}\mathrm{d}\mathrm{j}}}$. Here${\displaystyle B^{\mathrm{a}\mathrm{d}\mathrm{j}}}$ denotes theadjugate of${\displaystyle B}$ and is integral because${\displaystyle B}$ is integral. Therefore,Thus, if the matrix${\displaystyle A}$ of an ILP is totally unimodular, rather than use an ILP algorithm, the simplex method can be used to solve the LP relaxation and the solution will be integer.

When the matrix${\displaystyle A}$ is not totally unimodular, there are a variety of algorithms that can be used to solve integer linear programs exactly. One class of algorithms arecutting plane methods, which work by solving the LP relaxation and then adding linear constraints that drive the solution towards being integer without excluding any integer feasible points.

Another class of algorithms are variants of thebranch and bound method. For example, thebranch and cut method that combines both branch and bound and cutting plane methods. Branch and bound algorithms have a number of advantages over algorithms that only use cutting planes. One advantage is that the algorithms can be terminated early and as long as at least one integral solution has been found, a feasible, although not necessarily optimal, solution can be returned. Further, the solutions of the LP relaxations can be used to provide a worst-case estimate of how far from optimality the returned solution is. Finally, branch and bound methods can be used to return multiple optimal solutions.

Suppose${\displaystyle A}$ is anm-by-n integer matrix and${\displaystyle \mathbf{b}}$ is anm-by-1 integer vector. We focus on the feasibility problem, which is to decide whether there exists ann-by-1 vector${\displaystyle \mathbf{x}}$ satisfying${\displaystyle A\mathbf{x}\le \mathbf{b}}$.

LetV be the maximum absolute value of the coefficients in${\displaystyle A}$ and${\displaystyle \mathbf{b}}$. Ifn (the number of variables) is a fixed constant, then the feasibility problem can be solved in time polynomial inm and logV. This is trivial for the casen=1. The casen=2 was solved in 1981 byHerbert Scarf.^[13] The general case was solved in 1983 byHendrik Lenstra, combining ideas byLászló Lovász andPeter van Emde Boas.^[14]Doignon's theorem asserts that an integer program is feasible whenever every subset of${\displaystyle 2^n}$ constraints is feasible; a method combining this result with algorithms forLP-type problems can be used to solve integer programs in time that is linear in${\displaystyle m}$ andfixed-parameter tractable (FPT) in${\displaystyle n}$, but possiblydoubly exponential in${\displaystyle n}$, with no dependence on${\displaystyle V}$.^[15]

In the special case of 0-1 ILP, Lenstra's algorithm is equivalent to complete enumeration: the number of all possible solutions is fixed (2ⁿ), and checking the feasibility of each solution can be done in time poly(m, logV). In the general case, where each variable can be an arbitrary integer, complete enumeration is impossible. Here, Lenstra's algorithm uses ideas fromGeometry of numbers. It transforms the original problem into an equivalent one with the following property: either the existence of a solution${\displaystyle \mathbf{x}}$ is obvious, or the value of${\displaystyle x_n}$ (then-th variable) belongs to an interval whose length is bounded by a function ofn. In the latter case, the problem is reduced to a bounded number of lower-dimensional problems. The run-time complexity of the algorithm has been improved in several steps:

• The original algorithm of Lenstra^[14] had run-time${\displaystyle 2^{O(n^3)}\cdot (m\cdot \log V{)}^{O(1)}}$.
• Kannan^[16] presented an improved algorithm with run-time${\displaystyle n^{O(n)}\cdot (m\cdot \log V{)}^{O(1)}}$.^[17]
• Frank and Tardos^[18] presented an improved algorithm with run-time${\displaystyle n^{2.5n}\cdot 2^{O(n)}\cdot (m\cdot \log V{)}^{O(1)}}$.^{[19][20]: Prop.8}
• Dadush^[21] presented an improved algorithm with run-time${\displaystyle n^n\cdot 2^{O(n)}\cdot (m\cdot \log V{)}^{O(1)}}$.
• Reis and Rothvoss^[22] presented an improved algorithm with run-time${\displaystyle (\log n{)}^{O(n)}\cdot (m\cdot \log V{)}^{O(1)}}$.

These algorithms can also be used formixed integer linear programs (MILP) - programs in which some variables are integer and some variables are real.^[23] The original algorithm of Lenstra^{[14]: Sec.5} has run-time${\displaystyle 2^{O(n^3)}\cdot poly(d,L)}$, where n is the number of integer variables, d is the number of continuous variables, andL is the binary encoding size of the problem. Using techniques from later algorithms, the factor${\displaystyle 2^{O(n^3)}}$ can be improved to${\displaystyle 2^{O(n\log n)}}$ or to${\displaystyle n^n}$.^[23]

Since integer linear programming isNP-hard, many problem instances are intractable and so heuristic methods must be used instead. For example,tabu search can be used to search for solutions to ILPs.^[24] To use tabu search to solve ILPs, moves can be defined as incrementing or decrementing an integer constrained variable of a feasible solution while keeping all other integer-constrained variables constant. The unrestricted variables are then solved for. Short-term memory can consist of previously tried solutions while medium-term memory can consist of values for the integer constrained variables that have resulted in high objective values (assuming the ILP is a maximization problem). Finally, long-term memory can guide the search towards integer values that have not previously been tried.

Other heuristic methods that can be applied to ILPs include

• Hill climbing
• Simulated annealing
• Reactive search optimization
• Ant colony optimization
• Hopfield neural networks

There are also a variety of other problem-specific heuristics, such as thek-opt heuristic for the traveling salesman problem. A disadvantage of heuristic methods is that if they fail to find a solution, it cannot be determined whether it is because there is no feasible solution or whether the algorithm simply was unable to find one. Further, it is usually impossible to quantify how close to optimal a solution returned by these methods are.

It is often the case that the matrix${\displaystyle A}$ that defines the integer program issparse. In particular, this occurs when the matrix has ablock structure, which is the case in many applications. The sparsity of the matrix can be measured as follows. Thegraph of${\displaystyle A}$ has vertices corresponding to columns of${\displaystyle A}$, and two columns form an edge if${\displaystyle A}$ has a row where both columns have nonzero entries. Equivalently, the vertices correspond to variables, and two variables form an edge if they share an inequality. Thesparsity measure${\displaystyle d}$ of${\displaystyle A}$ is the minimum of thetree-depth of the graph of${\displaystyle A}$ and thetree-depth of the graph of the transpose of${\displaystyle A}$. Let${\displaystyle a}$ be thenumeric measure of${\displaystyle A}$ defined as the maximum absolute value of any entry of${\displaystyle A}$. Let${\displaystyle n}$ be the number of variables of the integer program. Then it was shown in 2018^[25] that integer programming can be solved instrongly polynomial andfixed-parameter tractable time parameterized by${\displaystyle a}$ and${\displaystyle d}$. That is, for some computable function${\displaystyle f}$ and some constant${\displaystyle k}$, integer programming can be solved in time${\displaystyle f(a,d)n^k}$. In particular, the time is independent of the right-hand side${\displaystyle b}$ and objective function${\displaystyle c}$. Moreover, in contrast to the classical result of Lenstra, where the number${\displaystyle n}$ of variables is a parameter, here the number${\displaystyle n}$ of variables is a variable part of the input.

• Constrained least squares
• Diophantine equation – Polynomial equation whose integer solutions are sought

• George L. Nemhauser; Laurence A. Wolsey (1988).Integer and combinatorial optimization. Wiley.ISBN 978-0-471-82819-8.
• Alexander Schrijver (1998).Theory of linear and integer programming. John Wiley and Sons.ISBN 978-0-471-98232-6.
• Laurence A. Wolsey (1998).Integer programming. Wiley.ISBN 978-0-471-28366-9.
• Dimitris Bertsimas; Robert Weismantel (2005).Optimization over integers. Dynamic Ideas.ISBN 978-0-9759146-2-5.
• John K. Karlof (2006).Integer programming: theory and practice. CRC Press.ISBN 978-0-8493-1914-3.
• H. Paul Williams (2009).Logic and Integer Programming. Springer.ISBN 978-0-387-92279-9.
• Michael Jünger; Thomas M. Liebling; Denis Naddef;George Nemhauser;William R. Pulleyblank; Gerhard Reinelt; Giovanni Rinaldi; Laurence A. Wolsey, eds. (2009).50 Years of Integer Programming 1958-2008: From the Early Years to the State-of-the-Art. Springer.ISBN 978-3-540-68274-5.
• Der-San Chen; Robert G. Batson; Yu Dang (2010).Applied Integer Programming: Modeling and Solution. John Wiley and Sons.ISBN 978-0-470-37306-4.
• Gerard Sierksma; Yori Zwols (2015).Linear and Integer Optimization: Theory and Practice. CRC Press.ISBN 978-1-498-71016-9.

• A Tutorial on Integer Programming
• ConferenceInteger Programming and Combinatorial Optimization, IPCO
• The Aussois Combinatorial Optimization Workshop

• Combinatorial optimization
• NP-complete problems

• Articles with short description
• Short description is different from Wikidata