Pagini recente » Cod sursa (job #1866407) | Cod sursa (job #2923427) | Cod sursa (job #699120) | Cod sursa (job #2559020) | Cod sursa (job #1198997)
#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;
}