Cod sursa(job #812984)

Utilizator Al3ks1002Alex Cociorva Al3ks1002 Data 14 noiembrie 2012 19:51:02
Problema Algoritmul lui Gauss Scor 20
Compilator cpp Status done
Runda Arhiva educationala Marime 0.98 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];
float 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("%f",&a[i][j]);
    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];
        for(k=j;k<=m+1;k++) a[i][k]/=v;
        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;
}