Cod sursa(job #1388791)

Utilizator robertstrecheStreche Robert robertstreche Data 15 martie 2015 18:39:22
Problema Algoritmul lui Gauss Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 1.05 kb
#include <cstdio>
#include <algorithm>

#define NMAX 305

using namespace std;

int n,m,k,kk;
float sol[NMAX],a[NMAX][NMAX];

int main()
{
    freopen("gauss.in","r",stdin);
    freopen("gauss.out","w",stdout);

    scanf("%d %d",&n,&m);m++;

    for (int i=1;i<=n;i++)
      for (int j=1;j<=m;j++)
        scanf("%f ",&a[i][j]);

    for (int i=1,j=1;i<=n && j<m;i++,j++)
     {
        for (k=i;k<=n;k++)
          if (a[k][j])break;
        if (k==n+1){i--;continue;}
        if (k!=i)
         for (int kk=j;kk<=m;kk++)
          swap(a[i][kk],a[k][kk]);
        for (k=j+1;k<=m;k++)
          a[i][k]/=a[i][j];
        a[i][j]=1;
        for (k=i+1;k<=n;k++)
         {for (kk=j+1;kk<=m;kk++)
           a[k][kk]-=a[i][kk]*a[k][j];
          a[k][j]=0;
         }
     }
    for (int i=m-1;i>=1;i--)
     {   sol[i]=a[i][m];
         for (int j=i+1;j<=m;j++)
           sol[i]-=sol[j]*a[i][j];
         sol[i]/=a[i][i];
     }
    for (int i=1;i<m;i++)
      printf("%f ",sol[i]);

    fclose(stdin);
    fclose(stdout);
}