Pagini recente » Cod sursa (job #2692150) | Cod sursa (job #1625049) | Cod sursa (job #477246) | Cod sursa (job #2322113) | Cod sursa (job #2719341)
#include <bits/stdc++.h>
#define dbg() cerr <<
#define name(x) (#x) << ": " << (x) << ' ' <<
// nu stiu ce sa i mai fac :(
using namespace std;
const long double EPS = 1e-3;
const int MAGIC = 1;
int Sgn(long double x) {
if (fabs(x) <= EPS) return 0;
return x < 0 ? -1 : +1;
}
vector<long double> Gauss(vector<vector<long double>> s) {
int n = s.size(), m = s[0].size();
int i = 0, j = 0;
while (i < n && j < m - 1) {
int k = i; while(k < n && Sgn(s[k][j]) == 0) ++k;
if (k == n) {
++j;
continue;
}
if (i != k) swap(s[i], s[k]);
for (k = 0; k < m; ++k) if (k != j) s[i][k] /= s[i][j];
s[i][j] = 1;
for (k = 0; k < n; ++k) if (i != k) {
auto coef = s[k][j];
if (Sgn(coef) != 0) {
for (int c = 0; c < m; ++c)
s[k][c] -= coef * s[i][c];
s[k][j] = 0;
}
}
++i, ++j;
}
vector<long double> x(n, -1);
for (i = 0; i < n; ++i) {
int pos = -1;
for (j = 0; j < m; ++j) {
if (Sgn(s[i][j]) != 0) {
pos = j;
break;
}
}
if (pos == m - 1) {
cout << "0\n";
exit(0);
}
if (pos != -1) x[pos] = s[i][m - 1];
}
return x;
}
int main() {
freopen("flux.in", "r", stdin);
freopen("flux.out", "w", stdout);
int n, m; cin >> n >> m;
vector<tuple<int, int, int>> edges(m);
vector<vector<long double>> s(n, vector<long double>(n + 1));
for (int i = 0; i < m; ++i) {
int a, b, c; cin >> a >> b >> c; --a, --b;
if (a > b) swap(a, b);
edges[i] = {a, b, c};
if (c == 0) continue;
if (a != 0) {
s[a][b] += 1; s[a][a] -= 1;
}
s[b][a] += 1; s[b][b] -= 1;
}
s[0][0] = 1; s[0][n] = 0; // dist[0] = 0
s[n - 1][n] = -MAGIC;
for (int i = 0; i < n; ++i) {
for (int j = 0; j <= n; ++j) {
cerr << setw(2) << s[i][j] << ' ';
}
cerr << endl;
}
auto dist = Gauss(move(s));
for (auto x : dist) dbg() name(x) endl;
long double ans = numeric_limits<long double>::max();
for (auto &e : edges) {
int a, b, c; tie(a, b, c) = e;
if (dist[a] < 0 || dist[b] < 0) continue;
long double flow = abs(dist[a] - dist[b]);
if (Sgn(flow) != 0) ans = min(ans, c / flow);
}
cout << fixed << setprecision(10) << ans * MAGIC << endl;
}