Cod sursa(job #1724532)

Utilizator andreiudilaUdila Andrei andreiudila Data 3 iulie 2016 14:23:17
Problema Algoritmul lui Gauss Scor 90
Compilator cpp Status done
Runda Arhiva educationala Marime 1.03 kb
#include <fstream>
#include <iomanip>
using namespace std;
ifstream fin("gauss.in");
ofstream fout("gauss.out");
int n,m,i,j,k;
double a[305][305];
double sol[305];
int main()
{
    fin>>n>>m;
    for(i=1;i<=n;++i)
        for(j=1;j<=m+1;++j)
            fin>>a[i][j];

   int ok=1;

   for(i=1;i<=n-1;++i)
   {
       j=i;
       while(j<=n && a[j][i]==0) ++j;

       if(j==n+1) {ok=0; break;}
       else if(j!=i)
       {
           for(k=1;k<=m+1;++k)
                swap(a[i][k],a[j][k]);
       }

       for(j=i+1;j<=n;++j)
       {
           double x=-a[j][i]/a[i][i];
           a[j][i]=0;

           for(k=i+1;k<=m+1;++k)
            a[j][k]+=a[i][k]*x;
       }

   }

   if(!ok) {fout<<"Imposibil"; return 0;}

   sol[m]=a[n][m+1]/a[n][m];
   for(i=m-1;i>=1;--i)
   {
       double suma=0;
       for(j=i+1;j<=m;++j)
        suma+=sol[j]*a[i][j];

       sol[i]=(a[i][m+1]-suma)/a[i][i];
   }

   for(i=1;i<=m;++i)
    fout<<setprecision(8)<<fixed<<sol[i]<<" ";
    return 0;
}