An example for the L A TEX package ORiONeng.sty

Similar documents
In Chapters 3 and 4 we introduced linear programming

A Brief Introduction to Proofs

Ch.03 Solving LP Models. Management Science / Prof. Bonghyun Ahn

LOWER BOUNDS FOR THE MAXIMUM NUMBER OF SOLUTIONS GENERATED BY THE SIMPLEX METHOD

Mathematics with Maple

Algebraic Simplex Active Learning Module 4

The simplex algorithm

Footnotes to Linear Algebra (MA 540 fall 2013), T. Goodwillie, Bases

Boolean Algebra CHAPTER 15

Math From Scratch Lesson 24: The Rational Numbers

OPTIMISATION 3: NOTES ON THE SIMPLEX ALGORITHM

Properties of a Simple Variant of Klee-Minty s LP and Their Proof

Using Simulation and Symbolic Computing in Linear Programming

The following techniques for methods of proofs are discussed in our text: - Vacuous proof - Trivial proof

Chapter 4 The Simplex Algorithm Part I

Bilinear Transformations via a Computer Algebra System

Introduction. How to use this book. Linear algebra. Mathematica. Mathematica cells

1.4 The Jacobian of a map

A CONSTRUCTION OF ARITHMETIC PROGRESSION-FREE SEQUENCES AND ITS ANALYSIS

Knapsack Feasibility as an Absolute Value Equation Solvable by Successive Linear Programming

Introduction to the Simplex Algorithm Active Learning Module 3

Introduction to integer programming II

PERIODIC POINTS OF THE FAMILY OF TENT MAPS

A Redundant Klee-Minty Construction with All the Redundant Constraints Touching the Feasible Region

Title of the Book. A. U. Thor XYZ University

Solutions of the diophantine equation 2 x + p y = z 2

CPSC 518 Introduction to Computer Algebra Schönhage and Strassen s Algorithm for Integer Multiplication

COMPUTATIONAL COMPLEXITY OF PARAMETRIC LINEAR PROGRAMMING +

Lecture 9 Tuesday, 4/20/10. Linear Programming

Simplex Method. Dr. rer.pol. Sudaryanto

A Type-Coercion Problem in Computer Algebra

Critical Reading of Optimization Methods for Logical Inference [1]

arxiv:math/ v1 [math.oa] 9 May 2005

Modified Berlekamp-Massey algorithm for approximating the k-error linear complexity of binary sequences

On the non-existence of maximal inference degrees for language identification

Pattern-Oriented Analysis and Design (POAD) Theory

Isoperimetric Sets of Integers

Math 200 A and B: Linear Algebra Spring Term 2007 Course Description

Identifying Graph Automorphisms Using Determining Sets

R u t c o r Research R e p o r t. Relations of Threshold and k-interval Boolean Functions. David Kronus a. RRR , April 2008

Ω R n is called the constraint set or feasible set. x 1

Michælmas 2012 Operations Research III/IV 1

Solving LFP by Converting it into a Single LP

TRANSPORTATION PROBLEMS

THe use of mixed data-type expressions, although not

Elementary Linear Algebra, Second Edition, by Spence, Insel, and Friedberg. ISBN Pearson Education, Inc., Upper Saddle River, NJ.

Limits: How to approach them?

Follow links Class Use and other Permissions. For more information, send to:

Linear Programming: Computer Solution and Sensitivity Analysis

ME DYNAMICAL SYSTEMS SPRING SEMESTER 2009

The spectra of super line multigraphs

OPRE 6201 : 3. Special Cases

Week 2. The Simplex method was developed by Dantzig in the late 40-ties.

This is logically equivalent to the conjunction of the positive assertion Minimal Arithmetic and Representability

A COMPUTER PROGRAM FOR THE DECONVOLUTION OF THERMOLUMINESCENCE GLOW CURVES K. S. Chung 1,, H. S. Choe 1, J. I. Lee 2, J. L. Kim 2 and S. Y.

Primitive sets in a lattice

The Simplex Algorithm

Tewodros Amdeberhan, Dante Manna and Victor H. Moll Department of Mathematics, Tulane University New Orleans, LA 70118

Description Logics: an Introductory Course on a Nice Family of Logics. Day 2: Tableau Algorithms. Uli Sattler

WeatherHawk Weather Station Protocol

Appendix B for The Evolution of Strategic Sophistication (Intended for Online Publication)

Calculus and Maximization I

Chapter One. The Real Number System

11 Division Mod n, Linear Integer Equations, Random Numbers, The Fundamental Theorem of Arithmetic

Important Properties of R

21 Markov Decision Processes

). In an old paper [11], I. N. Sanov

Two hours. Note that the last two pages contain inference rules for natural deduction UNIVERSITY OF MANCHESTER SCHOOL OF COMPUTER SCIENCE

1 The linear algebra of linear programs (March 15 and 22, 2015)

TIM 206 Lecture 3: The Simplex Method

Math 414, Fall 2016, Test I

A Simpler and Tighter Redundant Klee-Minty Construction

Mergible States in Large NFA

OPERATIONS RESEARCH. Michał Kulej. Business Information Systems

Introduction The Fibonacci numbers are traditionally described as a sequence F. = 0 and τ τ 1. =. To obtain the Fibonacci numbers as a

15-451/651: Design & Analysis of Algorithms September 13, 2018 Lecture #6: Streaming Algorithms last changed: August 30, 2018

Chapter 5 Linear Programming (LP)

On a correlation inequality of Farr

Introduction. An Introduction to Algorithms and Data Structures

The Simplex Method of Linear Programming

Fast Decoding Of Alternant Codes Using A Divison-Free Analog Of An Accelerated Berlekamp-Massey Algorithm

Math Models of OR: Handling Upper Bounds in Simplex

Lecture 5. Theorems of Alternatives and Self-Dual Embedding

ECOM Discrete Mathematics

Spring 2017 CO 250 Course Notes TABLE OF CONTENTS. richardwu.ca. CO 250 Course Notes. Introduction to Optimization

Math 309 Notes and Homework for Days 4-6

An Algorithm for Automatic Demonstration of Logical Theorems

LINEAR PROGRAMMING 2. In many business and policy making situations the following type of problem is encountered:

Laver Tables A Direct Approach

NUMBERS WITH INTEGER COMPLEXITY CLOSE TO THE LOWER BOUND

DO NOT COPY DO NOT COPY

Linear Algebra II. 2 Matrices. Notes 2 21st October Matrix algebra

The Laplace Rule of Succession Under A General Prior

IMA Preprint Series # 2176

OF ALL THE CHEMISTRY RELATED SOFTWARE

Continuity. Chapter 4

Absolute Value Programming

Pressure Swing Distillation with Aspen Plus V8.0

Myung-Hwan Kim and Byeong-Kweon Oh. Department of Mathematics, Seoul National University, Seoul , Korea

INDUCTION AND RECURSION

Transcription:

Operations Research Society of South Africa Submitted for publication in ORiON Operasionele Navorsingsvereniging van Suid-Afrika An example for the L A TEX package ORiONeng.sty Authors identities suppressed: Blind refereeing copy Abstract The format for papers published in ORiON is provided. Examples of how environments, such as tables, should be presented are also provided. The editorial format and standards of ORiON are maintained throughout this paper. The purpose is not to provide a mathematically correct or coherent content to this document, but rather to create a paper using specific examples of standard L A TEX elements. Key words: Style sheet, examples, L A TEX 1 Introduction This paper contains examples of several functions used in ORiON s L A TEX package. The purpose of this paper is not to provide readers with sensible content or proper logic. The theorems and definitions used in this paper may therefore appear not to make sense. The purpose is, however, not to provide the correct theorems and definitions, but rather to provide examples of these entities. The examples demonstrate how the standard environments and elements are created within the ORiON package. This paper also provides examples of the preferred format in which equations, references and other standard elements are handled in ORION. The examples provided in this paper are not exhaustive. Should an author wish to include specific or exceptional situations in a paper, the business manager of ORiON can be contacted for guidance. Frequently asked questions will be addressed as examples in further updates of this paper. There are various manuals available for L A TEX and TEX[, 5]. Should these not provide sufficient guidance, the internet can also be used as a source of information. Text is entered into the source document in the same manner as for standard L A TEX documents. The environments in Table 1 have already been defined. This document is funded by ORSSA. 1

Authors identities suppressed: Blind refereeing copy Text algorithm axiom corollary definition example heuristic lemma proposition theorem Table 1: Pre-defined environments in ORiONeng.sty To solve LPs with the help of the simplex algorithm, knowledge of linear algebra is required. Amongst other things, the concept of linearly dependent vectors is important. Definition 1 If a j with j = 1,, 3,..., r is a set of vectors, then the equation k 1 a 1 + k a +... + k r a r = 0 will have at least one solution, which is k 1 = k =... = k r = 0. If this is the only solution, the set, a j, is called linearly independent. If other solutions exist, the set is called linearly dependent. A set of column vectors, a j, with j = 1,, 3,..., k, is linear independent or linearly dependent. Assume the set of vectors is linearly dependent. Under these circumstances it is possible to reduce the number of positive variables step by step until the columns which correspond with the remaining positive variables are linearly independent. If a j, with j = 1,, 3,..., n, is linearly dependent, real numbers (not all equal to zero) exist, so that k a j = 0. (1) Equation (1) can be used to reduce an x j, say x r, in k x j a j = b, x j > 0 () to zero. Suppose any vector a r of the k vectors in equation (1), for which 0, can be expressed in terms of the remaining k 1 vectors. Thus a r = j r Substitute equation (3) into equation (), so that p ( x j x r a j. (3) ) a j = b. (4) This solution does not have more than k 1 non-zero variables. If a r is selected arbitrarily, there would be no assurance that some of the variables may not perhaps be negative. If a r is, however, selected carefully, the remaining k 1 variables will still be non-negative. Select a r so that x j x r 0 j = 1,..., p. (5)

ORiON: Example 3 For any j, where = 0, equation (5) will automatically hold. If 0, then equation (5), divided by, will give one of the following inequalities: x j From (6) and (7) follows that xr than any x j x r 0 if > 0, or (6) x j x r 0 if < 0. (7) can not be greater than any x j < 0. If we choose the vector a r so that > 0 and also not less { } x r xj = min, > 0, (8) j then all the variables in equation (4) will be non-negative. A feasible solution with no more then k 1 non-zero variables has been found. If this process of reduction is continued, the remaining x j s will eventually be linearly independent and the solution will be a basic feasible solution. The theorem which summarises the theory of the solution of LPs with the simplex algorithm, is the fundamental theorem of LP. Theorem (Fundamental theorem of linear programming (LP)) This theorem consists of two parts: 1. If a feasible solution exists for an LP, a basic feasible solution exists.. If an optimal solution exists for an LP, an optimal basic feasible solution exists. Only a partial proof for Theorem is presented. Proof: Assume that a feasible solution exists with k < n variables having, a positive value. Number these variables in such a way that the first k variables have a positive value. The solution can then be formulated as k x j a j = b with x j > 0 where j = 1,,..., k and x j = 0 for j = k +1,..., n. The vectors a j which are associated with the positive values can either be linearly independent or linearly dependent [6]. For a more detailed proof, see [4]. It is obviously much easier to solve LPs using software packages. A few of these programs are listed in Table.

4 Authors identities suppressed: Blind refereeing copy Software LINDO & LINGO QM for Windows Advantages Programmable Built-in mathematical functions User-friendly interface Interface with Excel 3 Algorithms Table : LP software with their respective advantages. The binary search algorithm starts its search in the centre of an alphabetical list of words. The algorithm compares the word in the middle of the list with the word that must be found (WRD) to test whether the middle word is, in fact, WRD. If this is not the case, it means that WRD is either after or before the word in the middle of the alphabetical list. If WRD comes alphabetically before the middle word, the word halfway between the first word and the middle word is tested to determine whether that word is WRD, while the other half (after the middle word) of the list is discarded. If WRD comes after the middle word, the word halfway between the middle and the last word is tested, while the first half of the list is discarded. This process of dividing the remaining list in half is repeated until the word is located or the list has been exhausted. This implies that only one half of the list has to be tested after the first step, one quarter of the list after the second step, one eighth of the list after the third step, etc. (The middle word is the word in the position in the list.) 1+n -th Algorithm 3 (Binary search algorithm) Find a word in an alphabetical list of words by comparing the middle word with the word that must be found and by dividing the list, if required, into sublists. 1. Let α and β be the first and last positions of the sublist. Initially α = 1 and β = n. Thus α 1 and β n. WRD the word to be located.. Repeat this step until α β, else go to 3. (a) Determine k, the middle position (between α and β). (b) If W (k) = WRD, stop. Output: k. α + β k (c) If WRD > W (k), then β k 1. Go to a. (d) If WRD < W (k), then α k + 1. Go to a. 3. Stop. Output: Not in the list. Let us examine the (worst case) complexity of the binary search algorithm. Note that step 1 and step 3 are both performed only once. In both of these steps, a constant number of calculations are performed. This means that the complexity of these steps is O(1). Let

ORiON: Example 5 B(n) be the maximum number of times that step must be repeated. After step has been performed for the first time, the remaining list will have a maximum of n elements to consider. This means that ( ) n B(n) = 1 + B. (9) We use induction to prove that B(n) 1 + log n (10) for each positive integer n. Equation (10) applies if n = 1. Accept that this result holds for an integer, k. We must then show that this inequality holds for k + 1. We therefore accept that B(k) 1 + log k. (11) From equation (9) we know that B(k + 1) = 1 + B, but from equation (11) it follows that B We also know that This means that log 1 + log. ( ) k + 1 log. ( ) k + 1 B(k + 1) + log + log (k + 1) log () 1 + log (k + 1), which gives the desired result that B(n) = O(log n). 4 Matrices The variable, x 3, enters the base with a value of 0 1 = 0, which is smaller than or equal to the upper bound. Thus a simplex iteration can be performed. The column below x 3 is given by 1 0 1 1 1 B1 1 p 3 = 0 1 0 1 = 1. 0 0 The elementary matrix is given by ( ) 1 1) ( E 1 = 1 ( 1) 0 1 0 0 1 0. 0 1

6 Authors identities suppressed: Blind refereeing copy We can then calculate B 1 = E 1 B 1 1 = 1 0 0 1 1 0 0 0 1 1 0 1 0 1 0 = 1 0 1 1 1 1. After the iteration, the respective values are given by 1 1 1 1 0 0 A = 1 0 1 0 1 0 B 1 = 1 1 0 0 0 1 Go to step 5. 1 0 1 1 1 1 b = c bv () = [ 1, 0, 1] c nbv () = [, 0, 0] x bv () = [x 3, s, x ] x nbv () = [ x 1, s 3, s 1 ]. Step 5: The current right hand sides are given by b() = B 1 b = 1 0 1 1 1 1 60 30 80 = 60 30 80 The basic variables x 3, s and x have the values 0, 50 and 80 respectively. None of these variables occur at a value above their upper bounds. Step 6: All the right hand sides are positive, which means that an optimal solution has been found. The optimal solution is given by x 1 = 80 x 1 = 80 0 = 80 x = 80 x = 80 80 = 0 x 3 = 0 s 1 = 0 s 1 = 0 0 = 0 s = 50 s 3 = 0. The optimal objection function value is 0 50 80. Objection function = c 1 x 1 + c x + c 3 x 3 = ( )(80) + ( 1)(0) + (1)(0) = 140. References [1] Chartrand G & Oellermann OR, 1993, Applied and Algorithmic Graph Theory, McGraw Hill, Singapore, pp 43 46

ORiON: Example 7 [] Doob M, 1996, TEX starting from 1, Springer-Verlag, New York. [3] Gass SI & Harris CM, 1996, Encyclopedia of Operations Research and Management Science, Kluwer Academic Publishers, Boston. [4] Hadley G, 196, Linear programming, Addison-Wesley Publishing Company Inc., Reading, Massachusetts. [5] Lamport L, 1994, L A TEX A document preparation system: User s guide and reference manual, Addison-Wesley Longman Inc., Reading, Massachusetts. [6] Nickolson WK, 1995, Linear algebra with applications, third edition, PWS Publishing Company, Boston.