[旧文章补档]BZOJ 1176 Mokia CDQ分治

CDQ分治论文例题。。
论文戳我

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstring>
#include <cmath>
#define MAXN 200005
#define N 2000005
#define LL long long
#define INF 1000000005
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
using namespace std;
const double eps = 1e-8;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,s,c[N],ans[MAXN];
inline int lowbit(int i){return i&(-i);}
struct t{int opt,to,add,x,y,f,pos;}p[MAXN],tmp[MAXN];
void add(int x,int v)
{
    for(int i=x;i<=n;i+=lowbit(i)) c[i]+=v;
}
int query(int x)
{
    int t=0;
    for(int i=x;i;i-=lowbit(i)) t+=c[i];
    return t;
}
void addquery(int x1,int y1,int x2,int y2)
{
    ans[0]++;
    p[++m].pos=m,p[m].to=ans[0],p[m].x=x1-1,p[m].y=y1-1,p[m].opt=0,p[m].f=1;
    p[++m].pos=m,p[m].to=ans[0],p[m].x=x2,p[m].y=y2,p[m].opt=0,p[m].f=1;
    p[++m].pos=m,p[m].to=ans[0],p[m].x=x1-1,p[m].y=y2,p[m].opt=0,p[m].f=-1;
    p[++m].pos=m,p[m].to=ans[0],p[m].x=x2,p[m].y=y1-1,p[m].opt=0,p[m].f=-1;
}
void solve(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(p[i].pos<=mid&&p[i].opt) add(p[i].y,p[i].add);
        if(p[i].pos>=mid+1&&(!p[i].opt)) ans[p[i].to]+=p[i].f*query(p[i].y);
    }
    for(int i=l;i<=r;i++)
        if(p[i].pos<=mid&&p[i].opt) add(p[i].y,-p[i].add);
    for(int i=l;i<=r;i++)
        if(p[i].pos<=mid) tmp[l1++]=p[i];
        else tmp[l2++]=p[i];
    for(int i=l;i<=r;i++) p[i]=tmp[i];
    solve(l,mid),solve(mid+1,r);
}
bool cmp(t a,t 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;
    else return a.x<b.x;
}
int main()
{
    s=read(),n=read();
    while(1)
    {
        int t=read();
        if(t==1)
            p[++m].pos=m,p[m].x=read(),p[m].y=read(),p[m].opt=1,p[m].add=read();
        else if(t==2)
        {
            int x1=read(),y1=read(),x2=read(),y2=read();
            addquery(x1,y1,x2,y2);
        }
        else break;
    }
    sort(p+1,p+1+m,cmp);
    solve(1,m);
    for(int i=1;i<=ans[0];i++)
        printf("%d\n",ans[i]);
    return 0;
}