BZOJ 3262 陌上花开

3262: 陌上花开

Time Limit: 20 Sec   Memory Limit: 256 MB
Submit: 1393   Solved: 626
[ 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

【题意如上】

【解题方法】典型三维偏序，CDQ板子题了。但是要注意记录某个数出现的次数。

【AC 代码】

//BZOJ 3262
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn=100010;
const int mod=1e9+7;
struct Point{
    int x,y,z;
    int cnt,sum;
    void read(){
        scanf("%d%d%d",&x,&y,&z);
        cnt=sum=0;
    }
    bool operator<(const Point &p)const{
        if(x!=p.x) return x<p.x;
        if(y!=p.y) return y<p.y;
        else return z<p.z;
    }
    bool operator==(const Point &p)const{
        return x==p.x&&y==p.y&&z==p.z;
    }
}p[maxn];
int lowbit(int x)
{
    return x&-x;
}
struct BIT{
    int n;
    int c[maxn];
    void init(int _n){
        n=_n;
        memset(c,0,sizeof(c));
    }
    void update(int i,int v){
        while(i<=n){
            c[i]+=v;
            i+=lowbit(i);
        }
    }
    int queryans(int i){
        int res=0;
        while(i){
            res+=c[i];
            i-=lowbit(i);
        }
        return res;
    }
}bit;
bool cmpy(Point a,Point b){
    return a.y<b.y;
}
void CDC(int l,int r)
{
    if(l==r) return ;
    int m=(l+r)/2;
    CDQ(l,m);
    CDQ(m+1,r);
    sort(p+l,p+m+1,cmpy);
    sort(p+m+1,p+r+1,cmpy);
    int j=l;
    for(int i=m+1; i<=r; i++){
        for(;j<=m&&p[j].y<=p[i].y; j++) bit.update(p[j].z,p[j].cnt);
        p[i].sum+=bit.queryans(p[i].z);
    }
    for(int i=l; i<j; i++) bit.update(p[i].z,-p[i].cnt);
}
int n,k,m;
int A[maxn];
int main()
{
    while(scanf("%d%d",&n,&k)!=EOF){
        bit.init(k);
        for(int i=1; i<=n; i++) p[i].read();
        sort(p+1,p+n+1);
        int m=0,j;
        for(int i=1; i<=n; i++){
            for(j=i; j<=n&&p[i]==p[j]; j++){
                ++p[i].cnt;
            }
            p[++m]=p[i];
            i=j-1;
        }
        CDQ(1,m);
        memset(A,0,sizeof(A));
        for(int i=1; i<=m; i++){
            A[p[i].sum+p[i].cnt]+=p[i].cnt;
        }
        for(int i=1; i<=n; i++){
            printf("%d\n",A[i]);
        }
    }
}