Cod sursa(job #1198997)

Utilizator gabrielinelusGabriel-Robert Inelus gabrielinelus Data 17 iunie 2014 21:09:07
Problema Drumuri minime Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.16 kb
#include <cstdio>
#include <queue>
#include <vector>
#include <algorithm>
#include <cstring>
#include <bitset>
#include <cmath>

#define Nmax 1605
#define MOD 104659
#define INF 0x3f3f3f3f
#define EPS1 0.0001
#define EPS2 0.00000001

using namespace std;

vector<pair<double,int> > G[Nmax];
vector<int> daddy[Nmax];
queue<int> Q;
bitset<Nmax> inQ;

int N,M;
double dist[Nmax];
int DP[Nmax];

void read()
{
    int a,b,c;
    scanf("%d%d",&N,&M);
    for(int i = 1; i <= M; ++i){
        scanf("%d%d%d",&a,&b,&c);
        G[a].push_back(make_pair(log2((double)c),b));
        G[b].push_back(make_pair(log2((double)c),a));
    }
}

void bellman_ford()
{
    Q.push(1);
    for(int i = 2; i <= N; ++i) dist[i] = INF;
    dist[1] = 0;
    int k;
    while(!Q.empty())
    {
        k = Q.front();Q.pop();inQ[k] = 0;
        for(vector<pair<double,int> >::iterator it = G[k].begin(); it != G[k].end(); ++it)
        {
            if(dist[it->second] > dist[k] + it->first + EPS2)
            {
                dist[it->second] = dist[k] + it->first;
                daddy[it->second].clear();
                daddy[it->second].push_back(k);
                inQ[it->second] = 1;
                Q.push(it->second);
                continue;
            }
            if(dist[it->second] - dist[k] - it->first + EPS2  >= 0 &&
               dist[it->second] - dist[k] - it->first  <= EPS1 )
                daddy[it->second].push_back(k);
        }
    }
}

int dynamic(int k)
{
    if(DP[k])return DP[k];
    int sum = 0;
    for(vector<int>::iterator it = daddy[k].begin(); it != daddy[k].end(); ++it)
        sum = (sum%MOD + dynamic(*it)%MOD) % MOD;
    DP[k] = sum;
    return DP[k] % MOD;
}
void prepare()
{
    for(int i = 1; i <= N; ++i)
    {
        vector<int>::iterator it = unique(daddy[i].begin(),daddy[i].end());
        daddy[i].resize(distance(daddy[i].begin(),it));
    }
}


int main()
{
    freopen("dmin.in","r",stdin);
    freopen("dmin.out","w",stdout);

    read();
    bellman_ford();
    DP[1] = 1;
    prepare();
    for(int i = 2; i <= N; ++i)
        printf("%d ",dynamic(i));

    return 0;
}