[bzoj5343][Ctsc2018]混合果汁——二分答案+主席树

题目描述:

小 R 热衷于做黑暗料理,尤其是混合果汁。

商店里有 nnn 种果汁,编号为 0,1,⋯,n−10,1,\cdots,n-10,1,,n1iii 号果汁的美味度是 did_idi,每升价格为 pip_ipi。小 R 在制作混合果汁时,还有一些特殊的规定,即在一瓶混合果汁中,iii 号果汁最多只能添加 lil_ili 升。

现在有 mmm 个小朋友过来找小 R 要混合果汁喝,他们都希望小 R 用商店里的果汁制作成一瓶混合果汁。其中,第 jjj 个小朋友希望他得到的混合果汁总价格不大于 gjg_jgj,体积不小于 LjL_jLj。在上述这些限制条件下,小朋友们还希望混合果汁的美味度尽可能地高,一瓶混合果汁的美味度等于所有参与混合的果汁的美味度的最小值。请你计算每个小朋友能喝到的最美味的混合果汁的美味度。

思路:

答案很明显具有单调性,考虑对于每一个询问二分答案。
二分一个最小值 d i d_i di后把所有美味值 ≥ d i \geq d_i di的果汁按照价格从小到大排序,然后贪心地一直选择,直到sum ≥ L i \geq L_i Li为之。
不难发现排序的过程复杂度较高,每一次是先按照 d i d_i di排序之后再按照 p i p_i pi排序,于是我们可以直接按照 d i d_i di的顺序以 p i p_i pi为下标构建一颗主席树,每一次直接在主席树上面二分查找就好了。

#include<bits/stdc++.h>
 
#define REP(i,a,b) for(int i=a,i##_end_=b;i<=i##_end_;++i)
#define DREP(i,a,b) for(int i=a,i##_end_=b;i>=i##_end_;--i)
#define MREP(i,x) for(int i=beg[x],v;v=to[i],i;i=las[i])
#define debug(x) cout<<#x<<"="<<x<<endl
#define fi first
#define se second
#define mk make_pair
#define pb push_back
typedef long long ll;
 
using namespace std;
 
void File(){
    freopen("bzoj5343.in","r",stdin);
    freopen("bzoj5343.out","w",stdout);
}
 
template<typename T>void read(T &_){
    T __=0,mul=1; char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')mul=-1;
        ch=getchar();
    }
    while(isdigit(ch))__=(__<<1)+(__<<3)+(ch^'0'),ch=getchar();
    _=__*mul;
}
 
const int maxn=1e5+10;
const int N=1e5;
int n,m,root[maxn];
struct node{
    ll d,p,w;
    bool operator < (const node & tt) const {
        if(d!=tt.d)return d>tt.d;
        return p<tt.p;
    }
}a[maxn];
 
struct Chairman_Tree{
#define mid ((l+r)>>1)
    int cnt;
    struct poi{
        int lc,rc;
        ll sum,cost;
    }b[maxn<<5];
    void insert(int &rt,int l,int r,ll p,ll x){
        int now=++cnt; b[now]=b[rt]; rt=now;
        if(l==r)b[rt].sum+=x,b[rt].cost=b[rt].sum*l;
        else{
            if(p<=mid)insert(b[rt].lc,l,mid,p,x);
            else insert(b[rt].rc,mid+1,r,p,x);
            b[rt].sum=b[b[rt].lc].sum+b[b[rt].rc].sum;
            b[rt].cost=b[b[rt].lc].cost+b[b[rt].rc].cost;
        }
    }
    bool judge(int rt,int l,int r,ll x,ll g){
        if(x>=0 && g<0)return false;
        if(x==0 && g>=0)return true;
        if(l==r)return b[rt].sum<x ? false : (g-x*l>=0);
        ll s=b[b[rt].lc].sum,c=b[b[rt].lc].cost;
        if(s<=x)return judge(b[rt].rc,mid+1,r,x-s,g-c);
        else return judge(b[rt].lc,l,mid,x,g);
    }
#undef mid
}T;
 
void init(){
    read(n),read(m);
    REP(i,1,n)read(a[i].d),read(a[i].p),read(a[i].w);
    sort(a+1,a+n+1);
    REP(i,1,n){
        root[i]=root[i-1];
        T.insert(root[i],1,N,a[i].p,a[i].w);
    }
    //REP(i,1,n)printf("%lld %lld %lld\n",a[i].d,a[i].p,a[i].w);
}
 
ll calc(ll x,ll g){
    int l=1,r=n;
    while(l<r){
        int mid=(l+r)>>1;
        if(T.judge(root[mid],1,N,x,g))r=mid;
        else l=mid+1;
    }
    if(!T.judge(root[l],1,N,x,g))return -1;
    else return a[l].d;
}
 
void work(){
    ll g,x;
    REP(i,1,m){
        read(g),read(x);
        printf("%lld\n",calc(x,g));
    }
}
 
int main(){
    //File();
    init();
    work();
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值