【BZOJ】【1013】【JSOI2008】球形空间产生器sphere

高斯消元

  高斯消元模板题

  说实话那个方程组(增广矩阵)怎么搞出来的我还没太弄明白……

技术分享
 1 /**************************************************************
 2     Problem: 1013
 3     User: Tunix
 4     Language: C++
 5     Result: Accepted
 6     Time:0 ms
 7     Memory:1272 kb
 8 ****************************************************************/
 9  
10 //BZOJ 1013
11 #include<cmath>
12 #include<cstdio>
13 #include<cstring>
14 #include<cstdlib>
15 #include<iostream>
16 #include<algorithm>
17 #define rep(i,n) for(int i=0;i<n;++i)
18 #define F(i,j,n) for(int i=j;i<=n;++i)
19 #define D(i,j,n) for(int i=j;i>=n;--i)
20 using namespace std;
21 void read(int &v){
22     v=0; int sign=1; char ch=getchar();
23     while(ch<0||ch>9){ if (ch==-) sign=-1; ch=getchar();}
24     while(ch>=0&&ch<=9){ v=v*10+ch-0; ch=getchar();}
25     v*=sign;
26 }
27 /******************tamplate*********************/
28 const int N=11;
29 typedef double Matrix[N][N];
30 void gauss(Matrix a,int n){
31     int i,j,k,r;
32     for(i=0;i<n;++i){
33         //选一行r并与第i行交换
34         r=i;
35         for(j=i+1;j<n;++j)
36             if (fabs(a[j][i]) > fabs(a[r][i])) r=j;
37         if (r!=i) for(j=0;j<=n;++j) swap(a[r][j],a[i][j]);
38          
39         //与第i+1~n行进行消元
40         for(k=i+1;k<n;++k){
41             double f=a[k][i]/a[i][i];//为了让a[k][i]=0,第i行所乘的倍数
42             for(j=i;j<=n;++j) a[k][j]-=f*a[i][j];
43         }
44     }
45     //回代过程
46     for(i=n-1;i>=0;i--){
47         for(j=i+1;j<n;j++)
48             a[i][n]-=a[j][n]*a[i][j];
49         a[i][n]/=a[i][i];
50     }
51 }
52 Matrix f;
53 int n;
54 int main(){
55     scanf("%d",&n);
56     F(i,0,n){
57         rep(j,n){
58             scanf("%lf",&f[i][j]);
59             f[i][n]-=f[i][j]*f[i][j];
60             f[i][j]*=-2.0;
61             if (i) f[i-1][j]-=f[i][j];
62         }
63         if (i) f[i-1][n]-=f[i][n];
64     }
65     gauss(f,n);
66     printf("%.3f",f[0][n]);
67     for(int i=1;i<n;++i) printf(" %.3f",f[i][n]);
68     return 0;
69 }
View Code

 

郑重声明:本站内容如果来自互联网及其他传播媒体,其版权均属原媒体及文章作者所有。转载目的在于传递更多信息及用于网络分享,并不代表本站赞同其观点和对其真实性负责,也不构成任何其他建议。