[BZOJ4004][JLOI2015]装备购买（线性基+贪心）

你是我命中意外，偶然闯入的小鹿

---

题目：

我是超链接

题解：

实际上可以发现线性基的过程和消上三角矩阵的过程有些相似，都是为每一位选定“关键向量”，并且保证关键向量下面的所有向量都没有这一位。那么可以发现，因为每一位都只有一个“关键向量”，那么如果从求出的线性基里面任意选择两个进行异或运算，设为ai和aj，不妨设i>j，那么由求解过程可得ai的第i位是它最高位的1，而aj的这一位是0。那么ai和aj异或就怎么也得不到0。推广到多个向量也是如此。也就是说，通过上面的过程，我们能够求出一组线性无关的向量，也就是给定向量的线性基。

注意到当这个过程结束以后，有些向量会被异或成0。那么这些向量一定可以被原向量组中其它某些向量线性表示，也就是说原向量组中存在一个组合使得它们的异或值为0。

而且找个很大的质数%一下比较好

代码：

#include <cstdio>
#include <algorithm>
#define LL long long
using namespace std;
const int mod=1e9+7;
struct hh{LL a[505];int val;}zb[505];
int b[505];
int cmp(hh a,hh b){return a.val<b.val;}
LL ksm(LL a,int k)
{
    LL ans=1;
    for (;k;k>>=1,a=a*a%mod)
      if (k&1) ans=ans*a%mod;
    return ans;
}
int main()
{
    int n,m;scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
      for (int j=1;j<=m;j++) scanf("%d",&zb[i].a[j]);
    for (int i=1;i<=n;i++) scanf("%d",&zb[i].val);   
    sort(zb+1,zb+n+1,cmp);
    int ans=0,cnt=0;
    for (int i=1;i<=n;i++)
      for (int j=1;j<=m;j++)
        if (zb[i].a[j])
        {
            if (!b[j])
            {
                b[j]=i;
                cnt++;ans+=zb[i].val;
                break;
            }
            else 
            {
                int t=zb[i].a[j]*ksm(zb[b[j]].a[j],mod-2)%mod;
                for (int k=j;k<=m;k++)
                  zb[i].a[k]=(zb[i].a[k]-t*zb[b[j]].a[k]%mod+mod)%mod;
            }
        }
    printf("%d %d",cnt,ans);
}