【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
3Description
有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
0
1
0
1
0
0
1
HINT

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

Source

看见三维偏序肯定想cdq分治
把询问重复的去掉之后开始分治
第一维排序第二维分治第三维建树状数组然后对每个询问求前缀和就是比他小的询问的数量.注意这个题对于相同的询问是要计入评级的.
因为给定的k很小还省去了离散化的功夫,直接建20W的Bit就行了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 100100
#define lowbit(x) (x&(-x))
using namespace std;
int n,k,cnt,top;
int c[MAXN<<1];
int ans[MAXN];
void in(int &x)
{
    char ch=getchar();x=0;
    while (!(ch>='0'&&ch<='9')) ch=getchar();
    while (ch>='0'&&ch<='9')    x=x*10+ch-'0',ch=getchar();
}
void add(int x,int delta)
{
    for (int i=x;i<=k;i+=lowbit(i)) c[i]+=delta;
}
int query(int x)
{
    int ret=0;
    for (int i=x;i;i-=lowbit(i))    ret+=c[i];
    return ret;
}
struct Query
{
    int a,b,c,sum,num;
    bool operator <(const Query& A)const
    {
        return b==A.b?c<A.c:b<A.b;
    }
    bool operator !=(const Query& A)const
    {
        return a!=A.a||b!=A.b||c!=A.c;
    }
}s[MAXN],q[MAXN];
bool cmp(Query a,Query b)
{
    if (a.a==b.a&&a.b==b.b) return a.c<b.c;
    if (a.a==b.a)   return a.b<b.b;
    return a.a<b.a;
}
void solve(int l,int r)
{
    int mid=(l+r)>>1;
    if (l==r)   return;
    solve(l,mid);solve(mid+1,r);
    sort(q+l,q+mid+1);sort(q+mid+1,q+r+1);
    int tp1=l;
    for (int i=mid+1;i<=r;i++)
    {
        while (tp1<=mid&&q[tp1].b<=q[i].b)  add(q[tp1].c,q[tp1].num),tp1++;
        q[i].sum+=query(q[i].c);
    }
    for (int i=l;i<tp1;i++) add(q[i].c,-q[i].num);
}
int main()
{
    in(n);in(k);
    for (int i=1;i<=n;i++)  in(s[i].a),in(s[i].b),in(s[i].c);
    sort(s+1,s+n+1,cmp);
    for (int i=1;i<=n;i++)
    {
        cnt++;
        if (s[i]!=s[i+1])   q[++top]=s[i],q[top].num=cnt,cnt=0;
    }
    solve(1,top);
    for (int i=1;i<=top;i++)    ans[q[i].sum+q[i].num-1]+=q[i].num;
    for (int i=0;i<n;i++)   printf("%d\n",ans[i]);
}