bzoj2683&&1176[Balkan2007]Mokia cdq分治

应该算是入门题,只不过我很久没做这种矩形题了,所以没有想到一个关键点导致我模型没转化过来。
还是像之前一样把询问和赋值多添加一个维,表示时间,然后添加很显然就不说了,每次直接bit加入就好,问题是查询。
一个子矩阵x1 y1 x2 y2 的答案等于ans[x2][y2]+ans[x1][y1]-ans[x1][y2]-ans[x2][y1]
注意一下边界就好了,然后直接cdq。
%GEOTCBRL
还有那个S是没有意义的,题目错了。

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=2e6+5;
typedef int ll;
int n;
struct node
{
    int id,x,y,v,t;
    node(){}
    node(int id,int x,int y,int v,int t):id(id),x(x),y(y),v(v),t(t){}
    bool operator < (const node a)const
    {
        if (x==a.x&&y==a.y)return t<a.t;
        if (x==a.x)return y<a.y;
        return x<a.x;
    }
}q[N],tmp[N];
inline int lowbit(int x)
{
    return x&(-x);
}
ll val[N];
inline void add(int x,int y)
{
    while (x<=n)
    {
        val[x]+=y;
        x+=lowbit(x);
    }
}
inline ll ask(int x)
{
    ll ret=0;
    while (x>0)
    {
        ret+=val[x];
        x-=lowbit(x);
    }
    return ret;
}
int tot,ques;
ll ans[N];
inline void ins(int x,int y,int v,int t)
{
    ++tot;
    q[tot]=node(tot,x,y,v,t);
}
inline void cdq(int l,int r)
{   
    if (l==r)return;
    int mid=(l+r)>>1;
    fo(i,l,r)
    {
        if (q[i].id<=mid&&!q[i].t)add(q[i].y,q[i].v);
        if (q[i].id>mid&&q[i].t)ans[q[i].t]+=ask(q[i].y)*q[i].v;
    }
    fo(i,l,r)if (q[i].id<=mid&&!q[i].t)add(q[i].y,-q[i].v);
    int l1=l,l2=mid+1;
    fo(i,l,r)
    if (q[i].id<=mid)tmp[l1++]=q[i];
    else tmp[l2++]=q[i];
    fo(i,l,r)q[i]=tmp[i];
    cdq(l,mid);
    cdq(mid+1,r);
}
int main()
{
    scanf("%d",&n);
    int opt;
    while (scanf("%d",&opt)&&opt!=3)
    {
        int x1,y1,x2,y2;
        scanf("%d%d%d",&x1,&y1,&x2);
        if (opt==1)ins(x1,y1,x2,0);
        else
        {
            ++ques;
            scanf("%d",&y2);
            ins(x1-1,y1-1,1,ques);
            ins(x1-1,y2,-1,ques);
            ins(x2,y1-1,-1,ques);
            ins(x2,y2,1,ques);
        }
    }
    sort(q+1,q+tot+1);
    cdq(1,tot);
    fo(i,1,ques)printf("%d\n",ans[i]);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值