Cod sursa(job #813023)

Utilizator Al3ks1002Alex Cociorva Al3ks1002 Data 14 noiembrie 2012 20:46:03
Problema Algoritmul lui Gauss Scor 60
Compilator cpp Status done
Runda Arhiva educationala Marime 1.03 kb
#include<cstdio>
#define eps 0.00001
#define zero(x) (x>=-eps && x<=eps)
#define difzero(x) (x<-eps || x>eps)
int n,m,i,j,k,p,L[301],C[301],x;
double a[301][301],v;
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("%d",&x);a[i][j]=(double)x;}
    for(j=1;j<=m;j++)
    {
        for(i=1;i<=n;i++)
            if(!L[i] && difzero(a[i][j])) break;
        L[i]=j;C[j]=i;
        v=a[i][j];
        if(v)for(k=j;k<=m+1;k++) a[i][k]/=v;
        else continue;
        for(k=1;k<=n;k++)
        {
            if(k==i || zero(a[k][j])) continue;
            v=a[k][j];
            for(p=j;p<=m+1;p++) a[k][p]-=a[i][p]*v;
        }
    }
    for(i=1;i<=n;i++)
        if(!L[i] && difzero(a[i][m+1]))
        {
            printf("Imposibil");
            return 0;
        }
    for(j=1;j<=m;j++)
    {
        if(C[j]) printf("%.10f ",a[C[j]][m+1]);
        else printf("0 ");
    }
    return 0;
}