网上的写法都什么辣鸡bzoj4004

网上说要什么拟阵证明巴拉巴拉。显然我是不会的，但是排个序我还是知道的。然后看了网上的程序一脸懵逼，讲道理高斯消元不是应该最外重循环是列数嘛。你这是怎么搞搞的，学长说这是把每行加进去判断之前的行能不能凑出这一行，这么讲好像也有些道理。但是这么丑的高斯消元我是拒绝的。。。所以我就想写个正常一点的。想想其实也不难。我只要不改变这些行的位置不就好了，不用swap操作就好了。写了一发交上去wa了我非常不服，然后网上标称一交也是wa，看来是精度问题，加了一下精度结果两种写法也都ac了，下面给出我的写法。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define eps 1e-8
#define maxn 510
#define Abs(x) ((x) > 0 ? (x) : (-(x)))
#define forup(i,a,b) for(int i=(a);i<=(b);i++)
#define fordown(i,a,b) for(int i=(a);i>=(b);i--)
using namespace std;

struct yts
{
long double b[maxn];
int val;
}a[maxn];

int n,ans,m,num;
int vis[maxn];

bool cmp(yts x,yts y)
{
return x.val<y.val;
}

int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
cin>>a[i].b[j];
for (int i=1;i<=n;i++) scanf("%d",&a[i].val);
sort(a+1,a+n+1,cmp);
/*for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
if (fabs(a[i].b[j])>eps)
{
if (!vis[j])
{
vis[j]=i;
ans+=a[i].val;num++;
break;
}
else
{
long double t=(long double)a[i].b[j]/(long double)a[vis[j]].b[j];
for (int k=j;k<=m;k++)
a[i].b[k]-=t*a[vis[j]].b[k];
}
}*/
forup(i,1,m)
{int b=0;
forup(j,1,n)
{ if(fabs(a[j].b[i])>eps&&!vis[j])
{ b=j;vis[j]=1;
break;
}
}
if(b!=0) {num++;ans=ans+a[b].val;}
else continue;
forup(j,1,n)
{if(j!=b&&fabs(a[j].b[i])>eps)
{ long double tmp=(long double)a[j].b[i]/(long double)a[b].b[i];
forup(k,1,m)
a[j].b[k]=a[j].b[k]-tmp*a[b].b[k];
}

}
}
printf("%d %d\n",num,ans);
return 0;
}