zoukankan      html  css  js  c++  java
  • 【BZOJ4004】【JLOI2015】装备购买

    Description

      
      脸哥最近在玩一款神奇的游戏,这个游戏里有 n 件装备,每件装备有 m 个属性,用向量zi(aj ,.....,am) 表示
      
    (1 <= i <= n; 1 <= j <= m),每个装备需要花费 ci,现在脸哥想买一些装备,但是脸哥很穷,所以总是盘算着
      
    怎样才能花尽量少的钱买尽量多的装备。对于脸哥来说,如果一件装备的属性能用购买的其他装备组合出(也就是
      
    说脸哥可以利用手上的这些装备组合出这件装备的效果),那么这件装备就没有买的必要了。严格的定义是,如果
      
    脸哥买了 zi1,.....zip这 p 件装备,那么对于任意待决定的 zh,不存在 b1,....,bp 使得 b1zi1 + ... + bpzi
      
    p = zh(b 是实数),那么脸哥就会买 zh,否则 zh 对脸哥就是无用的了,自然不必购买。举个例子,z1 =(1; 2;
      
    3);z2 =(3; 4; 5);zh =(2; 3; 4),b1 =1/2,b2 =1/2,就有 b1z1 + b2z2 = zh,那么如果脸哥买了 z1 和 z2
      
    就不会再买 zh 了。脸哥想要在买下最多数量的装备的情况下花最少的钱,你能帮他算一下吗?
      
      
      

    Input

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

    Output

      
      一行两个数,第一个数表示能够购买的最多装备数量,第二个数表示在购买最多数量的装备的情况下的最小花费
      

    Sample Input

      
      3 3
      1 2 3
      3 4 5
      2 3 4
      1 1 2
      

    Sample Output

      
      2 2
      
      
      

    Solution

      
      还是有关那个带权线性基的问题,但是我还未膜拜拟阵证明,所以就先写着吧。
      
      我们发现,这个新物品购买与否的判定机制,是看用已购买的物品能否线性表示出新物品。这不就是线性基的那一套插入、判定能否线性表示的做法吗?
      
      再考虑到那个贪心的方法——我们要使选中的物品价值尽可能小,那么依照套路我们按物品价值从小到大依次考虑,尝试加入线性基:能,那就买;中途消成0,就说明它可以被线性表示,不买。
      
      这题的实现不再能直接套用异或线性基的写法了,而是要用实数来模拟线性基。记线性基的向量为(b_1,b_2,...,b_m),其实每一个(b_i)类比于异或线性基中的每一个数,只不过异或线性基还把向量状压了。线性基用矩阵表示就长这样:

    [egin{bmatrix} b_{m,1}&b_{m,2} &b_{m,3} &... &b_{m,m}\ ...&...&...&...&...\ b_{3,1}&b_{3,2} &b_{3,3} &... &b_{3,m}\ b_{2,1}&b_{2,2} &b_{2,3} &... &b_{2,m}\ b_{1,1}&b_{1,2} &b_{1,3} &... &b_{1,m}\ end{bmatrix} ]

    ​  尝试加入一个向量(a=[a_1,a_2,...,a_m])时,我们从(b_m)开始,一行一行地循环到(b_1)。考虑(b_i)时,如果(a_i=0),continue;如果(b_{i,i}=0),我们把整个(b_i)赋值为(a);如果(b_{i,i} e0),我们用(b_i)消去(a_{i,i}),继续往下走。类比异或线性基的插入方法,就能很好地理解普通线性基的插入方法。
      
      由于BZOJ新加了三组毒瘤数据,因此实数要用long double保证精度。
      
       
      

    Code

      

    #include <cstdio>
    #include <algorithm>
    #include <cmath>
    using namespace std;
    typedef long double ld;
    const int N=505;
    const ld EPS=1e-8;
    int n,m;
    ld lb[N][N];
    struct Equip{
    	ld a[N];
    	int cost;
    }s[N];
    int use,sum;
    inline bool cmp(const Equip &x,const Equip &y){return x.cost<y.cost;}
    int main(){
    	scanf("%d%d",&n,&m);
    	for(int i=1;i<=n;i++)
    		for(int j=1,x;j<=m;j++) scanf("%llf",&s[i].a[j]);
    	for(int i=1;i<=n;i++) scanf("%d",&s[i].cost);
    	sort(s+1,s+1+n,cmp);
    	for(int x=1;x<=n;x++)
    		for(int i=1;i<=m;i++)
    			if(fabs(s[x].a[i])>=EPS){
    				if(fabs(lb[i][i])<EPS){
    					for(int j=i;j<=m;j++) lb[i][j]=s[x].a[j];
    					use++;
    					sum+=s[x].cost;
    					break;
    				}
    				for(int j=m;j>=i;j--)
    					s[x].a[j]-=lb[i][j]*(s[x].a[i]/lb[i][i]);
    			}
    	printf("%d %d
    ",use,sum);
    	return 0;
    }
    
  • 相关阅读:
    jmeter Transaction Controller、Throughput Controller 控制器
    图书系统的简易代码
    模板简单梳理
    自制过滤器
    自制URL转换器
    URL简单梳理
    window下的Django入门
    十一、常用内建模块
    九、进程与线程
    八、错误、调试与测试
  • 原文地址:https://www.cnblogs.com/RogerDTZ/p/9215566.html
Copyright © 2011-2022 走看看