高斯消元1

【训练题】高斯消元[1]

 

描述

 

　　贾老二是个品学兼优的好学生，但由于智商问题，算术学得不是很好，尤其是在解方程这个方面。虽然他解决 2x=2 这样的方程游刃有余，但是对于 {x+y=3 x-y=1} 这样的方程组就束手无策了。于是他要你来帮忙。前提是一次方程组且保证在integer的范围内可以处理所有问题。

输入

 

　　第一行一个数字N（1≤N≤100）表示要求的未知数的个数，同时也是所给的方程个数。　　第2到N+1行，每行N+1个数。前N个表示第1到N个未知数的系数。第N+1个数表示N个未知数乘以各自系数后的加和。（保证有唯一整数解）

输出

 

　　一行N个数，表示第1到N个未知数的值。

输入样例 1 

2
1 1 3
1 -1 1

输出样例 1

2 1

#include<cstdio>
#include<cmath>
#include<algorithm>
#define maxn 105
using namespace std;

int n;
double a[maxn],c[maxn][maxn];

int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n+1;j++) 
        scanf("%lf",&c[i][j]);
    }

    for(int i=1;i<=n;i++){
        for(int j=i;j<=n;j++){
            if(fabs(c[j][i])>1e-6){
                swap(c[i],c[j]);
                break;
            }
        }
        for(int j=1;j<=n;j++){
            if(i!=j){
                double rate=c[j][i]/c[i][i];
                for(int k=1;k<=n+1;k++){
                    c[j][k]-=c[i][k]*rate;
                }    
            }
        }
    }
    for(int i=1;i<=n;i++){
        a[i]=c[i][n+1]/c[i][i];
        a[i]=(a[i]*10+0.5)/10;
        printf("%.0lf ",a[i]);
    }
    return 0;
}