bzoj 1013: [JSOI2008]球形空间产生器sphere

最新推荐文章于 2019-09-16 09:45:00 发布

weixin_34004750

最新推荐文章于 2019-09-16 09:45:00 发布

阅读量54

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/xydddd/p/5226862.html

本文介绍了一种利用高斯消元法求解给定平面上n+1个点坐标对应的圆心坐标的算法实现。通过构造n个方程，采用C++语言实现了高斯消元的过程，并给出了具体的代码示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cmath>
 4 #define eps 1e-6
 5 #define M 15
 6 using namespace std;
 7 double a[M],b[M][M];
 8 int n;
 9 int main()
10 {
11     scanf("%d",&n);
12     for(int i=1;i<=n;i++)
13       scanf("%lf",&a[i]);
14     for(int i=1;i<=n;i++)
15       for(int j=1;j<=n;j++)
16         {
17             double t;
18             scanf("%lf",&t);
19             b[i][j]=2*(t-a[j]);
20             b[i][n+1]+=t*t-a[j]*a[j]; 
21         }
22     for(int i=1;i<=n;i++)
23       {
24         int to;
25         for(to=i;to<=n;to++)
26           if(b[to][i]!=0)
27             break;
28         if(i!=to)
29           for(int j=1;j<=n+1;j++)
30             swap(b[i][j],b[to][j]);
31         double t=b[i][i];
32         for(int j=i;j<=n+1;j++)
33           b[i][j]/=t;
34         for(int j=1;j<=n;j++)
35           if(j!=i)
36             {
37                 t=b[j][i];
38                 for(int k=1;k<=n+1;k++)
39                   b[j][k]-=t*b[i][k];
40             }
41       }
42     for(int i=1;i<n;i++)
43       printf("%.3lf ",b[i][n+1]);
44     printf("%.3lf\n",b[n][n+1]);
45     return 0;
46 }