题目大意:
n 件装备每件装备有m 件属性以及价格ci,属性可以用向量zi 表示,
现在想要购买尽量多的装备使得装备的属性向量线性无关,求最小代价
solution:
属性要求满足拟阵条件
贪心维护一个线性基即可
O(n^2*m)
拟阵这东西比较复杂···我们只需要知道它是满足贪心的就行了
线性基这东西也比较复杂···其实根本用不着,就是个线性表出,普通gauss就行
贪心的按照代价从小到大排序
然后枚举,枚举到i,如果i的所有属性都被消成0
则i可以被其他向量表出,否则将其加到答案里
然后把后面的都消掉
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define maxn 505
using namespace std;
int n,m,cnt,ans;
const double eps=1e-6;
struct qwqwq{
int c; long double a[maxn];//居然要用long double??????
}eq[maxn];
inline int rd(){
int x=0,f=1;char c=' ';
while(c<'0' || c>'9') {if(c=='-')f=-1;c=getchar();}
while(c<='9' && c>='0') x=x*10+c-'0',c=getchar();
return x*f;
}
bool cmp(qwqwq x,qwqwq y) {return x.c<y.c;}
void Gauss(){
for(int i=1;i<=n;i++){
bool flg=false;
for(int j=i;j<=m;j++) if(fabs(eq[i].a[j])>eps) {flg=true;break;}
if(!flg) continue;
cnt++; ans+=eq[i].c;
for(int j=i+1;j<=n;j++){
long double rate=eq[j].a[i]/eq[i].a[i];
for(int k=i;k<=m;k++) eq[j].a[k]-=rate*eq[i].a[k];
}
}
}
int main(){
n=rd(); m=rd();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++) eq[i].a[j]=rd();
for(int i=1;i<=n;i++) eq[i].c=rd();
sort(eq+1,eq+n+1,cmp);
Gauss();
printf("%d %d\n",cnt,ans);//输出格式不对wa了一次···
return 0;
}