Cod sursa(job #2192123)

Utilizator FrequeAlex Iordachescu Freque Data 4 aprilie 2018 19:58:50
Problema Algoritmul lui Gauss Scor 90
Compilator cpp Status done
Runda Arhiva educationala Marime 1.7 kb
#include <iostream>
#include <fstream>
#include <iomanip>

using namespace std;

ifstream fin("gauss.in");
ofstream fout("gauss.out");

const int NMAX = 300 + 5;
const int EPS = 0.0000001;

int n, m;
double x[NMAX];
double a[NMAX][NMAX];

void write()
{
    for (int i = 1; i <= m; ++i)
        fout << fixed << setprecision(8) << x[i] << " ";
}

void read()
{
    fin >> n >> m;
    for (int i = 1; i <= n; ++i)
        for (int j = 1; j <= m + 1; ++j)
            fin >> a[i][j];
}

int main()
{
    int i = 1, j = 1, k;

    read();

    while (i <= n && j <= m)
    {
        for (k = i; k <= n; ++k)
            if (a[k][j] < -EPS || a[k][j] > EPS)
                break;
        if (k == n + 1)
        {
            ++j;
            continue;
        }

        if (k != i)
            for (int l = 1; l <= m + 1; ++l)
                swap(a[i][l], a[k][l]);

        for (int l = j + 1; l <= m + 1; ++l)
            a[i][l] /= a[i][j];
        a[i][j] = 1;

        for (int u = i + 1; u <= n; ++u)
        {
            for (int l = j + 1; l <= m + 1; ++l)
                a[u][l] -= a[u][j] * a[i][l];
            a[u][j] = 0;
        }

        ++i;
        ++j;
    }

    for (int i = n; i; --i)
        for (int j = 1; j <= m + 1; ++j)
            if (a[i][j] < -EPS || a[i][j] > EPS)
            {
                if (j == m + 1)
                {
                    fout << "Imposibil\n";
                    return 0;
                }

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

    write();

    return 0;
}