【题解】洛谷P3265（同bzoj4004）[JSOI2015]装备购买 高斯消元

题目链接

题目描述

输入输出格式

输入格式：
第一行两个数 n;m。接下来 n 行，每行 m 个数，其中第 i 行描述装备 i 的各项属性值。接下来一行 n 个数，其中 ci 表示购买第 i 件装备的花费。

输出格式：
一行两个数，第一个数表示能够购买的最多装备数量，第二个数表示在购买最多数量的装备的情况下的最小花费

输入输出样例

输入样例#1：
3 3
1 2 3
3 4 5
2 3 4
1 1 2
输出样例#1：
2 2

说明

如题目中描述，选择装备 1 装备 2，装备 1 装备 3，装备 2 装备 3 均可，但选择装备 1 和装备 2 的花费最小，为 2。

对于 100% 的数据, 1 <= n;m <= 500; 0 <= aj <= 1000。

---

把n件装备看成n个长度为m的向量，求线性空间的基。
将属性看成系数矩阵，高斯消元求出该矩阵的秩。

#include<cstdio>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
#define _rep(i,a,b) for(int i=(a);i<=(b);i++)
typedef long double lb;
const double eps=1e-5;
const int N=510;
int n,m,cnt,cost;
struct node{
    lb buff[N];//属性 
    int c;//价格 
    bool operator<(const node&rhs)const{
    return c<rhs.c;}
}arm[N];
int z[N];//基底 
int main()
{
    //freopen("in.txt","r",stdin);
    std::ios::sync_with_stdio(false);//这玩意和freopen一起用就炸了 
    cin>>n>>m;
    _rep(i,1,n)_rep(j,1,m)cin>>arm[i].buff[j];
    _rep(i,1,n)cin>>arm[i].c;
    sort(arm+1,arm+n+1);
    _rep(i,1,n)_rep(j,1,m)if(fabs(arm[i].buff[j])>eps)
    {
        if(!z[j]){z[j]=i;cnt++;cost+=arm[i].c;break;}
        else
        {
            lb p=arm[i].buff[j]/arm[z[j]].buff[j];
            _rep(k,j,m)arm[i].buff[k]-=p*arm[z[j]].buff[k];
        }
    }
    printf("%d %d\n",cnt,cost);
}

总结

高斯消元求线性基