Pagini recente » Cod sursa (job #283868) | Cod sursa (job #2153766) | Cod sursa (job #230068) | Cod sursa (job #1012283) | Cod sursa (job #587659)
Cod sursa(job #587659)
#include <cstdio>
#include <algorithm>
#include <vector>
using namespace std;
const int N = 260;
int n;
double g[N][N];
vector <pair <int, int> > L[N];
void read() {
int m, x, y, c;
scanf("%d%d", &n, &m);
for (int i = 1; i <= m; ++ i) {
scanf("%d%d%d", &x, &y, &c);
L[x].push_back(make_pair(y, c));
L[y].push_back(make_pair(x, c));
}
}
void init_gauss() {
for (int i = 1; i < n; ++ i) {
g[i][i] = 1;
double c = (double)1 / L[i].size();
for (int j = 0; j < L[i].size(); ++ j) {
g[i][L[i][j].first] = - c;
g[i][0] += c * L[i][j].second;
}
}
g[n][n] = 1;
g[n][0] = 0;
}
void gauss() {
for (int j = 1; j <= n; ++ j) {
int i = 0;
for (i = j; i <= n; ++ i)
if (g[i][j])
break;
if (i != j)
for (int k = 0; k <= n; ++ k)
swap(g[i][k], g[j][k]);
for (i = 1; i <= n; ++ i)
if (i != j) {
double deinm = - g[i][j] / g[j][j];
for (int k = 1; k <= n; ++ k)
g[i][k] += deinm * g[j][k];
g[i][0] += deinm * g[j][0];
}
}
printf("%.4lf\n", g[1][0] / g[1][1]);
}
int main() {
freopen("tunel.in", "r", stdin);
freopen("tunel.out", "w", stdout);
read();
init_gauss();
gauss();
return 0;
}