Cod sursa(job #1877532)

Utilizator andreeainfo_dAndreea Dutulescu andreeainfo_d Data 13 februarie 2017 15:09:59
Problema Algoritmul lui Gauss Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 1.17 kb
#include <bits/stdc++.h>
using namespace std;
//#define EPS 0.000000001
int p[303],n,m,i,ii,j;
double a[303][303],sol[303],r, EPS=1e-8,pp;
double modul(double x)
{
    if(x<0)return -x;
    return x;
}
int main()
{
    freopen("gauss.in","r",stdin);
    freopen("gauss.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)
        for(j=1;j<=m+1;j++)
        scanf("%lf",&a[i][j]);
    for(i=1;i<=n;i++)
    {
        p[i]=0;
        for(j=1;j<=m+1;j++)
        {
            if(modul(a[i][j])>EPS)
            {
                p[i]=j;
                break;
            }
        }
        if(p[i]==m+1){printf("Imposibil\n");return 0;}
        if(p[i]!=0)
        {
            for(ii=1;ii<=n;ii++)
            {
                if(ii!=i&&modul(a[ii][p[i]])>EPS)
                {
                    r=a[ii][p[i]]/a[i][p[i]];
                    for(j=1;j<=m+1;j++)
                        a[ii][j]=a[ii][j]-a[i][j]*r;
                }
            }
        }
    }
    for(i=1;i<=n;i++)
    {
        if(p[i]>0)sol[p[i]]=a[i][n+1]/a[i][p[i]];
    }
    for(i=1;i<=m;i++)
    {
        printf("%.10lf ",sol[i]);
    }
    return 0;
}