Cod sursa(job #2301379)

Utilizator zdavid112zIon David-Gabriel zdavid112z Data 12 decembrie 2018 21:17:00
Problema Algoritmul lui Gauss Scor 40
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.12 kb
#include <bits/stdc++.h>
using namespace std;
#define EPS 1e-8
int p[303],n,m,i,k,j;
double a[303][303],sol[303],r,pp;
int main()
{
    freopen("gauss.in","r",stdin);
    freopen("gauss.out","w",stdout);
    scanf("%d%d",&m,&n);
    for(i=1; i<=m; i++)
        for(j=1; j<=n+1; j++)
            scanf("%lf",&a[i][j]);
    for(i=1; i<=m; i++)
    {
        p[i]=0;
        for(j=1; j<=n+1; j++)
        {
            if(abs(a[i][j])>EPS)
            {
                p[i]=j;
                break;
            }
        }
        if(p[i]==n+1)
        {
            printf("Imposibil");
            return 0;
        }
        if(p[i]==0)
            continue;

        for(j=1; j<=m; j++)
        {
            if(j!=i&&abs(a[j][p[i]])>EPS)
            {
                r=a[j][p[i]]/a[i][p[i]];
                for(k=1; k<=n+1; k++)
                    a[j][k]=a[j][k]-a[i][k]*r;
            }
        }

    }
    for(i=1; i<=m; i++)
    {
        if(p[i])
            sol[p[i]]=a[i][n+1]/a[i][p[i]];
    }
    for(i=1; i<=n; i++)
    {
        printf("%.8f ",sol[i]);
    }
    return 0;
}