Cod sursa(job #421036)

Utilizator ooctavTuchila Octavian ooctav Data 20 martie 2010 23:04:08
Problema Algoritmul Bellman-Ford Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.97 kb
#include <stdio.h>
#include <vector>
#include <set>
#include <algorithm>

using namespace std;

#define maxn 50010
#define maxm 250010
#define inf 1 << 30

struct muchie
{
    long a, b, c;
} e[maxm];//capetele si costul fiecarei muchii

long n, m, i, j, k, left, right, cost[maxn], f[maxn];
vector <long> v[maxn];
vector <pair<long, long> > h;

void citire()
{
	scanf("%d%d", &n, &m);
    for(i = 1; i <= m; i++)
    {
        scanf("%d%d%d", &e[i].a, &e[i].b, &e[i].c);
        v[e[i].a].push_back(i);
    }
}

void init()
{
    cost[1] = 0;
    for(int i = 2; i <= n; i++)
        cost[i] = inf;
	
    h.push_back(make_pair(0, 1));
    make_heap(h.begin(), h.end());
}

void solve()
{
    pair<long, long> per;
    long i, j, nod, poz;
    long long pas=0;
    while(h.size() && pas<=1LL*n*m)
    {
        pas++;
        pop_heap(h.begin(), h.end());
        per = h.back();
        h.pop_back();
        nod = per.second;
        f[nod] = 0;
        for(j = 0 ; j < v[nod].size() ; j++)
        {
            poz = v[nod][j];
            if(cost[e[poz].a] + e[poz].c<cost[e[poz].b])
            {
                cost[e[poz].b] = cost[e[poz].a] + e[poz].c;
                if(!f[e[poz].b])
                {
                    f[e[poz].b]=1;
                    h.push_back(make_pair(-cost[e[poz].b], e[poz].b));
                    push_heap(h.begin(), h.end());
                }
            }
        }
    }
}   

long check_negativ()
{
    long i;
    for(i = 1 ; i <= m ; i++)
        if(cost[e[i].a] + e[i].c < cost[e[i].b])
            return 1;
    return 0;
}         

void scrie()
{
	if(check_negativ())
    {
        printf("Ciclu negativ!\n");
        return ;
    }
    for(i = 2; i <= n; i++)
        printf("%d ", cost[i]);
}

int main()
{
    freopen("bellmanford.in", "r", stdin);
    freopen("bellmanford.out", "w", stdout);
	citire();
    init();
    solve();
	scrie();
    return 0;
}