Cod sursa(job #639116)

Utilizator swift90Ionut Bogdanescu swift90 Data 22 noiembrie 2011 15:21:51
Problema Tunelul groazei Scor 50
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.1 kb
#include<cstdio>
#include<fstream>
#include<vector>
using namespace std;
typedef pair<int,int> PII;
PII nr[260][260];
int N,M,grad[260];
double G[260][260];
void init(){
	for(int i=1;i<N;++i){
		G[i][i]=grad[i];
		for(int it=1;it<=N;++it){
			G[i][it]-=nr[i][it].first;
			G[i][N+1]+=nr[i][it].second;
		}
	}
	G[N][N]=1;
	G[N][N+1]=0;
}
void solve(){
	init();
	int i,j,k;
	for(i=1;i<=N;++i){
		for(j=i+1;j<=N;++j){
			for(k=N+1;k>=i;--k){
				if(G[j][i]){
					G[j][k]-=G[j][i]/G[i][i]*G[i][k];
				}
			}
		}
	}
	for(j=N;j;--j){
		G[j][N+1]/=G[j][j];
		G[j][j]=1;
		for(i=j-1;i;--i){
			G[i][N+1]-=G[i][j]*G[j][N+1];
			G[i][j]=0;
		}
	}
}
int main(){
	//freopen("tunel.in","r",stdin);
	ifstream f("tunel.in");
	freopen("tunel.out","w",stdout);
	int i,x,y,c;
	//scanf("%d%d",&N,&M);
	f>>N>>M;
	for(i=0;i<M;++i){
		//scanf("%d%d%d",&x,&y,&c);
		f>>x>>y>>c;
		++grad[x];
		++grad[y];
		++nr[x][y].first;
		nr[x][y].second+=c;
		++nr[y][x].first;
		nr[y][x].second+=c;
	}
	solve();
	printf("%.4lf\n",G[1][N+1]);
	
	fclose(stdin);
	fclose(stdout);
	return 0;
}