Cod sursa(job #661336)

Utilizator blustudioPaul Herman blustudio Data 14 ianuarie 2012 12:32:25
Problema Algoritmul lui Dijkstra Scor 60
Compilator cpp Status done
Runda Arhiva educationala Marime 2.49 kb
/*
 * Autor: Paul Herman
 * Problema: Dijkstra cu heap-uri
 * Data: 14.01.2012
 * Punctaj: -
 * Link: http://www.infoarena.ro/problema/dijkstra
 */
#include <fstream>
#include <vector>
#include <iostream>
using namespace std;

struct muchie
{
	int cost; //Costul pentru a ajunge de la A la B
	int b; //Nodul B
};

#define oo 0x3f3f3f

vector<muchie> vecini[50001]; //Nodurile grafului
int distanta[50001];
int n; //Numarul de noduri
int m; //Nr. de muchii
int heap[50001]; //Heapul ce pastreaza distanta nodurilor fata de sursa
int heap_size;

inline void hswap(int a, int b)
{
	int t = heap[a];
	heap[a] = heap[b];
	heap[b] = t;
}
inline void heap_push(int x)
{
	int pozitie = heap_size;
	int tata = (pozitie - 1) / 2;
	heap[pozitie] = x;
	heap_size++;
	while ((tata != pozitie) && (distanta[heap[tata]] > distanta[x]))
	{
		hswap(tata, pozitie);
		tata = (pozitie - 1) / 2;
	}
}
inline void heap_pop()
{
	heap_size--;
	swap(heap[0], heap[heap_size]);
	int modificat = 0;
	int fius, fiud, mic;
	do {
		fius = 2 * modificat + 1;
		fiud = fius + 1;
		mic = modificat;
		if ((fius < heap_size) && (distanta[heap[fius]] < distanta[heap[mic]]))
			mic = fius;
		if ((fiud < heap_size) && (distanta[heap[fiud]] < distanta[heap[mic]]))
			mic = fiud;
		if (modificat != mic)
			hswap(mic, modificat);
	} while (mic != modificat);
}
inline bool heap_empty()
{
	return heap_size == 0;
}
inline int heap_top()
{
	return heap[0];
}
inline void citire()
{
	int a;
	muchie b;
	ifstream fin("dijkstra.in");
	fin >> n >> m;
	for (int i = 1; i <= m ; i++)
	{
		fin >> a >> b.b >> b.cost;
		vecini[a].push_back(b);
	}
	fin.close();
}
inline void scriere()
{
	ofstream fout("dijkstra.out");
	for (int i = 2; i <= n; i++)
		if (distanta[i] < oo)
			fout << distanta[i] << " ";
		else
			fout << "0 ";
	fout.close();
}
inline void Dijkstra()
{
	int nod_curent = 1;
	distanta[1] = 0; //Costul pt. a ajunge la sursa e 0
	for (int i = 2; i <= n; i++)
		distanta[i] = oo; //Costul pt. a ajunge la restul nodurilor e infinit
	heap_push(1);
	while (heap_empty() == false)
	{
		//Cautam nodul nevizitat cel mai apropiat de sursa
		nod_curent = heap_top();
		heap_pop();
		for (int i = 0; i < vecini[nod_curent].size(); i++) //Parcurgem toti vecinii nodului curent
		{
			if ((distanta[nod_curent] + vecini[nod_curent][i].cost) < distanta[vecini[nod_curent][i].b])
			{
				distanta[vecini[nod_curent][i].b] = distanta[nod_curent] + vecini[nod_curent][i].cost;
				heap_push(vecini[nod_curent][i].b);
			}
		}		
	}
}
int main()
{
	citire();
	Dijkstra();
	scriere();
	return 0;
}