acwing 207 球形空间产生器 高斯消元

题目地址

题目指出有解。

设球心为(x1,x2,…,xn) ,则我们可以得到${\sum }_{j=0}^n(a_{i,j}-x_j{)}^2$ = C。

将其化简，让第一行减去第二行，依次类推。
可以得到${\sum }_{j=0}^n(a_{i,j}-a_{i+1,j})={\sum }_{j=0}^n(a_{i,j}^2-a_{i+1,j}^2)$

可以根据这个式子，把矩阵表示出来。
接下来处理矩阵，保证满足对角线上的元素不为0.

再消去其他方程x[i]的系数。

最后可以得到

$c_{1,1}$ * $x_1$ = b[1]
…

代码：

#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<iomanip>
#define int long long
#define sc scanf
#define pf printf
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
typedef unsigned long long LL;
 
const int INF = 0x3f3f3f3f;
const double eps = 1e-4;
const int mod = 10000;
const int N = 20;

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

signed main(){
//	IOS;
	#ifdef ddgo
		freopen("C:/Users/asus/Desktop/ddgoin.txt","r",stdin);
	#endif
	
	int n; cin>>n;
	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]) > 1e-8){//把矩阵处理为我们想要的，即让（i,i）位置为不为0的数 
				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 rate = c[j][i]/c[i][i];
			for(int k=1;k<=n;k++) c[j][k] -= c[i][k]*rate;
			b[j] -= b[i]*rate;
		}
	}
	for(int i=1;i<n;i++) cout<<fixed<<setprecision(3)<<b[i]/c[i][i]<<" ";
	cout<<fixed<<setprecision(3)<<b[n]/c[n][n]<<endl;
	
	return 0;
}