Cod sursa(job #2130653)

Utilizator Claudiu07Pana Claudiu Claudiu07 Data 13 februarie 2018 20:03:03
Problema Tunelul groazei Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.95 kb
#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;
}