Cod sursa(job #1420672)

Utilizator heracleRadu Muntean heracle Data 18 aprilie 2015 20:20:12
Problema Algoritmul lui Gauss Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.16 kb
#include <cstdio>

FILE* in=fopen("gauss.in","r");
FILE* out=fopen("gauss.out","w");

const int Q=308;

int n,m;

int v[Q][Q];

int sel[Q];

int main()
{
    fscanf(in,"%d%d",&n,&m);

    for(int i=1; i<=n; i++)
    {
        for(int j=1; j<=m+1; j++)
        {
            fscanf(in,"%d",&v[i][j]);
        }
    }

    int kept;

    for(int k=1; k<=m; k++)
    {
        int f;
        for(f=1; f<=n; f++)
        {
            if(sel[f]==0 && v[f][k]!=0)
            {
                sel[f]=1;

                for(int i=1; i<=n; i++)
                {
                    if(f==i)
                        continue;

                    kept=v[i][k];

                    for(int j=1; j<=m+1; j++)
                    {
                        v[i][j]=v[i][j]*v[f][k]-v[f][j]*kept;
                    }
                }

                break;
            }
        }
        if(f==n+1)
        {
            fprintf(out,"Imposibil");
            return 0;
        }

    }

    for(int j=1; j<=m; j++)
    {
        for(int i=1; i<=n; i++)
        {
            if(v[i][j]!=0)
                fprintf(out,"%.8lf ",(double)v[i][m+1]/v[i][j]);
        }
    }

    return 0;
}