bzoj2683（要改一点代码）&&bzoj1176: [Balkan2007]Mokia

仍然是一道cdq模版。。

那么对于一个询问，容斥一下分成四个，变成问(1,1)~(x,y)，那么对于x,y，修改只有x'<=x&&y'<=y，才对询问有影响，那么加上读入顺序，就是一个三维偏序了。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;

struct node
{
    int tp,x,y,v,be;//第一维时间已经有序
    //tp=0 表示修改 v为修改的值
    //tp=-1、1询问 v为所求值
}a[];int len;
void ins(int tp,int x,int y,int v,int be)
{
    len++;
    a[len].tp=tp;
    a[len].x=x;a[len].y=y;
    a[len].v=v;a[len].be=be;
}

//----------------------------------

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

//---------------树状数组----------------- 

node t[];
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].x<=a[j].x)
        {
            if(a[i].tp==)add(a[i].y,a[i].v);
            t[p++]=a[i++];
        }
        else
        {
            if(a[j].tp!=)a[j].v+=getsum(a[j].y);
            t[p++]=a[j++];
        }
    }
    while(i<=mid)
    {
        if(a[i].tp==)add(a[i].y,a[i].v);
        t[p++]=a[i++];
    }
    while(j<=r)
    {
        if(a[j].tp!=)a[j].v+=getsum(a[j].y);
        t[p++]=a[j++];
    }

    for(int i=l;i<=mid;i++)
        if(a[i].tp==)add(a[i].y,-a[i].v);

    for(int i=l;i<=r;i++)a[i]=t[i];
}

//-------------------cdq------------------------ 

int ans[],ansl;
int main()
{
    int S;
    scanf("%d%d",&S,&n);

    int op,x1,y1,x2,y2,v;len=;
    while(scanf("%d",&op)!=EOF)
    {
        if(op==)break;

        if(op==)
        {
            scanf("%d%d%d",&x1,&y1,&v);
            ins(,x1,y1,v,);
        }
        else
        {
            ansl++;
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);

            if(x1->&&y1->)ins(,x1-,y1-,,ansl);
            if(x1->)ins(-,x1-,y2,,ansl);
            if(y1->)ins(-,x2,y1-,,ansl);
            ins(,x2,y2,,ansl);
        }
    }

    cdq(,len);

    memset(ans,,sizeof(ans));
    for(int i=;i<=len;i++)
        if(a[i].tp!=)ans[a[i].be]+=a[i].v*a[i].tp;

    for(int i=;i<=ansl;i++)printf("%d\n",ans[i]);
    return ;
}