[JLOI2015]装备购买 （高斯消元）

Posted 2021-02-15 812-xiao-wen

[JLOI2015]装备购买

---

---

(solution:)

首先这道题的题面已经非常清晰的告诉我们这就是线性空间高斯消元的一道题（可以用某些装备来表示另一件装备，这已经不能再明显了），只是这道题要求我们求出花费最小的方案，这个我们可以直接贪心，如果有多个装备可以匹配某一个属性，那么我们选价值最小的那一个（这个可以用线性空间的定义证明），价值大的尽量往后再选（选到最后剩下一些价格大的不买即可）

---

---

(code:)

#include<iostream>
#include<cstdio>
#include<iomanip>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<ctime>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>

#define ll long long
#define db double
#define inf 0x7fffffff
#define rg register int

using namespace std;

const db cha=1e-4;

struct su{
    db a[505];
    int v;
    inline bool operator <(su y){return v<y.v;}
}g[505];

int n,m,tot,ans;
int b[505];

inline int qr(){
    char ch;
    while((ch=getchar())<'0'||ch>'9');
    int res=ch^48;
    while((ch=getchar())>='0'&&ch<='9')
        res=res*10+(ch^48);
    return res;
}

int main(){
    //freopen(".in","r",stdin);
    //freopen(".out","w",stdout);
    n=qr(); m=qr();
    for(rg i=1;i<=n;++i)
        for(rg j=1;j<=m;++j)
            g[i].a[j]=qr();
    for(rg i=1;i<=n;++i)
        g[i].v=qr();
    sort(g+1,g+n+1);
    for(rg i=1;i<=n;++i){
        for(rg j=1;j<=m;++j){
            if(fabs(g[i].a[j])>cha){
                if(!b[j]){
                    b[j]=i;++tot;
                    ans+=g[i].v;
                    break;
                }else{
                    db d=g[i].a[j]/g[b[j]].a[j];
                    for(rg k=j;k<=m;++k)
                        g[i].a[k]-=g[b[j]].a[k]*d;
                }
            }
        }
    }printf("%d %d
",tot,ans);
    return 0;
}