BZOJ 2683 简单题 cdq分治+树状数组

题意:链接

**方法:**cdq分治+树状数组

解析:

首先对于这道题，看了范围之后，二维的数据结构是显然不能过的，于是我们可能会考虑把一维排序之后另一位上数据结构什么的，然而cdq分治却能够很好的体现它的作用。

首先，对于每一个询问求和，显然是x在它左边的并且出现时间在它之前的所有的change对他可能会有影响。

我们按照x第一关键字,y第二关键字,操作第三关键字来排序所有的询问，然后在cdq的时候，每次递归处理左半区间，按照x动态的将y这一列的值加到树状数组里，来更新右半边的所有询问，注意这里的树状数组是需要清的，也就是每次cdq都是采用不同的树状数组。

另:这题神坑，数组开小不是RE而是WA

数组范围是要开到200000*4的，因为对于每个1操作，是一个操作，而2操作根据容斥原理，是4个操作。所以复杂度大概是多少呢O(nlog^2(n))cdq+树状数组。n<=800000

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 800010
#define Q q[tot]
using namespace std;
struct node
{
    int opt,x,y,A,no,belong;
}q[N],nq[N];
int n;
int jd,tot,T;
int ans[N];
int c[N];
int lowbit(int x){return x&(-x);}
int getsum(int x){int ret=0;while(x){ret+=c[x],x-=lowbit(x);}return ret;}
void update(int x,int v){while(x<=n){c[x]+=v;x+=lowbit(x);}}
int cmp(node a,node b)
{
    if(a.x==b.x&&a.y==b.y)return a.opt<b.opt;
    if(a.x==b.x)return a.y<b.y;
    return a.x<b.x;
}
void sov(int l,int r)
{
    if(l==r)return;
    int mid=(l+r)>>1,l1=l,l2=mid+1;
    for(int i=l;i<=r;i++)
    {
        if(q[i].no<=mid&&q[i].opt==1)update(q[i].y,q[i].A);
        if(q[i].no>mid&&q[i].opt==2)ans[q[i].belong]+=getsum(q[i].y);
    }
    for(int i=l;i<=r;i++)
    {
        if(q[i].no<=mid&&q[i].opt==1)update(q[i].y,-q[i].A);
    }
    l1=l,l2=mid+1;
    for(int i=l;i<=r;i++){if(q[i].no<=mid)nq[l1++]=q[i];else nq[l2++]=q[i];}
    for(int i=l;i<=r;i++)q[i]=nq[i];
    sov(l,mid),sov(mid+1,r);
}
int main()
{
    scanf("%d",&n);
    while(scanf("%d",&jd)&&jd^3)
    {
        if(jd==1){int x,y,a;scanf("%d%d%d",&x,&y,&a);q[++tot].opt=1,Q.x=x,Q.y=y,Q.A=a,Q.no=tot;}
        else
        {
            int x1,x2,y1,y2;
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            q[++tot].opt=2,Q.x=x1-1,Q.y=y1-1,Q.A=1,Q.no=tot,Q.belong=++T;
            q[++tot].opt=2,Q.x=x1-1,Q.y=y2,Q.A=-1,Q.no=tot,Q.belong=T;
            q[++tot].opt=2,Q.x=x2,Q.y=y1-1,Q.A=-1,Q.no=tot,Q.belong=T;
            q[++tot].opt=2,Q.x=x2,Q.y=y2,Q.A=1,Q.no=tot,Q.belong=T;
        }
    }
    sort(q+1,q+1+tot,cmp);
    sov(1,tot);
    for(int i=1;i<=T;i++)printf("%d\n",ans[i]);
}

版权声明：本文为博主原创文章，未经博主允许不得转载。