Open In App

Minimum Cost Path with Left, Right, Bottom and Up moves allowed

Improve
Improve
Like Article
Like
Save
Share
Report

Given a two-dimensional grid, each cell of which contains an integer cost which represents a cost to traverse through that cell, we need to find a path from the top left cell to the bottom right cell by which the total cost incurred is minimum.

Note: It is assumed that negative cost cycles do not exist in input matrix.

This problem is an extension of problem: Min Cost Path with right and bottom moves allowed.

In the previous problem only going right and the bottom was allowed but in this problem, we are allowed to go bottom, up, right and left i.e. in all 4 directions.

Examples:

A cost grid is given in below diagram, minimum 
cost to reach bottom right from top left 
is 327 (= 31 + 10 + 13 + 47 + 65 + 12 + 18 + 
6 + 33 + 11 + 20 + 41 + 20)

The chosen least cost path is shown in green.

It is not possible to solve this problem using dynamic programming similar to the previous problem because here current state depends not only on the right and bottom cells but also on the left and upper cells. We solve this problem using dijkstra’s algorithm. Each cell of the grid represents a vertex and neighbor cells adjacent vertices. We do not make an explicit graph from these cells instead we will use the matrix as it is in our Dijkstra’s algorithm. 

In the below code, Dijkstra’s algorithm’s implementation is used. The code implemented below is changed to cope with matrix represented implicit graph. Please also see use of dx and dy arrays in the below code, these arrays are taken for simplifying the process of visiting neighbor vertices of each cell.

Below is the implementation of the above approach:

C++




// C++ program to get least cost path in a grid from
// top-left to bottom-right
 
#include <bits/stdc++.h>
 
using namespace std;
 
#define ROW 5
#define COL 5
 
// structure for information of each cell
struct cell {
    int x, y;
    int distance;
    cell(int x, int y, int distance)
        : x(x)
        , y(y)
        , distance(distance)
    {
    }
};
 
// Utility method for comparing two cells
bool operator<(const cell& a, const cell& b)
{
    if (a.distance == b.distance) {
        if (a.x != b.x)
            return (a.x < b.x);
        else
            return (a.y < b.y);
    }
    return (a.distance < b.distance);
}
 
// Utility method to check whether a point is
// inside the grid or not
bool isInsideGrid(int i, int j)
{
    return (i >= 0 && i < ROW && j >= 0 && j < COL);
}
 
// Method returns minimum cost to reach bottom
// right from top left
int shortest(int grid[ROW][COL], int row, int col)
{
    int dis[row][col];
 
    // initializing distance array by INT_MAX
    for (int i = 0; i < row; i++)
        for (int j = 0; j < col; j++)
            dis[i][j] = INT_MAX;
 
    // direction arrays for simplification of getting
    // neighbour
    int dx[] = { -1, 0, 1, 0 };
    int dy[] = { 0, 1, 0, -1 };
 
    set<cell> st;
 
    // insert (0, 0) cell with 0 distance
    st.insert(cell(0, 0, 0));
 
    // initialize distance of (0, 0) with its grid value
    dis[0][0] = grid[0][0];
 
    // loop for standard dijkstra's algorithm
    while (!st.empty()) {
        // get the cell with minimum distance and delete
        // it from the set
        cell k = *st.begin();
        st.erase(st.begin());
 
        // looping through all neighbours
        for (int i = 0; i < 4; i++) {
            int x = k.x + dx[i];
            int y = k.y + dy[i];
 
            // if not inside boundary, ignore them
            if (!isInsideGrid(x, y))
                continue;
 
            // If distance from current cell is smaller,
            // then update distance of neighbour cell
            if (dis[x][y] > dis[k.x][k.y] + grid[x][y]) {
                // If cell is already there in set, then
                // remove its previous entry
                if (dis[x][y] != INT_MAX)
                    st.erase(
                        st.find(cell(x, y, dis[x][y])));
 
                // update the distance and insert new
                // updated cell in set
                dis[x][y] = dis[k.x][k.y] + grid[x][y];
                st.insert(cell(x, y, dis[x][y]));
            }
        }
    }
 
    // uncomment below code to print distance
    // of each cell from (0, 0)
    /*
    for (int i = 0; i < row; i++, cout << endl)
        for (int j = 0; j < col; j++)
            cout << dis[i][j] << " ";
    */
    // dis[row - 1][col - 1] will represent final
    // distance of bottom right cell from top left cell
    return dis[row - 1][col - 1];
}
 
// Driver code to test above methods
int main()
{
    int grid[ROW][COL]
        = { 31, 100, 65,  12,  18,  10, 13, 47,  157,
            6,  100, 113, 174, 11,  33, 88, 124, 41,
            20, 140, 99,  32,  111, 41, 20 };
 
    cout << shortest(grid, ROW, COL) << endl;
    return 0;
}


Java




// Java program to get least cost path
// in a grid from top-left to bottom-right
import java.io.*;
import java.util.*;
 
class GFG{
     
static int[] dx = { -1, 0, 1, 0 };
static int[] dy = { 0, 1, 0, -1 };
static int ROW = 5;
static int COL = 5;
 
// Custom class for representing
// row-index, column-index &
// distance of each cell
static class Cell
{
    int x;
    int y;
    int distance;
     
    Cell(int x, int y, int distance)
    {
        this.x = x;
        this.y = y;
        this.distance = distance;
    }
}
 
// Custom comparator for inserting cells
// into Priority Queue
static class distanceComparator
  implements Comparator<Cell>
{
    public int compare(Cell a, Cell b)
    {
        if (a.distance < b.distance)
        {
            return -1;
        }
        else if (a.distance > b.distance)
        {
            return 1;
        }
        else {return 0;}
    }
}
 
// Utility method to check whether current
// cell is inside grid or not
static boolean isInsideGrid(int i, int j)
{
    return (i >= 0 && i < ROW &&
            j >= 0 && j < COL);
}
 
// Method to return shortest path from
// top-corner to bottom-corner in 2D grid
static int shortestPath(int[][] grid, int row,
                                      int col)
{
    int[][] dist = new int[row][col];
     
    // Initializing distance array by INT_MAX
    for(int i = 0; i < row; i++)
    {
        for(int j = 0; j < col; j++)
        {
            dist[i][j] = Integer.MAX_VALUE;
        }
    }
     
    // Initialized source distance as
    // initial grid position value
    dist[0][0] = grid[0][0];
     
    PriorityQueue<Cell> pq = new PriorityQueue<Cell>(
                  row * col, new distanceComparator());
                   
    // Insert source cell to priority queue
    pq.add(new Cell(0, 0, dist[0][0]));
    while (!pq.isEmpty())
    {
        Cell curr = pq.poll();
        for(int i = 0; i < 4; i++)
        {
            int rows = curr.x + dx[i];
            int cols = curr.y + dy[i];
             
            if (isInsideGrid(rows, cols))
            {
                if (dist[rows][cols] >
                    dist[curr.x][curr.y] +
                    grid[rows][cols])
                {
                     
                    // If Cell is already been reached once,
                    // remove it from priority queue
                    if (dist[rows][cols] != Integer.MAX_VALUE)
                    {
                        Cell adj = new Cell(rows, cols,
                                       dist[rows][cols]);
                                        
                        pq.remove(adj);
                    }
                     
                    // Insert cell with updated distance
                    dist[rows][cols] = dist[curr.x][curr.y] +
                                       grid[rows][cols];
                                        
                    pq.add(new Cell(rows, cols,
                               dist[rows][cols]));
                }
            }
        }
    }
    return dist[row - 1][col - 1];
}
 
// Driver code
public static void main(String[] args)
throws IOException
{
    int[][] grid = { { 31, 100, 65, 12, 18 },
                     { 10, 13, 47, 157, 6 },
                     { 100, 113, 174, 11, 33 },
                     { 88, 124, 41, 20, 140 },
                     { 99, 32, 111, 41, 20 } };
                      
    System.out.println(shortestPath(grid, ROW, COL));
}
}
 
// This code is contributed by jigyansu


Python3




# Python program to get least cost path in a grid from
# top-left to bottom-right
from functools import cmp_to_key
 
ROW = 5
COL = 5
 
def mycmp(a,b):
     
    if (a.distance == b.distance):
        if (a.x != b.x):
            return (a.x - b.x)
        else:
            return (a.y - b.y)
    return (a.distance - b.distance)
 
# structure for information of each cell
class cell:
 
    def __init__(self,x, y, distance):
        self.x = x
        self.y = y
        self.distance = distance
 
# Utility method to check whether a point is
# inside the grid or not
def isInsideGrid(i, j):
    return (i >= 0 and i < ROW and j >= 0 and j < COL)
 
# Method returns minimum cost to reach bottom
# right from top left
def shortest(grid, row, col):
    dis = [[0 for i in range(col)]for j in range(row)]
 
    # initializing distance array by INT_MAX
    for i in range(row):
        for j in range(col):
            dis[i][j] = 1000000000
 
    # direction arrays for simplification of getting
    # neighbour
    dx = [-1, 0, 1, 0]
    dy = [0, 1, 0, -1]
 
    st = []
 
    # insert (0, 0) cell with 0 distance
    st.append(cell(0, 0, 0))
 
    # initialize distance of (0, 0) with its grid value
    dis[0][0] = grid[0][0]
 
    # loop for standard dijkstra's algorithm
    while (len(st)!=0):
 
        # get the cell with minimum distance and delete
        # it from the set
        k = st[0]
        st = st[1:]
 
        # looping through all neighbours
        for i in range(4):
 
            x = k.x + dx[i]
            y = k.y + dy[i]
 
            # if not inside boundary, ignore them
            if (isInsideGrid(x, y) == 0):
                continue
 
            # If distance from current cell is smaller, then
            # update distance of neighbour cell
            if (dis[x][y] > dis[k.x][k.y] + grid[x][y]):
                # update the distance and insert new updated
                # cell in set
                dis[x][y] = dis[k.x][k.y] + grid[x][y]
                st.append(cell(x, y, dis[x][y]))
 
        st.sort(key=cmp_to_key(mycmp))
 
    # uncomment below code to print distance
    # of each cell from (0, 0)
 
    # for i in range(row):
    #     for j in range(col):
    #         print(dis[i][j] ,end= " ")
    #     print()
 
    # dis[row - 1][col - 1] will represent final
    # distance of bottom right cell from top left cell
    return dis[row - 1][col - 1]
 
# Driver code to test above methods
 
grid = [[31, 100, 65, 12, 18], [10, 13, 47, 157, 6], [100, 113, 174, 11, 33], [88, 124, 41, 20, 140],[99, 32, 111, 41, 20]]
print(shortest(grid, ROW, COL))
 
# This code is contributed by shinjanpatra


C#




using System;
using System.Collections.Generic;
 
class GFG {
    static int[] dx = { -1, 0, 1, 0 };
    static int[] dy = { 0, 1, 0, -1 };
    static int ROW = 5;
    static int COL = 5;
 
    // Custom class for representing
    // row-index, column-index &
    // distance of each cell
    class Cell {
        public int x;
        public int y;
        public int distance;
 
        public Cell(int x, int y, int distance)
        {
            this.x = x;
            this.y = y;
            this.distance = distance;
        }
    }
 
    // Custom comparator for sorting cells
    // in ascending order of distance
    class DistanceComparer : IComparer<Cell> {
        public int Compare(Cell a, Cell b)
        {
            if (a.distance < b.distance) {
                return -1;
            }
            else if (a.distance > b.distance) {
                return 1;
            }
            else {
                return 0;
            }
        }
    }
 
    // Utility method to check whether current
    // cell is inside grid or not
    static bool IsInsideGrid(int i, int j)
    {
        return (i >= 0 && i < ROW && j >= 0 && j < COL);
    }
 
    // Method to return shortest path from
    // top-corner to bottom-corner in 2D grid
    static int ShortestPath(int[][] grid, int row, int col)
    {
        int[][] dist = new int[row][];
        for (int i = 0; i < row; i++) {
            dist[i] = new int[col];
            for (int j = 0; j < col; j++) {
                dist[i][j] = int.MaxValue;
            }
        }
 
        // Initialized source distance as
        // initial grid position value
        dist[0][0] = grid[0][0];
 
        List<Cell> pq = new List<Cell>();
        pq.Add(new Cell(0, 0, dist[0][0]));
        pq.Sort(new DistanceComparer());
 
        while (pq.Count > 0) {
            Cell curr = pq[0];
            pq.RemoveAt(0);
 
            for (int i = 0; i < 4; i++) {
                int rows = curr.x + dx[i];
                int cols = curr.y + dy[i];
 
                if (IsInsideGrid(rows, cols)) {
                    if (dist[rows][cols]
                        > dist[curr.x][curr.y]
                              + grid[rows][cols]) {
                        // If Cell is already been reached
                        // once, remove it from list
                        if (dist[rows][cols]
                            != int.MaxValue) {
                            Cell adj = new Cell(
                                rows, cols,
                                dist[rows][cols]);
                            pq.Remove(adj);
                        }
 
                        // Insert cell with updated distance
                        dist[rows][cols]
                            = dist[curr.x][curr.y]
                              + grid[rows][cols];
                        pq.Add(new Cell(rows, cols,
                                        dist[rows][cols]));
                        pq.Sort(new DistanceComparer());
                    }
                }
            }
        }
 
        return dist[row - 1][col - 1];
    }
 
    // Driver code
    static void Main(string[] args)
    {
        int[][] grid
            = { new int[] { 31, 100, 65, 12, 18 },
                new int[] { 10, 13, 47, 157, 6 },
                new int[] { 100, 113, 174, 11, 33 },
                new int[] { 88, 124, 41, 20, 140 },
                new int[] { 99, 32, 111, 41, 20 } };
 
        Console.WriteLine(ShortestPath(grid, ROW, COL));
    }
}
 
// This code is contributed by phasing17


Javascript




<script>
  
// Javascript program to get least cost path in a grid from
// top-left to bottom-right
var ROW = 5
var COL = 5
 
// structure for information of each cell
class cell
{
    constructor(x, y, distance)
    {
        this.x = x;
        this.y = y;
        this.distance = distance;
    }
};
 
 
// Utility method to check whether a point is
// inside the grid or not
function isInsideGrid(i, j)
{
    return (i >= 0 && i < ROW && j >= 0 && j < COL);
}
 
// Method returns minimum cost to reach bottom
// right from top left
function shortest(grid, row, col)
{
    var dis = Array.from(Array(row), ()=>Array(col).fill(0));
 
    // initializing distance array by INT_MAX
    for (var i = 0; i < row; i++)
        for (var j = 0; j < col; j++)
            dis[i][j] = 1000000000;
 
    // direction arrays for simplification of getting
    // neighbour
    var dx = [-1, 0, 1, 0];
    var dy = [0, 1, 0, -1];
 
    var st = [];
 
    // insert (0, 0) cell with 0 distance
    st.push(new cell(0, 0, 0));
 
    // initialize distance of (0, 0) with its grid value
    dis[0][0] = grid[0][0];
 
    // loop for standard dijkstra's algorithm
    while (st.length!=0)
    {
        // get the cell with minimum distance and delete
        // it from the set
        var k = st[0];
        st.shift();
 
        // looping through all neighbours
        for (var i = 0; i < 4; i++)
        {
            var x = k.x + dx[i];
            var y = k.y + dy[i];
 
            // if not inside boundary, ignore them
            if (!isInsideGrid(x, y))
                continue;
 
            // If distance from current cell is smaller, then
            // update distance of neighbour cell
            if (dis[x][y] > dis[k.x][k.y] + grid[x][y])
            {
                // update the distance and insert new updated
                // cell in set
                dis[x][y] = dis[k.x][k.y] + grid[x][y];
                st.push(new cell(x, y, dis[x][y]));
            }
        }
        st.sort((a,b)=>{
            if (a.distance == b.distance)
    {
        if (a.x != b.x)
            return (a.x - b.x);
        else
            return (a.y - b.y);
    }
    return (a.distance - b.distance);
        });
    }
 
    // uncomment below code to print distance
    // of each cell from (0, 0)
    /*
    for (int i = 0; i < row; i++, cout << endl)
        for (int j = 0; j < col; j++)
            cout << dis[i][j] << " ";
    */
    // dis[row - 1][col - 1] will represent final
    // distance of bottom right cell from top left cell
    return dis[row - 1][col - 1];
}
 
// Driver code to test above methods
var grid =
[
    [31, 100, 65, 12, 18],
    [10, 13, 47, 157, 6],
    [100, 113, 174, 11, 33],
    [88, 124, 41, 20, 140],
    [99, 32, 111, 41, 20]
];
document.write(shortest(grid, ROW, COL));
 
// This code is contributed by rutvik_56.
 
</script>


Output

327

Time Complexity: O(N2 log N), The Dijkstra’s algorithm used in the program has a time complexity of O(E log V) where E is the number of edges and V is the number of vertices. In this program, for each cell, we check its four neighbors, so the number of edges E in the worst case would be 4 times the number of cells (N^2) in the grid. Therefore, the time complexity of the program is O(4N^2logN^2), which simplifies to O(N^2 logN).
Auxiliary Space: O(N2), In this program, we are using a 2D array to store the distances from the top-left cell to each cell in the grid. The size of this array is N^2. We are also using a set to store the cells with their distances. In the worst case, all cells in the grid could be present in the set, so the size of the set could also be N^2. Therefore, the space complexity of the program is O(N^2).

 



Last Updated : 08 May, 2023
Like Article
Save Article
Previous
Next
Share your thoughts in the comments
Similar Reads