[KCOJ20170214]又一个背包

题目描述 Description

小W要去军训了！由于军训基地是封闭的，小W在军训期间将无法离开军训基地。所以他没有办法出去买他最爱吃的零食。万般无奈的小W只好事先买好他爱吃的零食，装在背包里带入军训基地。市场里的零食琳琅满目，纵然小W想把他们都带走，然而小W的背包只有有限的容量。带哪些零食好呢？小W给每种零食都评估了一个他的喜爱程度。另外，由于市场的零食都是散称售卖的，所以一种零食可以只买一部分带走。现在小W希望能挑选出最合适的购买方案，使得所购买的零食能装进背包且喜爱程度总分最大。因为零食的种类实在是太多了，所以小W找到了聪明的你，希望你能尽快（军训的车队即将出发，时间紧迫）告诉他购买方案。

输入描述 Input Description

第一行两个整数n和w，分别表示有n种零食种类和背包总容量。接下来n行，每行两个整数ci和vi，分别表示第i种食品占的体积和小W的喜爱程度。

输出描述 Output Description

一行一个数字，表示最优购买方案下，所能获得的最大喜爱程度总和。保留3位小数。

样例输入 Sample Input

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

样例输出 Sample Output

11.000

数据范围及提示 Data Size & Hint

n <= 2000000，w <= 2*10^9，0<=ci <= 100，0<=vi<=100

一道部分背包问题。这道题按性价比排序然后贪心取肯定是没有问题的，但是复杂度是O(n log n)有点大，像我们学校OJ对于n=2000000的肯定是跑不下来的。所以我们要O(n)做这道题。乍一想，发现完全没有思路，应该有一点灵感就是类似于O(n)求第K大数一样，一个用的是O(n log n)排序，一个是nth_element O(n)做的。于是我们就想分治做这道题，把问题分为更小的子问题。算出每个物品的性价比之后，我们nth_element求出中位数，这样这个数列的左边就变成比中位数小的，右边就变成比中位数大的。然后我们暴力扫一遍右边的价值和，判断一下，然后就好办了。复杂度分析的话，这个东西1+1/2+1/4+1/8+...+1/2^n是小于2的，所以类似的话，复杂度就是O（n）。本题还有一个坑点，就是空间有可能是0！

 #include<iostream>

 #include<algorithm>

 #include<cstdio>

 #include<cstdlib>

 #include<cstring>

 #include<queue>

 using namespace std;

 typedef long long LL;

 inline int read()

 {

     int x=,f=;char c=getchar();

     while(!isdigit(c)){if(c=='-')f=-;c=getchar();}

     while(isdigit(c)){x=x*+c-'';c=getchar();}

     return x*f;

 }

 const int maxn=;

 struct Item

 {

     double w,v,s;

     bool operator < (const Item &t)const  {return s<t.s;}

 }a[maxn];

 int n,ans;

 double V;

 double solve(int l,int r,double val)//val是当前背包容量

 {

     if(l>r)return ;

     if(l==r)

     {

         if(val<a[l].v)return val*a[l].s;

         else return a[l].w;

     }

     int mid=(l+r)/;

     double sum=,sum2=;

     nth_element(a+l,a+mid,a+r+);

     for(int i=mid+;i<=r;i++)sum+=a[i].w,sum2+=a[i].v;

     //printf("%d %d %d %lf %lf %lf\n",l,r,mid,sum,sum2,val);

     if(val>sum2)return sum+solve(l,mid,val-sum2);

     if(val==sum2)return sum;

     if(val<sum2)return solve(mid+,r,val);

 }

 int main()

 {

     n=read();V=(double)read();

     int i=;

     while(i<=n)

     {

          a[i].v=(double)read();a[i].w=(double)read();

          if(a[i].v==)ans+=a[i].w,n--;

          else a[i].s=a[i].w/a[i].v,i++;

     }

     printf("%.3f\n",ans+solve(,n,V));

     return ;

 }

巴特西

[KCOJ20170214]又一个背包

最新文章

热门文章