bzoj4004 [JLOI2015]装备购买（贪心+高斯消元）

求权值和最小的极大线性无关集。
我们贪心地按权值从小往大加，如果线性相关了就删掉。
由于bzoj貌似卡精度，于是可以写模意义下的…
$O(nm^2)$

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define N 510
#define mod 1000000007
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,a[N][N],c[N],ans=0,b[N],aa[N],cnt=0;
inline bool cmp(int x,int y){return c[x]<c[y];}
inline int ksm(int x,int k){
    int res=1;for(;k;k>>=1,x=(ll)x*x%mod) if(k&1) res=(ll)res*x%mod;return res;
}
inline void dec(int &x,int y){x-=y;if(x<0) x+=mod;}
inline bool ins(int x){
    for(int i=1;i<=m;++i){
        if(!a[x][i]) continue;
        if(!b[i]){b[i]=x;return 1;}
        int t=(ll)a[x][i]*ksm(a[b[i]][i],mod-2)%mod;
        for(int j=i;j<=m;++j) dec(a[x][j],(ll)a[b[i]][j]*t%mod);
    }return 0;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j) a[i][j]=read();
    for(int i=1;i<=n;++i) aa[i]=i,c[i]=read();
    sort(aa+1,aa+n+1,cmp);
    for(int i=1;i<=n;++i) if(ins(aa[i])) ans+=c[aa[i]],++cnt;
    printf("%d %d\n",cnt,ans);
    return 0;
}