bzoj 1176 [Balkan2007]Mokia 【CDQ分治】

W过大，很难在线维护，考虑离线算法

给每个操作加一个时间属性n，显然，对于n=i的询问，对它有影响的修改只在n<i中，所以可以CDQ（因为是按时间序读进来的，所以不用排序了

对于统计矩形和，可以使用二维前缀和的思想，即只需要统计四个点即可

这样就转化成了三维偏序问题，只是询问和修改要分开处理。

（初始值不用管，输出的时候加上就行

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=2000005;
int n,s,qu,tot,ans[N],t[N],p[N],cnt,q[N];
struct qwe
{
    int o,x,y,v,n;
    qwe(){}
    qwe(int a1,int a2,int a3,int a4,int a5)
    {
        o=a1,x=a2,y=a3,v=a4,n=a5;
    }
}a[N];
bool cmp(const qwe &a,const qwe &b)
{
    return (a.x==b.x&&a.y==b.y&&a.o<b.o)||(a.x==b.x&&a.y<b.y)||(a.x<b.x);
}
int read()
{
    int r=0,f=1;
    char p=getchar();
    while(p>'9'||p<'0')
    {
        if(p=='-')
            f=-1;
        p=getchar();
    }
    while(p>='0'&&p<='9')
    {
        r=r*10+p-48;
        p=getchar();
    }
    return r*f;
}
int lb(int x)
{
    return x&(-x);
}
void update(int x,int v)
{
    for(int i=x;i<=n;i+=lb(i))
        t[i]+=v;
}
int ques(int x)
{
    int r=0;
    for(int i=x;i>=1;i-=lb(i))
        r+=t[i];
    return r;
}
void cdq(int l,int r)
{
    if(l==r)
        return;
    int mid=(l+r)/2,tmp=0;
    cdq(l,mid);
    cdq(mid+1,r);
    sort(a+l,a+mid+1,cmp);
    sort(a+mid+1,a+r+1,cmp);
    int i=l,j=mid+1;
    while(j<=r)
    {
        for(;a[i].o==2&&i<=mid;i++);
        for(;a[j].o==1&&j<=r;j++);
        if(i<=mid&&a[i].x<=a[j].x)
            update(a[i].y,a[i].v),i++,tmp=i-1;
        else if(j<=r)
            ans[a[j].n]+=ques(a[j].y),j++;
    }
    for(int t=l;t<=tmp;t++)
        if(a[t].o==1)
            update(a[t].y,-a[t].v);
}
int main()
{
    s=read();n=read();
    while(1)
    {
        qu=read();
        if(qu==3)
            break;
        if(qu==1)
        {
            int x=read(),y=read(),z=read();
            a[++tot]=qwe(1,x,y,z,tot);
        }
        else
        {
            int x1=read(),y1=read(),x2=read(),y2=read();
            p[++cnt]=tot;
            q[cnt]=(x2-x1+1)*(y2-y1+1)*s;
            a[++tot]=qwe(2,x1-1,y1-1,0,tot);
            a[++tot]=qwe(2,x2,y2,0,tot);
            a[++tot]=qwe(2,x1-1,y2,0,tot);
            a[++tot]=qwe(2,x2,y1-1,0,tot);
        }
    }
    cdq(1,tot);
    for(int i=1;i<=cnt;i++)
        printf("%d
",q[i]+ans[p[i]+1]+ans[p[i]+2]-ans[p[i]+3]-ans[p[i]+4]);
    return 0;
}