next up previous
Next: References Up: No Title Previous: Integer Linear Programming Model

Mixed Integer Programming-Network Model

Given a directed graph, G = (V,A) (see Bondy and Murty [2], and West [9]), the uncapacitated fixed charge network flow problem, (UF), is the following mixed-integer program

eqnarray1094

where is the single source vertex for all flow, is a collection of sink vertices with rational demand , denotes the flow in arc at variable cost . Arcs of also carry fixed cost with indicating that arc (i,j) is selected. In general we assume that all costs are nonnegative. Figure 1(a) shows a simple instance. gif

  

The (UF) model is an important one in discrete optimization because a wide variety of applied problems can be effectively modeled as uncapacitated fixed charge network flow problems. These problems include lot-sizing (planning production setups) problems, designing single commodity utility and logistics networks, and planning for warehouse and distribution systems.

Problem (UF) is NP-Hard for because it contains as a special case the hard Steiner Tree Problem (Garey and Johnson [4]). (There are polynomial algorithms for . We denote by the linear programming relaxation of any problem (P) and by its optimal value.) Thus, one way of obtaining approximate (UF) solutions is by solving its linear programming relaxation, and rounding ``up'' to one all fractional in the relaxation optimum. Unfortunately, that LP-relaxation usually provides a weak approximation to its corresponding discrete problem, i.e. often . For example, Figures 1 (b) and (e) show versus .

Now, we present an approach to obtain better LP-relaxation values as approximations to .

A well known structural description of the optimal solutions in the directed graph associated with (UF) can be stated as follows:

Rem1138

Thus there is no loss of generality in thinking of the flow in any arc in terms of the sinks for which it is bound. We will term solutions conforming to Remark 1 as extreme.

This observation led Rardin and Choe [7] to propose the multicommodity extended formulation, (MCE), which is obtained from (UF) by artificially extending the single commodity in (UF) to multiple commodity networks, one for each sink t. New variables, , represent the fraction of flow destined for sink t passing through arc (i,j).

  eqnarray1146

It is easy to see is optimal in (UF) if and only if there are w's so that is optimal in (MCE). The point of this extended formulation is to make optimal values of the y variables in LP-relaxations more accurately reflect their integer values. In each optimal will equal the fraction of total demand passing through arc (i,j). In , will equal the maximum of the fraction of particular demands satisfied through arc (i,j). For example in Figure gif(b) in because only 2/5 of total demand passes through arc (1,2). Corresponding results in Figure gif(c) have , the maximum of and . These tighter constraints on y-variables produce a relaxation optimal value of versus .

Computational testing (see for example Balakrishnan, Magnanti and Wong [1]) confirms that values are usually much sharper (closer to the solution value of the discrete model) than those of . However, need not be exact (e.g. in Figure gif), and there are instances in which there is a relatively large gap between and .

We investigate new commodity family extended formulations (CFE) that seek to improve on . (See Ng and Rardin [6]). This approach is to introduce a family tex2html_wrap_inline1784 of commodity variables

displaymath1385

Using these notions we are finally ready to develop a rudimentary statement of the commodity family extended formulation associated with a given commodity family tex2html_wrap_inline1784 and representation function tex2html_wrap_inline1794 .

    eqnarray1207

Constraints (gif) assume flow originates at source s for all commodities , the representation of the whole sink set. Demands are fulfilled at by singleton commodities . If tex2html_wrap_inline1784 contains more than those singleton I, this implies flows must split along the way from source to sink into commodities for smaller and smaller demand subsets.

Variables in the above formulation provide the required commodity interfaces. In particular,

displaymath1422

Each is a subset of T, so that the multicommodity model (MCE) is the special case where .

Figure gif(d) shows the result is an LP-relaxation strictly improving on . (The formulations of (CFE) and its relaxations are in the next section.) For this example, the stronger (CFE) form achieves integer optimal value .

Clearly not all subsets of T can be included in the tex2html_wrap_inline1784 of any practical (CFE) because their number grows exponentially with |T|. Still, there are an enormous variety of incomplete tex2html_wrap_inline1784 's of manageable size.

Ng and Rardin [6] described some properties of ``good'' families, most of which are trivial in the (MCE) setting but subtle and complex in richer commodity families. They also characterized which families should be considered. For each such family, they derived elegant axioms characterizing when a tex2html_wrap_inline1784 possesses the desired property.


next up previous
Next: References Up: No Title Previous: Integer Linear Programming Model

Peh H. Ng
Thu Mar 27 16:07:36 CST 1997