题意:
给定一个n*n的螺旋矩阵，每个方格的权值为当前数字各位之和，现在给出m个宝殿的坐标。再给出q个询问，询问给出的矩形中宝典的权值之和

解析:
这道题最难的地方就是给定坐标求出螺旋矩阵对应的权值

给出对应代码，细节大家自己看

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;
}

第二个难点就在于N的范围太大，我们不能用二维数组来求前缀和(会炸空间)，所以我们要给m个宝殿坐标和q个询问坐标进行离线处理。
对于求二维前缀和的问题，我们把s[x2][y2]-s[x1-1][y2]-s[x2][y1-1]+s[x1-1][y1-1]拆解成4个，用flag来标记哪些需要±
之后我们可以运用扫描线的方法，先按照x排序，然后对y进行映射，用树状数组维护映射后y上的值。用一条竖直的扫描线不断向右扫，遇到一个询问点就查询一次。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2e6+10;
ll c[N*2],cnt;
ll ans[N];
ll t,n,m,p;
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;
}


struct node
{
    
	ll x,y,flag,pos;
	bool operator<(const node &W){
    
		if(x!=W.x) return x<W.x;
		if(y!=W.y) return y<W.y;
		return pos<W.pos;
	}
}a[N*2];
ll lowbit(ll x)
{
    
	return x&(-x);
}
void add(int i,ll val,int n)
{
    
	while(i<=n)
	{
    
		c[i]+=val;
		i+=lowbit(i);
	}
}
ll sum(int i,int n)
{
    
	ll res=0;
	while(i)
	{
    
		res+=c[i];
		i-=lowbit(i);
	}
	return res;
}
int main()
{
    
	scanf("%lld",&t);
	while(t--)
	{
    
		scanf("%lld %lld %lld",&n,&m,&p);
		memset(c,0,sizeof c);
		memset(ans,0,sizeof ans);
		cnt=0;
		for(int i=1;i<=m;i++)
		{
    
			ll x,y;
			scanf("%lld %lld",&x,&y);
			ll tmp=getval(x,y,n);
			ll ans=0;
			while(tmp)
			{
    
				ans+=tmp%10;
				tmp/=10;
			}
			a[++cnt]={
    x,y,ans,i};
		}
		for(int i=m+1;i<=m+p;i++)
		{
    
			ll x1,x2,y1,y2;
			scanf("%lld %lld %lld %lld",&x1,&y1,&x2,&y2);
			a[++cnt]={
    x2,y2,1,i};
			a[++cnt]={
    x1-1,y2,-1,i};
			a[++cnt]={
    x2,y1-1,-1,i};
			a[++cnt]={
    x1-1,y1-1,1,i};
		}
		sort(a+1,a+1+cnt);
		for(int i=1;i<=cnt;i++)
		{
    
			if(a[i].pos<=m) add(a[i].y,a[i].flag,n);
			else
			{
    
				ans[a[i].pos-m]+=sum(a[i].y,n)*a[i].flag;
			}
		}
		for(int i=1;i<=p;i++) printf("%lld\n",ans[i]);
	}
}