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

时间:2023-03-09 13:39:33
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[],ba[];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=;a[len].f=;
}
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=;i<=n;i++)
scanf("%d%d%d",&x,&y,&z),ins(x,y,z);
sort(a+,a+len+,cmp);
n=;
for(int i=;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[];
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=;
while(x)
{
ans+=s[x];
x-=lowbit(x);
}
return ans;
} void cdq(int l,int r)
{
if(l==r)return ;
int mid=(l+r)/;
cdq(l,mid);cdq(mid+,r); int i=l,j=mid+,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[];
int main()
{
scanf("%d%d",&n,&m);
qtt(); memset(s,,sizeof(s));
cdq(,n); memset(ans,,sizeof(ans));
for(int i=;i<=n;i++)ans[a[i].f+a[i].tot-]+=a[i].tot;
for(int i=;i<len;i++)printf("%d\n",ans[i]);
return ;
}