bzoj 3262(cdq分治)

3262: 陌上花开

Time Limit: 20 Sec   Memory Limit: 256 MB
Submit: 1414   Solved: 634
[ 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分治。步骤对第一维进行排序，然后进行分治，

分别处理两边，然后处理右边对左边的影响，用树状数组维护。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
int n,k,cnt;
struct ss
 {
    int x,y,z,dui,size;
 }q[110000],g[110000];
int zhi[210000];
int ans[210000];
int zhong[210000];
 
inline int read()
{
    char y; int x=0,f=1; y=getchar();
    while (y<'0' || y>'9') {if (y=='-')f=-1; y=getchar();}
    while (y>='0' && y<='9') {x=x*10+int(y)-48; y=getchar();}
    return x*f;
}
 
bool cog(ss xg,ss yg)
 {
    if (xg.y!=yg.y)
     return xg.y<yg.y;
    return xg.z<yg.z;
 }
 
void change(int o,int sug)
 {
    while (o<=k)
     {
        zhi[o]+=sug; o+=o&-o;
     }
 }
 
int query(int o) 
{
    int sug=0;
    while (o>0)
     {
        sug+=zhi[o]; o-=o&-o;
     }
    return sug;
}
 
void solve(int l,int r)
 {
    if (l==r) return;
    int mid=(l+r)/2;
    solve(l,mid); solve(mid+1,r);
    sort(g+l,g+mid+1,cog); sort(g+mid+1,g+r+1,cog);
    int tail=l; 
    for (int i=mid+1;i<=r;++i)
     {
      while (tail<mid+1 && g[i].y>=g[tail].y)
       {
        change(g[tail].z,g[tail].size);
        ++tail;
       }
      ans[g[i].dui]+=query(g[i].z);
     }
    for (int i=l;i<=tail-1;++i)
     {
        change(g[i].z,-g[i].size);
     }
 }
 
bool cmp(ss xg,ss yg)
 {
    if (xg.x!=yg.x)
     return xg.x<yg.x;
    if (xg.y!=yg.y)
     return xg.y<yg.y;
    return xg.z<yg.z;     
 }
 
int main()
{
    n=read(); k=read();
    for (int i=1;i<=n;++i)
     {
      q[i].x=read(); q[i].y=read(); q[i].z=read(); q[i].dui=i;
     }
    sort(q+1,q+n+1,cmp); cnt=0;
    for (int i=1;i<=n;++i)
     {
      if (q[i].x!=q[i-1].x || q[i].y!=q[i-1].y || q[i].z!=q[i-1].z) ++cnt,g[cnt]=q[i],g[cnt].dui=cnt;
++g[cnt].size;
}
    solve(1,cnt);
    for (int i=1;i<=n;++i)
     {
        zhong[ans[g[i].dui]+g[i].size-1]+=g[i].size;
     }
    for (int i=0;i<=n-1;++i)
     {
        printf("%d\n",zhong[i]);
     }
}