高斯消元

#include<bits/stdc++.h>
using namespace std;
double m[111][111];
double ans[111];
double eps=1e-8;
int main(){
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            scanf("%lf",&m[i][j]);
    for(int i=1;i<=n;i++){
        int maxx=i;//找出第一列最大的系数。
        for(int j=i+1;j<=n;j++)
            if(fabs(m[maxx][i])<fabs(m[j][i]))
                maxx=j;
        if(fabs(m[maxx][i])<eps){
            printf("No Solution");//无解。
            return 0;
        }
        if(i!=maxx) 
            swap(m[i],m[maxx]); //交换第一行和系数最大行。
        double div=m[i][i];//化简这一横行系数，让除0外第一个为1。
        for(int j=i;j<=n+1;j++)
            m[i][j]/=div;
        for(int j=i+1;j<=n;j++){
            div=m[j][i];
            for(int k=i;k<=n+1;k++)
                m[j][k]-=m[i][k]*div;
        }
    }
    ans[n]=m[n][n+1];
    for(int i=n-1;i>=1;i--){
        ans[i]=m[i][n+1];
        for(int j=i+1;j<=n;j++)
            ans[i]-=(m[i][j]*ans[j]);
    }
    for(int i=1;i<=n;i++)
        printf("%.2lf\n",ans[i]);
    // while(1);
    return 0;
}