Cod sursa(job #1880403)

Utilizator FloareaCucura Floarea Ludovica Floarea Data 15 februarie 2017 18:52:46
Problema Drumuri minime Scor 100
Compilator cpp Status done
Runda becreative17 Marime 1.59 kb
#include<fstream>
#include<vector>
#include<cmath>
#include<queue>
#define MOD 104659
#define EPS 0.0000000001
#define INF 20000000000.0
using namespace std;
ifstream fin("dmin.in");
ofstream fout("dmin.out");
int n, m, a, b, c;
vector< pair<int, double> > v[1505];
double d[1505];
int cnt[1505];
priority_queue< pair< double, int >, vector< pair< double, int > >, greater< pair< double, int > > > q;

double modul( double x ){
    if( x > 0 ){
        return x;
    }
    return -x;
}

int main(){
    fin >> n >> m;
    for( int i = 1; i <= m; i++ ){
        fin >> a >> b >> c;
        v[a].push_back( make_pair( b, log(c) ) );
        v[b].push_back( make_pair( a, log(c) ) );
    }
    for( int i = 2; i <= n; i++ ){
        d[i] = INF;
    }
    cnt[1] = 1;
    q.push( make_pair( 0, 1 ) );
    while( !q.empty() ){
        int nod = q.top().second;

        if( q.top().first - d[nod] >= EPS ){
            q.pop();
            continue;
        }

        for( int i = 0; i < v[nod].size(); i++ ){
            int vecin = v[nod][i].first;
            double cost = v[nod][i].second + d[nod];
            if( modul( d[vecin] - cost ) <= EPS ){
                cnt[vecin] = ( cnt[vecin] + cnt[nod] ) % MOD;
            }else{
                if( cost < d[vecin] ){
                    d[vecin] = cost;
                    cnt[vecin] = cnt[nod];
                    q.push( make_pair( d[vecin], vecin ) );
                }
            }
        }

        q.pop();
    }
    for( int i = 2; i <= n; i++ ){
        fout << cnt[i] << " ";
    }
    return 0;
}