Cod sursa(job #920972)

Utilizator razvan.popaPopa Razvan razvan.popa Data 20 martie 2013 18:26:02
Problema Drumuri minime Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.57 kb
#include<cstdio>
#include<cmath>
#include<vector>
#include<queue>
#include<algorithm>
#define INF (1 << 30)
#define pb push_back
#define mkp make_pair
#define nxt (*it).first
#define cst (*it).second
#define FOR(i,a,b)\
   for(int i=a; i<=b; ++i)
#define ALL(g)\
   for(typeof(g.begin()) it=g.begin(); it!=g.end(); ++it)
#define infile "dmin.in"
#define outfile "dmin.out"
#define nMax 1505
#define MOD 104659
#define EPS 0.000000001
using namespace std;

priority_queue < pair < double, int > > PQ;

vector < pair < int, double > > v[nMax];

double D[nMax];

int NR[nMax];

int N, M;

void read(){
    freopen(infile, "r", stdin);

    scanf("%d %d", &N, &M);

    int x, y, c;
    while(M--){
        scanf("%d %d %d", &x, &y, &c);

        v[x].pb(mkp(y,log(c)));
        v[y].pb(mkp(x,log(c)));
    }

    fclose(stdin);
}


void Dijkstra(){
    FOR(i,2,N)
       D[i] = INF;

    NR[1] = 1;
    PQ.push(mkp(0.0, 1));

    while(!PQ.empty()){
        double dx = -PQ.top().first;
        int x = PQ.top().second;
        PQ.pop();

        if(dx > D[x])
           continue;

        ALL(v[x])
           if(dx + cst < D[nxt]){
               D[nxt] = dx + cst;
               NR[nxt] = NR[x];

               PQ.push(mkp(-D[nxt], nxt));
           }
           else if(abs(D[nxt] - dx - cst) < EPS)
               NR[nxt] = (0LL + NR[nxt] + NR[x]) % MOD;
    }
}

void print(){
    freopen(outfile, "w", stdout);

    FOR(i,2,N)
       printf("%d ", NR[i]);

    fclose(stdout);
}

int main(){
    read();
    Dijkstra();
    print();

    return 0;
}