bzoj3262

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

分析：

一眼就是三维偏序~~
首先按x值排序，将l~r分成l~mid和mid+1~r此时左边对右边有贡献。
再将l~mid和mid+1~r区间分别按y值排序。用两个标记$l_1$和$l_2$分别从l和mid+1开始扫，不断将l~mid中的z值加入树状数组中，当$a[l_1].y>a[l_2].y$时统计答案并将$l_2$向右移。

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cctype>
using namespace std;
typedef long long LL;
const int M=262144;
inline int rd(){
    int x;
    scanf("%d",&x);
    return x;
}

int tcas;
int mx;
int n;
int ans[M],dp[M];
int tms[M];
int c[M];
int T;

struct node{
    int x,y,z,id,dp;
    void input(int ii){x=rd(),y=rd(),z=rd(),id=ii;dp=0;}
}a[M],t[M];

bool operator ==(node x,node y){
    return x.x==y.x&&x.y==y.y&&x.z==y.z;
}

bool cmp(node x,node y){
    if(x.x==y.x&&x.y==y.y&&x.z==y.z) return x.id<y.id;
    if(x.x==y.x&&x.y==y.y)return x.z<y.z;
    if(x.x==y.x)return x.y<y.y;
    return x.x<y.x;
}

bool cmpy(node x,node y){
    if(x.y==y.y&&x.z==y.z) return x.id<y.id;
    if(x.y==y.y)return x.z<y.z;
    return x.y<y.y;
}

inline int lb(int x){return x&(-x);}

void ins(int x,int d){
    for(;x<=mx;x+=lb(x))c[x]+=d;
}

int get(int x){
    int res=0;
    for(;x>0;x-=lb(x))res+=c[x];
    return res;
}

void calc(int l,int r){
    if(l==r)return;
    int i,j,mid=l+r>>1;
    for (i=l;i<=r;i++)
        t[i]=a[i];
    sort(a+l,a+mid+1,cmpy);
    sort(a+mid+1,a+r+1,cmpy);
    for(i=l,j=mid+1;j<=r;j++){
        for(;a[i].y<=a[j].y&&i<=mid;i++)
            ins(a[i].z,1);
        dp[a[j].id]+=get(a[j].z);
    }
    for (j=l;j<i;j++) ins(a[j].z,-1);
    for (i=l;i<=r;i++)
        a[i]=t[i];
    calc(l,mid);
    calc(mid+1,r);
}

int siz[M];

int main(){
    freopen("a.in","r",stdin);
    int i;
    tcas=1;
    while(tcas--){
        n=rd(),mx=rd();
        for(i=1;i<=n;i++) a[i].input(i);
        sort(a+1,a+n+1,cmp);
        calc(1,n);
        sort(a+1,a+n+1,cmp);
        for(i=n-1;i>0;i--) if(a[i]==a[i+1]) dp[i]=dp[i+1];
        for(i=1;i<=n;i++) ans[a[i].id]=dp[i];
        for (i=1;i<=n;i++) siz[ans[i]]++;
        for (i=0;i<n;i++) printf("%d\n",siz[i]);
    }
    return 0;
}