Cod sursa(job #3319316)

Utilizator nistor_dora_valentinaNistor Dora Valentina nistor_dora_valentina Data 31 octombrie 2025 16:55:08
Problema Algoritmul lui Gauss Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.64 kb
#include <fstream>
#include <iomanip>
#include <cmath>

#define INF 0.0000001
using namespace std;
ifstream fin("gauss.in");
ofstream fout("gauss.out");
int n, m, i, j;
double a[305][305], x[305];
void swap_matrice(int x, int y)
{
    for(int i=1; i<=m+1; i++)
        swap(a[x][i], a[y][i]);
}
int main()
{
    fin>>n>>m;
    for(i=1; i<=n; i++)
        for(j=1; j<=m+1; j++)
            fin>>a[i][j];
    for(i=1; i<=n; i++)
    {
        for(int j=i; j<=n; j++)
            if(fabs(a[j][i])>INF)
            {
                swap_matrice(i, j);
                break;
            }
        if(fabs(a[i][i])>INF)
        {
            for(j=i+1; j<=n; j++)
            {
                double coef=a[j][i]/a[i][i];
                for(int k=i; k<=m+1; k++)
                    a[j][k]-=a[i][k]*coef;
            }
        }
    }
    //daca coef mei sunt 0, dar a[..][m+1]!=0
    for(int i=n; i>=1; i--)
    {
        bool ok=false;
        for(j=1; j<=m; j++)
            if(fabs(a[i][j])>INF)
            {
                ok=true;
                break;
            }
        if(ok==false && fabs(a[i][m+1])>INF)
        {
            fout<<"Imposibil";
            return 0;
        }
    }
    for(i=n; i>=1; i--)
    {
        for(j=1; j<=m; j++)
        {
            if(fabs(a[i][j])>INF)
            {
                x[j]=a[i][m+1];
                for(int k=j+1; k<=m; k++)
                    x[j]-=x[k]*a[i][k];
                x[j]/=a[i][j];
                break;
            }
        }
    }
    for(i=1; i<=m; i++)
        fout<<fixed<<setprecision(8)<<x[i]<<" ";
    return 0;
}