Given a cost matrix cost[][] and a position (m, n) in cost[][], write a function that returns cost of minimum cost path to reach (m, n) from (0, 0). Each cell of the matrix represents a cost to traverse through that cell. Total cost of a path to reach (m, n) is sum of all the costs on that path (including both source and destination). You can only traverse down, right and diagonally lower cells from a given cell, i.e., from a given cell (i, j), cells (i+1, j), (i, j+1) and (i+1, j+1) can be traversed. You may assume that all costs are positive integers.
For example, in the following figure, what is the minimum cost path to (2, 2)?
The path with minimum cost is highlighted in the following figure. The path is (0, 0) –> (0, 1) –> (1, 2) –> (2, 2). The cost of the path is 8 (1 + 2 + 2 + 3).
1) Optimal Substructure
The path to reach (m, n) must be through one of the 3 cells: (m-1, n-1) or (m-1, n) or (m, n-1). So minimum cost to reach (m, n) can be written as “minimum of the 3 cells plus cost[m][n]”.
minCost(m, n) = min (minCost(m-1, n-1), minCost(m-1, n), minCost(m, n-1)) + cost[m][n]
2) Overlapping Subproblems
Following is simple recursive implementation of the MCP (Minimum Cost Path) problem. The implementation simply follows the recursive structure mentioned above.
/* A Naive recursive implementation of MCP(Minimum Cost Path) problem */
#include<stdio.h>
#include<limits.h>
#define R 3
#define C 3
int
min(
int
x,
int
y,
int
z);
/* Returns cost of minimum cost path from (0,0) to (m, n) in mat[R][C]*/
int
minCost(
int
cost[R][C],
int
m,
int
n)
{
if
(n < 0 || m < 0)
return
INT_MAX;
else
if
(m == 0 && n == 0)
return
cost[m][n];
else
return
cost[m][n] + min( minCost(cost, m-1, n-1),
minCost(cost, m-1, n),
minCost(cost, m, n-1) );
}
/* A utility function that returns minimum of 3 integers */
int
min(
int
x,
int
y,
int
z)
{
if
(x < y)
return
(x < z)? x : z;
else
return
(y < z)? y : z;
}
/* Driver program to test above functions */
int
main()
{
int
cost[R][C] = { {1, 2, 3},
{4, 8, 2},
{1, 5, 3} };
printf
(
" %d "
, minCost(cost, 2, 2));
return
0;
}
|
It should be noted that the above function computes the same subproblems again and again. See the following recursion tree, there are many nodes which apear more than once. Time complexity of this naive recursive solution is exponential and it is terribly slow.
mC refers to minCost() mC(2, 2) / | \ / | \ mC(1, 1) mC(1, 2) mC(2, 1) / | \ / | \ / | \ / | \ / | \ / | \ mC(0,0) mC(0,1) mC(1,0) mC(0,1) mC(0,2) mC(1,1) mC(1,0) mC(1,1) mC(2,0)
So the MCP problem has both properties (see this and this) of a dynamic programming problem. Like other typical Dynamic Programming(DP) problems, recomputations of same subproblems can be avoided by constructing a temporary array tc[][] in bottom up manner.
#include <stdio.h>
#include <stdlib.h>
int three_min(int a,int b,int c)
{
int max = a;
if(b > max)
max = b;
if(c > max)
max = c;
return max;
}
int min_cost_path(int** cost, int m,int n)
{
int **dp = new int*[m];
for(int i = 0; i < m;i++)
dp[i] = new int[n];
dp[0][0] = cost[0][0];
for(int i = 0; i < m; i++)
dp[i][0] += cost[i][0];
for(int j = 0; j < n; j++)
dp[0][j] += cost[0][j];
for(int i = 1; i < m; i++)
{
for(int j = 1; j < n; j++)
dp[i][j] = three_min(dp[i-1][j],dp[i][j-1],dp[i-1][j-1]) + cost[i][j];
}
int result = dp[m-1][n-1];
for(int i = 0; i < m;i++)
delete [] dp[i];
delete [] dp;
return result;
}
int main()
{
int cost[3][3] = {
{1,2,3},
{4,8,2},
{1,5,3}
};
printf(" %d \n", min_cost_path((int **)cost,3,3));
return 0;
}