Cod sursa(job #820965)

Utilizator geniucosOncescu Costin geniucos Data 21 noiembrie 2012 14:41:34
Problema Algoritmul lui Gauss Scor 20
Compilator cpp Status done
Runda Arhiva educationala Marime 1.19 kb
#include<cstdio>
using namespace std;
int k,i,n,m,j,nr;
double eps,s,a[308][308],x[308],b[308][308];
double mod(double x)
{
    if(x<0) return -x;
    return x;
}
int main()
{
freopen("gauss.in","r",stdin);
freopen("gauss.out","w",stdout);
scanf("%d",&n);
scanf("%d",&m);
for(i=1;i<=n;i++)
    for(j=1;j<=m+1;j++)
    {
        scanf("%lf",&a[i][j]);
        b[i][j]=a[i][j];
    }
for(i=1;i<=n;i++)
{
    if(a[i][i]==0)
    {
        printf("Imposibil\n");
        return 0;
    }
    for(j=m+1;j>=i;j--)
        a[i][j]=(double)a[i][j]/a[i][i];
    for(k=i+1;k<=n;k++)
    {
        for(j=m+1;j>=i;j--)
            a[k][j]=(double)a[k][j]-a[k][i]*a[i][j];
    }
}
for(j=1;j<=m;j++)
    if(a[n][j]==0) nr++;
/*if(nr!=m-1)
{
    printf("Imposibil\n");
    return 0;
}*/
x[n]=(double)a[n][m+1];
for(i=n-1;i>=1;i--)
{
    s=0;
    for(j=i+1;j<=m;j++)
        s=s+a[i][j]*x[j];
    x[i]=a[i][m+1]-s;
}
eps=0.0000001;
for(i=1;i<=n;i++)
{
    s=0;
    for(j=1;j<=m;j++)
        s=(double)s+b[i][j]*x[j];
    if(mod(s-b[i][m+1])>eps)
    {
        printf("Imposibil\n");
        return 0;
    }
}
for(i=1;i<=m;i++)
    printf("%.6lf ",x[i]);
return 0;
}