传送门
纪念我第一个CDQ分治代码(当然调了很久)还有这个很像某手机品牌的题目名称
二位数点但是差值太大树套树MLE+TLE
但是离线很重要,这个时候CDQ分治派上用场了(%%%陈丹琦dalao发明了这个东西)
我们先把所有操作按照x坐标先y坐标后排序(双关键字保平安)
然后区间l,r取mid,计算l~mid的修改对mid+1~r的查询的影响
我的做法是先把编号小于等于mid的换到左边去,然后直接找右边的查询计算左边贡献,最后向下递归
求贡献的话我们用一个树状数组维护一下合法的值(最后别忘记清空,几乎很多CDQ的题都要套上一两个数据结构吧。。。)
这是CDQ的主要思想,这题求区间矩阵和我们可以用矩阵前缀和的方式来实现
也就是说询问操作一分四然后分别操作即可
其实这个向下递归有点奇怪,我把它放前面就会出问题(照理来说应该不会有问题把。。。)希望大神解答>_<
(其实本蒟蒻对于CDQ的理解不深慢慢长进)
贴一下我拙劣的代码
#include<bits/stdc++.h>
using namespace std;
struct ppap{
int b,j,x,y,v,k;
}c[1000001],d[1000001];
int f[2000001],ans[1000001],S,n,cnt,q;
inline bool cmp(ppap a,ppap b){return a.x==b.x?a.y<b.y:a.x<b.x;}
inline int lowbit(int x){return x&-x;}
inline void add(int x,int v){for(int i=x;i<=n;i+=lowbit(i))f[i]+=v;}
inline int sum(int x){int ans=0;for(int i=x;i;i-=lowbit(i))ans+=f[i];return ans;}
void cdq(int l,int r){
if(l==r)return;
int mid=l+r>>1,dd1=l-1,dd2=mid;
for(int i=l;i<=r;i++)if(c[i].b<=mid)d[++dd1]=c[i];
else d[++dd2]=c[i];
for(int i=l;i<=r;i++)c[i]=d[i];
int cc=l;
for(int i=mid+1;i<=r;i++)if(c[i].j==2){
for(;cc<=mid&&c[cc].x<=c[i].x;cc++)if(c[cc].j==1)add(c[cc].y,c[cc].v);
ans[c[i].k]+=sum(c[i].y)*c[i].v;
}
for(int i=l;i<cc;i++)if(c[i].j==1)add(c[i].y,-c[i].v);
cdq(l,mid);cdq(mid+1,r);
}
int main()
{
scanf("%d%d",&S,&n);int k;
while(1){
scanf("%d",&k);if(k==3)break;
if(k==1){
int x,y,z;scanf("%d%d%d",&x,&y,&z);
c[++cnt].b=cnt;c[cnt].j=1;c[cnt].x=x;c[cnt].y=y;c[cnt].v=z;c[cnt].k=0;
}else{
q++;int x,y,xx,yy;scanf("%d%d%d%d",&x,&y,&xx,&yy);
// ans[q]=S*(xx-x)*(yy-y);
//还有这题的S并没有用。。。上面的预处理也就不用啦
c[++cnt].b=cnt;c[cnt].j=2;c[cnt].x=xx;c[cnt].y=yy;c[cnt].v=1;c[cnt].k=q;
c[++cnt].b=cnt;c[cnt].j=2;c[cnt].x=x-1;c[cnt].y=y-1;c[cnt].v=1;c[cnt].k=q;
c[++cnt].b=cnt;c[cnt].j=2;c[cnt].x=x-1;c[cnt].y=yy;c[cnt].v=-1;c[cnt].k=q;
c[++cnt].b=cnt;c[cnt].j=2;c[cnt].x=xx;c[cnt].y=y-1;c[cnt].v=-1;c[cnt].k=q;
}
}
sort(c+1,c+cnt+1,cmp);
cdq(1,cnt);
for(int i=1;i<=q;i++)printf("%d\n",ans[i]);
return 0;
}
呵呵,我是不会告诉你这题是双倍经验的。。。QAQ
双倍经验传送门(BZOJ2683非权限号的就别想啦。。。)