bzoj 3262 陌上花开



3262: 陌上花开

Time Limit: 20 Sec  Memory Limit:
256 MB
Submit: 2142  Solved: 956

[Submit][Status][Discuss]

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分治还是半懂不懂...(智商感人qwqqq)

对于三个属性，y排序，x做CDQ分治，z做树状数组处理。

具体实现中，我们在处理[l,r]时，把整个区间扫一遍，然后把x<=mid的放在左边，否则放在右边，由于y已经排好序，所以左边区间y依旧单增，右边也是。

所以对于右边的每一个花，考虑左半边的贡献，就要用一个指针线性扫描一下，如果左边指针l1对应的y<=右边指针l2对应的y，那么可以把l1.z加入树状数组中统计方案。

注意tim数组，由于树状数组要储存很多次分治信息，一直memset很慢，所以对每次分治打一个标记来清除树状数组中的信息。

不得不说好神啊qwqqq

【代码】

//bzoj 3262 陌上花开
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=200005;
int n,m,T,tot;
int c[mxn],ans[mxn],tim[mxn];
struct node
{
int x,y,z,p,cnt; //p:评级
}a[mxn],tmp[mxn];
inline bool operator < (node x,node y)
{
if(x.x!=y.x) return x.x<y.x;
if(x.y!=y.y) return x.y<y.y;
return x.z<y.z;
}
inline bool comp(node x,node y)
{
if(x.y!=y.y) return x.y<y.y;
if(x.z!=y.z) return x.z<y.z;
return x.x<y.x;
}
inline int lowbit(int x) {return x&-x;}
inline void add(int x,int y)
{
for(int i=x;i<=m;i+=lowbit(i))
if(tim[i]!=T) tim[i]=T,c[i]=y;
else tim[i]=T,c[i]+=y;
}
inline int getsum(int x)
{
int sum=0;
for(int i=x;i>=1;i-=lowbit(i))
if(tim[i]==T) sum+=c[i];
return sum;
}
inline void CDQ(int l,int r)
{
int i,j,mid=l+r>>1;
if(l==r)
{
a[l].p+=a[l].cnt-1;
return;
}
int l1=l,l2=mid+1;
fo(i,l,r)
{
if(a[i].x<=mid) tmp[l1++]=a[i];
else tmp[l2++]=a[i];
}
memcpy(a+l,tmp+l,sizeof(a[0])*(r-l+1));
j=l,T++;
fo(i,mid+1,r)
{
for(j;j<=mid && a[j].y<=a[i].y;j++)
add(a[j].z,a[j].cnt);
a[i].p+=getsum(a[i].z);
}
CDQ(l,mid),CDQ(mid+1,r);
}
int main()
{
int i,j;
scanf("%d%d",&n,&m);
fo(i,1,n)
scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z),a[i].cnt=1;
sort(a+1,a+n+1);
fo(i,1,n)
{
if(i==1 || a[i-1]<a[i])
a[++tot]=a[i];
else
a[tot].cnt++;
}
fo(i,1,tot) a[i].x=i;
sort(a+1,a+tot+1,comp);
CDQ(1,tot);
fo(i,1,tot)
ans[a[i].p]+=a[i].cnt;
fo(i,0,n-1)
printf("%d\n",ans[i]);
return 0;
}