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

高斯消元

　　高斯消元模板题

 /**************************************************************
     Problem: 1013
     User: Tunix
     Language: C++
     Result: Accepted
     Time:0 ms
     Memory:1272 kb
 ****************************************************************/

 //BZOJ 1013
 #include<cmath>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 using namespace std;
 void read(int &v){
     v=; int sign=; char ch=getchar();
     while(ch<''||ch>''){ if (ch=='-') sign=-; ch=getchar();}
     while(ch>=''&&ch<=''){ v=v*+ch-''; ch=getchar();}
     v*=sign;
 }
 /******************tamplate*********************/
 const int N=;
 typedef double Matrix[N][N];
 void gauss(Matrix a,int n){
     int i,j,k,r;
     for(i=;i<n;++i){
         //选一行r并与第i行交换
         r=i;
         for(j=i+;j<n;++j)
             if (fabs(a[j][i]) > fabs(a[r][i])) r=j;
         if (r!=i) for(j=;j<=n;++j) swap(a[r][j],a[i][j]);

         //与第i+1～n行进行消元
         for(k=i+;k<n;++k){
             double f=a[k][i]/a[i][i];//为了让a[k][i]=0,第i行所乘的倍数
             for(j=i;j<=n;++j) a[k][j]-=f*a[i][j];
         }
     }
     //回代过程
     for(i=n-;i>=;i--){
         for(j=i+;j<n;j++)
             a[i][n]-=a[j][n]*a[i][j];
         a[i][n]/=a[i][i];
     }
 }
 Matrix f;
 int n;
 int main(){
     scanf("%d",&n);
     F(i,,n){
         rep(j,n){
             scanf("%lf",&f[i][j]);
             f[i][n]-=f[i][j]*f[i][j];
             f[i][j]*=-2.0;
             if (i) f[i-][j]-=f[i][j];
         }
         if (i) f[i-][n]-=f[i][n];
     }
     gauss(f,n);
     printf("%.3f",f[][n]);
     for(int i=;i<n;++i) printf(" %.3f",f[i][n]);
     return ;
 }