2019ICPC南京网络赛A The beautiful values of the palace （排序+树状数组 or cdq分治 不带修改矩形内和）

大致题意

给一个矩阵的生成方式，这个矩阵为 nxn (n是奇数，1<=n<=1000000)，现给出m个位置，填上对应的矩阵位置的值（这个值是矩阵中原数的各位数字之和），p次询问，询问矩形 x1,y1,x2,y2，的范围内的点的权值和。

思路

首先总结出矩阵任意坐标的值的常数级别求法，剩下的就是一个求解平面上有些点，点有些权值，然后询问某个矩形内的权值和的问题，而且不带修改，是静态的。
由于对算法理解还不够深入，只想到直接当作三维偏序cdq分治。首先是操作看作一维，按时间顺序，前面的操作影响后面，另外 x 坐标看作一维，把查询拆成4个问题，每次查询的其实是 （0，0) 和 （x,y) 所围的矩形的价值和，所以 x 看作第二维，小的x会对后面的范围产生贡献。然后第3维 y 用树状数组维护价值的前缀和。

但是实际上，这道题目不需要修改，所有的赋值操作都在询问操作前面，所以第二维的分治过程其实就是在按照x排序，因为每个询问只和比自己小的赋值操作有关，所以分治的最终结果就是将所有插入操作的x 比 询问操作的x 小的贡献都给询问操作算上，所以直接按x排序后，从前往后用树状数组维护，遇到插入操作，就插入值，遇到询问就求和。

下面贴一下cdq分治的代码。树状数组的比较好写

代码

#include<bits/stdc++.h>
using namespace std;
#define maxn 2000005
#define maxm 200005
#define ll long long int
#define INF 0x3f3f3f3f
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
#define mem(a) memset(a,0,sizeof(a))
#define sqr(x) (x*x)
#define inf (ll)2e18+1
#define mod 998244353
#define ls x<<1
#define rs x<<1|1
int read(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
     while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
     return f*x;
}
int T,n,m,p,cnt,cntq;
ll ans[maxn];
struct node{int x,y,z,type;
    bool operator < (const node &rid)const{
        if(x!=rid.x)return x < rid.x;
        else if(y!=rid.y)return y < rid.y;
        else return type < rid.type;
    }
}q[maxn],b[maxn];
ll c[maxn];
inline void add(int x,int y){
    while(x<maxn){
        c[x]+=y;
        x+=x&(-x);
    }
}
inline ll ask(int x){
    ll res=0;
    while(x){
        res+=c[x];
        x-=x&(-x);
    }
    return res;
}
inline ll getval(ll x, ll y, ll n) {
 ll t = min(min(x, y), min(n - x + 1, n - y + 1));
 ll ta = 4 * (t - 1) * (n - t + 1);
 if (x == n - t + 1) ta += n - t - y + 2;//在所在圈层的右侧边
 else if (x == t) ta += 2 * n - 5 * t + y + 3;//左侧边
 else if (y == t) ta += 2 * n - 3 * t - x + 3;//下侧边
 else ta += 3 * n - 7 * t + x + 4;//上侧边
 return ta;
}
inline void cdq(int l,int r){
    if(l==r)return ;
    int mid=(l+r)>>1;
    cdq(l,mid);cdq(mid+1,r);
    int p1=l,p2=mid+1,t=l;
    while(p1<=mid&&p2<=r){
        if(q[p1]<q[p2]){
            if(q[p1].type==0)add(q[p1].y,q[p1].z);
            b[t++]=q[p1++];
        }
        else {
            if(q[p2].type==1)ans[q[p2].z]-=ask(q[p2].y);
            else if(q[p2].type==2)ans[q[p2].z]+=ask(q[p2].y);
            b[t++]=q[p2++];
        }
    }
    while(p1<=mid){
        if(q[p1].type==0)add(q[p1].y,q[p1].z);
        b[t++]=q[p1++];
    }
    while(p2<=r){
        if(q[p2].type==1)ans[q[p2].z]-=ask(q[p2].y);
        else if(q[p2].type==2)ans[q[p2].z]+=ask(q[p2].y);
        b[t++]=q[p2++];
    }
    inc(i,l,mid)if(q[i].type==0)add(q[i].y,-q[i].z);
    inc(i,l,r)q[i]=b[i];
}
int main()
{
    T=read();
    while(T--){
        cnt=0;
        //mem(c);
        n=read();m=read();p=read();
        int x,y,z,r;
        inc(i,1,m){
            x=read();y=read();
            ll tmp=getval(x,y,n);
            z=0;
            while(tmp){
                z+=tmp%10;
                tmp/=10;
            }
            q[++cnt]={x,y,z,0};
        }
        inc(i,1,p){
            x=read();y=read();z=read();r=read();
            q[++cnt]={x-1,y-1,i,2};
            q[++cnt]={x-1,r,i,1};
            q[++cnt]={z,y-1,i,1};
            q[++cnt]={z,r,i,2};
        }
        cdq(1,cnt);
        inc(i,1,p)printf("%lld\n",ans[i]);
        inc(i,1,p)ans[i]=0;
    }
    return 0;
}