BZOJ 1176: [Balkan2007]Mokia&&2683: 简单题|cdq分治

最新推荐文章于 2019-07-15 14:13:24 发布

ws_yzy

最新推荐文章于 2019-07-15 14:13:24 发布

阅读量699

点赞数

分类专栏：分治文章标签： cdq分治

本文链接：https://blog.csdn.net/ws_yzy/article/details/50783200

版权

分治专栏收录该内容

2 篇文章 0 订阅

订阅专栏

标准cdq分治，网上的题解已经烂大街了，不多赘述了233333
$ps：$ 有什么问题可以留言

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
#include<set>
#include<map>
#define lowbit(x) (x&(-x))
#define N 2000500
using namespace std;
struct W{int x,y,v,id,op;}q[N/10],q1[N/10];
int pos[N],tr[N],ans[N];
int tot,n,cnt,S;
int sc()
{
    int i=0,f=1; char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i*f;
}
void change(int x,int f)
{
    for(;x<=n;x+=lowbit(x))tr[x]+=f;
}
int ask(int x)
{
    int ans=0;
    for(;x;x-=lowbit(x))ans+=tr[x];
    return ans;
}
void cdq_solve(int l,int r)
{
    if(l==r)return;
    int mid=l+r>>1;
    cdq_solve(l,mid),cdq_solve(mid+1,r);
    int i=l,j=mid+1,last=0;
    while(j<=r)
    {
        while(q[i].op==2&&i<=mid)i++;
        while(q[j].op==1&&j<=r)j++;
        if(i<=mid&&q[i].x<=q[j].x) change(q[i].y,q[i].v),last=i++;
        else if(j<=r) ans[q[j].id]+=ask(q[j].y),j++;
    }
    for(int i=l;i<=last;i++)if(q[i].op==1)change(q[i].y,-q[i].v);
    if(l==1&&r==cnt)return;
    int l1=l,l2=mid+1;
    for(int i=l;i<=r;i++)
        if(l2>r||(l1<=mid&&q[l1].x<=q[l2].x))q1[i]=q[l1++];else q1[i]=q[l2++];
    for(int i=l;i<=r;i++)q[i]=q1[i];
}
void pr(int x)
{
    printf("%d %d\n",x,ans[x]);
}
int main()
{
    S=sc(),n=sc();
    for(;;)
    {
        int x=sc(); 
        if(x==1)
            q[++cnt].op=1,q[cnt].x=sc(),q[cnt].y=sc(),q[cnt].v=sc();
        else if(x==2)
        {
            int x1=sc(),y1=sc(),x2=sc(),y2=sc();pos[++tot]=cnt;ans[cnt+1]=(x2-x1+1)*(y2-y1+1)*S;
            q[++cnt].op=2,q[cnt].id=cnt,q[cnt].x=x2,q[cnt].y=y2;
            q[++cnt].op=2,q[cnt].id=cnt,q[cnt].x=x1-1,q[cnt].y=y1-1;
            q[++cnt].op=2,q[cnt].id=cnt,q[cnt].x=x1-1,q[cnt].y=y2;
            q[++cnt].op=2,q[cnt].id=cnt,q[cnt].x=x2,q[cnt].y=y1-1;
        }
        else break;
    }
    cdq_solve(1,cnt);
    for(int i=1,x=pos[i];i<=tot;x=pos[++i])
        printf("%d\n",ans[x+1]+ans[x+2]-ans[x+3]-ans[x+4]);
    return 0;
}