Cod sursa(job #988219)

Utilizator DaNutZ2UuUUBB Bora Dan DaNutZ2UuU Data 22 august 2013 12:27:51
Problema Drumuri minime Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.37 kb
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<deque>
#include<vector>
#define nmax 1501
#define E 0.0000000001
#define mod 104659
#define inf 2000000000
using namespace std;
int n,i,j,k,m,x,y,rez[nmax];
double dist[nmax];
vector<pair<int,double> >v[nmax];
deque<int>c;
bool viz[nmax];
  
int main()
{
    freopen("dmin.in","r",stdin);
    freopen("dmin.out","w",stdout);
    scanf("%d %d",&n,&m);
    for (i=1;i<=m;i++)
    {
        double cost;
        scanf("%d %d %lf",&x,&y,&cost);
        cost=log(cost);
        v[x].push_back(make_pair(y,cost));
        v[y].push_back(make_pair(x,cost));
    }
    for(i=2;i<=n;i++) dist[i]=inf;
    c.push_back(1);rez[1]=1;viz[1]=true;
    while (!c.empty())
    {
        x=c.front();c.pop_front();viz[x]=false;
        for (j=0;j<v[x].size();j++)
        {
            y=v[x][j].first;
            double cost=v[x][j].second;
            if (E<dist[y]-dist[x]-cost)
            {
                dist[y]=dist[x]+cost;
                rez[y]=rez[x];
                if (!viz[y]) c.push_back(y),viz[y]=true;
            }else
                if (E>fabs(dist[x]+cost-dist[y]))
                {
                    rez[y]+=rez[x];
                    if (rez[y]>=mod) rez[y]-=mod;
                }
        }
    }
    for (i=2;i<=n;i++) printf("%d ",rez[i]);
    return 0;
}