【金牌导航】【洛谷 P4305】【高斯消元】球形空间产生器
题目
解题思路
因为所有点到球心的距离相等
得:(a-x1)2+(b-x2)2+(c-x3)2=(a1-x1)2+(b1-x2)2+(c1-x3)2
将式子拆开得到 2*(a-a1+b-b1+c-c1)=a2-a12+b2-b12+c2-c12
然后运用高斯消元求得方程解即可
代码
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n;
double a[15][15],f[15][15];
int main()
{
scanf("%d",&n);
for (int i=1;i<=n+1;i++)
{
for (int j=1;j<=n;j++)
scanf("%lf",&a[i][j]);
if (i>1)
{
double js=0;
for (int j=1;j<=n;j++)
{
f[i-1][j]=2*(a[i][j]-a[i-1][j]);
js+=a[i][j]*a[i][j]-a[i-1][j]*a[i-1][j];
}
f[i-1][n+1]=js;
}
}
for (int i=1;i<=n;i++)
{
int w=i;
for (int j=i+1;j<=n;j++)
if (fabs(f[j][i])>fabs(f[w][i])) w=j;
if (w!=i)
for (int j=i;j<=n+1;j++)
swap(f[i][j],f[w][j]);
for (int j=1;j<=n;j++)
if (i!=j)
{
double tmp=1.0*f[j][i]/f[i][i];
for (int k=1;k<=n+1;k++)
f[j][k]-=tmp*f[i][k];
}
}
for (int i=1;i<=n;i++)
printf("%.3lf ",f[i][n+1]/f[i][i]);
return 0;
}