网上的写法都什么辣鸡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;
}