|
|
| Line 1: |
Line 1: |
| {{For|the programming paradigm|Dynamic programming language}}
| | I'm Hannah and I live in a seaside city in northern Austria, Kohfidisch. I'm 29 and I'm will soon finish my study at Anthropology and Sociology.<br><br>Look into my web page: [https://www.youtube.com/watch?v=rG24yqlTLfg stone floor care in collin county] |
| In [[mathematics]], [[computer science]], [[economics]], and [[bioinformatics]], '''dynamic programming''' is a method for solving complex problems by breaking them down into simpler subproblems. It is applicable to problems exhibiting the properties of [[overlapping subproblem]]s<ref>S. Dasgupta, C.H. Papadimitriou, and U.V. Vazirani, ''''Algorithms'''', p 173, available at http://www.cs.berkeley.edu/~vazirani/algorithms.html</ref> and [[optimal substructure]] (described below). When applicable, the method takes far less time than naive methods that don't take advantage of the subproblem overlap (like [[depth-first search]]).
| |
| | |
| <!-- Note to copyeditors: the term used here is "memoization", NOT "memorization". -->
| |
| The idea behind dynamic programming is quite simple. In general, to solve a given problem, we need to solve different parts of the problem (subproblems), then combine the solutions of the subproblems to reach an overall solution. Often when using a more naive method, many of the subproblems are generated and solved many times. The dynamic programming approach seeks to solve each subproblem only once, thus reducing the number of computations: once the solution to a given subproblem has been computed, it is stored or "[[memoization|memo-ized]]": the next time the same solution is needed, it is simply looked up. This approach is especially useful when the number of repeating subproblems [[exponential growth|grows exponentially]] as a function of the size of the input.
| |
| | |
| Dynamic programming algorithms are used for optimization (for example, finding the shortest path between two points, or the fastest way to multiply many matrices). A dynamic programming algorithm will examine all possible ways to solve the problem and will pick the best solution. Therefore, we can roughly think of dynamic programming as an intelligent, brute-force method that enables us to go through all possible solutions to pick the best one. If the scope of the problem is such that going through all possible solutions is possible and fast enough, dynamic programming guarantees finding the optimal solution. The alternatives are many, such as using a greedy algorithm, which picks the best possible choice "at any possible branch in the road". While a greedy algorithm does not guarantee the optimal solution, it is faster. Fortunately, some greedy algorithms (such as minimum spanning trees) are proven to lead to the optimal solution.
| |
| | |
| For example, let's say that you have to get from point A to point B as fast as possible, in a given city, during rush hour. A dynamic programming algorithm will look into the entire traffic report, looking into all possible combinations of roads you might take, and will only then tell you which way is the fastest. Of course, you might have to wait for a while until the algorithm finishes, and only then can you start driving. The path you will take will be the fastest one (assuming that nothing changed in the external environment). On the other hand, a greedy algorithm will start you driving immediately and will pick the road that looks the fastest at every intersection. As you can imagine, this strategy might not lead to the fastest arrival time, since you might take some "easy" streets and then find yourself hopelessly stuck in a traffic jam.
| |
| | |
| Sometimes, applying memoization to a naive basic recursive solution already results in an optimal dynamic programming solution, however many problems require more sophisticated dynamic programming algorithms. Some of these may be recursive as well but parametrized differently from the naive solution. Others can be more complicated and cannot be implemented as a recursive function with memoization. Examples of these are the two solutions to the Egg Dropping puzzle below.
| |
| | |
| ==History==
| |
| The term ''dynamic programming'' was originally used in the 1940s by [[Richard Bellman]] to describe the process of solving problems where one needs to find the best decisions one after another. By 1953, he refined this to the modern meaning, referring specifically to nesting smaller decision problems inside larger decisions,<ref>http://www.wu-wien.ac.at/usr/h99c/h9951826/bellman_dynprog.pdf</ref> and the field was thereafter recognized by the [[IEEE]] as a [[systems analysis]] and [[engineering]] topic. Bellman's contribution is remembered in the name of the [[Bellman equation]], a central result of dynamic programming which restates an optimization problem in [[Recursion (computer science)|recursive]] form.
| |
| | |
| Bellman explains the reasoning behind the term ''dynamic programming'' in his autobiography, ''Eye of the Hurricane: An Autobiography'' (1984). He explains:
| |
| | |
| :"I spent the Fall quarter (of 1950) at RAND. My first task was to find a name for multistage decision processes. An interesting question is, Where did the name, dynamic programming, come from? The 1950s were not good years for mathematical research. We had a very interesting gentleman in Washington named [[Charles Erwin Wilson|Wilson]]. He was Secretary of Defense, and he actually had a pathological fear and hatred of the word research. I’m not using the term lightly; I’m using it precisely. His face would suffuse, he would turn red, and he would get violent if people used the term research in his presence. You can imagine how he felt, then, about the term mathematical. The RAND Corporation was employed by the Air Force, and the Air Force had Wilson as its boss, essentially. Hence, I felt I had to do something to shield Wilson and the Air Force from the fact that I was really doing mathematics inside the RAND Corporation. What title, what name, could I choose? In the first place I was interested in planning, in decision making, in thinking. But planning, is not a good word for various reasons. I decided therefore to use the word “programming”. I wanted to get across the idea that this was dynamic, this was multistage, this was time-varying I thought, lets kill two birds with one stone. Lets take a word that has an absolutely precise meaning, namely dynamic, in the classical physical sense. It also has a very interesting property as an adjective, and that is its impossible to use the word dynamic in a pejorative sense. Try thinking of some combination that will possibly give it a pejorative meaning. It's impossible. Thus, I thought dynamic programming was a good name. It was something not even a Congressman could object to. So I used it as an umbrella for my activities."
| |
| | |
| The word ''dynamic'' was chosen by Bellman to capture the time-varying aspect of the problems, and because it sounded impressive.<ref name="Eddy">Eddy, S. R., What is dynamic programming?, Nature Biotechnology, 22, 909–910 (2004).</ref> The word ''programming'' referred to the use of the method to find an optimal ''program'', in the sense of a military schedule for training or logistics. This usage is the same as that in the phrases ''[[linear programming]]'' and ''mathematical programming'', a synonym for [[mathematical optimization]].<ref>Nocedal, J.; Wright, S. J.: Numerical Optimization, page 9, Springer, 2006..</ref>
| |
| | |
| ==Overview==
| |
| [[Image:Shortest path optimal substructure.svg|200px|right|thumb|'''Figure 1.''' Finding the shortest path in a graph using optimal substructure; a straight line indicates a single edge; a wavy line indicates a shortest path between the two vertices it connects (other nodes on these paths are not shown); the bold line is the overall shortest path from start to goal.]]
| |
| Dynamic programming is both a mathematical optimization method and a computer programming method. In both contexts it refers to simplifying a complicated problem by breaking it down into simpler subproblems in a [[Recursion|recursive]] manner. While some decision problems cannot be taken apart this way, decisions that span several points in time do often break apart recursively; Bellman called this the "[[Bellman equation#Bellman's Principle of Optimality|Principle of Optimality]]". Likewise, in computer science, a problem that can be broken down recursively is said to have [[optimal substructure]].
| |
| | |
| If subproblems can be nested recursively inside larger problems, so that dynamic programming methods are applicable, then there is a relation between the value of the larger problem and the values of the subproblems.<ref>Cormen, T. H.; Leiserson, C. E.; Rivest, R. L.; Stein, C. (2001), Introduction to Algorithms (2nd ed.), MIT Press & McGraw–Hill, ISBN 0-262-03293-7 . pp. 327–8.</ref> In the optimization literature this relationship is called the [[Bellman equation]].
| |
| | |
| ===Dynamic programming in mathematical optimization===
| |
| In terms of [[mathematical optimization]], dynamic programming usually refers to simplifying a decision by breaking it down into a sequence of decision steps over time. This is done by defining a sequence of '''value functions''' ''V''<sub>1</sub>, ''V''<sub>2</sub>, ..., ''V''<sub>''n''</sub>, with an argument ''y'' representing the '''[[State variable|state]]''' of the system at times ''i'' from 1 to ''n''. The definition of ''V''<sub>''n''</sub>(''y'') is the value obtained in state ''y'' at the last time ''n''. The values ''V''<sub>''i''</sub> at earlier times ''i'' = ''n'' −1, ''n'' − 2, ..., 2, 1 can be found by working backwards, using a [[Recursion|recursive]] relationship called the [[Bellman equation]]. For ''i'' = 2, ..., ''n'', ''V''<sub>''i''−1</sub> at any state ''y'' is calculated from ''V''<sub>''i''</sub> by maximizing a simple function (usually the sum) of the gain from a decision at time ''i'' − 1 and the function ''V''<sub>''i''</sub> at the new state of the system if this decision is made. Since ''V''<sub>''i''</sub> has already been calculated for the needed states, the above operation yields ''V''<sub>''i''−1</sub> for those states. Finally, ''V''<sub>1</sub> at the initial state of the system is the value of the optimal solution. The optimal values of the decision variables can be recovered, one by one, by tracking back the calculations already performed.
| |
| | |
| ===Dynamic programming in bioinformatics===
| |
| Dynamic programming is widely used in bioinformatics for the tasks such as [[sequence alignment]], protein folding, RNA structure prediction and protein-DNA binding. First dynamic programming algorithms for protein-DNA binding were developed in the 1970s independently by [[Charles DeLisi]] in USA and [[Georgii Gurskii]] and [[Alexander Zasedatelev]] in USSR (DeLisi, Biopolymers, 1974, Volume 13, Issue 7, pages 1511–1512, July 1974; Gurskiĭ GV, Zasedatelev AS, Biofizika, 1978 Sep-Oct;23(5):932-46). Recently these algorithms have become very popular in bioinformatics and computational biology, particularly in the studies of [[nucleosome]] positioning and [[transcription factor]] binding (Granek and Clarke 2005; Nechipurenko et al. 2005; Hermsen et al. 2006; Segal et al. 2006; Segal et al. 2008; He et al. 2009; Laurila et al. 2009; Morozov et al. 2009; Wasson and Hartemink 2009; He et al. 2010; Hermsen et al. 2010; Teif and Rippe 2012).
| |
| | |
| ===Dynamic programming in computer programming===
| |
| There are two key attributes that a problem must have in order for dynamic programming to be applicable: [[optimal substructure]] and [[overlapping subproblem]]s. If a problem can be solved by combining optimal solutions to ''non-overlapping'' subproblems, the strategy is called "[[Divide and conquer algorithm|divide and conquer]]" instead. This is why [[mergesort]] and [[quicksort]] are not classified as dynamic programming problems.
| |
| | |
| ''Optimal substructure'' means that the solution to a given optimization problem can be obtained by the combination of optimal solutions to its subproblems. Consequently, the first step towards devising a dynamic programming solution is to check whether the problem exhibits such optimal substructure. Such optimal substructures are usually described by means of [[recursion]]. For example, given a graph ''G=(V,E)'', the shortest path ''p'' from a vertex ''u'' to a vertex ''v'' exhibits optimal substructure: take any intermediate vertex ''w'' on this shortest path ''p''. If ''p'' is truly the shortest path, then it can be split into subpaths ''p<sub>1</sub>'' from ''u'' to ''w'' and ''p<sub>2</sub>'' from ''w'' to ''v'' such that these, in turn, are indeed the shortest paths between the corresponding vertices (by the simple cut-and-paste argument described in ''[[Introduction to Algorithms]]''). Hence, one can easily formulate the solution for finding shortest paths in a recursive manner, which is what the [[Bellman–Ford algorithm]] or the [[Floyd–Warshall algorithm]] does.
| |
| | |
| ''Overlapping'' subproblems means that the space of subproblems must be small, that is, any recursive algorithm solving the problem should solve the same subproblems over and over, rather than generating new subproblems. For example, consider the recursive formulation for generating the Fibonacci series: ''F''<sub>''i''</sub> = ''F''<sub>''i''−1</sub> + ''F''<sub>''i''−2</sub>, with base case ''F''<sub>1</sub> = ''F''<sub>2</sub> = 1. Then ''F''<sub>43</sub> = ''F''<sub>42</sub> + ''F''<sub>41</sub>, and ''F''<sub>42</sub> = ''F''<sub>41</sub> + ''F''<sub>40</sub>. Now ''F''<sub>41</sub> is being solved in the recursive subtrees of both ''F''<sub>43</sub> as well as ''F''<sub>42</sub>. Even though the total number of subproblems is actually small (only 43 of them), we end up solving the same problems over and over if we adopt a naive recursive solution such as this. Dynamic programming takes account of this fact and solves each subproblem only once.
| |
| | |
| [[Image:Fibonacci dynamic programming.svg|thumb|108px|'''Figure 2.''' The subproblem graph for the Fibonacci sequence. The fact that it is not a [[tree structure|tree]] indicates overlapping subproblems.]]
| |
| | |
| This can be achieved in either of two ways:{{Citation needed|date=May 2009}}
| |
| | |
| * ''[[top-down|Top-down approach]]'': This is the direct fall-out of the recursive formulation of any problem. If the solution to any problem can be formulated recursively using the solution to its subproblems, and if its subproblems are overlapping, then one can easily [[memoization|memoize]] or store the solutions to the subproblems in a table. Whenever we attempt to solve a new subproblem, we first check the table to see if it is already solved. If a solution has been recorded, we can use it directly, otherwise we solve the subproblem and add its solution to the table.
| |
| | |
| * ''[[Top-down and bottom-up design|Bottom-up approach]]'': Once we formulate the solution to a problem recursively as in terms of its subproblems, we can try reformulating the problem in a bottom-up fashion: try solving the subproblems first and use their solutions to build-on and arrive at solutions to bigger subproblems. This is also usually done in a tabular form by iteratively generating solutions to bigger and bigger subproblems by using the solutions to small subproblems. For example, if we already know the values of ''F''<sub>41</sub> and ''F''<sub>40</sub>, we can directly calculate the value of ''F''<sub>42</sub>.
| |
| | |
| Some [[programming language]]s can automatically [[memoization|memoize]] the result of a function call with a particular set of arguments, in order to speed up [[call-by-name]] evaluation (this mechanism is referred to as ''[[call-by-need]]''). Some languages make it possible portably (e.g. [[Scheme (programming language)|Scheme]], [[Common Lisp]] or [[Perl]]), some need special extensions (e.g. [[C++]], see<ref>[http://www.apl.jhu.edu/~paulmac/c++-memoization.html Automated Memoization in C]. Apl.jhu.edu (1998-08-21). Retrieved on 2013-08-10.</ref>). Some languages have automatic [[memoization]] <!-- still not a typo for "memor-" --> built in, such as tabled [[Prolog]] and [[J (programming language)|J]], which supports memoization with the ''M.'' adverb.<ref>{{cite web|title=M. Memo|url=http://www.jsoftware.com/help/dictionary/dmcapdot.htm|work=J Vocabulary|publisher=J Software|accessdate=28 October 2011}}</ref> In any case, this is only possible for a [[referential transparency (computer science)|referentially transparent]] function.
| |
| | |
| == Example: Mathematical optimization ==
| |
| | |
| ===Optimal consumption and saving===
| |
| A mathematical optimization problem that is often used in teaching dynamic programming to economists (because it can be solved by hand<ref>Stokey et al., 1989, Chap. 1</ref>) concerns a consumer who lives over the periods <math>t=0,1,2,\ldots,T</math> and must decide how much to consume and how much to save in each period.
| |
| | |
| Let <math>c_t</math> be consumption in period <math>t</math>, and assume consumption yields [[utility]] <math>u(c_t)=\ln(c_t)</math> as long as the consumer lives. Assume the consumer is impatient, so that he [[discounting|discounts]] future utility by a factor <math>b</math> each period, where <math>0<b<1</math>. Let <math>k_t</math> be [[capital (economics)|capital]] in period <math>t</math>. Assume initial capital is a given amount <math>k_0>0</math>, and suppose that this period's capital and consumption determine next period's capital as <math>k_{t+1}=Ak^a_t - c_t</math>, where <math>A</math> is a positive constant and <math>0<a<1</math>. Assume capital cannot be negative. Then the consumer's decision problem can be written as follows:
| |
| | |
| : <math>\max \sum_{t=0}^T b^t \ln(c_t)</math> subject to <math>k_{t+1}=Ak^a_t - c_t \geq 0</math> for all <math>t=0,1,2,\ldots,T</math>
| |
| | |
| Written this way, the problem looks complicated, because it involves solving for all the choice variables <math>c_0, c_1, c_2, \ldots , c_T</math>. (Note that <math>k_0</math> is not a choice variable—the consumer's initial capital is taken as given.)
| |
| | |
| The dynamic programming approach to solving this problem involves breaking it apart into a sequence of smaller decisions. To do so, we define a sequence of ''value functions'' <math>V_t(k)</math>, for <math>t=0,1,2,\ldots,T,T+1</math> which represent the value of having any amount of capital <math>k</math> at each time <math>t</math>. Note that <math>V_{T+1}(k)=0</math>, that is, there is (by assumption) no utility from having capital after death.
| |
| | |
| The value of any quantity of capital at any previous time can be calculated by [[backward induction]] using the [[Bellman equation]]. In this problem, for each <math>t=0,1,2,\ldots,T</math>, the Bellman equation is
| |
| | |
| : <math>V_t(k_t) \, = \, \max \left( \ln(c_t) + b V_{t+1}(k_{t+1}) \right)\text{ subject to }k_{t+1}=Ak^a_t - c_t \geq 0</math>
| |
| | |
| This problem is much simpler than the one we wrote down before, because it involves only two decision variables, <math>c_t</math> and <math>k_{t+1}</math>. Intuitively, instead of choosing his whole lifetime plan at birth, the consumer can take things one step at a time. At time <math>t</math>, his current capital <math>k_t</math> is given, and he only needs to choose current consumption <math>c_t</math> and saving <math>k_{t+1}</math>.
| |
| | |
| To actually solve this problem, we work backwards. For simplicity, the current level of capital is denoted as <math>k</math>. <math>V_{T+1}(k)</math> is already known, so using the Bellman equation once we can calculate <math>V_T(k)</math>, and so on until we get to <math>V_0(k)</math>, which is the ''value'' of the initial decision problem for the whole lifetime. In other words, once we know <math>V_{T-j+1}(k)</math>, we can calculate <math>V_{T-j}(k)</math>, which is the maximum of <math>\ln(c_{T-j}) + b V_{T-j+1}(Ak^a-c_{T-j})</math>, where <math>c_{T-j}</math> is the choice variable and <math>Ak^a-c_{T-j} \ge 0</math>.
| |
| | |
| Working backwards, it can be shown that the value function at time <math>t=T-j</math> is
| |
| | |
| : <math>V_{T-j}(k) \, = \, a \sum_{i=0}^j a^ib^i \ln k + v_{T-j}</math>
| |
| | |
| where each <math>v_{T-j}</math> is a constant, and the optimal amount to consume at time <math>t=T-j</math> is
| |
| | |
| : <math>c_{T-j}(k) \, = \, \frac{1}{\sum_{i=0}^j a^ib^i} Ak^a</math>
| |
| | |
| which can be simplified to
| |
| | |
| : <math>c_{T}(k) \, = \, Ak^a</math>, and <math>c_{T-1}(k) \, = \, \frac{Ak^a}{1+ab}</math>, and <math>c_{T-2}(k) \, = \, \frac{Ak^a}{1+ab+a^2b^2}</math>, etc.
| |
| | |
| We see that it is optimal to consume a larger fraction of current wealth as one gets older, finally consuming all remaining wealth in period <math>T</math>, the last period of life.
| |
| | |
| ==Examples: Computer algorithms==
| |
| | |
| === Dijkstra's algorithm for the shortest path problem ===
| |
| | |
| From a dynamic programming point of view, [[Dijkstra's algorithm]] for the [[shortest path problem]] is a successive approximation scheme that solves the dynamic programming functional equation for the shortest path problem by the '''Reaching''' method.<ref name=sniedovich_06>{{Citation | last = Sniedovich | first = M. | title = Dijkstra’s algorithm revisited: the dynamic programming connexion | journal = Journal of Control and Cybernetics | volume = 35 | issue = 3 | pages = 599–620 | year = 2006 | url = http://matwbn.icm.edu.pl/ksiazki/cc/cc35/cc3536.pdf | format = PDF | postscript = .}} [http://www.ifors.ms.unimelb.edu.au/tutorial/dijkstra_new/index.html Online version of the paper with interactive computational modules.]</ref><ref name=denardo_03>{{Citation | last = Denardo | first = E.V. | title = Dynamic Programming: Models and Applications | publisher = [[Dover Publications]] | location = Mineola, NY | year = 2003 | isbn = 978-0-486-42810-9}}</ref><ref name=sniedovich_10>{{Citation | last = Sniedovich | first = M. | title = Dynamic Programming: Foundations and Principles | publisher = [[Taylor & Francis]] | year = 2010 | isbn = 978-0-8247-4099-3 }}</ref>
| |
| | |
| In fact, Dijkstra's explanation of the logic behind the algorithm,<ref>{{harvnb|Dijkstra|1959|p=270}}</ref> namely
| |
| {{quote|
| |
| '''Problem 2.''' Find the path of minimum total length between two given nodes <math>P</math> and <math>Q</math>.
| |
| | |
| We use the fact that, if <math>R</math> is a node on the minimal path from <math>P</math> to <math>Q</math>, knowledge of the latter implies the knowledge of the minimal path from <math>P</math> to <math>R</math>.
| |
| }}
| |
| | |
| is a paraphrasing of [[Richard Bellman|Bellman's]] famous [[Principle of Optimality]] in the context of the [[shortest path problem]].
| |
| | |
| === Fibonacci sequence ===
| |
| {{unreferenced section|date=May 2013}}
| |
| Here is a naïve implementation of a function finding the ''n''th member of the [[Fibonacci sequence]], based directly on the mathematical definition:
| |
| | |
| '''function''' fib(n)
| |
| '''if''' n = 0 '''return''' 0
| |
| '''if''' n = 1 '''return''' 1
| |
| '''return''' fib(n − 1) + fib(n − 2)
| |
| | |
| Notice that if we call, say, <code>fib(5)</code>, we produce a call tree that calls the function on the same value many different times:
| |
| # <code>fib(5)</code>
| |
| # <code>fib(4) + fib(3)</code>
| |
| # <code>(fib(3) + fib(2)) + (fib(2) + fib(1))</code>
| |
| # <code>((fib(2) + fib(1)) + (fib(1) + fib(0))) + ((fib(1) + fib(0)) + fib(1))</code>
| |
| # <code>(((fib(1) + fib(0)) + fib(1)) + (fib(1) + fib(0))) + ((fib(1) + fib(0)) + fib(1))</code>
| |
| | |
| In particular, <code>fib(2)</code> was calculated three times from scratch. In larger examples, many more values of <code>fib</code>, or ''subproblems'', are recalculated, leading to an exponential time algorithm.
| |
| | |
| Now, suppose we have a simple [[Associative array|map]] object, ''m'', which maps each value of <code>fib</code> that has already been calculated to its result, and we modify our function to use it and update it. The resulting function requires only [[Big-O notation|O]](''n'') time instead of exponential time (but requires [[Big-O notation|O]](''n'') space):
| |
| | |
| '''var''' m := '''''map'''''(0 → 0, 1 → 1)
| |
| '''function''' fib(n)
| |
| '''if ''key''''' n '''is not in ''map''''' m
| |
| m[n] := fib(n − 1) + fib(n − 2)
| |
| '''return''' m[n]
| |
| | |
| This technique of saving values that have already been calculated is called ''[[memoization]]''; <!-- Yes, memoization, not memorization. Not a typo. --> this is the top-down approach, since we first break the problem into subproblems and then calculate and store values.
| |
| | |
| In the '''bottom-up''' approach, we calculate the smaller values of <code>fib</code> first, then build larger values from them. This method also uses O(''n'') time since it contains a loop that repeats n − 1 times, but it only takes constant (O(1)) space, in contrast to the top-down approach which requires O(''n'') space to store the map.
| |
| | |
| '''function''' fib(n)
| |
| '''if''' n = 0
| |
| '''return''' 0
| |
| '''else'''
| |
| '''var''' previousFib := 0, currentFib := 1
| |
| '''repeat''' n − 1 '''times''' ''// loop is skipped if n = 1''
| |
| '''var''' newFib := previousFib + currentFib
| |
| previousFib := currentFib
| |
| currentFib := newFib
| |
| '''return''' currentFib
| |
| | |
| In both examples, we only calculate <code>fib(2)</code> one time, and then use it to calculate both <code>fib(4)</code> and <code>fib(3)</code>, instead of computing it every time either of them is evaluated.
| |
| | |
| Note that the above method actually takes <math>\Omega(n^2)</math> time for large n because addition of two integers with <math>\Omega(n)</math> bits each takes <math>\Omega(n)</math> time. (The ''n''<sup>th</sup> fibonacci number has <math>\Omega(n)</math> bits.) Also, there is a closed form for the Fibonacci sequence, [[Jacques Philippe Marie Binet#Binet's Fibonacci number formula|known as Binet's formula]], from which the <math>n</math>-th term can be [[Computational complexity of mathematical operations|computed]] in approximately <math>O(n(\log n)^2)</math> time, which is more efficient than the above dynamic programming technique. However, the simple recurrence directly gives [[Fibonacci sequence#Matrix form|the matrix form]] that leads to an approximately <math>O(n\log n)</math> algorithm by fast matrix exponentiation.
| |
| | |
| === A type of balanced 0–1 matrix ===
| |
| {{unreferenced section|date=May 2013}}
| |
| Consider the problem of assigning values, either zero or one, to the positions of an {{math|<var>n</var> × <var>n</var>}} matrix, with {{math|<var>n</var>}} even, so that each row and each column contains exactly {{math|<var>n</var> / 2}} zeros and {{math|<var>n</var> / 2}} ones. We ask how many different assignments there are for a given <math>n</math>. For example, when {{math|<var>n</var> {{=}} 4}}, four possible solutions are
| |
| | |
| :<math>\begin{bmatrix}
| |
| 0 & 1 & 0 & 1 \\
| |
| 1 & 0 & 1 & 0 \\
| |
| 0 & 1 & 0 & 1 \\
| |
| 1 & 0 & 1 & 0
| |
| \end{bmatrix} \text{ and } \begin{bmatrix}
| |
| 0 & 0 & 1 & 1 \\
| |
| 0 & 0 & 1 & 1 \\
| |
| 1 & 1 & 0 & 0 \\
| |
| 1 & 1 & 0 & 0
| |
| \end{bmatrix} \text{ and } \begin{bmatrix}
| |
| 1 & 1 & 0 & 0 \\
| |
| 0 & 0 & 1 & 1 \\
| |
| 1 & 1 & 0 & 0 \\
| |
| 0 & 0 & 1 & 1
| |
| \end{bmatrix} \text{ and } \begin{bmatrix}
| |
| 1 & 0 & 0 & 1 \\
| |
| 0 & 1 & 1 & 0 \\
| |
| 0 & 1 & 1 & 0 \\
| |
| 1 & 0 & 0 & 1
| |
| \end{bmatrix}.</math>
| |
| | |
| There are at least three possible approaches: [[Brute-force search|brute force]], [[backtracking]], and dynamic programming.
| |
| | |
| Brute force consists of checking all assignments of zeros and ones and counting those that have balanced rows and columns (<math>n/2</math> zeros and <math>n/2</math> ones). As there are <math>\tbinom{n}{n/2}^n</math> possible assignments, this strategy is not practical except maybe up to <math>n=6</math>.
| |
| | |
| Backtracking for this problem consists of choosing some order of the matrix elements and recursively placing ones or zeros, while checking that in every row and column the number of elements that have not been assigned plus the number of ones or zeros are both at least ''n / 2''. While more sophisticated than brute force, this approach will visit every solution once, making it impractical for ''n'' larger than six, since the number of solutions is already 116,963,796,250 for ''n = 10'', as we shall see.
| |
| | |
| Dynamic programming makes it possible to count the number of solutions without visiting them all. Imagine backtracking values for the first row – what information would we require about the remaining rows, in order to be able to accurately count the solutions obtained for each first row value? We consider {{math|<var>k</var> × <var>n</var>}} boards, where {{math|1 ≤ <var>k</var> ≤ <var>n</var>}}, whose <math>k</math> rows contain <math>n/2</math> zeros and <math>n/2</math> ones. The function ''f'' to which [[memoization]] is applied maps vectors of ''n'' pairs of integers to the number of admissible boards (solutions). There is one pair for each column, and its two components indicate respectively the number of zeros and ones that have yet to be placed in that column. We seek the value of <math> f((n/2, n/2), (n/2, n/2), \ldots (n/2, n/2)) </math> (<math>n</math> arguments or one vector of <math>n</math> elements). The process of subproblem creation involves iterating over every one of <math>\tbinom{n}{n/2}</math> possible assignments for the top row of the board, and going through every column, subtracting one from the appropriate element of the pair for that column, depending on whether the assignment for the top row contained a zero or a one at that position. If any one of the results is negative, then the assignment is invalid and does not contribute to the set of solutions (recursion stops). Otherwise, we have an assignment for the top row of the {{math|<var>k</var> × <var>n</var>}} board and recursively compute the number of solutions to the remaining {{math|(<var>k</var> − 1) × <var>n</var>}} board, adding the numbers of solutions for every admissible assignment of the top row and returning the sum, which is being memoized. The base case is the trivial subproblem, which occurs for a {{math|1 × <var>n</var>}} board. The number of solutions for this board is either zero or one, depending on whether the vector is a permutation of {{math|<var>n</var> / 2}} <math>(0, 1)</math> and {{math|<var>n</var> / 2}} <math>(1, 0)</math> pairs or not.
| |
| | |
| For example, in the first two boards shown above the sequences of vectors would be
| |
| <PRE WIDTH=80>
| |
| ((2, 2) (2, 2) (2, 2) (2, 2)) ((2, 2) (2, 2) (2, 2) (2, 2)) k = 4
| |
| 0 1 0 1 0 0 1 1
| |
| | |
| ((1, 2) (2, 1) (1, 2) (2, 1)) ((1, 2) (1, 2) (2, 1) (2, 1)) k = 3
| |
| 1 0 1 0 0 0 1 1
| |
| | |
| ((1, 1) (1, 1) (1, 1) (1, 1)) ((0, 2) (0, 2) (2, 0) (2, 0)) k = 2
| |
| 0 1 0 1 1 1 0 0
| |
| | |
| ((0, 1) (1, 0) (0, 1) (1, 0)) ((0, 1) (0, 1) (1, 0) (1, 0)) k = 1
| |
| 1 0 1 0 1 1 0 0
| |
| | |
| ((0, 0) (0, 0) (0, 0) (0, 0)) ((0, 0) (0, 0), (0, 0) (0, 0))
| |
| </PRE>
| |
| | |
| The number of solutions {{OEIS|id=A058527}} is
| |
| :<math> 1,\, 2,\, 90,\, 297200,\, 116963796250,\, 6736218287430460752, \ldots </math>
| |
| | |
| Links to the MAPLE implementation of the dynamic programming approach may be found among the [[Dynamic programming#External links|external links]].
| |
| | |
| === Checkerboard ===
| |
| {{unreferenced section|date=May 2013}}
| |
| Consider a [[checkerboard]] with ''n'' × ''n'' squares and a cost-function ''c''(''i'', ''j'') which returns a cost associated with square ''i'', ''j'' (''i'' being the row, ''j'' being the column). For instance (on a 5 × 5 checkerboard),
| |
| | |
| {| class="wikitable" style="text-align:center"
| |
| |-
| |
| ! 5
| |
| | 6 || 7 || 4 || 7 || 8
| |
| |-
| |
| ! 4
| |
| | 7 || 6 || 1 || 1 || 4
| |
| |-
| |
| ! 3
| |
| | 3 || 5 || 7 || 8 || 2
| |
| |-
| |
| ! 2
| |
| | – || 6 || 7 || 0 || –
| |
| |-
| |
| ! 1
| |
| | – || – || *5* || – || –
| |
| |-
| |
| !width="15"| !! style="width:15px;"|1 !! style="width:15px;"|2 !! style="width:15px;"|3 !! style="width:15px;"|4 !! style="width:15px;"|5
| |
| |}
| |
| Thus ''c''(1, 3) = 5
| |
| | |
| Let us say you had a checker that could start at any square on the first rank (i.e., row) and you wanted to know the shortest path (sum of the costs of the visited squares are at a minimum) to get to the last rank, assuming the checker could move only diagonally left forward, diagonally right forward, or straight forward. That is, a checker on (1,3) can move to (2,2), (2,3) or (2,4).
| |
| {| class="wikitable" style="text-align:center"
| |
| |-
| |
| ! 5
| |
| | || || || ||
| |
| |-
| |
| ! 4
| |
| | || || || ||
| |
| |-
| |
| ! 3
| |
| | || || || ||
| |
| |-
| |
| ! 2
| |
| | || x || x || x ||
| |
| |-
| |
| ! 1
| |
| | || || o || ||
| |
| |-
| |
| !width="15"| !! style="width:15px;"|1 !! style="width:15px;"|2 !! style="width:15px;"|3 !! style="width:15px;"|4 !! style="width:15px;"|5
| |
| |}
| |
| | |
| This problem exhibits '''optimal substructure'''. That is, the solution to the entire problem relies on solutions to subproblems. Let us define a function ''q''(''i'', ''j'') as
| |
| | |
| :''q''(''i'', ''j'') = the minimum cost to reach square (''i'', ''j'').
| |
| | |
| If we can find the values of this function for all the squares at rank ''n'', we pick the minimum and follow that path backwards to get the shortest path.
| |
| | |
| Note that ''q''(''i'', ''j'') is equal to the minimum cost to get to any of the three squares below it (since those are the only squares that can reach it) plus ''c''(''i'', ''j''). For instance:
| |
| | |
| {| class="wikitable" style="text-align:center"
| |
| |-
| |
| ! 5
| |
| | || || || ||
| |
| |-
| |
| ! 4
| |
| | || || A || ||
| |
| |-
| |
| ! 3
| |
| | || B || C || D ||
| |
| |-
| |
| ! 2
| |
| | || || || ||
| |
| |-
| |
| ! 1
| |
| | || || || ||
| |
| |-
| |
| !width="15"| !! style="width:15px;"|1 !! style="width:15px;"|2 !! style="width:15px;"|3 !! style="width:15px;"|4 !! style="width:15px;"|5
| |
| |}
| |
| | |
| : <math>q(A) = \min(q(B),q(C),q(D))+c(A) \, </math>
| |
| | |
| Now, let us define ''q''(''i'', ''j'') in somewhat more general terms:
| |
| | |
| : <math>q(i,j)=\begin{cases} \infty & j < 1 \text{ or }j > n \\ c(i, j) & i = 1 \\ \min(q(i-1, j-1), q(i-1, j), q(i-1, j+1)) + c(i,j) & \text{otherwise.}\end{cases}</math>
| |
| | |
| The first line of this equation is there to make the recursive property simpler (when dealing with the edges, so we need only one recursion). The second line says what happens in the last rank, to provide a base case. The third line, the recursion, is the important part. It is similar to the A,B,C,D example. From this definition we can make a straightforward recursive code for ''q''(''i'', ''j''). In the following pseudocode, ''n'' is the size of the board, <code>c(i, j)</code> is the cost-function, and <code>min()</code> returns the minimum of a number of values:
| |
| | |
| '''function''' minCost(i, j)
| |
| '''if''' j < 1 '''or''' j > n
| |
| '''return''' infinity
| |
| '''else if''' i = 1
| |
| '''return''' c(i, j)
| |
| '''else'''
| |
| '''return''' '''min'''( minCost(i-1, j-1), minCost(i-1, j), minCost(i-1, j+1) ) + c(i, j)
| |
| | |
| It should be noted that this function only computes the path-cost, not the actual path. We will get to the path soon. This, like the Fibonacci-numbers example, is horribly slow since it wastes time recomputing the same shortest paths over and over. However, we can compute it much faster in a bottom-up fashion if we store path-costs in a two-dimensional array <code>q[i, j]</code> rather than using a function. This avoids recomputation; before computing the cost of a path, we check the array <code>q[i, j]</code> to see if the path cost is already there.
| |
| | |
| We also need to know what the actual shortest path is. To do this, we use another array <code>p[i, j]</code>, a ''predecessor array''. This array implicitly stores the path to any square ''s'' by storing the previous node on the shortest path to ''s'', i.e. the predecessor. To reconstruct the path, we lookup the predecessor of ''s'', then the predecessor of that square, then the predecessor of that square, and so on, until we reach the starting square. Consider the following code:
| |
| | |
| '''function''' computeShortestPathArrays()
| |
| '''for''' x '''from''' 1 '''to''' n
| |
| q[1, x] := c(1, x)
| |
| '''for''' y '''from''' 1 '''to''' n
| |
| q[y, 0] := infinity
| |
| q[y, n + 1] := infinity
| |
| '''for''' y '''from''' 2 '''to''' n
| |
| '''for''' x '''from''' 1 '''to''' n
| |
| m := min(q[y-1, x-1], q[y-1, x], q[y-1, x+1])
| |
| q[y, x] := m + c(y, x)
| |
| '''if''' m = q[y-1, x-1]
| |
| p[y, x] := -1
| |
| '''else if''' m = q[y-1, x]
| |
| p[y, x] := 0
| |
| '''else'''
| |
| p[y, x] := 1
| |
| | |
| Now the rest is a simple matter of finding the minimum and printing it.
| |
| | |
| '''function''' computeShortestPath()
| |
| computeShortestPathArrays()
| |
| minIndex := 1
| |
| min := q[n, 1]
| |
| '''for''' i '''from''' 2 '''to''' n
| |
| '''if''' q[n, i] < min
| |
| minIndex := i
| |
| min := q[n, i]
| |
| printPath(n, minIndex)
| |
| | |
| '''function''' printPath(y, x)
| |
| '''print'''(x)
| |
| '''print'''("<-")
| |
| '''if''' y = 2
| |
| '''print'''(x + p[y, x])
| |
| '''else'''
| |
| printPath(y-1, x + p[y, x])
| |
| | |
| === Sequence alignment ===
| |
| | |
| In [[genetics]], [[sequence alignment]] is an important application where dynamic programming is essential.<ref name="Eddy"/> Typically, the problem consists of transforming one sequence into another using edit operations that replace, insert, or remove an element. Each operation has an associated cost, and the goal is to find the [[edit distance|sequence of edits with the lowest total cost]].
| |
| | |
| The problem can be stated naturally as a recursion, a sequence A is optimally edited into a sequence B by either:
| |
| | |
| # inserting the first character of B, and performing an optimal alignment of A and the tail of B
| |
| # deleting the first character of A, and performing the optimal alignment of the tail of A and B
| |
| # replacing the first character of A with the first character of B, and performing optimal alignments of the tails of A and B.
| |
| | |
| The partial alignments can be tabulated in a matrix, where cell (i,j) contains the cost of the optimal alignment of A[1..i] to B[1..j]. The cost in cell (i,j) can be calculated by adding the cost of the relevant operations to the cost of its neighboring cells, and selecting the optimum.
| |
| | |
| Different variants exist, see [[Smith–Waterman algorithm]] and [[Needleman–Wunsch algorithm]].
| |
| | |
| === Tower of Hanoi puzzle ===
| |
| [[Image:Tower of Hanoi.jpeg|300px|thumb|A model set of the Towers of Hanoi (with 8 disks)]]
| |
| [[Image:Tower of Hanoi 4.gif|300px|thumb|An animated solution of the '''Tower of Hanoi''' puzzle for ''T(4,3)''.]]
| |
| The '''[[Tower of Hanoi]]''' or '''Towers of [[Hanoi]]''' is a [[mathematical game]] or [[puzzle]]. It consists of three rods, and a number of disks of different sizes which can slide onto any rod. The puzzle starts with the disks in a neat stack in ascending order of size on one rod, the smallest at the top, thus making a conical shape.
| |
| | |
| The objective of the puzzle is to move the entire stack to another rod, obeying the following rules:
| |
| | |
| * Only one disk may be moved at a time.
| |
| * Each move consists of taking the upper disk from one of the rods and sliding it onto another rod, on top of the other disks that may already be present on that rod.
| |
| * No disk may be placed on top of a smaller disk.
| |
| | |
| The dynamic programming solution consists of solving the [[Bellman equation|functional equation]]
| |
| | |
| : S(n,h,t) = S(n-1,h, not(h,t)) ; S(1,h,t) ; S(n-1,not(h,t),t)
| |
| | |
| where n denotes the number of disks to be moved, h denotes the home rod, t denotes the target rod, not(h,t) denotes the third rod (neither h nor t), ";" denotes concatenation, and
| |
| | |
| : S(n, h, t) := solution to a problem consisting of n disks that are to be moved from rod h to rod t.
| |
| | |
| Note that for n=1 the problem is trivial, namely S(1,h,t) = "move a disk from rod h to rod t" (there is only one disk left).
| |
| | |
| The number of moves required by this solution is 2<sup>''n''</sup> − 1. If the objective is to '''maximize''' the number of moves (without cycling) then the dynamic programming [[Bellman equation|functional equation]] is slightly more complicated and 3<sup>''n''</sup> − 1 moves are required.<ref>{{Citation |author=Moshe Sniedovich |title= OR/MS Games: 2. The Towers of Hanoi Problem, |journal=INFORMS Transactions on Education |volume=3 |issue=1 |year=2002 |pages=34–51 |url=http://archive.ite.journal.informs.org/Vol3No1/Sniedovich/ |postscript=.}}</ref>
| |
| | |
| === Egg dropping puzzle ===
| |
| The following is a description of the instance of this famous [[puzzle]] involving n=2 eggs and a building with H=36 floors:<ref>Konhauser J.D.E., Velleman, D., and Wagon, S. (1996). [http://www.cambridge.org/uk/catalogue/catalogue.asp?isbn=9780883853252 Which way did the Bicycle Go?] Dolciani Mathematical Expositions – No 18. [[The Mathematical Association of America]].</ref>
| |
| :Suppose that we wish to know which stories in a 36-story building are safe to drop eggs from, and which will cause the eggs to break on landing (using [[U.S. English]] terminology, in which the first floor is at ground level). We make a few assumptions:
| |
| :* An egg that survives a fall can be used again.
| |
| :* A broken egg must be discarded.
| |
| :* The effect of a fall is the same for all eggs.
| |
| :* If an egg breaks when dropped, then it would break if dropped from a higher window.
| |
| :* If an egg survives a fall, then it would survive a shorter fall.
| |
| :* It is not ruled out that the first-floor windows break eggs, nor is it ruled out that the 36th-floor windows do not cause an egg to break.
| |
| | |
| : If only one egg is available and we wish to be sure of obtaining the right result, the experiment can be carried out in only one way. Drop the egg from the first-floor window; if it survives, drop it from the second-floor window. Continue upward until it breaks. In the worst case, this method may require 36 droppings. Suppose 2 eggs are available. What is the least number of egg-droppings that is guaranteed to work in all cases?
| |
| | |
| To derive a dynamic programming [[Bellman equation|functional equation]] for this puzzle, let the '''state''' of the dynamic programming model be a pair s = (n,k), where
| |
| | |
| : ''n'' = number of test eggs available, ''n'' = 0, 1, 2, 3, ..., ''N'' − 1.
| |
| : ''k'' = number of (consecutive) floors yet to be tested, ''k'' = 0, 1, 2, ..., ''H'' − 1.
| |
| | |
| For instance, ''s'' = (2,6) indicates that two test eggs are available and 6 (consecutive) floors are yet to be tested. The initial state of the process is ''s'' = (''N'',''H'') where ''N'' denotes the number of test eggs available at the commencement of the experiment. The process terminates either when there are no more test eggs (''n'' = 0) or when ''k'' = 0, whichever occurs first. If termination occurs at state ''s'' = (0,''k'') and ''k'' > 0, then the test failed.
| |
| | |
| Now, let
| |
| | |
| : ''W''(''n'',''k'') = minimum number of trials required to identify the value of the critical floor under the worst case scenario given that the process is in state ''s'' = (''n'',''k'').
| |
| | |
| Then it can be shown that<ref name="sniedovich_03">Sniedovich, M. (2003). [http://archive.ite.journal.informs.org/Vol4No1/Sniedovich/index.php The joy of egg-dropping in Braunschweig and Hong Kong]. INFORMS Transactions on Education, 4(1) 48–64.</ref>
| |
| | |
| : ''W''(''n'',''k'') = 1 + min{max(''W''(''n'' − 1, ''x'' − 1), ''W''(''n'',''k'' − ''x'')): ''x'' = 1, 2, ..., ''k'' }
| |
| | |
| with ''W''(''n'',1) = 1 for all ''n'' > 0 and ''W''(1,''k'') = ''k'' for all ''k''. It is easy to solve this equation iteratively by systematically increasing the values of ''n'' and ''k''.
| |
| | |
| An [http://archive.ite.journal.informs.org/Vol4No1/Sniedovich/index.php interactive online facility] is available for experimentation with this model as well as with other versions of this puzzle (e.g. when the objective is to minimize the '''expected value''' of the number of trials.)<ref name="sniedovich_03" />
| |
| | |
| ====Faster DP solution using a different parametrization====
| |
| {{unreferenced section|date=May 2013}}
| |
| Notice that the above solution takes <math>O( n k^2 )</math> time. This can be improved to <math>O( n k \log k )</math> time by binary searching on the optimal <math>x</math> in the above recurrence, since <math>W(n-1,x-1)</math> is increasing in <math>x</math> while <math>W(n,k-x)</math> is decreasing in <math>x</math>, thus a local minimum of <math>\max(W(n-1,x-1),W(n,k-x))</math> is a global minimum. However, there is a much faster solution that involves a different parametrization of the problem:
| |
| | |
| Let <math>k</math> be the total number of floors such that the eggs break when dropped from the <math>k</math>th floor (The example above is equivalent to taking <math>k=37</math>).
| |
| | |
| Let <math>m</math> be the minimum floor from which the egg must be dropped to be broken.
| |
| | |
| Let <math>f(t,n)</math> be the maximum number of values of <math>m</math> that are distinguishable using <math>t</math> tries and <math>n</math> eggs.
| |
| | |
| Then <math>f(t,0) = f(0,n) = 1</math> for all <math>t,n \geq 0</math>.
| |
| | |
| Let <math>x</math> be the floor from which the first egg is dropped in the optimal strategy.
| |
| | |
| If the first egg broke, <math>m</math> is from <math>1</math> to <math>x</math> and distinguishable using at most <math>t-1</math> tries and <math>n-1</math> eggs.
| |
| | |
| If the first egg did not break, <math>m</math> is from <math>x+1</math> to <math>k</math> and distinguishable using <math>t-1</math> tries and <math>n</math> eggs.
| |
| | |
| Therefore <math>f(t,n) = f(t-1,n-1) + f(t-1,n)</math>.
| |
| | |
| Then the problem is equivalent to finding the minimum <math>x</math> such that <math>f(x,n) \geq k</math>.
| |
| | |
| To do so, we could compute <math>{ f(t,i) : 0 \leq i \leq n }</math> in order of increasing <math>t</math>, which would take <math>O( n x )</math> time.
| |
| | |
| Thus, if we separately handle the case of <math>n=1</math>, the algorithm would take <math>O( n \sqrt{k} )</math> time.
| |
| | |
| But the recurrence relation can in fact be solved, giving <math>f(t,n) = \sum_{i=0}^{n}{ \binom{t}{i} }</math>, which can be computed in <math>O(n)</math> time using the identity <math>\binom{t}{i+1} = \binom{t}{i} \frac{t-i}{i+1}</math> for all <math>i \geq 0</math>.
| |
| | |
| Since <math>f(t,n) \leq f(t+1,n)</math> for all <math>t \geq 0</math>, we can binary search on <math>t</math> to find <math>x</math>, giving an <math>O( n \log k )</math> algorithm.
| |
| | |
| === Matrix chain multiplication ===
| |
| {{unreferenced section|date=May 2013}}
| |
| <!--Show how the placement of parentheses affects the number of scalar multiplications required when multiplying a bunch of matrices.
| |
| Show how to write a dynamic program to calculate the optimal parentheses placement.
| |
| This is such a long example that it might be better to make it its own article.-->
| |
| | |
| Matrix chain multiplication is a well known example that demonstrates utility of dynamic programming. For example, engineering applications often have to multiply a chain of matrices. It is not surprising to find matrices of large dimensions, for example 100×100. Therefore, our task is to multiply matrices A1, A2, .... An. As we know from basic linear algebra, matrix multiplication is not commutative, but is associative; and we can multiply only two matrices at a time. So, we can multiply this chain of matrices in many different ways, for example:
| |
| | |
| : ((A<sub>1</sub> × A<sub>2</sub>) × A<sub>3</sub>) × ... A<sub>n</sub>
| |
| | |
| : A<sub>1</sub>×(((A<sub>2</sub>×A<sub>3</sub>)× ... ) × A<sub>n</sub>)
| |
| | |
| : (A<sub>1</sub> × A<sub>2</sub>) × (A<sub>3</sub> × ... A<sub>n</sub>)
| |
| | |
| and so on. There are numerous ways to multiply this chain of matrices. They will all produce the same final result, however they will take more or less time to compute, based on which particular matrices are multiplied. If matrix A has dimensions m×n and matrix B has dimensions n×q, then matrix C=A×B will have dimensions m×q, and will require m*n*q scalar multiplications (using a simplistic matrix multiplication algorithm for purposes of illustration).
| |
| | |
| For example, let us multiply matrices A, B and C. Let us assume that their dimensions are m×n, n×p, and p×s, respectively. Matrix A×B×C will be of size m×s and can be calculated in two ways shown below:
| |
| | |
| # Ax(B×C) This order of matrix multiplication will require nps + mns scalar multiplications.
| |
| # (A×B)×C This order of matrix multiplication will require mnp + mps scalar calculations.
| |
| | |
| Let us assume that m = 10, n = 100, p = 10 and s = 1000. So, the first way to multiply the chain will require 1,000,000 + 1,000,000 calculations. The second way will require only 10,000+100,000 calculations. Obviously, the second way is faster, and we should multiply the matrices using that arrangement of parenthesis.
| |
| | |
| Therefore, our conclusion is that the order of parenthesis matters, and that our task is to find the optimal order of parenthesis.
| |
| | |
| At this point, we have several choices, one of which is to design a dynamic programming algorithm that will split the problem into overlapping problems and calculate the optimal arrangement of parenthesis. The dynamic programming solution is presented below.
| |
| | |
| Let's call m[i,j] the minimum number of scalar multiplications needed to multiply a chain of matrices from matrix i to matrix j (i.e. A<sub>i</sub> × .... × A<sub>j</sub>, i.e. i<=j). We split the chain at some matrix k, such that i <= k < j, and try to find out which combination produces minimum m[i,j].
| |
| | |
| The formula is:
| |
| '''if''' i = j, m[i,j]= 0
| |
| '''if''' i < j, m[i,j]= min over all possible values of k ({{tmath|m[i,k]+m[k+1,j] + p_{{(}}i-1{{)}}*p_k*p_j}})
| |
| where ''k'' is changed from ''i'' to ''j'' − 1.
| |
| | |
| *{{tmath|p_{{(}}i-1{{)}}}} is the row dimension of matrix i,
| |
| *{{tmath|p_k}} is the column dimension of matrix k,
| |
| *{{tmath|p_j}} is the column dimension of matrix j.
| |
| | |
| This formula can be coded as shown below, where input parameter "chain" is the chain of matrices, i.e. {{tmath|A_1, A_2, ... A_n}}:
| |
| | |
| '''function''' OptimalMatrixChainParenthesis(chain)
| |
| n = length(chain)
| |
| '''for''' i = 1, n
| |
| m[i,i] = 0 ''//since it takes no calculations to multiply one matrix''
| |
| '''for''' len = 2, n
| |
| '''for''' i = 1, n - len + 1
| |
| '''for''' j = i, len -1
| |
| m[i,j] = infinity ''//so that the first calculation updates ''
| |
| '''for''' k = i, j-1
| |
| {{tmath|1=q = m[i, k] + m[k+1, j] + p_{{(}}i-1{{)}}*p_k*p_j}}
| |
| '''if''' q < m[i, j] ''// the new order of parenthesis is better than what we had''
| |
| m[i, j] = q ''//update''
| |
| s[i, j] = k ''//record which k to split on, i.e. where to place the parenthesis''
| |
| | |
| So far, we have calculated values for all possible {{math|''m''[''i'', ''j'']}}, the minimum number of calculations to multiply a chain from matrix ''i'' to matrix ''j'', and we have recorded the corresponding "split point"{{math|''s''[''i'', ''j'']}}. For example, if we are multiplying chain A<sub>1</sub>×A<sub>2</sub>×A<sub>3</sub>×A<sub>4</sub>, and it turns out that {{math|1=''m''[1, 3] = 100}} and {{math|1=''s''[1, 3] = 2}}, that means that the optimal placement of parenthesis for matrices 1 to 3 is (A1×A2)×A3 and to multiply those matrices will require 100 scalar calculation.
| |
| | |
| This algorithm will produce "tables" ''m''[, ] and ''s''[, ] that will have entries for all possible values of i and j. The final solution for the entire chain is m[1, n], with corresponding split at s[1, n]. Unraveling the solution will be recursive, starting from the top and continuing until we reach the base case, i.e. multiplication of single matrices.
| |
| | |
| Therefore, the next step is to actually split the chain, i.e. to place the parenthesis where they (optimally) belong. For this purpose we could use the following algorithm:
| |
| | |
| '''function''' PrintOptimalParenthesis(s, i, j)
| |
| '''if''' i = j
| |
| print "A"i
| |
| '''else'''
| |
| print "(" PrintOptimalParenthesis(s, i, s[i, j]) PrintOptimalParenthesis(s, s[i, j] + 1, j) ")"
| |
|
| |
| Of course, this algorithm is not useful for actual multiplication. This algorithm is just a user-friendly way to see what the result looks like.
| |
| | |
| To actually multiply the matrices using the proper splits, we need the following algorithm:
| |
| <source lang="javascript" enclose="div">
| |
| function MatrixChainMultiply(chain from 1 to n) //returns the final matrix, i.e. A1×A2×... ×An
| |
| OptimalMatrixChainParenthesis(chain from 1 to n) //this will produce s[ . ] and m[ . ] "tables"
| |
| OptimalMatrixMultiplication(s, chain from 1 to n) //actually multiply
| |
| | |
| function OptimalMatrixMultiplication(s, i, j) //returns the result of multiplying a chain of matrices from Ai to Aj in optimal way
| |
| if i < j
| |
| //keep on splitting the chain and multiplying the matrices in left and right sides
| |
| LeftSide = OptimalMatrixMultiplication(s, i, s[i, j])
| |
| RightSide = OptimalMatrixMultiplication(s, s[i, j] + 1, j)
| |
| return MatrixMultiply(LeftSide, RightSide)
| |
| else if i = j
| |
| return Ai //matrix at position i
| |
| else
| |
| print "error, i <= j must hold"
| |
| | |
| function MatrixMultiply(A, B) //function that multiplies two matrices
| |
| if columns(A) = rows(B)
| |
| for i = 1, rows(A)
| |
| for j = 1, columns(B)
| |
| C[i, j] = 0
| |
| for k = 1, columns(A)
| |
| C[i, j] = C[i, j] + A[i, k]*B[k, j]
| |
| return C
| |
| else
| |
| print "error, incompatible dimensions."
| |
| </source>
| |
| | |
| == Algorithms that use dynamic programming ==
| |
| {{unreferenced section|date=May 2013}}
| |
| * Recurrent solutions to [[lattice models]] for protein-DNA binding
| |
| * [[Backward induction]] as a solution method for finite-horizon [[discrete-time]] dynamic optimization problems
| |
| * [[Method of undetermined coefficients]] can be used to solve the [[Bellman equation]] in infinite-horizon, discrete-time, [[discounting|discounted]], [[Time-invariant system|time-invariant]] dynamic optimization problems
| |
| * Many [[String (computer science)|string]] algorithms including [[longest common subsequence problem|longest common subsequence]], [[Longest increasing subsequence problem|longest increasing subsequence]], [[Longest common substring problem|longest common substring]], [[Levenshtein distance]] (edit distance)
| |
| * Many algorithmic problems on [[undirected graph|graphs]] can be solved efficiently for graphs of bounded [[treewidth]] or bounded [[clique-width]] by using dynamic programming on a [[tree decomposition]] of the graph.
| |
| * The [[CYK algorithm|Cocke–Younger–Kasami (CYK) algorithm]] which determines whether and how a given string can be generated by a given [[context-free grammar]]
| |
| * [[Word wrap|Knuth's word wrapping algorithm]] that minimizes raggedness when word wrapping text
| |
| * The use of [[transposition table]]s and [[refutation table]]s in [[computer chess]]
| |
| * The [[Viterbi algorithm]] (used for [[hidden Markov model]]s)
| |
| * The [[Earley algorithm]] (a type of [[chart parser]])
| |
| * The [[Needleman–Wunsch algorithm|Needleman–Wunsch]] and other algorithms used in [[bioinformatics]], including [[sequence alignment]], [[structural alignment]], [[RNA structure|RNA structure prediction]]
| |
| * [[Floyd–Warshall algorithm|Floyd's all-pairs shortest path algorithm]]
| |
| * Optimizing the order for [[chain matrix multiplication]]
| |
| * [[Pseudo-polynomial time]] algorithms for the [[subset sum problem|subset sum]] and [[Knapsack problem|knapsack]] and [[Partition problem|partition]] problems
| |
| * The [[dynamic time warping]] algorithm for computing the global distance between two time series
| |
| * The [[Patricia Selinger|Selinger]] (a.k.a. [[IBM System R|System R]]) algorithm for relational database query optimization
| |
| * [[De Boor algorithm]] for evaluating B-spline curves
| |
| * [[Duckworth–Lewis method]] for resolving the problem when games of cricket are interrupted
| |
| * The value iteration method for solving [[Markov decision process]]es
| |
| * Some graphic image edge following selection methods such as the "magnet" selection tool in [[Photoshop]]
| |
| * Some methods for solving [[interval scheduling]] problems
| |
| * Some methods for solving [[word wrap]] problems
| |
| * Some methods for solving the [[travelling salesman problem]], either exactly (in [[exponential time]]) or approximately (e.g. via the [[bitonic tour]])
| |
| * [[Recursive least squares]] method
| |
| * [[Beat (music)|Beat]] tracking in [[music information retrieval]]
| |
| * Adaptive-critic training strategy for [[artificial neural networks]]
| |
| * Stereo algorithms for solving the [[correspondence problem]] used in stereo vision
| |
| * [[Seam carving]] (content aware image resizing)
| |
| * The [[Bellman–Ford algorithm]] for finding the shortest distance in a graph
| |
| * Some approximate solution methods for the [[linear search problem]]
| |
| * Kadane's algorithm for the [[maximum subarray problem]]
| |
| | |
| == See also ==
| |
| {{Portal|Computer science}}
| |
| * [[Convexity in economics]]
| |
| * [[Greedy algorithm]]
| |
| * [[Non-convexity (economics)]]
| |
| * [[Stochastic programming]]
| |
| | |
| ==References==
| |
| {{Reflist|30em}}
| |
| | |
| == Further reading ==
| |
| *{{citation|last1=Adda|first1=Jerome|last2=Cooper|first2=Russell|year=2003|url=http://www.eco.utexas.edu/~cooper/dynprog/dynprog1.html|title=Dynamic Economics|publisher=MIT Press}}. An accessible introduction to dynamic programming in economics. The link contains sample programs.
| |
| *{{citation|first=Richard|last=Bellman|authorlink=Richard Bellman|title=The theory of dynamic programming|journal=[[Bulletin of the American Mathematical Society]]|year=1954|volume=60|pages=503–516|doi=10.1090/S0002-9904-1954-09848-8|mr=0067459|issue=6}}. Includes an extensive bibliography of the literature in the area, up to the year 1954.
| |
| *{{citation|first=Richard|last=Bellman|authorlink=Richard Bellman|year=1957|title=Dynamic Programming|publisher=Princeton University Press}}. Dover paperback edition (2003), ISBN 0-486-42809-5.
| |
| *{{citation|last=Bertsekas|first=D. P.|year=2000|title=Dynamic Programming and Optimal Control|edition=2nd|publisher=Athena Scientific|isbn=1-886529-09-4}}. In two volumes.
| |
| *{{citation|last1=Cormen|first1=Thomas H.|author1-link=Thomas H. Cormen|last2=Leiserson|first2=Charles E.|author2-link=Charles E. Leiserson|last3=Rivest|first3=Ronald L.|author3-link=Ronald L. Rivest|last4=Stein|first4=Clifford|author4-link=Clifford Stein|year=2001|title=[[Introduction to Algorithms]]|edition=2nd|publisher=MIT Press & McGraw–Hill|isbn=0-262-03293-7}}. Especially pp. 323–69.
| |
| *{{citation|first1=Stuart E.|last1=Dreyfus|first2=Averill M.|last2=Law|year=1977|title=The art and theory of dynamic programming|publisher=Academic Press|isbn=978-0-12-221860-6}}.
| |
| *{{citation|last1=Giegerich|first1=R.|last2=Meyer|first2=C.|last3=Steffen|first3=P.|year=2004|url=http://bibiserv.techfak.uni-bielefeld.de/adp/ps/GIE-MEY-STE-2004.pdf|title=A Discipline of Dynamic Programming over Sequence Data|journal=Science of Computer Programming|volume=51|pages=215–263|issue=3|doi=10.1016/j.scico.2003.12.005}}.
| |
| *{{citation|first=Sean|last=Meyn|url= https://netfiles.uiuc.edu/meyn/www/spm_files/CTCN/CTCN.html|title=Control Techniques for Complex Networks|publisher=Cambridge University Press|year=2007|isbn=978-0-521-88441-9}}.
| |
| *S. S. Sritharan, "Dynamic Programming of the Navier-Stokes Equations," in Systems and Control Letters, Vol. 16, No. 4, 1991, pp. 299–307.
| |
| *{{citation|last1=Stokey|first1=Nancy|author1-link=Nancy Stokey|last2=Lucas|first2=Robert E.|author2-link=Robert E. Lucas|last3=Prescott|first3=Edward|author3-link=Edward Prescott|year=1989|title=Recursive Methods in Economic Dynamics|publisher=Harvard Univ. Press|isbn=978-0-674-75096-8}}.
| |
| | |
| ==External links==
| |
| * [http://20bits.com/articles/introduction-to-dynamic-programming/ An Introduction to Dynamic Programming]
| |
| * [http://www.dyna.org Dyna], a declarative programming language for dynamic programming algorithms
| |
| * Wagner, David B., 1995, "[http://citeseer.ist.psu.edu/268391.html Dynamic Programming.]" An introductory article on dynamic programming in [[Mathematica]].
| |
| * [http://www.cse.ohio-state.edu/~gurari/course/cis680/cis680Ch21.html Ohio State University: CIS 680: class notes on dynamic programming], by Eitan M. Gurari
| |
| * [http://mat.gsia.cmu.edu/classes/dynamic/dynamic.html A Tutorial on Dynamic programming]
| |
| * [http://ocw.mit.edu/OcwWeb/Electrical-Engineering-and-Computer-Science/6-046JFall-2005/VideoLectures/detail/embed15.htm MIT course on algorithms] – Includes a video lecture on DP along with lecture notes
| |
| * [http://www.csse.monash.edu.au/~lloyd/tildeAlgDS/Dynamic More DP Notes]
| |
| *King, Ian, 2002 (1987), "[http://researchspace.auckland.ac.nz/bitstream/handle/2292/190/230.pdf A Simple Introduction to Dynamic Programming in Macroeconomic Models.]" An introduction to dynamic programming as an important tool in economic theory.
| |
| * [http://www.topcoder.com/tc?module=Static&d1=tutorials&d2=dynProg Dynamic Programming: from novice to advanced] A TopCoder.com article by Dumitru on Dynamic Programming
| |
| * [http://bibiserv.techfak.uni-bielefeld.de/adp/ Algebraic Dynamic Programming] – a formalized framework for dynamic programming, including an [http://bibiserv.techfak.uni-bielefeld.de/dpcourse entry-level course] to DP, University of Bielefeld
| |
| * Dreyfus, Stuart, "[http://www.eng.tau.ac.il/~ami/cd/or50/1526-5463-2002-50-01-0048.pdf Richard Bellman on the birth of Dynamic Programming.]"
| |
| * [http://www.avatar.se/lectures/molbioinfo2001/dynprog/dynamic.html Dynamic programming tutorial]
| |
| *[http://www.cambridge.org/resources/0521882672/7934_kaeslin_dynpro_new.pdf A Gentle Introduction to Dynamic Programming and the Viterbi Algorithm]
| |
| * Tabled Prolog [http://www.probp.com BProlog] and [http://xsb.sourceforge.net/ XSB]
| |
| * [http://ifors.org/tutorial/category/dynamic-programming/ Online interactive dynamic programming modules] including, shortest path, traveling salesman, knapsack, false coin, egg dropping, bridge and torch, replacement, chained matrix products, and critical path problem.
| |
| | |
| {{optimization algorithms|combinatorial|state=expanded}}
| |
| | |
| {{DEFAULTSORT:Dynamic Programming}}
| |
| [[Category:Dynamic programming| ]]
| |
| [[Category:Mathematical optimization]]
| |
| [[Category:Optimization algorithms and methods]]
| |
| [[Category:Equations]]
| |
| [[Category:Systems engineering]]
| |
| [[Category:Optimal control]]
| |
| [[Category:Operations research]]
| |