码迷,mamicode.com
首页 > Web开发 > 详细

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

时间:2015-01-24 17:18:33      阅读:144      评论:0      收藏:0      [点我收藏+]

标签:

高斯消元

  高斯消元模板题

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

技术分享
 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

 

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

标签:

原文地址:http://www.cnblogs.com/Tunix/p/4246082.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!