题目链接:点击打开链接
题目大意:
就是一个w*w的矩阵,多次操作,每次将某个位置的值+val,询问某个矩阵当前所有值的和。
解题思路:
刚开始没看题目范围直接准备裸敲二维线段树的,结果队友提醒才发现w太大写不了,之后想了起来之前做的一道类似的题目,最后我是用线段树动态开点过的,然而这道题就算动态开点也是恐怖,然后就想到当时我查那道题题解的时候查到过cdq分治,然后凉了,当时没学cdq去用动态开点过了。。。赛后补题的时候发现这果然是cdq分治的一道论文题。在此附上论文链接:点击打开链接
论文作者的简写就是cdq,膜拜= = 。
再谈下cdq分治的算法,讲道理学的时候很痛苦,搞不懂这个算法到底是什么,后来照着代码手推了一遍,才发现cdq的奥妙所在,同时发现之前查的资料说的都挺清晰的,这里我也只能大概的说下cdq的算法,个人感觉是小范围的离线操作。cdq在于先把问题分解为小问题,就跟归并排序很像,之后合并问题。对于区间slove(l,r),我们会利用(l,mid)里的更新操作去计算它对(mid+1,r)里面查询操作的贡献,实现一个个人感觉是小范围离散的操作。不过最好推荐还是稍微手推一下,还是比较容易理解的。
Ac代码:
#include<bits/stdc++.h>
#define lson rt<<1
#define rson rt<<1|1
using namespace std;
typedef long long ll;
const int maxn=4e5+5;
const int INF=1e9+7;
const int mod=998244353;
int n,w,cnt,tot,bit[maxn],ans[maxn];
struct node //记录一些需要的值
{
int x,y,v,tp,id;
node(int xx=0,int yy=0,int vv=0,int t=0,int i=0){x=xx;y=yy;v=vv;tp=t;id=i;}
bool operator<(const node &p) const
{
return x<p.x;
}
}q[maxn];
int lowbit(int x) { return x&(-x);}
void add(int x,int val)
{
while(x<=w)
{
bit[x]+=val;
x+=lowbit(x);
}
}
int sum(int x)
{
int s=0;
while(x)
{
s+=bit[x];
x-=lowbit(x);
}
return s;
}
void cdq(int l,int r)
{
if(l==r) return ;
int mid=(l+r)>>1;
cdq(l,mid),cdq(mid+1,r);
sort(q+l,q+mid+1); //分别排序两个区间以进行离散操作
sort(q+mid+1,q+r+1); //分开排序保证左区间的更新在右区间的查询之前,
int j=l;
for(int i=mid+1;i<=r;i++) //利用离线的想法对右区间的查询更新贡献
{
while(j<=mid&&q[j].x<=q[i].x)
{
if(q[j].tp==1) add(q[j].y,q[j].v);
j++;
}
if(q[i].tp==2) ans[q[i].id]+=q[i].v*sum(q[i].y);
}
for(int i=l;i<j;i++) //将更新操作取消 便于下次合并
if(q[i].tp==1) add(q[i].y,-q[i].v);
}
int main()
{
scanf("%d%d",&n,&w);
{
cnt=0,tot=0; //cnt记录一共的操作次数 tot记录一共的查询操作
while(n--)
{
int flag,a,b,c,d;
scanf("%d",&flag);
if(flag==1)
{
scanf("%d%d%d",&a,&b,&c);
q[++cnt]=node(a,b,c,1,0); //标记当前为什么操作
}
if(flag==2)
{
tot++;
scanf("%d%d%d%d",&a,&b,&c,&d);
q[++cnt]=node(a-1,b-1,1,2,tot); //这里将查询分为四个小部分
q[++cnt]=node(c,d,1,2,tot);
q[++cnt]=node(a-1,d,-1,2,tot);
q[++cnt]=node(c,b-1,-1,2,tot);
}
}
cdq(1,cnt); //对所有操作次数进行cdq分治
for(int i=1;i<=tot;i++)
printf("%d\n",ans[i]);
}
//system("pause");
}