Description
有一个球形空间产生器能够在
Input
第一行是一个整数
Output
有且只有一行,依次给出球心的
Sample Input
2
0.0 0.0
-1.0 1.0
1.0 0.0
Sample Output
0.500 1.500
HINT
提示:
给出两个定义:
1、 球心:到球面上任意一点距离都相等的点。
2、 距离:设两个
Source
思路
高斯消元模板。
设球心
代码
#include <cstdio>
#include <cmath>
#include <algorithm>
const int maxn=10;
const double eps=0.0001;
int n;
struct matrix
{
double a[maxn+1][maxn+2];
inline int gauss()
{
for(register int i=1; i<n; ++i)
{
int t=i;
for(register int j=i+1; j<=n; ++j)
{
if(fabs(a[t][i])<fabs(a[j][i]))
{
t=j;
}
}
if(fabs(a[t][i])<eps)
{
return 1;
}
for(register int j=i; j<=n+1; ++j)
{
std::swap(a[t][j],a[i][j]);
}
for(register int j=i+1; j<=n; ++j)
{
if(!a[j][i])
{
continue;
}
double f=a[j][i]/a[i][i];
for(register int k=i; k<=n+1; ++k)
{
a[j][k]-=f*a[i][k];
}
}
}
for(register int i=n; i; --i)
{
a[i][n+1]=a[i][n+1]/a[i][i];
a[i][i]=1;
for(register int j=1; j<i; ++j)
{
a[j][n+1]-=a[j][i]*a[i][n+1];
}
}
return 0;
}
inline int print_ans()
{
for(register int i=1; i<n; ++i)
{
printf("%.3lf ",a[i][n+1]);
}
printf("%.3lf",a[n][n+1]);
return 0;
}
};
matrix t;
double point[maxn+2][maxn+1];
int main()
{
scanf("%d",&n);
for(register int i=1; i<=n+1; ++i)
{
for(register int j=1; j<=n; ++j)
{
scanf("%lf",&point[i][j]);
}
}
for(register int i=1; i<=n; ++i)
{
double sum=0;
for(register int j=1; j<=n; ++j)
{
t.a[i][j]=2*(point[i][j]-point[i+1][j]);
sum+=point[i][j]*point[i][j]-point[i+1][j]*point[i+1][j];
}
t.a[i][n+1]=sum;
}
t.gauss();
t.print_ans();
return 0;
}