[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>?rx:(-rx);
 }

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

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