第十四届华中科技大学程序设计竞赛H:Magnificent Tree(cdq分治)

题目链接：点击打开链接

题目大意：

就是一个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");
}