洛谷P3389 【模板】高斯消元法

高斯—约当消元法：

#include<bits/stdc++.h>
using namespace std;
double a[105][105];
double eps=1e-7;
int main(){
    int n;scanf("%d",&n);
    for(int i=1;i<=n;i++) 
        for(int j=1;j<=n+1;j++) scanf("%lf",&a[i][j]);
    for(int i=1;i<=n;i++){//枚举列 
        int max=i;
        for(int j=i+1;j<=n;j++)  
            if(fabs(a[j][i])>fabs(a[max][i])) max=j;//选取该列最大系数，真实目的是选一个非0系数 
        for(int j=1;j<=n+1;j++) swap(a[i][j],a[max][j]);//移到前面
        if(fabs(a[i][i])<eps){//对角线上的主元系数为0，说明没有唯一解 
            puts("No Solution");
            return 0;
        } 
        for(int j=n+1;j>=1;j--) a[i][j]=a[i][j]/a[i][i];//把这一行的主元系数变为1
        for(int j=1;j<=n;j++){//消去主元所在列的其他行的主元 
             if(j!=i){
                double temp=a[j][i]/a[i][i];
                for(int k=1;k<=n+1;k++) a[j][k]-=a[i][k]*temp;
            }
        } 
    }
    for(int i=1;i<=n;i++) printf("%0.2lf\n",a[i][n+1]);
    return 0;
}