bzoj3262: 陌上花开(cdq分治+树状数组)

3262: 陌上花开

题目：传送门 

---

题解：

　　 %%%cdq分治 

　　 很强大的一个暴力...感觉比分块高级多了

　　 这道题目就是一个十分经典的三维偏序的例题：

　　 一维直接暴力排序x

　　 二维用csq维护y

 　　三维用树状数组来搞

　　 最后ans处理答案，注意：全部值相等，相互之间也算自己更加漂亮

---

代码：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
struct flower
{
    int x,y,z,f,tot;
}a[210000],ba[210000];int n,m,len;
void ins(int x,int y,int z)
{
    len++;
    a[len].x=x;a[len].y=y;a[len].z=z;
    a[len].tot=1;a[len].f=0;
}
bool cmp(flower n1,flower n2)
{
    if(n1.x!=n2.x)return n1.x<n2.x;
    if(n1.y!=n2.y)return n1.y<n2.y;
    return n1.z<n2.z;
}
void qtt()
{
    int x,y,z;
    for(int i=1;i<=n;i++)
        scanf("%d%d%d",&x,&y,&z),ins(x,y,z);
    sort(a+1,a+len+1,cmp);
    n=1;
    for(int i=2;i<=len;i++)
        if(a[n].x==a[i].x && a[n].y==a[i].y && a[n].z==a[i].z)
            a[n].tot++;
        else 
            a[++n]=a[i];
}

int s[210000];
int lowbit(int x){return x&-x;}
void add(int x,int k)
{
    while(x<=m)
    {
        s[x]+=k;
        x+=lowbit(x);
    }
}
int getsum(int x)
{
    int ans=0;
    while(x)
    {
        ans+=s[x];
        x-=lowbit(x);
    }
    return ans;
}

void cdq(int l,int r)
{
    if(l==r)return ;
    int mid=(l+r)/2;
    cdq(l,mid);cdq(mid+1,r);
    
    int i=l,j=mid+1,p=l;
    while(i<=mid && j<=r)
    {
        if(a[i].y<=a[j].y)add(a[i].z,a[i].tot),ba[p++]=a[i++];
        else a[j].f+=getsum(a[j].z),ba[p++]=a[j++];
    }
    while(i<=mid)add(a[i].z,a[i].tot),ba[p++]=a[i++];
    while(j<=r)a[j].f+=getsum(a[j].z),ba[p++]=a[j++];
    
    for(int i=l;i<=mid;i++)add(a[i].z,-a[i].tot);
    
    for(int i=l;i<=r;i++)a[i]=ba[i];
}
int ans[110000];
int main()
{
    scanf("%d%d",&n,&m);
    qtt();
    
    memset(s,0,sizeof(s));
    cdq(1,n);
    
    memset(ans,0,sizeof(ans));
    for(int i=1;i<=n;i++)ans[a[i].f+a[i].tot-1]+=a[i].tot;
    for(int i=0;i<len;i++)printf("%d
",ans[i]);
    return 0;
}