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

以下这个好像叫高斯约旦消元法，没有回代

https://www.luogu.org/blog/37781/solution-p3389

#include<cstdio>
#include<algorithm>
#include<cmath>
#define eps 1e-7
using namespace std;
double a[][];
int n;
void swap_line(int x,int y)
{
    for(int i=;i<=n+;i++)    swap(a[x][i],a[y][i]);
}
int dcmp(double x)
{
    if(fabs(x)<eps)    return ;
    else    return x>?:-;
}
int main()
{
    int i,j,k,nx;double nmax,t;
    scanf("%d",&n);
    for(i=;i<=n;i++)
        for(j=;j<=n+;j++)
            scanf("%lf",&a[i][j]);
    for(i=;i<=n;i++)
    {
        nmax=;
        for(j=i;j<=n;j++)
            if(fabs(a[j][i])>fabs(nmax))
            {
                nmax=a[j][i];
                nx=j;
            }
        if(dcmp(nmax)==)
        {
            puts("No Solution");
            return ;
        }
        swap_line(nx,i);
        for(j=;j<=n+;j++)    a[i][j]/=nmax;
        for(j=;j<=n;j++)
            if(j!=i)
            {
                t=a[j][i]/a[i][i];//可以直接t=a[j][i];因为a[i][i]一定是1
                for(k=;k<=n+;k++)
                    a[j][k]-=t*a[i][k];
            }
    }
    for(i=;i<=n;i++)    printf("%.2lf\n",a[i][n+]);
    return ;
}

还可以卡常

#include<cstdio>
#include<algorithm>
#include<cmath>
#define eps 1e-7
using namespace std;
double a[][];
int n;
void swap_line(int x,int y)
{
    for(int i=;i<=n+;i++)    swap(a[x][i],a[y][i]);
}
int dcmp(double x)
{
    if(fabs(x)<eps)    return ;
    else    return x>?:-;
}
int main()
{
    int i,j,k,nx;double nmax;
    scanf("%d",&n);
    for(i=;i<=n;i++)
        for(j=;j<=n+;j++)
            scanf("%lf",&a[i][j]);
    for(i=;i<=n;i++)
    {
        nmax=;
        for(j=i;j<=n;j++)
            if(fabs(a[j][i])>fabs(nmax))
            {
                nmax=a[j][i];
                nx=j;
            }
        if(dcmp(nmax)==)
        {
            puts("No Solution");
            return ;
        }
        swap_line(nx,i);
        for(j=i+;j<=n+;j++)    a[i][j]/=nmax;
        for(j=;j<=n;j++)
            if(j!=i)
                for(k=i+;k<=n+;k++)
                    a[j][k]-=a[j][i]*a[i][k];
    }
    for(i=;i<=n;i++)    printf("%.2lf\n",a[i][n+]);
    return ;
}