Cod sursa(job #2920986)

Utilizator Savu_Stefan_CatalinSavu Stefan Catalin Savu_Stefan_Catalin Data 26 august 2022 20:35:46
Problema Tunelul groazei Scor 80
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.23 kb
#include <fstream>
#include <iomanip>
using namespace std;
ifstream in("tunel.in");
ofstream out("tunel.out");
double ep,a[260][260],z,ras,ma2[260],ma3[260],ma[260][260];
int n,m,x,y,i,j,k,t;
double ze(double a)
{
    if (a>ep||a<-ep) return 1;
    return 0;
}
int main()
{
    in>>n>>m;
    ep=0.000000001;
    for(i=1;i<=m;++i)
    {
        in>>x>>y>>z;
        ma[x][y]++;
        ma[y][x]++;
        ma2[x]+=z;
        ma2[y]+=z;
        ma3[x]++;
        ma3[y]++;
    }
    a[1][n]=-1;
    for (i=1;i<=n-1;++i)
        for (j=1;j<=n-1;++j)
    {
        if (i==j) {a[i][j]=-1;}
        else {a[i][j]=ma[i][j]/ma3[j];}
    }
    i=j=1;
    while (i<n)
    {
        for (k=i;k<n;++k)
            if (ze(a[k][j])) break;
        for (t=1;t<=n;++t)
            swap(a[i][t],a[k][t]);
        for (k=j+1;k<=n;++k)
            a[i][k]/=a[i][j];
        a[i][j]=1;
        for (k=1;k<n;++k)
            if (k!=i)
        {double ss=a[k][j];
            for (t=j+1;t<=n;++t)
            a[k][t]-=ss*a[i][t];
            a[k][j]=0;
        }
        i++; j++;
    }
    for (i=1;i<n;++i)
    {
        ras=ras+ma2[i]/ma3[i]*a[i][n];
    }
    out<<fixed<<setprecision(3)<<ras;
    return 0;
}
///val[i]=