Cod sursa(job #1948293)

Utilizator caprariuapCaprariu Alin-Paul caprariuap Data 31 martie 2017 22:52:45
Problema Algoritmul lui Gauss Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 1.52 kb
#include <fstream>
#include <iomanip>
#define eps 0.000001

using namespace std;

ifstream fin("gauss.in");
ofstream fout("gauss.out");

int n,m,i,j,k;
long double a[305][305],aux,ans[305];

int main()
{
    fin >> n >> m;
    for (i=1; i<=n; i++)
        for (j=1; j<=m+1; j++)
            fin >> a[i][j];
    i=1;
    j=1;
    while (i<=n&&j<=m)
    {
        for (k=i; k<=n; k++)
            if (a[k][j]<eps||a[k][j]>-eps)
                break;
        if (k==n+1)
        {
            j++;
            continue;
        }
        if (k!=i)
            for (int h=1; h<=m+1; h++)
            {
                long double aux=a[i][h];
                a[i][h]=a[k][h];
                a[k][h]=aux;
            }
        for (int h=j+1; h<=m+1; h++)
            a[i][h]/=a[i][j];
        a[i][j]=1;
        for (int g=i+1; g<=n; g++)
        {
            for (int h=j+1; h<=m+1; h++)
                a[g][h]=a[g][h]-a[i][h]*a[g][j];
            a[g][j]=0;
        }
        i++;
        j++;
    }
    for (i=n; i; i--)
        for (j=1; j<=m; j++)
            if (a[i][j]>eps||a[i][j]<-eps)
            {
                if (j==m+1)
                {
                    fout << "Imposibil\n";
                    return 0;
                }
                ans[j]=a[i][m+1];
                for (int g=j+1; g<=m; g++)
                    ans[j]=ans[j]-a[i][g]*ans[g];
                break;
            }
    for (i=1; i<=m; i++)
    fout << fixed << setprecision(10) << ans[i] << ' ';
}