【BZOJ 3262】陌上开花 CDQ分治

最新推荐文章于 2020-07-12 14:47:32 发布

pbihao

最新推荐文章于 2020-07-12 14:47:32 发布

阅读量299

点赞数

CC 4.0 BY-SA版权

分类专栏： ac之路 bzoj CDQ分治

本文链接：https://blog.youkuaiyun.com/pbihao/article/details/54909834

ac之路同时被 3 个专栏收录

303 篇文章

订阅专栏

bzoj

236 篇文章

订阅专栏

CDQ分治

3 篇文章

订阅专栏

本文介绍了一种基于离线算法解决区间查询问题的方法。通过使用树状数组进行区间更新和查询操作，配合分治思想优化算法效率。适用于处理大量查询请求并能够有效地给出答案。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

不想说了，感觉自己是傻逼

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define maxn 200021
using namespace std;
int n,c[maxn],K,ans[maxn],z[maxn],cc[maxn],bl[maxn];
void update(int x,int add){
	while(x<=K){c[x]+=add;x+=x&-x;}
}
int query(int x){
	int ans=0;
	while(x>0){ans+=c[x];x-=x&-x;}
	return ans;
}
struct ques{
	int x,y,z,id;
	bool operator<(const ques& b)const{
		if(z!=b.z)return z<b.z;
		if(x!=b.x)return x<b.x;
		return y<b.y;
	}
}q[maxn],p[maxn];
bool cmp(const ques& a,const ques& b){
	if(a.x!=b.x)return a.x<b.x;
	if(a.y!=b.y)return a.y<b.y;
	return a.z<b.z;
}

void solve(int l,int r){
	if(l>=r)return;
	for(int i=l;i<=r;i++)z[i]=q[i].x;
	sort(z+l,z+r+1);
	int mid=l+r>>1,x=z[mid];
	
	int t1=l,t2=mid+1,cnt=0,z=0;
	for(int i=l;i<=r;i++)if(q[i].x<x)z++;
	for(int i=l;i<=r;i++){
		if(q[i].x<x)p[t1]=q[i],bl[i]=t1++,z--;
		else if(q[i].x==x&&z+t1<=mid)p[t1]=q[i],bl[i]=t1++;
		else p[t2]=q[i],bl[i]=t2++;
	}
	
	for(int j,i=l;i<=r;i=j){
		cnt=0;
		for(j=i+1;j<=r;j++){
			if(q[j].x!=q[j-1].x||q[j].y!=q[j-1].y||q[j].z!=q[j-1].z)break;
			else if(q[j].x==x&&bl[j]>mid)cnt++;
		}if(j>i+1&&q[i].x==x&&bl[i]>mid)cnt++;
		for(int k=i;k<j;k++)if(q[k].x==x&&bl[k]<=mid)ans[q[k].id]+=cnt;
	}
	
	for(int i=l;i<=r;i++){
		if(bl[i]<=mid)update(q[i].y,1);
		else ans[q[i].id]+=query(q[i].y);
	}
	for(int i=l;i<=r;i++){if(bl[i]<=mid)update(q[i].y,-1);}
	
	for(int i=l;i<=r;i++){
		if(bl[i]>mid&&q[i].x==x)update(q[i].y,1);
		else if(q[i].x==x)ans[q[i].id]+=query(q[i].y);
	}
	for(int i=l;i<=r;i++){if(bl[i]>mid&&q[i].x==x)update(q[i].y,-1);}
	for(int i=l;i<=r;i++)q[i]=p[i];
	
	solve(l,mid),solve(mid+1,r);
}

int main(){
	scanf("%d%d",&n,&K);
	for(int i=1;i<=n;i++)scanf("%d%d%d",&q[i].x,&q[i].y,&q[i].z),q[i].id=i;
	sort(q+1,q+1+n);
	solve(1,n);
	for(int i=1;i<=n;i++)cc[ans[i]]++;
	for(int i=0;i<n;i++)printf("%d\n",cc[i]);
	return 0;
}