球形空间产生器[BZOJ1013]

最新推荐文章于 2021-02-12 11:32:18 发布

zsyz_lb2003

最新推荐文章于 2021-02-12 11:32:18 发布

阅读量100

点赞数 1

分类专栏：蓝书做题记录蓝书数论 caioj数论高斯消元

本文链接：https://blog.csdn.net/zyszlb2003/article/details/89530541

版权

蓝书做题记录同时被 3 个专栏收录

66 篇文章 2 订阅

订阅专栏

caioj数论

54 篇文章 0 订阅

订阅专栏

蓝书数论

29 篇文章 1 订阅

订阅专栏

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

传送门

思路

因为一个球体的球面上的所有点到圆心的距离相等，因此只需求出一个点( $x_1,x_2,x_3,......,x_n$ )，使得:

$\large\sum_{j=1}^{n}(a_{i,j}-x_j)^2=C (i\in[1,n+1])$

其中 $C$ 为常量，球面上第 $i$ 个已知点的坐标为( $a_{i,1},a_{i,2},......,a_{i,n}$ )。

如何将这个 $n + 1$ 个 $n$ 元二次方程组转化为线性方程组呢？观察到，方程右边都有 $C$

我们不妨两两做差，将 $C$ 消掉。

得到:
$\large\sum_{j=1}^{n}(a_{i,j}^2-a_{i+1,j}^2-2x_j(a_{i,j}-a_{i+1,j}))=0(i\in[1,n])$

将变量放左，常量放右，有:
$\large\sum_{j=1}^n2(a_{i,j}-a_{i+1,j})x_j=\sum_{j=1}^n(a_{i,j}^2-a_{i+1,j}^2)(i\in[1,n])$

这就是一个线性方程组了。

$\begin{bmatrix}2(a_{1,1}-a_{2,1})&2(a_{1,2}-a_{2,2})&...&2(a_{1,n}-a_{2,n})&\large\sum_{j=1}^n(a_{1,j}^2-a_{2,j}^2)\\2(a_{2,1}-a_{3,1})&2(a_{2,2}-a_{3,2})&...&2(a_{2,n}-a_{3,n})&\large\sum_{j=1}^n(a_{2,j}^2-a_{3,j}^2)\\\vdots&\vdots&\ddots&\vdots&\vdots\\2(a_{n,1}-a_{n+1,1})&2(a_{n,2}-a_{n+1,2})&...&2(a_{n,n}-a_{n+1,n})&\large\sum_{j=1}^n(a_{n,j}^2-a_{n+1,j}^2)\end{bmatrix}$

高斯消元求 $x_j$ 即可。

AC code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#define eps 1e-10
using namespace std;
double a[20][20],c[20][20],b[20];int n;
void gauss()
{
	for(int i=1;i<=n;i++)
	{
		int x=i;
		while(fabs(c[x][i])<eps&&x<=n)x++;
		for(int j=1;j<=n;j++)swap(c[x][j],c[i][j]);swap(b[x],b[i]);
		for(int k=1;k<=n;k++)
		{
			if(k==i)continue;
			b[k]-=c[k][i]/c[i][i]*b[i];
			for(int j=n;j>=i;j--)c[k][j]-=c[k][i]/c[i][i]*c[i][j];
		}
	}
	for(int i=1;i<n;i++)printf("%.3lf ",b[i]/c[i][i]);
	printf("%.3lf\n",b[n]/c[n][n]);
}
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]);
	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];
		}
	gauss();
	return 0;
}