#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;
}
ZOJ 3645 BiliBili (高斯消元)
最新推荐文章于 2017-07-10 23:51:14 发布