BZOJ1013(JSOI2008)[球形空间产生器sphere]--高斯消元

最新推荐文章于 2021-01-03 16:35:55 发布

Greninja_Wu

最新推荐文章于 2021-01-03 16:35:55 发布

阅读量713

点赞数

分类专栏： -------------题解--------------- BZOJ题解 -------------数论--------------- 高斯消元 BZOJ题解 By Greninja 部分常见刷题网站的题解 By Greninja 文章标签： BZOJ 高斯消元

本文链接：https://blog.csdn.net/CHNWJD/article/details/78057138

版权

-------------题解--------------- 同时被 3 个专栏收录

119 篇文章 0 订阅

订阅专栏

部分常见刷题网站的题解 By Greninja

118 篇文章 1 订阅

订阅专栏

BZOJ题解

90 篇文章 0 订阅

订阅专栏

【题目链接】
bzoj1013

【解题报告】

设球心坐标为( $O_1$ , $O_2$ ,…, $O_n$ )。
第i个点坐标为( $X_{(i,1)}$ , $X_{(i,2)}$ ,…, $X_{(i,n)}$ )
由题意得
$\sum(X_{(1,i)}^2-2*X_{(1,i)}*O_i+O_i^2)=\sum(X_{(2,i)}^2-2*X_{(2,i)}*O_i+O_i^2)=...=\sum(X_{(n+1,i)}^2-2*X_{(n+1,i)}*O_i+O_i^2)$
所以 $\sum(2*(X_{(2,i)}-X_{(1,i)})*O_i)=\sum(X_{(2,i)}^2-X_{(1,i)}^2)$
$\sum(2*(X_{(3,i)}-X_{(1,i)})*O_i)=\sum(X_{(3,i)}^2-X_{(1,i)}^2)$
…
$\sum(2*(X_{(n+1,i)}-X_{(1,i)})*O_i)=\sum(X_{(n+1,i)}^2-X_{(1,i)}^2)$
就可以用高斯消元解线性方程了

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=15;
const double P=1e-10;
int n;
double a[maxn][maxn],M[maxn][maxn],ans[maxn];
inline int Read()
{
    int res=0;
    char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') res=res*10+ch-48,ch=getchar();
    return res;
}
bool Fc(double x) {if (x<=P&&x>=-P) return 1; return 0;}
int main()
{
    freopen("1013.in","r",stdin);
    freopen("1013.out","w",stdout);
    n=Read();
    for (int i=1; i<=n+1; i++)
     for (int j=1; j<=n; j++) scanf("%lf",&a[i][j]);
    for (int i=2; i<=n+1; i++)
     for (int j=1; j<=n; j++)
      M[i-1][j]=(a[i][j]-a[1][j])*2,M[i-1][n+1]+=a[i][j]*a[i][j]-a[1][j]*a[1][j];
    for (int i=1; i<=n; i++)
    {
        int where;
        for (int j=i; j<=n; j++)
         if (!Fc(M[j][i])) {where=j; break;}
        swap(M[where],M[i]);
        for (int j=i+1; j<=n; j++) {double num=M[j][i]/M[i][i]; for (int k=i; k<=n+1; k++) M[j][k]-=num*M[i][k];}
    }                
    for (int i=n; i; i--)   
    {
        ans[i]=M[i][n+1];
        for (int j=i+1; j<=n; j++) ans[i]-=M[i][j]*ans[j];
        ans[i]/=M[i][i];
    }
    for (int i=1; i<n; i++) printf("%.3lf ",ans[i]); printf("%.3lf\n",ans[n]);
    return 0;
}