【金牌导航】【洛谷 P4305】【高斯消元】球形空间产生器

【金牌导航】【洛谷 P4305】【高斯消元】球形空间产生器

题目

---

解题思路

因为所有点到球心的距离相等
得：(a-x1)²+(b-x2)²+(c-x3)²=(a1-x1)²+(b1-x2)²+(c1-x3)²
将式子拆开得到 2*(a-a1+b-b1+c-c1)=a²-a1²+b²-b1²+c²-c1²
然后运用高斯消元求得方程解即可

---

代码

#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;
}