【CF799E】Aquarium decoration

题意:有n个物品,小A和小B各自喜欢其中的某些物品,一件物品可能既被小A喜欢又被小B喜欢,也可能既不被小A喜欢又不被小B喜欢。每个物品都有一个价格$c_i$,让你选出其中的m件物品,满足小A和小B都至少喜欢其中的k件,且总价格最小。

$n\le 2\times 10^5,c_i\le 10^9$

题解:我们把物品分成四部分:0,A,B,AB,并分别排序,然后我们枚举AB中选了i个。接着我们可以二分:在A,B,0中选择的价格最大的物品的价格mid,假设A中有a个数小于等于mid,则我们要从中选max(a,k)个,同理B中选max(b,k)个,C中选m-i-max(a,k)-max(b,k)个。发现这个显然是可以二分的。如果把二分放到线段树上的话,则时间复杂度是$O(n\log n)$。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define lson x<<1
#define rson x<<1|1
using namespace std;
typedef long long ll;
const int maxn=200010;
ll ans;
int n,m,k,na,nb,nc,nd,N;
int A[maxn],B[maxn],C[maxn],D[maxn],v[maxn],tag[maxn],p[maxn];
int ca[maxn<<2],cb[maxn<<2],cd[maxn<<2],ta[maxn],tb[maxn],td[maxn];
ll sa[maxn],sb[maxn],sd[maxn],ref[maxn]; bool cmp(const int &a,const int &b)
{
return v[a]<v[b];
}
void insert(int l,int r,int x,int a,int b,int c)
{
if(c==0) ca[x]=b;
if(c==1) cb[x]=b;
if(c==2) cd[x]=b;
if(l==r) return ;
int mid=(l+r)>>1;
if(a<=mid) insert(l,mid,lson,a,b,c);
else insert(mid+1,r,rson,a,b,c);
}
ll query(int l,int r,int x,int a)
{
if(l==r)
return sa[max(ca[x],a)]+sb[max(cb[x],a)]+sd[cd[x]]-ref[l]*(max(ca[x],a)+max(cb[x],a)+cd[x]+k-a-m);
int mid=(l+r)>>1,t=max(ca[lson],a)+max(cb[lson],a)+cd[lson]+k-a;
if(t>=m) return query(l,mid,lson,a);
return query(mid+1,r,rson,a);
}
void build(int l,int r,int x)
{
if(l==r)
{
ca[x]=ta[l],cb[x]=tb[l],cd[x]=td[l];
return ;
}
int mid=(l+r)>>1;
build(l,mid,lson),build(mid+1,r,rson);
ca[x]=ca[rson],cb[x]=cb[rson],cd[x]=cd[rson];
}
inline int rd()
{
int ret=0,f=1; char gc=getchar();
while(gc<'0'||gc>'9') {if(gc=='-') f=-f; gc=getchar();}
while(gc>='0'&&gc<='9') ret=ret*10+(gc^'0'),gc=getchar();
return ret*f;
}
int main()
{
n=rd(),m=rd(),k=rd();
int i,a,b;
for(i=1;i<=n;i++) v[i]=rd(),p[i]=i;
sort(p+1,p+n+1,cmp);
for(i=1;i<=n;i++)
{
if(i==1||v[p[i]]>ref[N]) ref[++N]=v[p[i]];
v[p[i]]=N;
}
a=rd();
for(i=1;i<=a;i++) b=rd(),tag[b]|=1;
a=rd();
for(i=1;i<=a;i++) b=rd(),tag[b]|=2;
for(i=1;i<=n;i++)
{
if(tag[i]==0) D[++nd]=ref[v[i]],td[v[i]]++;
if(tag[i]==1) A[++na]=ref[v[i]],ta[v[i]]++;
if(tag[i]==2) B[++nb]=ref[v[i]],tb[v[i]]++;
if(tag[i]==3) C[++nc]=ref[v[i]];
}
sort(A+1,A+na+1),sort(B+1,B+nb+1),sort(C+1,C+nc+1),sort(D+1,D+nd+1);
for(i=1;i<=N;i++) ta[i]+=ta[i-1],tb[i]+=tb[i-1],td[i]+=td[i-1];
for(i=1;i<=na;i++) sa[i]=sa[i-1]+A[i];
for(i=1;i<=nb;i++) sb[i]=sb[i-1]+B[i];
for(i=1;i<=nd;i++) sd[i]=sd[i-1]+D[i];
ll sum=0; ans=1ll<<60;
build(0,N,1);
for(i=0;i<=min(nc,m);i++)
{
sum+=C[i];
if(min(na,nb)>=k-i&&2*k-i<=m&&i+na+nb+nd>=m)
ans=min(ans,sum+query(0,N,1,k-i));
}
if(ans==1ll<<60) puts("-1");
else printf("%lld",ans);
return 0;
}

最新文章

  1. jsp:中文乱码解决
  2. gameui-for-phaser-js更新到最新版本
  3. Java多线程时内存模型
  4. poj 3254(状态压缩基础题)
  5. 导入flash 逐帧动画
  6. CodeForces 609C Load Balancing
  7. Scrapy爬取西刺代理ip流程
  8. FFmpeg的H.264解码器源代码简单分析:环路滤波(Loop Filter)部分
  9. SpringBoot自动配置源码调试
  10. Android横竖屏切换的生命周期
  11. PHP isset 和 array_key_exists 对比
  12. 在IIS上启用Gzip压缩(HTTP压缩)
  13. 导出可运行jar包
  14. 设计模式-装饰者模式(Decorator Pattern)
  15. 如何在PHP项目中使用phinx进行数据迁移和建表
  16. hive 锁
  17. Chrome 插件推荐
  18. Visual Studio 下nuget命令的使用
  19. centos7嵌入式环境搭建
  20. Linux du查询文件大小

热门文章

  1. Thinkphp5笔记六:公共模块common的使用
  2. pytest集成Allure Report
  3. RMAN:简单的duplicate创建新数据库
  4. kettle教程二
  5. 关于UIGestureRecognizerState
  6. ios学习之UIWebView网页视图
  7. JS有趣的单线程
  8. jinja2主要语法
  9. 嵌入式系统之ubootENV环境变量
  10. ISD9160学习笔记05_ISD9160语音识别代码分析