Pagini recente » Cod sursa (job #1204210) | Cod sursa (job #708792) | Cod sursa (job #91429) | Cod sursa (job #2106201) | Cod sursa (job #2130653)
#include <iostream>
#include <fstream>
#include <iomanip>
#include <cmath>
using namespace std;
ifstream f("tunel.in");
ofstream g("tunel.out");
const double eps=1E-8;
double m[256][257],x[257];
int N,M;
void gauss()
{
for(int i=1; i<N; i++)
{
for(int j=i+1; j<=N+1; j++)
m[i][j]/=m[i][i];
m[i][i]=1;
for(int j=i+1; j<N; j++)
{
for(int k=i+1; k<=N+1; k++)
m[j][k]-=m[i][k]*m[j][i];
m[j][i]=0;
}
}
}
void sol()
{
for(int i=N-1; i>0; i--)
{
x[i]=m[i][N+1];
for(int j=i+1; j<N; j++)
x[i]-=m[i][j]*x[j];
}
}
int main()
{
int a,b,c;
f>>N>>M;
for(int i=1; i<=M; i++)
{
f>>a>>b>>c;
m[a][a]++; m[a][b]--; m[a][N+1]+=c;
m[b][b]++; m[b][a]--; m[b][N+1]+=c;
}
gauss();
sol();
g<<fixed<<setprecision(3)<<x[1]<<'\n';
return 0;
}