[CTSC2018]混合果汁

一、题目

点此看题

二、解法

每个询问的二分方式是一样的，所以可以用整体二分。

维护询问$[ql,qr]$的答案是$[l,r]$，我们二分一个答案$mid$（美味度的下标，从大到小），然后把$[1,mid]$加入权值线段树，看取够$L$需要花多少钱，和拥有的钱相比较，看分到左边还是右边。

我们不能把$[1,mid]$直接加入权值线段树，考虑每次指针会从一个边界开始，移动到$mid$，然后再移动回去，不难发现这样的复杂度是$O(n\log n)$，所以总复杂度$O(n{\log }^{2}n)$

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
#define int long long
const int M = 100005;
int read()
{	
	int x=0,flag=1;char c;
	while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
	while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return x*flag;
}
int n,m,up,cur,ans[M],s1[4*M],s2[4*M];
struct node
{
	int d,p,l;
	bool operator < (const node &b) const
	{
		return d>b.d;
	}
}a[M];
struct data
{
	int g,l,id;
}p[M],L[M],R[M];
void pushup(int i)
{
	s1[i]=s1[i<<1]+s1[i<<1|1];
	s2[i]=s2[i<<1]+s2[i<<1|1];
}
void ins(int i,int l,int r,int id,int v)
{
	if(l==r)
	{
		s1[i]+=v;s2[i]+=l*v; 
		return ;
	}
	int mid=(l+r)>>1;
	if(mid>=id) ins(i<<1,l,mid,id,v);
	else ins(i<<1|1,mid+1,r,id,v);
	pushup(i);
}
int ask(int i,int l,int r,int v)
{
	if(!v) return 0;
	if(l==r) return l*v;
	int mid=(l+r)>>1;
	if(s1[i<<1]>=v) return ask(i<<1,l,mid,v);
	return s2[i<<1]+ask(i<<1|1,mid+1,r,v-s1[i<<1]);
}
void cdq(int l,int r,int ql,int qr)
{
	if(l>r || ql>qr) return ;
	if(l==r)
	{
		for(int i=ql;i<=qr;i++)
			ans[p[i].id]=a[l].d;
		return ;
	}
	int mid=(l+r)>>1,tl=0,tr=0;
	while(cur<mid) cur++,ins(1,1,up,a[cur].p,a[cur].l);
	while(cur>mid) ins(1,1,up,a[cur].p,-a[cur].l),cur--;
	for(int i=ql;i<=qr;i++)
	{
		if(p[i].l>s1[1]) R[++tr]=p[i];
		else if(ask(1,1,up,p[i].l)<=p[i].g) L[++tl]=p[i];
		else R[++tr]=p[i];
	}
	for(int i=1;i<=tl;i++) p[ql+i-1]=L[i];
	for(int i=1;i<=tr;i++) p[ql+tl+i-1]=R[i];
	cdq(l,mid,ql,ql+tl-1);
	cdq(mid+1,r,ql+tl,qr); 
}
signed main()
{
	n=read();m=read();
	for(int i=1;i<=n;i++)
	{
		a[i]=node{read(),read(),read()};
		up=max(up,a[i].p);
	}
	a[++n]=node{-1,0,1e18};
	sort(a+1,a+1+n);
	for(int i=1;i<=m;i++)
		p[i]=data{read(),read(),i};
	cdq(1,n,1,m);
	for(int i=1;i<=m;i++)
		printf("%lld\n",ans[i]);
}