大致题意
给一个矩阵的生成方式,这个矩阵为 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;
}