Cod sursa(job #1784710)

Utilizator team_nameUPB Banu Popa Visan team_name Data 20 octombrie 2016 13:44:05
Problema Algoritmul lui Gauss Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.51 kb
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <algorithm>
using namespace std;

const int nmax = 310;
const double eps = 1e-12;

int n, m;
double a[nmax][nmax], sol[nmax];

void gauss()
{
    int ec = 1, nec = 1;
    while (ec <= n && nec <= m)
    {
        int k = ec;
        while (k <= n && fabs(a[k][nec]) < eps)
            k ++;
        if (k > n)
        {
            nec ++;
            continue;
        }
        for (int j = 1; j <= m + 1; ++ j)
            swap(a[ec][j], a[k][j]);

        for (int i = 1; i <= n; ++ i)
            if (i != ec)
            {
                double f = a[i][nec] / a[ec][nec];
                for (int j = 1; j <= m + 1; ++ j)
                    a[i][j] -= a[ec][j] * f;
            }
        ec ++;
        nec ++;
    }
    for (int i = 1; i <= n; ++ i)
    {
        int j;
        for (j = 1; j <= m + 1; ++ j)
            if (fabs(a[i][j]) > eps)
                break;
        if (j == m + 1)
        {
            printf("Imposibil\n");
            exit(0);
        }
    }

    for (int i = 1; i <= n; ++ i)
        sol[i] = a[i][m + 1] / a[i][i];

    for (int i = 1; i <= m; ++ i)
        printf("%.8f ", sol[i]);
}

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

    //freopen("in", "r", stdin);

    scanf("%i %i", &n, &m);

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

    gauss();
}