BZOJ 3262 陌上花开

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;}