ZOJ 3645 BiliBili （高斯消元）

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <iostream>
#include <queue>
#include <stack>
using namespace std;
double re[13][13];
double ans[13];
const double EPS = 1e-12;
bool zero(double x){return -EPS<x&&x<EPS;}
void xch(int a,int b)
{
    for(int i=0;i<12;i++)
    swap(re[a][i],re[b][i]);
}
void gs()
{
    for(int i=0;i<11;i++)
    {
        int a=i;
        for(int j=i+1;j<11;j++)
        if(fabs(re[j][i])>fabs(re[a][i])) a = i;
        if(a!=i)
        xch(a,i);
        for(int j=i+1;j<11;j++)
        {
            double tmp = re[j][i]/re[i][i];
            for(int k=i;k<12;k++)
            re[j][k] = re[j][k] - re[i][k]*tmp;
        }
    }
    for(int i=10;i>=0;i--)
    {
        double tmp = re[i][11];
        for(int j=i+1;j<11;j++)
        tmp -= ans[j]*re[i][j];
        ans[i] = tmp/re[i][i];
    }
    for(int i=0;i<11;i++)
    {
        if(i) printf(" ");
         printf("%.2lf",ans[i]);
    }printf("\n");

}
int main()
{
    freopen("in.txt","r",stdin);
    int cas;scanf("%d",&cas);
    while(cas--)
    {
        for(int i=0;i<12;i++)
        for(int j=0;j<12;j++)
            scanf("%lf",&re[i][j]);
        for(int i=0;i<11;i++)
        {
            re[i][11] = re[11][11]*re[11][11] - re[i][11]*re[i][11];
            for(int j=0;j<11;j++)
            {
                re[i][11] += re[i][j]*re[i][j]-re[11][j]*re[11][j];
                re[i][j] = 2*(re[i][j]-re[11][j]);
            }
        }
        gs();
    }
    return 0;
}