Pagini recente » Cod sursa (job #62140) | Cod sursa (job #718923) | Cod sursa (job #1842007) | Cod sursa (job #685753) | Cod sursa (job #587202)
Cod sursa(job #587202)
#include <iostream>
#include <vector>
#include <string>
using namespace std;
#define NM 105
#define MM 5005
#define PB push_back
#define EPS 0.0000000001
double matrix[NM][NM];
int N, M, muchii[MM][3];
vector <int> G[NM];
bool viz[NM];
void df (int nod)
{
viz[nod] = 1;
for (int i = 0; i < G[nod].size(); ++i)
{
int nnod = G[nod][i];
if (viz[nnod]) continue;
df (nnod);
}
}
int check_connectivity()
{
df (1);
for (int i = 1; i <= N; ++i) if (!viz[i]) return 0;
return 1;
}
double ab (double val)
{
if (val < 0) return -val;
return val;
}
int iszero (double val)
{
if (ab(val) < EPS) return 1;
return 0;
}
void swap_rows (int r1, int r2)
{
for (int i = 1; i <= N+1; ++i)
{
double aux = matrix[r1][i];
matrix[r1][i] = matrix[r2][i];
matrix[r2][i] = aux;
}
}
void substract_rows (int r1, int r2, double val)
{
for (int i = 1; i <= N+1; ++i) matrix[r1][i] -= matrix[r2][i] * val;
}
vector <double> solve_system ()
{
vector <double> ans;
for (int i = 1; i < N; ++i)
{
if (iszero(matrix[i][i]))
for (int j = i + 1; j <= N; ++j)
if (!iszero(matrix[j][i]))
{
swap_rows(i, j);
break;
}
for (int j = i + 1; j <= N; ++j) substract_rows (j, i, matrix[j][i]/matrix[i][i]);
}
for (int i = N; i >= 1; --i)
{
for (int j = 1; j < i; ++j) substract_rows(j, i, matrix[j][i]/matrix[i][i]);
ans.push_back(matrix[i][N+1]/matrix[i][i]);
}
int s = 0, e = N - 1;
while (s < e)
{
double aux = ans[s];
ans[s] = ans[e];
ans[e] = aux;
++s, --e;
}
return ans;
}
int main()
{
double X[NM], smax = 2000000000;
memset (matrix, 0, sizeof(matrix));
freopen ("flux.in", "r", stdin);
freopen ("flux.out", "w", stdout);
scanf ("%d %d", &N, &M);
for (int i = 1; i <= M; ++i)
{
int a, b, c;
scanf ("%d %d %d", &a, &b, &c);
G[a].PB(b);
G[b].PB(a);
muchii[i][0] = a;
muchii[i][1] = b;
muchii[i][2] = c;
matrix[a-1][b-1] += 1;
matrix[a-1][a-1] -= 1;
matrix[b-1][a-1] += 1;
matrix[b-1][b-1] -= 1;
}
if (!check_connectivity())
{
printf ("0.000000");
return 0;
}
matrix[N-1][N] = -1;
--N;
vector <double> values = solve_system();
++N;
X[1] = 0;
for (int i = 2; i <= N; ++i) X[i] = values[i-2];
for (int i = 1; i <= M; ++i)
{
int a = muchii[i][0];
int b = muchii[i][1];
int c = muchii[i][2];
double cc = ab(X[b]-X[a]);
smax = min (smax, (double)c/cc);
}
printf ("%lf", smax);
return 0;
}