Optimization Methods: Integer Programming Integer Linear Programming 1. Module 7 Lecture Notes 1. Integer Linear Programming

Similar documents
GSLM Operations Research II Fall 13/14

Complex Numbers. Now we also saw that if a and b were both positive then ab = a b. For a second let s forget that restriction and do the following.

LECTURE NOTES Duality Theory, Sensitivity Analysis, and Parametric Programming

2x x l. Module 3: Element Properties Lecture 4: Lagrange and Serendipity Elements

R s s f. m y s. SPH3UW Unit 7.3 Spherical Concave Mirrors Page 1 of 12. Notes

5 The Primal-Dual Method

A NEW APPROACH FOR SOLVING LINEAR FUZZY FRACTIONAL TRANSPORTATION PROBLEM

LOOP ANALYSIS. The second systematic technique to determine all currents and voltages in a circuit

Support Vector Machines

Outline. Third Programming Project Two-Dimensional Arrays. Files You Can Download. Exercise 8 Linear Regression. General Regression

Classification / Regression Support Vector Machines

Intra-Parametric Analysis of a Fuzzy MOLP

Solving two-person zero-sum game by Matlab

Solutions to Programming Assignment Five Interpolation and Numerical Differentiation

Analysis of Continuous Beams in General

SENSITIVITY ANALYSIS IN LINEAR PROGRAMMING USING A CALCULATOR

For instance, ; the five basic number-sets are increasingly more n A B & B A A = B (1)

Multiple Instance Learning via Multiple Kernel Learning *

Assignment # 2. Farrukh Jabeen Algorithms 510 Assignment #2 Due Date: June 15, 2009.

Mathematics 256 a course in differential equations for engineering students

A Facet Generation Procedure. for solving 0/1 integer programs

A new paradigm of fuzzy control point in space curve

Accounting for the Use of Different Length Scale Factors in x, y and z Directions

Multicast Tree Rearrangement to Recover Node Failures. in Overlay Multicast Networks

Introduction to Geometrical Optics - a 2D ray tracing Excel model for spherical mirrors - Part 2

Improving Low Density Parity Check Codes Over the Erasure Channel. The Nelder Mead Downhill Simplex Method. Scott Stransky

Sum of Linear and Fractional Multiobjective Programming Problem under Fuzzy Rules Constraints

3D vector computer graphics

The Greedy Method. Outline and Reading. Change Money Problem. Greedy Algorithms. Applications of the Greedy Strategy. The Greedy Method Technique

Support Vector Machines. CS534 - Machine Learning

6.854 Advanced Algorithms Petar Maymounkov Problem Set 11 (November 23, 2005) With: Benjamin Rossman, Oren Weimann, and Pouya Kheradpour

i v v 6 i 2 i 3 v + (1) (2) (3) (4) (5) Substituting (4) and (5) into (3) (6) = 2 (7) (5) and (6) (8) (4) and (6) ˆ

Unbalanced Transportation Problems in Fuzzy Environment using Centroid Ranking Technique

Loop Transformations, Dependences, and Parallelization

Generating Fuzzy Term Sets for Software Project Attributes using and Real Coded Genetic Algorithms

Smoothing Spline ANOVA for variable screening

Proposed Simplex Method For Fuzzy Linear Programming With Fuzziness at the Right Hand Side

MODULE - 9 LECTURE NOTES 1 FUZZY OPTIMIZATION

Heuristic Methods for Locating Emergency Facilities

The Codesign Challenge

Solitary and Traveling Wave Solutions to a Model. of Long Range Diffusion Involving Flux with. Stability Analysis

UNIT 2 : INEQUALITIES AND CONVEX SETS

On Some Entertaining Applications of the Concept of Set in Computer Science Course

On-line Scheduling Algorithm with Precedence Constraint in Embeded Real-time System

A New Approach For the Ranking of Fuzzy Sets With Different Heights

Introduction. Leslie Lamports Time, Clocks & the Ordering of Events in a Distributed System. Overview. Introduction Concepts: Time

ON SOME ENTERTAINING APPLICATIONS OF THE CONCEPT OF SET IN COMPUTER SCIENCE COURSE

Lecture 5: Multilayer Perceptrons

An Adaptive Hyperplane Approach for Multiple Objective Optimization Problems with Complex Constraints

Machine Learning 9. week

Today Using Fourier-Motzkin elimination for code generation Using Fourier-Motzkin elimination for determining schedule constraints

INF 4300 Support Vector Machine Classifiers (SVM) Anne Solberg

Reading. 14. Subdivision curves. Recommended:

An Application of Network Simplex Method for Minimum Cost Flow Problems

Type-2 Fuzzy Non-uniform Rational B-spline Model with Type-2 Fuzzy Data

Algorithm To Convert A Decimal To A Fraction

Meta-heuristics for Multidimensional Knapsack Problems

Vertex Degrees in Planar Maps

An Optimal Algorithm for Prufer Codes *

Parallelism for Nested Loops with Non-uniform and Flow Dependences

TN348: Openlab Module - Colocalization

A Simple and Efficient Goal Programming Model for Computing of Fuzzy Linear Regression Parameters with Considering Outliers

Design for Reliability: Case Studies in Manufacturing Process Synthesis

An Efficient Fault-Tolerant Multi-Bus Data Scheduling Algorithm Based on Replication and Deallocation

Compiler Design. Spring Register Allocation. Sample Exercises and Solutions. Prof. Pedro C. Diniz

User Behavior Recognition based on Clustering for the Smart Home

Intro. Iterators. 1. Access

Review of approximation techniques

Parallel Branch and Bound Algorithm - A comparison between serial, OpenMP and MPI implementations

Lecture #15 Lecture Notes

Loop Permutation. Loop Transformations for Parallelism & Locality. Legality of Loop Interchange. Loop Interchange (cont)

X- Chart Using ANOM Approach

NAG Fortran Library Chapter Introduction. G10 Smoothing in Statistics

Data Representation in Digital Design, a Single Conversion Equation and a Formal Languages Approach

CMPS 10 Introduction to Computer Science Lecture Notes

Quality Improvement Algorithm for Tetrahedral Mesh Based on Optimal Delaunay Triangulation

Kinematics of pantograph masts

Topology Design using LS-TaSC Version 2 and LS-DYNA

Course Introduction. Algorithm 8/31/2017. COSC 320 Advanced Data Structures and Algorithms. COSC 320 Advanced Data Structures and Algorithms

ECONOMICS 452* -- Stata 11 Tutorial 6. Stata 11 Tutorial 6. TOPIC: Representing Multi-Category Categorical Variables with Dummy Variable Regressors

A CLASS OF TRANSFORMED EFFICIENT RATIO ESTIMATORS OF FINITE POPULATION MEAN. Department of Statistics, Islamia College, Peshawar, Pakistan 2

Mesh Editing in ROI with Dual Laplacian

S1 Note. Basis functions.

Load Balancing for Hex-Cell Interconnection Network

y and the total sum of

Large Margin Nearest Neighbor Classifiers

Harmonic Coordinates for Character Articulation PIXAR

AMath 483/583 Lecture 21 May 13, Notes: Notes: Jacobi iteration. Notes: Jacobi with OpenMP coarse grain

Multicriteria Decision Making

Kent State University CS 4/ Design and Analysis of Algorithms. Dept. of Math & Computer Science LECT-16. Dynamic Programming

MIXED INTEGER-DISCRETE-CONTINUOUS OPTIMIZATION BY DIFFERENTIAL EVOLUTION Part 1: the optimization method

AP PHYSICS B 2008 SCORING GUIDELINES

Programming in Fortran 90 : 2017/2018

Edge Detection in Noisy Images Using the Support Vector Machines

Esc101 Lecture 1 st April, 2008 Generating Permutation

CHAPTER 4. Applications of Boolean Algebra/ Minterm and Maxterm Expansions

Machine Learning. Support Vector Machines. (contains material adapted from talks by Constantin F. Aliferis & Ioannis Tsamardinos, and Martin Law)

Pattern Classification of Back-Propagation Algorithm Using Exclusive Connecting Network

User Authentication Based On Behavioral Mouse Dynamics Biometrics

Parallel Numerics. 1 Preconditioning & Iterative Solvers (From 2016)

Transcription:

Optzaton Methods: Integer Prograng Integer Lnear Prograng Module Lecture Notes Integer Lnear Prograng Introducton In all the prevous lectures n lnear prograng dscussed so far, the desgn varables consdered are supposed to take any real value. However n practcal probles lke nzaton of labor needed n a proect, t akes lttle sense n assgnng a value lke 5.6 to the nuber of labourers. In stuatons lke ths, one natural dea for obtanng an nteger soluton s to gnore the nteger constrants and use any of the technques prevously dscussed and then round-off the soluton to the nearest nteger value. However, there are several fundaental probles n usng ths approach:. The rounded-off solutons ay not be feasble.. The obectve functon value gven by the rounded-off solutons (even f soe are feasble) ay not be the optal one. 3. Even f soe of the rounded-off solutons are optal, checkng all the rounded-off solutons s coputatonally epensve ( for an n varable proble) n possble round-off values to be consdered Types of Integer Prograng When all the varables n an optzaton proble are restrcted to take only nteger values, t s called an all nteger prograng proble. When the varables are restrcted to take only dscrete values, the proble s called a dscrete prograng proble. When only soe varable values are restrcted to take nteger or dscrete, t s called ed nteger or dscrete prograng proble. When the varables are constraned to take values of ether zero or, then the proble s called zero one prograng proble. Integer Lnear Prograng Integer Lnear Prograng (ILP) s an etenson of lnear prograng, wth an addtonal restrcton that the varables should be nteger valued. The standard for of an ILP s of the for, D Nagesh Kuar, IISc, Bangalore ML

Optzaton Methods: Integer Prograng Integer Lnear Prograng a subect to c T X AX b X 0 X ust be nteger valued The assocated lnear progra droppng the nteger restrctons s called lnear relaaton LR. Thus, LR s less constraned than ILP. If the obectve functon coeffcents are nteger, then for nzaton, the optal obectve for ILP s greater than or equal to the rounded-off value of the optal obectve for LR. For azaton, the optal obectve for ILP s less than or equal to the rounded-off value of the optal obectve for LR. For a nzaton ILP, the optal obectve value for LR s less than or equal to the optal obectve for ILP and for a azaton ILP, the optal obectve value for LR s greater than or equal to that of ILP. If LR s nfeasble, then ILP s also nfeasble. Also, f LR s optzed by nteger varables, then that soluton s feasble and optal for IP. A ost popular ethod used for solvng all-nteger and ed-nteger lnear prograng probles s the cuttng plane ethod by Goory (Goory, 95). Goory s Cuttng Plane Method for All Integer Prograng Consder the followng optzaton proble. Maze subect to Z = 3 3, + 9 + 0 6 45 and are ntegers The graphcal soluton for the lnear relaaton of ths proble s shown below. D Nagesh Kuar, IISc, Bangalore ML

Optzaton Methods: Integer Prograng Integer Lnear Prograng 3 6 5 A 4 3 3 + 9 45 B ( 4 5,3 3 ) Z = 4 6 0 D C 0 3 4 5 6 It can be seen that the soluton s = 5, 4 = 3 3 and the optal value of Z = 4. The feasble solutons accountng the nteger constrants are shown by red dots. These ponts are called nteger lattce ponts. The orgnal feasble regon s reduced to a new feasble regon by ncludng soe addtonal constrants such that an etree pont of the new feasble regon becoes an optal soluton after accountng for the nteger constrants. The graphcal soluton for the eaple prevously dscussed takng and as ntegers are shown below. Two addtonal constrants (MN and OP) are ncluded so that the orgnal feasble regon ABCD s reduced to a new feasble regon AEFGCD. Thus the soluton for ths ILP s =, 3 and the optal value s Z = 5. 4 = D Nagesh Kuar, IISc, Bangalore ML

Optzaton Methods: Integer Prograng Integer Lnear Prograng 4 Addtonal constrants M 6 5 A O 4 3 E F (4,3) B ( 4 5,3 3 ) Z = 4 G Z = 5 P D C N 0 0 3 4 5 6 Goary proposed a systeatc ethod to develop these addtonal constrants known as Goory constrants. Generaton of Goory Constrants: Let the fnal tableau of an LP proble consst of n basc varables (orgnal varables) and non basc varables (slack varables) as shown n the table below. The basc varables are represented as (=,,,n) and the non basc varables are represented as y (=,,,). D Nagesh Kuar, IISc, Bangalore ML

Optzaton Methods: Integer Prograng Integer Lnear Prograng 5 Bass Z Table Varables n y y y y br Z 0 0 0 0 c c c c b 0 0 0 0 c c c c b 0 0 0 0 c c c c b 0 0 0 0 c 3 c 3 c 3 c 3 b n 0 0 0 0 c 4 c 4 c 4 c 4 b n Choose any basc varable wth the hghest fractonal value. If there s a te between two basc varables, arbtrarly choose any of the as th. Then fro the equaton of table, = b c = y..() Epress both b and c as an nteger value plus a fractonal part. b c = b = c + β + α...()...(3) where b, c denote the nteger part and β, α denote the fractonal part. β wll be a strctly postve fracton ( 0 < β < ) and α s a non-negatve fracton ( 0 < ) Substtutng equatons () and (3) n (), equaton () can be wrtten as = α. β α y = b c y.(4) = For all the varables and y to be ntegers, the rght hand sde of equaton (4) should be an nteger. β α y = nteger (5) = D Nagesh Kuar, IISc, Bangalore ML

Optzaton Methods: Integer Prograng Integer Lnear Prograng 6 Snce α are non-negatve ntegers and y are non-negatve ntegers, the ter α wll always be a non-negatve nuber. Thus we have, β < α y β.(6) = Hence the constrant can be epressed as β α y 0.() = = y By ntroducng a slack varable s (whch should also be an nteger), the Goory constrant can be wrtten as s = α y = β.(8) General procedure for solvng ILP:. Solve the gven proble as an ordnary LP proble neglectng the nteger constrants. If the optu values of the varables are ntegers tself, then there s nothng ore to be done.. If any of the basc varables has fractonal values, ntroduce the Goory constrants as dscussed n the prevous secton. Insert a new row wth the coeffcents of ths constrant, to the fnal tableau of the ordnary LP proble (Table ). 3. Solve ths by applyng the dual sple ethod. Snce the value of y = 0 n Table, the Goory constrant equaton becoes s = β whch s a negatve value and thus nfeasble. Dual sple ethod s used to obtan a new optal soluton that satsfes the Goory constrant. 4. Check whether the new soluton s all-nteger or not. If all values are not ntegers, then a new Goory constrant s developed fro the new sple tableau and the dual sple ethod s appled agan. 5. Ths process s contnued untl an optal nteger soluton s obtaned or t shows that the proble has no feasble nteger soluton. D Nagesh Kuar, IISc, Bangalore ML

Optzaton Methods: Integer Prograng Integer Lnear Prograng Thus, the fundaental dea behnd cuttng planes s to add constrants to a lnear progra untl the optal basc feasble soluton takes on nteger values. Goory cuts have the property that they can be generated for any nteger progra, but has the dsadvantage that the nuber of constrants generated can be enorous dependng upon the nuber of varables. D Nagesh Kuar, IISc, Bangalore ML