高斯消元

高斯消元最朴素的板子......

求解线性方程组。

先转换为增广矩阵。

从上到下一行一行消元呗。

每次挑个首列最大的行消。

求出来可能不是整数，要用double。

具体看代码。

最后求出来x_n之后回带求出x₁、x₂......x_n-1。

给个洛谷的模板题链接：传送门

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

int n;
double c[105][105];
double x[105];
double eps=1e-7;

double v(double rx)
{
    return rx>0?rx:-rx;
}

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++)
    {
        int p=i;
        for(int j=i+1;j<=n;j++)
            if(v(c[j][i])>v(c[p][i]))p=j;
        if(v(c[p][i])<eps)
        {
            printf("No Solution");
            return 0;
        }
        if(i!=p)
        {
            for(int j=i;j<=n+1;j++)
                swap(c[i][j],c[p][j]);
        }
        double div=c[i][i];
        for(int j=i;j<=n+1;j++)
            c[i][j]/=div;
        for(int j=i+1;j<=n;j++)
        {
            div=c[j][i];
            for(int k=i;k<=n+1;k++)
                c[j][k]-=c[i][k]*div;
        }
    }
    x[n]=c[n][n+1];
    for(int i=n-1;i;i--)
    {
        x[i]=c[i][n+1];
        for(int j=i+1;j<=n;j++)
            x[i]-=c[i][j]*x[j];
    }
    for(int i=1;i<=n;i++)
        printf("%.2lf
",x[i]);
    return 0;
}

这样就OK啦。