Cod sursa(job #1724540)

Utilizator andreiudilaUdila Andrei andreiudila Data 3 iulie 2016 14:41:47
Problema Algoritmul lui Gauss Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.26 kb
#include <fstream>
#include <iomanip>
#include <algorithm>
using namespace std;
ifstream fin("gauss.in");
ofstream fout("gauss.out");
int n,m,i,j,k;
double eps=0.000000001;
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 && fabs(a[j][i])<eps) ++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;
       }

   }
   int nr=0;

   for(i=1;i<=n;++i)
    {
        nr=0;

    for(j=1;j<=m+1;++j)
       if(fabs(a[i][j])>eps) nr++;

       if(nr==1&& fabs(a[i][m+1])>eps) ok=0;
   }


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



   sol[m]=a[m][m+1]/a[m][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;
}