Synchronized Base Stock Policies

Size: px
Start display at page:

Download "Synchronized Base Stock Policies"

Transcription

1 Synchronized Base Stock Policies Ton de Kok Department of Technology Management Technische Universiteit Eindhoven

2 Contents Value Network Model Network transformation principles Divergent structures and allocation rules Optimality conditions on base stock levels Computational aspects Comparison against LP rolling schedule Managerial insights Conclusions and further research 2

3 Value Network Model 3

4 4

5 Multi-item multi-echelon inventory systems Arbitrary BOM no alternative recipees, no multi-sourcing Stochastic end-item demand I.i.d. periodic Linear holding costs Customer service objectives End-item linear shortage costs End-item service level constraints Item order lead times Constant (planned) Unsatisfied end-item demand is backordered 5

6 Network transformation principles 6

7 Basic principle: Divergent node structures At time 0 item(s) with longest cumulative lead time has (have) been ordered At time t items m C 0 have been ordered, of which some at time t At time t+l another set of items C is ordered End-items in P(C) use items in C ( ) PC ( ) PC ( ) 0 c L Immediate excess coverage P( C 0 ) C 0 C 0 ( ) P( C ) P C ( ) P C 0 0 C 0 C C 0 C P(C P(C ) = ) P(C P(C 0 ) 0 ) 7

8 Product structure requirements for SBSpolicies Acyclic BOM If two components i and i 2 are common to items j and j 2 then a a ij 2 i j aij = = c a i j ii 2 Condition allows for elimination of i 2 with respect to (j,j 2 ) a D + a D ij j ij j 2 2 a D + a D i j j i j j = c ii 2 mutually non-zero elements of rows i and i 2 are dependent 8

9 Elimination of different consumption ratios Create versions of items, whereby sum of consumption rates of versions satisfy the original product structure constraint P={3} 3 P={3} P={3,4} C={,2,5} C={3} C={2} P={4} 4 P={4} C={,2} C={4} 9

10 Elimination of cycles in BOM If item i is used at n different levels in the BOM structure, then create n different 3 versions of item i. P={} P={4} P={4} P={4} P={4} 5 P={4} C={9}={7} C={8}={6} C={2,3} 3 3 C={2,3,7} 2 C={2,3,6} P={5} C={4} P={5} C={,2} C={5} 3 3 0

11 Natural hierarchy BOM and lead time structures determine a natural hierarchy defined by divergent structurs Divergent structures relate sets of end-items to sets of items Divergent structures determine forecasting requirements An order of a set of items associated with a node can be released based on the forecast of demand for the set of end-items associated with this node over its cumulative lead time (plus one period) Safety stock requirements over the cumulative lead time

12 Dead stocks average stock in divergent structures [ ] ˆ PC ( ) dead i C i E X = E X + E X { Ci C} Due to synchronization mechanism not all combinations of average stocks of items are possible If average inventories are taken from real life data and SBS policies are used, dead stock may result net stock of an item cannot get below some positive number Computation of performance for given average stock levels requires determination of dead stocks there is no unique solution to this minimize value of dead stock 2

13 Divergent structures and allocation rules 3

14 ME divergent networks W i s u p p li 0 U i i k E i c u st o m er e r s s 4

15 Definitions h i p k U i W i E i Y i X i I ki added value created at (holding item i cost) at item i penalty cost per unit short of end-item k set of predecessors of i, i.e. all child items set of all items in the echelon of i set of all end-items in the echelon of i echelon inventory position of i net stock of item i net stock of end-item k in a system with root node i 5

16 Base stock policies and linear allocation rules Shortfall of i ( ] ( ) Y () t = S q Z ( t L) + D t L, t j j j i i i i j i ( ) S = S + q i j j j j V End-item demand during lead time of i + Parameter that determines stock kept at i on behalf of j j-specific eliminates the effect of zero-cost successors on non-zero-cost successors: no stock of item i ( ) [ ] = ( 0, ] E X q E Z D L + i j j i i i j Vi 6

17 Allocation fractions Eppen and Schrage (98): equal stockout probability q j σ j = σ m V i m Van der Heijden et al (997): minimal imbalance σ µ 2 2 j j q j = σ m 2 µ m j Vi j Vi 7

18 Optimization problem Given allocation fractions, determine base stock levels for all items Expressions for costs and performance are similar to expressions for serial systems Imbalance should be sufficiently low to ensure validity of analytical results in case imbalance too high, base stock levels should be such that average stocks are increased at upstream stages with high imbalance requires experimental research 8

19 Optimality conditions on base stock levels 9

20 Optimal order-up-to-policies for divergent MEIS: Generalized Newsvendor equations theorem for periodic review systems without setup costs Optimal order-up-to-policies and allocation functions under balance assumption satisfy generalized Newsvendor equations: Probability of non-stockout at downstream stockpoint k in the divergent subsystem with stockpoint i as most upstream stage equals or equivalently { 0} P I ki = p k + j U k k j j U i h p + h + h j p + h = p + h + h P I j Ui j Uk k j k k j ki k { 0} 20

21 Sample path condition on control policies Assume that each node i is controlled according to a policy, which is defined by a parameter ξ i and possibly other parameters and functions Sample path condition If ξ j is increased by E j ε for all j W j, and there is upstream availability then a one-time additional flow of ε units is created towards each k E j 2

22 Main theorem If sample path condition holds then i I p + h = p + h + h P I k j k k j ki k E j U \ W k E j U i k i i k { 0} Proof based on echelon costs and similar to proof of single-item Newsvendor equation Relationship can be used as heuristic 22

23 Computational aspects 23

24 Finite horizon ruin probabilities Non-stockout probabilities P{I ki 0} can be written as finite horizon ruin probabilities j P{ Iki 0 } = P Xk ξ j, j =,..., i k = Expressions can be accurately approximated recursively j Gi( ξ,..., ξi) = P Xk ξ j, j =,..., i k = 24

25 Recursive computations Define random variables Y i { } = { } PY x P X x i Gi( ξ,..., ξi, x) =, Gi ( ξ,..., ξi ) 0, =2,...,. { } P Y x x i N Theorem { } { i + i, i ξi } PY { ξ } P X Y x Y PY x =, i= 2,..., N i i i 25

26 Two-moment recursion Fit mixture of Erlang distributions on E[Y i ] and σ 2 (Y i ) [ i] [ i] [ i i ξi ] ( ) ( ) ( ) EY = E X + EY Y, i= 2,..., N- σ Y = σ X + σ Y Y ξ, i = 2,..., N - i i i i i Compute P{I ki 0} = G i (ξ,..., ξ ) ι recursively G( ξ) = P{ Y ξ} G ( ξ,..., ξ ) = P{ Y ξ } G ( ξ,..., ξ ), i = 2,..., N i i i i i i 26

27 Computational efficiency Solving cost-balance equations or generalized Newsvendor equations is equivalent to solving a one-dimensional bisection scheme to find a target value in (0,) Finding base stock levels for N-item divergent structure requires solution of N bisection schemes As efficient as solving N single-echelon systems Finding base stock levels for N-item general systems requires solution of at most 3N bisection schemes 27

28 Comparison against LP rolling schedule 28

29 LP-based rolling schedule approach Define objective function linear holding and penalty costs incorporate safety stocks by goal programming formulation Feasibility constraints inventory balance equations at item level material availability constraints w.r.t. dependent requirements No lot sizing constraints Solve LP problem Implement immediate decisions 29

30 Experimental comparison of SCOP paradigms: LP vs. SBS for uncapacitated systems Example product structure 4 final products E[D] = 00 P = 95% added value $0 assembly lead time 4 specific items value $0 2 semi-specific items value $30 common item value $50 30

31 An -item system L s 5 L sc 9 L f (c,α ) L s 6 L f 2 (c 2,α 2 ) L c L s 7 L sc 0 L f 3 (c 3,α 3 ) L s 8 L f 4 (c 4,α 4 ) 3

32 Results simulation study 2 Supply Chain Inventory capital Customer service c i (L f, L s,l sc,l c ) SBS ana SBS sim LP % P,SBS P,LP 0.25 (,,2,4) % 95% 95% 0.25 (,4,2,) % 95% 95% 0.25 (,,4,2) % 95% 95% 0.5 (,,2,4) % 95% 95% 0.5 (,4,2,) % 95% 95% 0.5 (,,4,2) % 95% 95% (,,2,4) % 95% 95% (,4,2,) % 95% 95% (,,4,2) % 95% 95% 2 (,,2,4) % 94% 95% 2 (,4,2,) % 95% 95% 2 (,,4,2) % 94% 95% 32

33 Impact of final product cost on safety stocks, LP (c 2, c 2 2, c 3 2, c 4 2 ) (h, h 2, h 3, h 4 ) (L f,l s,l sc,l c ) LP-based safety stocks P P 2 ss ss 2 ss 3 ss 4 (0.25, 0.25, 0.25, 0.25) (20, 20, 20, 0) (,,4,2) (0.25, 0.25, 0.25, 0.25) (20, 20, 20, 0) (,,2,4) (0.50, 0.50, 0.50, 0.50) (20, 20, 20, 0) (,,2,4) (0.50, 0.50, 0.50, 0.50) (20, 20, 20, 0) (,,4,2) (,,, ) (20, 20, 20, 0) (,,2,4) (,,, ) (20, 20, 20, 0) (,,4,2) (2, 2, 2, 2) (20, 20, 20, 0) (,,2,4) (2, 2, 2, 2) (20, 20, 0, 0) (,,4,2) (0.25, 0.25, 0.50, 0.50) (0, 0, 0, 0) (,,4,2)

34 Impact of final product cost on safety stocks, SBS (c 2, c 2 2, c 3 2, c 4 2 ) (h, h 2, h 3, h 4 ) (L f,l s,l sc,l c ) SBS-based safety stocks P P 2 ss Ss 2 ss 3 ss 4 (0.25, 0.25, 0.25, 0.25) (20, 20, 20, 0) (,,4,2) (0.25, 0.25, 0.25, 0.25) (20, 20, 20, 0) (,,2,4) (0.50, 0.50, 0.50, 0.50) (20, 20, 20, 0) (,,2,4) (0.50, 0.50, 0.50, 0.50) (20, 20, 20, 0) (,,4,2) (,,, ) (20, 20, 20, 0) (,,2,4) (,,, ) (20, 20, 20, 0) (,,4,2) (2, 2, 2, 2) (20, 20, 20, 0) (,,2,4) (2, 2, 2, 2) (20, 20, 0, 0) (,,4,2) (0.25, 0.25, 0.50, 0.50) (0, 0, 0, 0) (,,4,2)

35 Managerial insights Most of inventory capital tied up downstream informally speaking optimal upstream stock never exceeds 0% of overall supply chain inventory capital Flexibility only pays off downstream lead time reduction replenishment frequency Focus on efficiency upstream Postponement should come almost for free Shortages and overages should be shared among parent items Moving CODP upstream yields substantial savings Control parameters, i.e. safety stocks, depend on the operational planning logic applied 35

36 Inventory capital allocation CW lead time inventory investm ent 00% 90% 80% 70% 60% 50% 40% 30% 20% 0% 0% 0% 20% 30% 40% 50% 60% 70% 80% 90% 00% value at CW central local 36

37 Inventory capital allocation CW lead time 5 inventory investment 00% 90% 80% 70% 60% 50% 40% 30% 20% 0% 0% 0% 20% 30% 40% 50% 60% 70% 80% 90% 00% value at CW central local 37

38 Inventory capital allocation CW lead time 0 inventory investment 00% 90% 80% 70% 60% 50% 40% 30% 20% 0% 0% 0% 20% 30% 40% 50% 60% 70% 80% 90% 00% value at CW central local 38

39 Lot sizing and allocation rules R i review period of item i Ri rj =, rj, j Vi R V it θ jt j set of successors of item i that order at time t the last ordering moment of successor j at or before time t Allocation policy should be such that net stock of i non-negative Linear allocation policies defined before in general do not satisfy this constraint 39

40 Lot sizing and allocation rules ( ) X () t = q W() t i j j i j V i + Then we can recursively define the random variables W(t), ( ] W( L) = Z (0) + D 0, L i i i i Wt () = Zi( t Li) + Di( t Lt i, ] Dj( θ jt, t qw j ( θ jt) q j j Vit j Vit j V it The allocation policy follows from ( ) Y () t = S q W() t, if j orders at time t j j j j + 40

41 Conclusions SBS policies allow for control of general multiitem multi-echelon networks Close-to-optimal SBS policies can be efficiently computed SBS policies perform well compared with standard rolling scheduling SBS policies can be implemented for operational control requires forecast-dependent base stock levels 4

42 Challenges for further research Extension of SBS with (period order) lot sizing Incorporation of multi-level forecasting future detailed forecasts are conditional on today's aggregate forecast multi-level forecast structure should be aligned with value network structure Finite capacity Should finite capacity be modelled explicitly through constraints or explicitly through model of transformation process or implicitly through planned lead time, i.e. loosely coupling SCOPlevel and shopfloor control level Empirical research to test various modelling options 42

Newsvendor Bounds and Heuristics for Series Systems

Newsvendor Bounds and Heuristics for Series Systems Newsvendor Bounds and Heuristics for Series Systems August 6, 2002 Abstract We propose a heuristic for a multi-stage serial system based on solving a single newsvendor problem per stage. The heuristic

More information

Column Generation II : Application in Distribution Network Design

Column Generation II : Application in Distribution Network Design Column Generation II : Application in Distribution Network Design Teo Chung-Piaw (NUS) 27 Feb 2003, Singapore 1 Supply Chain Challenges 1.1 Introduction Network of facilities: procurement of materials,

More information

Rollout Algorithms for Discrete Optimization: A Survey

Rollout Algorithms for Discrete Optimization: A Survey Rollout Algorithms for Discrete Optimization: A Survey by Dimitri P. Bertsekas Massachusetts Institute of Technology Cambridge, MA 02139 dimitrib@mit.edu August 2010 Abstract This chapter discusses rollout

More information

Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)

Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers) Heuristics for multi-item two-echelon spare parts inventory control problem with batch ordering in the central warehouse Topan, E.; Bayindir, Z.P.; Tan, T. Published: 01/01/2010 Document Version Publisher

More information

Mixture Models and the EM Algorithm

Mixture Models and the EM Algorithm Mixture Models and the EM Algorithm Padhraic Smyth, Department of Computer Science University of California, Irvine c 2017 1 Finite Mixture Models Say we have a data set D = {x 1,..., x N } where x i is

More information

10. Network dimensioning

10. Network dimensioning Partly based on slide material by Samuli Aalto and Jorma Virtamo ELEC-C7210 Modeling and analysis of communication networks 1 Contents Introduction Parameters: topology, routing and traffic Dimensioning

More information

FLEXIBLE ASSEMBLY SYSTEMS

FLEXIBLE ASSEMBLY SYSTEMS FLEXIBLE ASSEMBLY SYSTEMS Job Shop and Flexible Assembly Job Shop Each job has an unique identity Make to order, low volume environment Possibly complicated route through system Very difficult Flexible

More information

Computing intersections in a set of line segments: the Bentley-Ottmann algorithm

Computing intersections in a set of line segments: the Bentley-Ottmann algorithm Computing intersections in a set of line segments: the Bentley-Ottmann algorithm Michiel Smid October 14, 2003 1 Introduction In these notes, we introduce a powerful technique for solving geometric problems.

More information

Note Set 4: Finite Mixture Models and the EM Algorithm

Note Set 4: Finite Mixture Models and the EM Algorithm Note Set 4: Finite Mixture Models and the EM Algorithm Padhraic Smyth, Department of Computer Science University of California, Irvine Finite Mixture Models A finite mixture model with K components, for

More information

Cost Optimization in the (S 1, S) Lost Sales Inventory Model with Multiple Demand Classes

Cost Optimization in the (S 1, S) Lost Sales Inventory Model with Multiple Demand Classes Cost Optimization in the (S 1, S) Lost Sales Inventory Model with Multiple Demand Classes A.A. Kranenburg, G.J. van Houtum Department of Technology Management, Technische Universiteit Eindhoven, Eindhoven,

More information

AM 221: Advanced Optimization Spring 2016

AM 221: Advanced Optimization Spring 2016 AM 221: Advanced Optimization Spring 2016 Prof Yaron Singer Lecture 3 February 1st 1 Overview In our previous lecture we presented fundamental results from convex analysis and in particular the separating

More information

Tuesday, April 10. The Network Simplex Method for Solving the Minimum Cost Flow Problem

Tuesday, April 10. The Network Simplex Method for Solving the Minimum Cost Flow Problem . Tuesday, April The Network Simplex Method for Solving the Minimum Cost Flow Problem Quotes of the day I think that I shall never see A poem lovely as a tree. -- Joyce Kilmer Knowing trees, I understand

More information

BUFFER STOCKS IN KANBAN CONTROLLED (TRADITIONAL) UNSATURATED MULTI-STAGE PRODUCTION SYSTEM

BUFFER STOCKS IN KANBAN CONTROLLED (TRADITIONAL) UNSATURATED MULTI-STAGE PRODUCTION SYSTEM VOL. 3, NO., FEBRUARY 008 ISSN 89-6608 006-008 Asian Research Publishing Network (ARPN). All rights reserved. BUFFER STOCKS IN KANBAN CONTROLLED (TRADITIONAL) UNSATURATED MULTI-STAGE PRODUCTION SYSTEM

More information

4 LINEAR PROGRAMMING (LP) E. Amaldi Fondamenti di R.O. Politecnico di Milano 1

4 LINEAR PROGRAMMING (LP) E. Amaldi Fondamenti di R.O. Politecnico di Milano 1 4 LINEAR PROGRAMMING (LP) E. Amaldi Fondamenti di R.O. Politecnico di Milano 1 Mathematical programming (optimization) problem: min f (x) s.t. x X R n set of feasible solutions with linear objective function

More information

Maximizing Network Utilization with Max-Min Fairness in Wireless Sensor Networks

Maximizing Network Utilization with Max-Min Fairness in Wireless Sensor Networks 1 Maximizing Network Utilization with Max-Min Fairness in Wireless Sensor Networks Avinash Sridharan and Bhaskar Krishnamachari {asridhar,bkrishna}@usc.edu Department of Electrical Engineering University

More information

Cost Optimization in the (S 1, S) Lost Sales Inventory Model with Multiple Demand Classes

Cost Optimization in the (S 1, S) Lost Sales Inventory Model with Multiple Demand Classes Cost Optimization in the (S 1, S) Lost Sales Inventory Model with Multiple Demand Classes A.A. Kranenburg, G.J. van Houtum Department of Technology Management, Technische Universiteit Eindhoven, Eindhoven,

More information

The Multi-Item Capacitated Lot-Sizing Problem With Safety Stocks and Demand Shortage Costs

The Multi-Item Capacitated Lot-Sizing Problem With Safety Stocks and Demand Shortage Costs The Multi-Item Capacitated Lot-Sizing Problem With Safety Stocks and Demand Shortage Costs Nabil Absi 1, Safia Kedad-Sidhoum 2 1 Ecole des Mines de Saint-Etienne, CMP - Site Georges Charpak - 880 Avenue

More information

Job Management & Network Flows

Job Management & Network Flows Job Management & Network Flows How to solve job management problems with networks flow techniques Francesca Armandillo BMI paper Vre Universiteit Faculty of Sciences Business Mathematics and Informatics

More information

Foundations of Computer Science Spring Mathematical Preliminaries

Foundations of Computer Science Spring Mathematical Preliminaries Foundations of Computer Science Spring 2017 Equivalence Relation, Recursive Definition, and Mathematical Induction Mathematical Preliminaries Mohammad Ashiqur Rahman Department of Computer Science College

More information

Outline. Petri nets. Introduction Examples Properties Analysis techniques. 1 EE249Fall04

Outline. Petri nets. Introduction Examples Properties Analysis techniques. 1 EE249Fall04 Outline Petri nets Introduction Examples Properties Analysis techniques 1 Petri Nets (PNs) Model introduced by C.A. Petri in 1962 Ph.D. Thesis: Communication with Automata Applications: distributed computing,

More information

A Heuristic Approach to the Design of Kanban Systems

A Heuristic Approach to the Design of Kanban Systems A Heuristic Approach to the Design of Kanban Systems Chuda Basnet Department of Management Systems University of Waikato, Hamilton Abstract Kanbans are often used to communicate replenishment requirements

More information

MAXIMAL FLOW THROUGH A NETWORK

MAXIMAL FLOW THROUGH A NETWORK MAXIMAL FLOW THROUGH A NETWORK L. R. FORD, JR. AND D. R. FULKERSON Introduction. The problem discussed in this paper was formulated by T. Harris as follows: "Consider a rail network connecting two cities

More information

UNIT 4 Branch and Bound

UNIT 4 Branch and Bound UNIT 4 Branch and Bound General method: Branch and Bound is another method to systematically search a solution space. Just like backtracking, we will use bounding functions to avoid generating subtrees

More information

d) Consider the single item inventory model. D = 20000/year, C0 = Rs 1000/order, C = Rs 200/item, i = 20%, TC at EOQ is a) 1000 b) 15000

d) Consider the single item inventory model. D = 20000/year, C0 = Rs 1000/order, C = Rs 200/item, i = 20%, TC at EOQ is a) 1000 b) 15000 Assignment 12 1. A model of the form Y = a + bt was fit with data for four periods. The two equations are 90 = 4a + 10b and 240 = 10a + 30b. The forecast for the fifth period is a) 120 b) 30 c) 160 d)

More information

15-780: MarkovDecisionProcesses

15-780: MarkovDecisionProcesses 15-780: MarkovDecisionProcesses J. Zico Kolter Feburary 29, 2016 1 Outline Introduction Formal definition Value iteration Policy iteration Linear programming for MDPs 2 1988 Judea Pearl publishes Probabilistic

More information

Distributed non-convex optimization

Distributed non-convex optimization Distributed non-convex optimization Behrouz Touri Assistant Professor Department of Electrical and Computer Engineering University of California San Diego/University of Colorado Boulder AFOSR Computational

More information

Mixed-Integer Optimization for the Combined capacitated Facility Location-Routing Problem

Mixed-Integer Optimization for the Combined capacitated Facility Location-Routing Problem Mixed-Integer Optimization for the Combined capacitated Facility Location-Routing Problem Dimitri Papadimitriou 1, Didier Colle 2, Piet Demeester 2 dimitri.papadimitriou@nokia.com, didier.colle@ugent.be,

More information

IJSER

IJSER International Journal of Scientific & Engineering Research, Volume 4, Issue 1, October-213 1399 EFFECT OF KANBANS IN THE PERFORMANCE OF, AND E G.G Sastry 1, Dr. Mukesh Saxena 2, Dr. Rajnish Garg 3 Abstract

More information

Shortest-route formulation of mixed-model assembly line balancing problem

Shortest-route formulation of mixed-model assembly line balancing problem European Journal of Operational Research 116 (1999) 194±204 Theory and Methodology Shortest-route formulation of mixed-model assembly line balancing problem Erdal Erel a, *, Hadi Gokcen b a Faculty of

More information

Balancing the Migration of Virtual Network Functions with Replications in Data Centers

Balancing the Migration of Virtual Network Functions with Replications in Data Centers Balancing the Migration of Virtual Network Functions with Replications in Data Centers Francisco Carpio, Admela Jukan and Rastin Pries Technische Universität Braunschweig, Germany, Nokia Bell Labs, Germany

More information

Summary of Raptor Codes

Summary of Raptor Codes Summary of Raptor Codes Tracey Ho October 29, 2003 1 Introduction This summary gives an overview of Raptor Codes, the latest class of codes proposed for reliable multicast in the Digital Fountain model.

More information

Integer Programming ISE 418. Lecture 7. Dr. Ted Ralphs

Integer Programming ISE 418. Lecture 7. Dr. Ted Ralphs Integer Programming ISE 418 Lecture 7 Dr. Ted Ralphs ISE 418 Lecture 7 1 Reading for This Lecture Nemhauser and Wolsey Sections II.3.1, II.3.6, II.4.1, II.4.2, II.5.4 Wolsey Chapter 7 CCZ Chapter 1 Constraint

More information

11. APPROXIMATION ALGORITHMS

11. APPROXIMATION ALGORITHMS 11. APPROXIMATION ALGORITHMS load balancing center selection pricing method: vertex cover LP rounding: vertex cover generalized load balancing knapsack problem Lecture slides by Kevin Wayne Copyright 2005

More information

Probabilistic Modeling of Leach Protocol and Computing Sensor Energy Consumption Rate in Sensor Networks

Probabilistic Modeling of Leach Protocol and Computing Sensor Energy Consumption Rate in Sensor Networks Probabilistic Modeling of Leach Protocol and Computing Sensor Energy Consumption Rate in Sensor Networks Dezhen Song CS Department, Texas A&M University Technical Report: TR 2005-2-2 Email: dzsong@cs.tamu.edu

More information

Approximate Linear Programming for Average-Cost Dynamic Programming

Approximate Linear Programming for Average-Cost Dynamic Programming Approximate Linear Programming for Average-Cost Dynamic Programming Daniela Pucci de Farias IBM Almaden Research Center 65 Harry Road, San Jose, CA 51 pucci@mitedu Benjamin Van Roy Department of Management

More information

CS839: Probabilistic Graphical Models. Lecture 10: Learning with Partially Observed Data. Theo Rekatsinas

CS839: Probabilistic Graphical Models. Lecture 10: Learning with Partially Observed Data. Theo Rekatsinas CS839: Probabilistic Graphical Models Lecture 10: Learning with Partially Observed Data Theo Rekatsinas 1 Partially Observed GMs Speech recognition 2 Partially Observed GMs Evolution 3 Partially Observed

More information

1. What do you get as the integer and noninteger parts if you factor this as we did with our cutting planes:

1. What do you get as the integer and noninteger parts if you factor this as we did with our cutting planes: 1. What do you get as the integer and noninteger parts if you factor this as we did with our cutting planes: x 1 = -8.75-3.1 x 2 + 4.2 x 3 + 7 x 5-8 x 6 2. What constraint is added? 3. What do you need

More information

Multithreaded Algorithms Part 1. Dept. of Computer Science & Eng University of Moratuwa

Multithreaded Algorithms Part 1. Dept. of Computer Science & Eng University of Moratuwa CS4460 Advanced d Algorithms Batch 08, L4S2 Lecture 11 Multithreaded Algorithms Part 1 N. H. N. D. de Silva Dept. of Computer Science & Eng University of Moratuwa Announcements Last topic discussed is

More information

Job-shop scheduling with limited capacity buffers

Job-shop scheduling with limited capacity buffers Job-shop scheduling with limited capacity buffers Peter Brucker, Silvia Heitmann University of Osnabrück, Department of Mathematics/Informatics Albrechtstr. 28, D-49069 Osnabrück, Germany {peter,sheitman}@mathematik.uni-osnabrueck.de

More information

11. APPROXIMATION ALGORITHMS

11. APPROXIMATION ALGORITHMS 11. APPROXIMATION ALGORITHMS load balancing center selection pricing method: vertex cover LP rounding: vertex cover generalized load balancing knapsack problem Lecture slides by Kevin Wayne Copyright 2005

More information

2.3.4 Optimal paths in directed acyclic graphs

2.3.4 Optimal paths in directed acyclic graphs .3.4 Optimal paths in directed acyclic graphs Definition: A directed graph G = (N, A) is acyclic if it contains no circuits. A directed acyclic graph is referred to as DAG. circuit Problem Given a directed

More information

Approximation Algorithms: The Primal-Dual Method. My T. Thai

Approximation Algorithms: The Primal-Dual Method. My T. Thai Approximation Algorithms: The Primal-Dual Method My T. Thai 1 Overview of the Primal-Dual Method Consider the following primal program, called P: min st n c j x j j=1 n a ij x j b i j=1 x j 0 Then the

More information

Design of Flexible Assembly Line to Minimize Equipment Cost

Design of Flexible Assembly Line to Minimize Equipment Cost Design of Flexible Assembly Line to Minimize Equipment Cost Joseph Bukchin Department of Industrial Engineering Faculty of Engineering, Tel-Aviv University, Tel-Aviv 69978 ISRAEL Tel: 972-3-640794; Fax:

More information

Quality-Assured Cloud Bandwidth Auto-Scaling for Video-on-Demand Applications

Quality-Assured Cloud Bandwidth Auto-Scaling for Video-on-Demand Applications Quality-Assured Cloud Bandwidth Auto-Scaling for Video-on-Demand Applications Di Niu, Hong Xu, Baochun Li University of Toronto Shuqiao Zhao UUSee, Inc., Beijing, China 1 Applications in the Cloud WWW

More information

Optimization. Using Analytic Solver Platform REVIEW BASED ON MANAGEMENT SCIENCE

Optimization. Using Analytic Solver Platform REVIEW BASED ON MANAGEMENT SCIENCE Optimization Using Analytic Solver Platform REVIEW BASED ON MANAGEMENT SCIENCE What We ll Cover Today Introduction Frontline Systems Session Ι beta training program goals Overview of Analytic Solver Platform

More information

Measurement Based Routing Strategies on Overlay Architectures

Measurement Based Routing Strategies on Overlay Architectures Measurement Based Routing Strategies on Overlay Architectures Student: Tuna Güven Faculty: Bobby Bhattacharjee, Richard J. La, and Mark A. Shayman LTS Review February 15 th, 2005 Outline Measurement-Based

More information

7.3 Spanning trees Spanning trees [ ] 61

7.3 Spanning trees Spanning trees [ ] 61 7.3. Spanning trees [161211-1348 ] 61 7.3 Spanning trees We know that trees are connected graphs with the minimal number of edges. Hence trees become very useful in applications where our goal is to connect

More information

A Path Decomposition Approach for Computing Blocking Probabilities in Wavelength-Routing Networks

A Path Decomposition Approach for Computing Blocking Probabilities in Wavelength-Routing Networks IEEE/ACM TRANSACTIONS ON NETWORKING, VOL. 8, NO. 6, DECEMBER 2000 747 A Path Decomposition Approach for Computing Blocking Probabilities in Wavelength-Routing Networks Yuhong Zhu, George N. Rouskas, Member,

More information

Lower-Bounded Facility Location

Lower-Bounded Facility Location Lower-Bounded Facility Location Zoya Svitkina Abstract We study the lower-bounded facility location problem, which generalizes the classical uncapacitated facility location problem in that it comes with

More information

1 The range query problem

1 The range query problem CS268: Geometric Algorithms Handout #12 Design and Analysis Original Handout #12 Stanford University Thursday, 19 May 1994 Original Lecture #12: Thursday, May 19, 1994 Topics: Range Searching with Partition

More information

Robust time-varying shortest path with arbitrary waiting time at vertices

Robust time-varying shortest path with arbitrary waiting time at vertices Croatian Operational Research Review 525 CRORR 8(2017), 525 56 Robust time-varying shortest path with arbitrary waiting time at vertices Gholamhassan Shirdel 1, and Hassan Rezapour 1 1 Department of Mathematics,

More information

EECS 571 Principles of Real-Time Embedded Systems. Lecture Note #8: Task Assignment and Scheduling on Multiprocessor Systems

EECS 571 Principles of Real-Time Embedded Systems. Lecture Note #8: Task Assignment and Scheduling on Multiprocessor Systems EECS 571 Principles of Real-Time Embedded Systems Lecture Note #8: Task Assignment and Scheduling on Multiprocessor Systems Kang G. Shin EECS Department University of Michigan What Have We Done So Far?

More information

Department of Computer Science and Engineering Analysis and Design of Algorithm (CS-4004) Subject Notes

Department of Computer Science and Engineering Analysis and Design of Algorithm (CS-4004) Subject Notes Page no: Department of Computer Science and Engineering Analysis and Design of Algorithm (CS-00) Subject Notes Unit- Greedy Technique. Introduction: Greedy is the most straight forward design technique.

More information

Lecture notes on Transportation and Assignment Problem (BBE (H) QTM paper of Delhi University)

Lecture notes on Transportation and Assignment Problem (BBE (H) QTM paper of Delhi University) Transportation and Assignment Problems The transportation model is a special class of linear programs. It received this name because many of its applications involve determining how to optimally transport

More information

1 Proximity via Graph Spanners

1 Proximity via Graph Spanners CS273: Algorithms for Structure Handout # 11 and Motion in Biology Stanford University Tuesday, 4 May 2003 Lecture #11: 4 May 2004 Topics: Proximity via Graph Spanners Geometric Models of Molecules, I

More information

Effect of memory latency

Effect of memory latency CACHE AWARENESS Effect of memory latency Consider a processor operating at 1 GHz (1 ns clock) connected to a DRAM with a latency of 100 ns. Assume that the processor has two ALU units and it is capable

More information

Petri Nets ee249 Fall 2000

Petri Nets ee249 Fall 2000 Petri Nets ee249 Fall 2000 Marco Sgroi Most slides borrowed from Luciano Lavagno s lecture ee249 (1998) 1 Models Of Computation for reactive systems Main MOCs: Communicating Finite State Machines Dataflow

More information

Approximating the Virtual Network Embedding Problem: Theory and Practice

Approximating the Virtual Network Embedding Problem: Theory and Practice Approximating the Virtual Network Embedding Problem: Theory and Practice 2 5 AC 3 B 2 2 2 D 2 0 3 23rd International Symposium on Mathematical Programming 208 Bordeaux, France Matthias Rost Technische

More information

Load Balancing for Problems with Good Bisectors, and Applications in Finite Element Simulations

Load Balancing for Problems with Good Bisectors, and Applications in Finite Element Simulations Load Balancing for Problems with Good Bisectors, and Applications in Finite Element Simulations Stefan Bischof, Ralf Ebner, and Thomas Erlebach Institut für Informatik Technische Universität München D-80290

More information

Advanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs

Advanced Operations Research Techniques IE316. Quiz 1 Review. Dr. Ted Ralphs Advanced Operations Research Techniques IE316 Quiz 1 Review Dr. Ted Ralphs IE316 Quiz 1 Review 1 Reading for The Quiz Material covered in detail in lecture. 1.1, 1.4, 2.1-2.6, 3.1-3.3, 3.5 Background material

More information

Advanced Synchronization and Deadlock

Advanced Synchronization and Deadlock Advanced Synchronization and Deadlock A house of cards? Locks + CV/signal a great way to regulate access to a single shared object......but general multi-threaded programs touch multiple shared objects

More information

Lecture 009 (November 25, 2009) Suppose that I is a small category, and let s(r Mod) I be the category of I-diagrams in simplicial modules.

Lecture 009 (November 25, 2009) Suppose that I is a small category, and let s(r Mod) I be the category of I-diagrams in simplicial modules. Lecture 009 (November 25, 2009) 20 Homotopy colimits Suppose that I is a small category, and let s(r Mod) I be the category of I-diagrams in simplicial modules. The objects of this category are the functors

More information

Machine Learning

Machine Learning Machine Learning 10-601 Tom M. Mitchell Machine Learning Department Carnegie Mellon University March 4, 2015 Today: Graphical models Bayes Nets: EM Mixture of Gaussian clustering Learning Bayes Net structure

More information

Computer Vision Group Prof. Daniel Cremers. 4. Probabilistic Graphical Models Directed Models

Computer Vision Group Prof. Daniel Cremers. 4. Probabilistic Graphical Models Directed Models Prof. Daniel Cremers 4. Probabilistic Graphical Models Directed Models The Bayes Filter (Rep.) (Bayes) (Markov) (Tot. prob.) (Markov) (Markov) 2 Graphical Representation (Rep.) We can describe the overall

More information

Generalized Network Flow Programming

Generalized Network Flow Programming Appendix C Page Generalized Network Flow Programming This chapter adapts the bounded variable primal simplex method to the generalized minimum cost flow problem. Generalized networks are far more useful

More information

Analysis of Algorithms

Analysis of Algorithms Lab Manual Analysis of Algorithms S.E. Computer Sem IV Index Table S r. N o. Tit le of Pr ogr amm ing Assignm ents 1. Recursive Binar y Sear ch 2. Recursive Mer ge Sort 3. Heap Sort 4. Random ized Q uick

More information

Distributed Resource Allocation for Stream Data Processing

Distributed Resource Allocation for Stream Data Processing Distributed Resource Allocation for Stream Data Processing Ao Tang,ZhenLiu,CathyXia,andLiZhang California Institute of Technology, Department of Electrical Engineering, Pasadena, CA 95. IBM T.J. Watson

More information

11. APPROXIMATION ALGORITHMS

11. APPROXIMATION ALGORITHMS Coping with NP-completeness 11. APPROXIMATION ALGORITHMS load balancing center selection pricing method: weighted vertex cover LP rounding: weighted vertex cover generalized load balancing knapsack problem

More information

Optimal Control of a Production-Inventory System with both Backorders and Lost Sales

Optimal Control of a Production-Inventory System with both Backorders and Lost Sales Optimal Control of a Production-Inventory System with both Backorders and Lost Sales Saif Benjaafar Mohsen ElHafsi 2 Tingliang Huang 3 Industrial & Systems Engineering, Department of Mechanical Engineering,

More information

Evolutionary tree reconstruction (Chapter 10)

Evolutionary tree reconstruction (Chapter 10) Evolutionary tree reconstruction (Chapter 10) Early Evolutionary Studies Anatomical features were the dominant criteria used to derive evolutionary relationships between species since Darwin till early

More information

Apprenticeship Learning for Reinforcement Learning. with application to RC helicopter flight Ritwik Anand, Nick Haliday, Audrey Huang

Apprenticeship Learning for Reinforcement Learning. with application to RC helicopter flight Ritwik Anand, Nick Haliday, Audrey Huang Apprenticeship Learning for Reinforcement Learning with application to RC helicopter flight Ritwik Anand, Nick Haliday, Audrey Huang Table of Contents Introduction Theory Autonomous helicopter control

More information

(2,4) Trees Goodrich, Tamassia. (2,4) Trees 1

(2,4) Trees Goodrich, Tamassia. (2,4) Trees 1 (2,4) Trees 9 2 5 7 10 14 (2,4) Trees 1 Multi-Way Search Tree ( 9.4.1) A multi-way search tree is an ordered tree such that Each internal node has at least two children and stores d 1 key-element items

More information

Introduction to Complex Analysis

Introduction to Complex Analysis Introduction to Complex Analysis George Voutsadakis 1 1 Mathematics and Computer Science Lake Superior State University LSSU Math 413 George Voutsadakis (LSSU) Complex Analysis October 2014 1 / 50 Outline

More information

PRODUCTION SYSTEMS ENGINEERING:

PRODUCTION SYSTEMS ENGINEERING: PRODUCTION SYSTEMS ENGINEERING: Optimality through Improvability Semyon M. Meerkov EECS Department University of Michigan Ann Arbor, MI 48109-2122, USA YEQT-IV (Young European Queueing Theorists Workshop)

More information

ONLINE SUPPLEMENT TO: OPTIMIZING STRATEGIC SAFETY STOCK PLACEMENT IN GENERAL ACYCLIC NETWORKS

ONLINE SUPPLEMENT TO: OPTIMIZING STRATEGIC SAFETY STOCK PLACEMENT IN GENERAL ACYCLIC NETWORKS ONLINE SUPPLEMENT TO: OPTIMIZING STRATEGIC SAFETY STOCK PLACEMENT IN GENERAL ACYCLIC NETWORKS SALAL HUMAIR School of Science and Engineering, Lahore University of Management Sciences Lahore 54792, Pakistan,

More information

Reinforcement Learning: A brief introduction. Mihaela van der Schaar

Reinforcement Learning: A brief introduction. Mihaela van der Schaar Reinforcement Learning: A brief introduction Mihaela van der Schaar Outline Optimal Decisions & Optimal Forecasts Markov Decision Processes (MDPs) States, actions, rewards and value functions Dynamic Programming

More information

Cumulative Reputation Systems for Peer-to-Peer Content Distribution

Cumulative Reputation Systems for Peer-to-Peer Content Distribution Cumulative Reputation Systems for Peer-to-Peer Content Distribution B. Mortazavi and G. Kesidis Pennsylvania State University also with Verizon Wireless CISS Princeton March 23, 2006 1 Outline P2P CDNs

More information

AN APPROXIMATE INVENTORY MODEL BASED ON DIMENSIONAL ANALYSIS. Victoria University, Wellington, New Zealand

AN APPROXIMATE INVENTORY MODEL BASED ON DIMENSIONAL ANALYSIS. Victoria University, Wellington, New Zealand AN APPROXIMATE INVENTORY MODEL BASED ON DIMENSIONAL ANALYSIS by G. A. VIGNAUX and Sudha JAIN Victoria University, Wellington, New Zealand Published in Asia-Pacific Journal of Operational Research, Vol

More information

56:272 Integer Programming & Network Flows Final Examination -- December 14, 1998

56:272 Integer Programming & Network Flows Final Examination -- December 14, 1998 56:272 Integer Programming & Network Flows Final Examination -- December 14, 1998 Part A: Answer any four of the five problems. (15 points each) 1. Transportation problem 2. Integer LP Model Formulation

More information

An Optimal Algorithm for the Indirect Covering Subtree Problem

An Optimal Algorithm for the Indirect Covering Subtree Problem An Optimal Algorithm for the Indirect Covering Subtree Problem Joachim Spoerhase Lehrstuhl für Informatik I, Universität Würzburg, Am Hubland, 97074 Würzburg, Germany Key words: Graph algorithm, coverage,

More information

8.1. Optimal Binary Search Trees:

8.1. Optimal Binary Search Trees: DATA STRUCTERS WITH C 10CS35 UNIT 8 : EFFICIENT BINARY SEARCH TREES 8.1. Optimal Binary Search Trees: An optimal binary search tree is a binary search tree for which the nodes are arranged on levels such

More information

1 The Permanent Income Hypothesis

1 The Permanent Income Hypothesis The Permanent Income Hypothesis. A two-period model Consider a two-period model where households choose consumption ( 2 ) to solve + 2 max log + log 2 { 2 } µ + + where isthediscountfactor, theinterestrate.

More information

Clustering. SC4/SM4 Data Mining and Machine Learning, Hilary Term 2017 Dino Sejdinovic

Clustering. SC4/SM4 Data Mining and Machine Learning, Hilary Term 2017 Dino Sejdinovic Clustering SC4/SM4 Data Mining and Machine Learning, Hilary Term 2017 Dino Sejdinovic Clustering is one of the fundamental and ubiquitous tasks in exploratory data analysis a first intuition about the

More information

Branch-and-bound: an example

Branch-and-bound: an example Branch-and-bound: an example Giovanni Righini Università degli Studi di Milano Operations Research Complements The Linear Ordering Problem The Linear Ordering Problem (LOP) is an N P-hard combinatorial

More information

Approximation Algorithms

Approximation Algorithms Approximation Algorithms Given an NP-hard problem, what should be done? Theory says you're unlikely to find a poly-time algorithm. Must sacrifice one of three desired features. Solve problem to optimality.

More information

Reducing Clocks in Timed Automata while Preserving Bisimulation

Reducing Clocks in Timed Automata while Preserving Bisimulation Reducing Clocks in Timed Automata while Preserving Bisimulation Shibashis Guha Chinmay Narayan S. Arun-Kumar Indian Institute of Technology Delhi {shibashis, chinmay, sak}@cse.iitd.ac.in arxiv:1404.6613v2

More information

Strategic Design of Robust Global Supply Chains under Uncertainty

Strategic Design of Robust Global Supply Chains under Uncertainty Strategic Design of Robust Global Supply Chains under Uncertainty Marc Goetschalckx, Shabbir Ahmed, Alexander Shapiro, Tjendera Santoso, and Gonzalo Cordova School of Industrial and Systems Engineering,

More information

Optimized Choice of Parameters in interior-point methods for linear programming. Luiz-Rafael dos Santos

Optimized Choice of Parameters in interior-point methods for linear programming. Luiz-Rafael dos Santos Optimized Choice of Parameters in interior-point methods for linear programming Luiz-Rafael dos Santos Joint work with F. Villas-Bôas, A. Oliveira and C. Perin 1st Brazilian Workshop on Interior Point

More information

Chapter S:V. V. Formal Properties of A*

Chapter S:V. V. Formal Properties of A* Chapter S:V V. Formal Properties of A* Properties of Search Space Graphs Auxiliary Concepts Roadmap Completeness of A* Admissibility of A* Efficiency of A* Monotone Heuristic Functions S:V-1 Formal Properties

More information

Simulation. Lecture O1 Optimization: Linear Programming. Saeed Bastani April 2016

Simulation. Lecture O1 Optimization: Linear Programming. Saeed Bastani April 2016 Simulation Lecture O Optimization: Linear Programming Saeed Bastani April 06 Outline of the course Linear Programming ( lecture) Integer Programming ( lecture) Heuristics and Metaheursitics (3 lectures)

More information

Horn Formulae. CS124 Course Notes 8 Spring 2018

Horn Formulae. CS124 Course Notes 8 Spring 2018 CS124 Course Notes 8 Spring 2018 In today s lecture we will be looking a bit more closely at the Greedy approach to designing algorithms. As we will see, sometimes it works, and sometimes even when it

More information

Computer Vision Group Prof. Daniel Cremers. 4. Probabilistic Graphical Models Directed Models

Computer Vision Group Prof. Daniel Cremers. 4. Probabilistic Graphical Models Directed Models Prof. Daniel Cremers 4. Probabilistic Graphical Models Directed Models The Bayes Filter (Rep.) (Bayes) (Markov) (Tot. prob.) (Markov) (Markov) 2 Graphical Representation (Rep.) We can describe the overall

More information

Introduction to Linear Programming

Introduction to Linear Programming Introduction to Linear Programming Eric Feron (updated Sommer Gentry) (updated by Paul Robertson) 16.410/16.413 Historical aspects Examples of Linear programs Historical contributor: G. Dantzig, late 1940

More information

Machine Learning for Video-Based Rendering

Machine Learning for Video-Based Rendering Machine Learning for Video-Based Rendering Arno Schadl arno@schoedl. org Irfan Essa irjan@cc.gatech.edu Georgia Institute of Technology GVU Center / College of Computing Atlanta, GA 30332-0280, USA. Abstract

More information

Lecture: Simulation. of Manufacturing Systems. Sivakumar AI. Simulation. SMA6304 M2 ---Factory Planning and scheduling. Simulation - A Predictive Tool

Lecture: Simulation. of Manufacturing Systems. Sivakumar AI. Simulation. SMA6304 M2 ---Factory Planning and scheduling. Simulation - A Predictive Tool SMA6304 M2 ---Factory Planning and scheduling Lecture Discrete Event of Manufacturing Systems Simulation Sivakumar AI Lecture: 12 copyright 2002 Sivakumar 1 Simulation Simulation - A Predictive Tool Next

More information

MIP-based heuristics for multi-item capacitated lot-sizing problem. with setup times and shortage costs

MIP-based heuristics for multi-item capacitated lot-sizing problem. with setup times and shortage costs MIP-based heuristics for multi-item capacitated lot-sizing problem with setup times and shortage costs Nabil Absi 1,2, Safia Kedad-Sidhoum 1 1 Laboratoire LIP6, 4 place Jussieu, 75 252 Paris Cedex 05,

More information

Quota Management Overview

Quota Management Overview CHAPTER 1 Quota Management Overview Published: May 27, 2013, Introduction This chapter provides an overview of the Quota Manager. The chapter contains the following sections: Information About the Quota

More information

Homework 4: Clustering, Recommenders, Dim. Reduction, ML and Graph Mining (due November 19 th, 2014, 2:30pm, in class hard-copy please)

Homework 4: Clustering, Recommenders, Dim. Reduction, ML and Graph Mining (due November 19 th, 2014, 2:30pm, in class hard-copy please) Virginia Tech. Computer Science CS 5614 (Big) Data Management Systems Fall 2014, Prakash Homework 4: Clustering, Recommenders, Dim. Reduction, ML and Graph Mining (due November 19 th, 2014, 2:30pm, in

More information

Conflict Graphs for Combinatorial Optimization Problems

Conflict Graphs for Combinatorial Optimization Problems Conflict Graphs for Combinatorial Optimization Problems Ulrich Pferschy joint work with Andreas Darmann and Joachim Schauer University of Graz, Austria Introduction Combinatorial Optimization Problem CO

More information

Timing-Anomaly Free Dynamic Scheduling of Task-Based Parallel Applications

Timing-Anomaly Free Dynamic Scheduling of Task-Based Parallel Applications Timing-Anomaly Free Dynamic Scheduling of Task-Based Parallel Applications Petros Voudouris, Per Stenström, Risat Pathan Chalmers University of Technology, Sweden {petrosv, per.stenstrom, risat}@chalmers.se

More information