[SDOI2006] 线性方程组

洛谷 P2455 传送门

刚开始写了个消成上三角的，结果狂wa。

后来经过研究发现，消成上三角那种不能直接判断无解或无穷多解，需要其它的操作。

所以干脆学了个消成对角线的，写了一发A了。

其实两种消元的方法都差不多，就当存个板子吧。

#include<cstdio>

void swap(double &x,double &y)
{
    double t=x;x=y;y=t;
}

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

int n;
double a[55][55];
double x[55];
const double eps=0.00001;

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