[JSOI2008]球形空间产生器 高斯消元

njH2Q
评测方式 云端评测
标签 各省省选2008江苏
难度 省选/NOI-
时空限制 1000ms / 128MB

题目描述
有一个球形空间产生器能够在 nn 维空间中产生一个坚硬的球体。现在，你被困在了这个 nn 维球体中，你只知道球面上 n+1n+1 个点的坐标，你需要以最快的速度确定这个 nn 维球体的球心坐标，以便于摧毁这个球形空间产生器。

输入输出格式
输入格式：
第一行是一个整数 nn (1<=N=10)(1<=N=10) 。接下来的 n+1n+1 行，每行有 nn 个实数，表示球面上一点的 nn 维坐标。每一个实数精确到小数点后 66 位，且其绝对值都不超过 2000020000 。

输出格式：
有且只有一行，依次给出球心的 nn 维坐标（ nn 个实数），两个实数之间用一个空格隔开。每个实数精确到小数点后 33 位。数据保证有解。你的答案必须和标准输出一模一样才能够得分。

输入输出样例
输入样例#1：
2
0.0 0.0
-1.0 1.0
1.0 0.0
输出样例#1：
0.500 1.500

省选难度的，（其实还好

高斯消元即可；

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<map>
#include<set>
#include<vector>
#include<queue>
#include<string>
#include<bitset>
#include<ctime>
#include<deque>
#include<stack>
#include<functional>
#include<sstream>
using namespace std;
#define maxn 500005
#define inf 0x3f3f3f3f
#define INF 0x7fffffff
typedef long long  ll;
typedef unsigned long long ull;
#define ms(x) memset(x,0,sizeof(x))
const long long int mod = 1e9 + 7;
#define eps 1e-7

inline int read()
{
	int x = 0, k = 1; char c = getchar();
	while (c < '0' || c > '9') { if (c == '-')k = -1; c = getchar(); }
	while (c >= '0' && c <= '9')x = (x << 3) + (x << 1) + (c ^ 48), c = getchar();
	return x * k;
}

int quickpow(int a,int b)
{
    int ans=1;
    while(b){
        if(b&1)ans=ans*a;
        b>>=1;
        a=a*a;
    }
    return ans;
}

double a[20][20],b[20],c[20][20];

int n;


int main()
{
	//ios::sync_with_stdio(false);
    n=read();
    for(int  i=1;i<=n+1;i++)
        for(int j=1;j<=n;j++)cin>>a[i][j];
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
                
            c[i][j]=2*(a[i][j]-a[i+1][j]);
            
            b[i]+=a[i][j]*a[i][j]-a[i+1][j]*a[i+1][j];
        }
    }

    for(int i=1;i<=n;i++){
        for(int j=i;j<=n;j++){
            if(fabs(c[j][i])>eps){
                for(int k=1;k<=n;k++)swap(c[i][k],c[j][k]);
                swap(b[i],b[j]);
            }
        }
        for(int j=1;j<=n;j++){
            if(i==j)continue;
            double rte=c[j][i]*1.0/c[i][i];
            for(int k=i;k<=n;k++)c[j][k]-=c[i][k]*rte;
            b[j]-=b[i]*rte;
        }
    }
    for(int i=1;i<n;i++)printf("%.3f ",1.0*b[i]/c[i][i]);
    printf("%.3f\n",b[n]/c[n][n]);
}