BZOJ 3262 陌上花开_bzoj陌上花开-CSDN博客

本文链接：https://blog.csdn.net/SenyeLicone/article/details/58582151

本文介绍了一种处理大规模三维数据集的有效算法——CDQ分治算法，通过使用树状数组进行维护，解决了花的美丽程度评级问题。该算法首先对花的属性进行排序，然后运用分治思想在不同维度上递归处理，最终统计出每个级别的花的数量。

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

Description

有n朵花，每朵花有三个属性：花形(s)、颜色(c)、气味(m)，又三个整数表示。现要对每朵花评级，一朵花的级别是它拥有的美丽能超过的花的数量。定义一朵花A比另一朵花B要美丽，当且仅当Sa>=Sb,Ca>=Cb,Ma>=Mb。显然，两朵花可能有同样的属性。需要统计出评出每个等级的花的数量。

Input

第一行为N,K (1 <= N <= 100,000, 1 <= K <= 200,000 ), 分别表示花的数量和最大属性值。

以下N行，每行三个整数si, ci, mi (1 <= si, ci, mi <= K)，表示第i朵花的属性

Output

包含N行，分别表示评级为0...N-1的每级花的数量。

Sample Input

10 3
3 3 3
2 3 3
2 3 1
3 1 1
3 1 2
1 3 1
1 1 2
1 2 2
1 3 2
1 2 1

Sample Output

3
1
3
0
1
0
1
0
0
1

HINT

1 <= N <= 100,000, 1 <= K <= 200,000

Source

树套树 CDQ分治

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

CDQ分治~

三维，x排序，yCDQ分治，z树状数组维护。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

int n,m,tot,t,c[200001],tim[200001],ans[100001];

struct node{
	int x,y,z,cnt,ans;
}a[100001],b[100001];

bool operator < (node u,node v)
{
	if(u.x==v.x) return u.y==v.y ? u.z<v.z:u.y<v.y;
	return u.x<v.x;
}

bool cmp(node u,node v)
{
	if(u.y==v.y) return u.z==v.z ? u.x<v.x:u.z<v.z;
	return u.y<v.y;
}

int read()
{
	int totnum=0,f=1;char ch=getchar();
	while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0' && ch<='9') {totnum=(totnum<<1)+(totnum<<3)+ch-'0';ch=getchar();}
	return totnum*f;
}

void add(int u,int v)
{
	for(;u<=m;u+=u&(-u))
	{
		if(tim[u]!=t) c[u]=0;
		tim[u]=t;c[u]+=v;
	}
}

int cal(int u)
{
	int now=0;
	for(;u;u-=u&(-u)) if(tim[u]==t) now+=c[u];
	return now;
}

void findd(int l,int r)
{
	int mid=(l+r)>>1,l1=l,l2=mid+1,j;
	if(l==r)
	{
		a[l].ans+=a[l].cnt-1;return;
	}
	for(int i=l;i<=r;i++)
	  if(a[i].x<=mid) b[l1++]=a[i];
	  else b[l2++]=a[i];
	memcpy(a+l,b+l,sizeof(a[0])*(r-l+1));
	findd(l,mid);t++;j=l;
	for(int i=mid+1;i<=r;i++)
	{
		for(;j<=mid && a[j].y<=a[i].y;j++) add(a[j].z,a[j].cnt);
		a[i].ans+=cal(a[i].z);
	}
	findd(mid+1,r);
	l1=l;l2=mid+1;
	for(int i=l;i<=r;i++)
	  if((cmp(a[l1],a[l2]) || l2>r) && l1<=mid) b[i]=a[l1++];
	  else b[i]=a[l2++];
	memcpy(a+l,b+l,sizeof(a[0])*(r-l+1));
}

int main()
{
	n=read();m=read();
	for(int i=1;i<=n;i++) a[i].x=read(),a[i].y=read(),a[i].z=read(),a[i].cnt++;
	sort(a+1,a+n+1);
	for(int i=1;i<=n;i++)
	  if(i==1 || a[i-1]<a[i]) a[++tot]=a[i];
	  else a[tot].cnt++;
	for(int i=1;i<=tot;i++) a[i].x=i;
	sort(a+1,a+tot+1,cmp);
	findd(1,tot);
	for(int i=1;i<=tot;i++) ans[a[i].ans]+=a[i].cnt;
	for(int i=0;i<n;i++) printf("%d\n",ans[i]);
	return 0;
}