【BZOJ4004】【JLOI2015】装备购买 线性基

链接：

#include <stdio.h>
int main()
{
    puts("转载请注明出处[vmurder]谢谢");
    puts("网址：blog.youkuaiyun.com/vmurder/article/details/45171081");
}

题解：

基本裸题，不写题解了。看这篇博客吧。
【BZOJ3105】【cqoi2013】新Nim游戏 线性基神马的
http://blog.youkuaiyun.com/vmurder/article/details/42560691
这篇博客没写拟阵部分，只写了线性基部分，刚刚好。
哎考试时时间都喂T2了。这种裸题竟然没去AC它哭死啊啊啊啊。

对了，据BK说这道题可能卡精度，然后就需要一些黑科技。
但是因为我根本没有加特技，裸写一发直接过了，就不介绍了。

代码：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 505
#define eps 1e-5
#define inf 0x3f3f3f3f
using namespace std;
struct Eli
{
    int val;
    double x[N];
    void read(int m){for(int i=1;i<=m;i++)scanf("%lf",&x[i]);}
    void v(){scanf("%d",&val);}
    bool operator < (const Eli &A)const
    {return val<A.val;}
}a[N];
int n,m;
int ins[N],num,ans;
void getins()
{
    int i,j,k;
    double rate;
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++)if(fabs(a[i].x[j])>eps)
        {
            if(ins[j])
            {
                rate=a[i].x[j]/a[ins[j]].x[j];
                for(k=j;k<=m;k++)a[i].x[k]-=rate*a[ins[j]].x[k];
            }
            else {ins[j]=i,ans+=a[i].val,num++;break;}
        }
    }
}
int main()
{
    freopen("test.in","r",stdin);

    int i,j,k;

    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)a[i].read(m);
    for(i=1;i<=n;i++)a[i].v();
    sort(a+1,a+n+1);
    getins();
    printf("%d %d\n",num,ans);
    return 0;
}