csp模拟 骨粉【二分答案】【主席树】

传送门

存在离线做法但我不会。

因为数据范围可以看出是个二分答案。

对于每次mid，验证是O(n)的，复杂度总共n^2logn。

我们从简化验证的复杂度入手。

因为农作物会自然生长，所以实际上验证的高度是mid+t（询问的t）

令v=mid+t；

对于每个作物，要骨粉次，我们来简化这个式子。

我们可以简单地通过二分计算出前面的和，对于余数，我们用数据结构维护。

使用主席树可以做到区间查询。

完了。

#include<bits/stdc++.h>
using namespace std;
#define in read()
#define int long long
int in{
	int cnt=0,f=1;char ch=0;
	while(!isdigit(ch)){
		ch=getchar();if(ch=='-')f=-1;
	}
	while(isdigit(ch)){
		cnt=cnt*10+ch-48;ch=getchar();
	}return cnt*f;
} 
const int N=100003,M=3200003;
int L[M],R[M],size[M],rt[N],sumsuf[N],a[N],b[N],bcnt;
int n,m,x;int id;
int build(int l,int r){
	int root=++id;
	if(l<r){
		int mid=(l+r)>>1;
		L[root]=build(l,mid);R[root]=build(mid+1,r);
	}return root;
}
int insert(int pre,int l,int r,int key){
	int root=++id;
	L[root]=L[pre];R[root]=R[pre];size[root]=size[pre]+1;
	if(l<r){
		int mid=(l+r)>>1;
		if(key<=mid)L[root]=insert(L[pre],l,mid,key);
		else R[root]=insert(R[pre],mid+1,r,key);
	}return root;
}
int query(int u,int l,int r,int ql,int qr){
	if(!u)return 0;
	if(r<ql||l>qr)return 0;
	if(ql<=l&&r<=qr)return size[u];int mid=(l+r)>>1;
	return query(L[u],l,mid,ql,qr)+query(R[u],mid+1,r,ql,qr);
}
bool check(int height,int use){
	int l=1,r=n;
//	cout<<"height "<<height<<endl;
	while(l<r){
		int mid=(l+r)>>1;
		if(height>=a[mid])l=mid+1;
		else r=mid;
	}
	int pos=l;//cout<<"pos"<<" "<<pos<<endl;
	if(height>=a[n]&&pos==n)return true;
	int gu=height/x;
	int ans=sumsuf[pos]-(n-pos+1)*gu;int res=height%x;
	if(ans<=use)return true;if(res==0)return false;
	l=0,r=bcnt;
	while(l<r){
		int mid=(l+r+1)>>1;
		if(res>=b[mid])l=mid;
		else r=mid-1;
	}
	if(l){int x=query(rt[pos],1,bcnt,1,l);ans-=x;}
//	cout<<"x "<<x<<l<<" "<<bcnt<<endl;
	return ans<=use;
}int maxx;
signed main(){
	n=in;m=in;x=in;
	for(register int i=1;i<=n;i++)maxx=max(maxx,a[i]=in);sort(a+1,a+n+1);
	for(register int i=n;i>=1;i--)sumsuf[i]=(sumsuf[i+1]+(a[i]-1)/x+1);
	for(register int i=1;i<=n;i++)b[++bcnt]=a[i]%x;
	sort(b+1,b+bcnt+1);bcnt=unique(b+1,b+bcnt+1)-b-1;
	rt[n+1]=build(1,bcnt);
	for(register int i=n;i>=1;i--){
		int gu=a[i]%x;if(!gu)continue;
		gu=lower_bound(b+1,b+bcnt+1,a[i]%x)-b;
//		cout<<gu<<" ";
		rt[i]=insert(rt[i+1],1,bcnt,gu);
	}//cout<<endl;
	while(m--){
		int que=in;
		int l=0,r=maxx;
		while(l<r){
			int mid=(l+r)>>1;
			if(check(mid+que,que))r=mid;
			else l=mid+1;
		}cout<<l<<'\n';
	}
	return 0;
}