[BZOJ1013][JSOI2008]球形空间产生器sphere（高斯消元）

最新推荐文章于 2018-05-15 13:50:51 发布

Clove_unique

最新推荐文章于 2018-05-15 13:50:51 发布

阅读量674

点赞数

分类专栏：题解省选高斯消元文章标签： BZOJ JSOI 数论

本文链接：https://blog.csdn.net/Clove_unique/article/details/51226479

版权

题解同时被 3 个专栏收录

1097 篇文章 4 订阅

订阅专栏

省选

211 篇文章 0 订阅

订阅专栏

高斯消元

23 篇文章 0 订阅

订阅专栏

题目描述

传送门

题解

按照题意列出方程组
$(x_1-x)^2+(y_1-y)^2+(z_1-z)^2=r^2$
$(x_2-x)^2+(y_2-y)^2+(z_2-z)^2=r^2$
…
$(x_{n+1}-x)^2+(y_{n+1}-y)^2+(z_{n+1}-z)^2=r^2$
各项展开
$x_1^2+x^2+y_1^2+y^2+z_1^2+z^2=2x_1x+2y_1y+2z_1z+r^2$
$x_2^2+x^2+y_2^2+y^2+z_2^2+z^2=2x_2x+2y_2y+2z_2z+r^2$
…
$x_{n+1}^2+x^2+y_{n+1}^2+y^2+z_{n+1}^2+z^2=2x_{n+1}x+2y_{n+1}y+2z_{n+1}z+r^2$
用每个方程与第一个方程左右分别做差就可以把所有未知量的二次项全部消去
$x_2^2-x_1^2+y_2^2-y_1^2+z_2^2-z_1^2=2(x_2-x_1)x+2(y_2-y_1)y+2(z_2-z_1)z$
$x_3^2-x_1^2+y_3^2-y_1^2+z_3^2-z_1^2=2(x_3-x_1)x+2(y_3-y_1)y+2(z_3-z_1)z$
…
$x_{n+1}^2-x_1^2+y_{n+1}^2-y_1^2+z_{n+1}^2-z_1^2=2(x_{n+1}-x_1)x+2(y_{n+1}-y_1)y+2(z_{n+1}-z_1)z$
得到n个线性方程
然后用高斯消元求解即可

代码

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define N 20

const double eps=1e-6;
int n;
double x[N],a[N][N],b[N],ans[N];

void gauss()
{
    for (int i=1;i<=n;++i)
    {
        int num=i;
        for (int j=i+1;j<=n;++j)
            if (fabs(a[j][i])>fabs(a[num][i]))
                num=j;
        for (int j=1;j<=n;++j) swap(a[i][j],a[num][j]);
        swap(b[i],b[num]);
        for (int j=i+1;j<=n;++j)
        {
            if (fabs(a[j][i])<=eps) continue;
            double t=a[j][i]/a[i][i];
            for (int k=1;k<=n;++k)
                a[j][k]-=a[i][k]*t;
            b[j]-=b[i]*t;
        }
    }
    for (int i=n;i>=1;--i)
    {
        ans[i]=b[i]/a[i][i];
        for (int j=i;j>=1;--j)
            b[j]-=ans[i]*a[j][i];
    }
}
int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;++i) scanf("%lf",&x[i]);
    for (int i=1;i<=n;++i)
    {
        for (int j=1;j<=n;++j)
        {
            scanf("%lf",&x[0]);
            a[i][j]=2*(x[0]-x[j]);
            b[i]+=x[0]*x[0]-x[j]*x[j];
        }
    }
    gauss();
    for (int i=1;i<=n;++i) printf("%.3lf%c",ans[i]," \n"[i==n]);
}