[JLOI2015]装备购买

版权声明：欢迎转载（标记出处），写得不好还请多指教 https://blog.csdn.net/quan_tum/article/details/82668638

传送门

此题卡精度。

按权值排序从小到大插入，若存在向量线性相关，则不插入。

用类似线性基的方法维护向量集。

#include<bits/stdc++.h>
using namespace std;
#define ld long double
#define getchar()(p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
inline int read(){
    int x=0,f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=getchar()) x=(x+(x<<2)<<1)+c-48;
    return x*f;
}
const ld eps=1e-6;
struct A{ld p[505];int w;}a[505];
bool cmp(A x,A y){return x.w<y.w;}
int n,m,ele,ans;ld p[505][505];
int main(){
    n=read(),m=read();
    for(int i=1;i<=n;++i) for(int j=1;j<=m;++j) a[i].p[j]=read();
    for(int i=1;i<=n;++i) a[i].w=read();
    sort(a+1,a+n+1,cmp);
    for(int t=1;t<=n;++t) for(int i=1;i<=m;++i){
        if(fabs(a[t].p[i])<eps) continue;
        if(fabs(p[i][i])<eps){
            for(int j=i;j<=n;++j) p[i][j]=a[t].p[j];
            ++ele;ans+=a[t].w;break;
        }for(int j=m;j>=i;--j) a[t].p[j]-=p[i][j]*a[t].p[i]/p[i][i];
    }return !printf("%d %d",ele,ans);
}