Classic Dijkstra

Dijkstra invented one of the most classic algorithms to find the min/max path in a graph - we call it today simply Dijkstra's Algorithm.

This problem requires a classic use of Dijkstra's Algorithm: https://leetcode.com/problems/path-with-maximum-probability/

1514. Path with Maximum Probability
Medium

You are given an undirected weighted graph of n nodes (0-indexed), represented by an edge list where edges[i] = [a, b] is an undirected edge connecting the nodes a and b with a probability of success of traversing that edge succProb[i].

Given two nodes start and end, find the path with the maximum probability of success to go from start to end and return its success probability.

If there is no path from start to end, return 0. Your answer will be accepted if it differs from the correct answer by at most 1e-5.

Constraints:

  • 2 <= n <= 10^4
  • 0 <= start, end < n
  • start != end
  • 0 <= a, b < n
  • a != b
  • 0 <= succProb.length == edges.length <= 2*10^4
  • 0 <= succProb[i] <= 1
  • There is at most one edge between every two nodes.
Accepted
6,747
Submissions
20,137

In a nutshell, Dijkstra's can be explained in the following pseudo-code:

A) Build or use a Priority Queue
B) For the problem here, the priority will be a double and it will be in descend order
C) Enqueue the start index, with a probability 1
D) While the queue isn't empty
E) Dequeue, process (check for the end) and mark as visited
F) Get the neighbors
G) For each one, check if it isn't visited and add to the queue with the proper probability
H) Notice that on (G), you do not mark as visited, meaning that the queue can contain multiple instances of the same node

Works well, code is below, cheers, ACC.



public class Solution
{
public double MaxProbability(int n, int[][] edges, double[] succProb, int start, int end)
{
Hashtable graph = new Hashtable();

for (int i = 0; i < edges.Length; i++)
{
if (!graph.ContainsKey(edges[i][0])) graph.Add(edges[i][0], new Hashtable());
Hashtable map = (Hashtable)graph[edges[i][0]];
if (!map.ContainsKey(edges[i][1])) map.Add(edges[i][1], succProb[i]);

if (!graph.ContainsKey(edges[i][1])) graph.Add(edges[i][1], new Hashtable());
map = (Hashtable)graph[edges[i][1]];
if (!map.ContainsKey(edges[i][0])) map.Add(edges[i][0], succProb[i]);
}

if (!graph.ContainsKey(start)) return 0;

Hashtable visited = new Hashtable();
PriorityQueue pQueue = new PriorityQueue();
pQueue.Enqueue(start, 1.0);

while (pQueue.Count > 0)
{
double probability = 0;
int index = (int)pQueue.Dequeue(out probability);

if (index == end) return probability;
if (!visited.Contains(index)) visited.Add(index, true);

if (graph.ContainsKey(index))
{
Hashtable map = (Hashtable)graph[index];
foreach (int neighbor in map.Keys)
{
if (!visited.ContainsKey(neighbor))
{
pQueue.Enqueue(neighbor, probability * (double)map[neighbor]);
}
}
}
}

return 0;
}

/*
public static void Main()
{
Solution sol = new Solution();
}
*/
}


public class PriorityQueue
{
public struct HeapEntry
{
private object item;
private double priority;
public HeapEntry(object item, double priority)
{
this.item = item;
this.priority = priority;
}
public object Item
{
get
{
return item;
}
}
public double Priority
{
get
{
return priority;
}
}
}

private int count;
private int capacity;
private HeapEntry[] heap;

public int Count
{
get
{
return this.count;
}
}

public PriorityQueue()
{
capacity = 1000000;
heap = new HeapEntry[capacity];
}

public object Dequeue(out double priority)
{
priority = heap[0].Priority;
object result = heap[0].Item;
count--;
trickleDown(0, heap[count]);
return result;
}

public void Enqueue(object item, double priority)
{
count++;
bubbleUp(count - 1, new HeapEntry(item, priority));
}

private void bubbleUp(int index, HeapEntry he)
{
int parent = (index - 1) / 2;
// note: (index > 0) means there is a parent
while ((index > 0) && (heap[parent].Priority < he.Priority))
{
heap[index] = heap[parent];
index = parent;
parent = (index - 1) / 2;
}
heap[index] = he;
}

private void trickleDown(int index, HeapEntry he)
{
int child = (index * 2) + 1;
while (child < count)
{
if (((child + 1) < count) &&
(heap[child].Priority < heap[child + 1].Priority))
{
child++;
}
heap[index] = heap[child];
index = child;
child = (index * 2) + 1;
}
bubbleUp(index, he);
}
}

Comments

Post a Comment

[フレーム]

Popular posts from this blog

Quasi FSM (Finite State Machine) problem + Vibe

Not really an FSM problem since the state isn't changing, it is just defined by the current input. Simply following the instructions should do it. Using VSCode IDE you can also engage the help of Cline or Copilot for a combo of coding and vibe coding, see below screenshot. Cheers, ACC. Process String with Special Operations I - LeetCode You are given a string  s  consisting of lowercase English letters and the special characters:  * ,  # , and  % . Build a new string  result  by processing  s  according to the following rules from left to right: If the letter is a  lowercase  English letter append it to  result . A  '*'   removes  the last character from  result , if it exists. A  '#'   duplicates  the current  result  and  appends  it to itself. A  '%'   reverses  the current  result . Return the final string  result  after processing all char...

Shortest Bridge – A BFS Story (with a Twist)

Here's another one from the Google 30 Days challenge on LeetCode — 934. Shortest Bridge . The goal? Given a 2D binary grid where two islands (groups of 1s) are separated by water (0s), flip the fewest number of 0s to 1s to connect them. Easy to describe. Sneaky to implement well. 🧭 My Approach My solution follows a two-phase Breadth-First Search (BFS) strategy: Find and mark one island : I start by scanning the grid until I find the first 1 , then use BFS to mark all connected land cells as 2 . I store their positions for later use. Bridge-building BFS : For each cell in the marked island, I run a BFS looking for the second island. Each BFS stops as soon as it hits a cell with value 1 . The minimum distance across all these searches gives the shortest bridge. πŸ” Code Snippet Here's the core logic simplified: public int ShortestBridge(int[][] grid) { // 1. Mark one island as '2' and gather its coordinates List<int> island = FindAndMark...

Classic Dynamic Programming IX

A bit of vibe code together with OpenAI O3. I asked O3 to just generate the sieve due to laziness. Sieve is used to calculate the first M primes (when I was using Miller-Rabin, was giving me TLE). The DP follows from that in a straightforward way: calculate the numbers from i..n-1, then n follows by calculating the min over all M primes. Notice that I made use of Goldbach's Conjecture as a way to optimize the code too. Goldbach's Conjecture estates that any even number greater than 2 is the sum of 2 primes. The conjecture is applied in the highlighted line. Cheers, ACC. PS: the prompt for the sieve was the following, again using Open AI O3 Advanced Reasoning: " give me a sieve to find the first M prime numbers in C#. The code should produce a List<int> with the first M primes " Minimum Number of Primes to Sum to Target - LeetCode You are given two integers  n  and  m . You have to select a multiset of  prime numbers  from the  first   m  pri...