题目描述

解析

差点活活恶心死
搬砖题
（其实细节没有那么多，还是代码能力太差）

利用矩阵的二维差分
加上树状数组搞一搞
就完事了（我实在不想再写了 ）

洛谷的双倍经验没有那么好水
MLE！
这题洛谷的空间限制非常苛刻
几乎没有多余的空间

以后要增强对空间复杂度的敏感性

代码

（看我的print就能看出我调了多久）

#include<bits/stdc++.h>
using namespace std;
const int N=2050;
#define ll long long
int d[N][N];
ll id[N][N],jd[N][N],ijd[N][N];
int n,m;void add_d(int x,int y,ll v){
//	printf("\nx=%d y=%d v=%lld\n",x,y,v);for(int i=x;i<=n;i+=i&-i){for(int j=y;j<=m;j+=j&-j){d[i][j]+=v;
//			printf("x=%d y=%d\n",i,j);}}
}
ll ask_d(int x,int y){
//	printf("\nask:\nx=%d y=%d\n",x,y);ll ans=0;for(int i=x;i;i-=i&-i){for(int j=y;j;j-=j&-j){ans+=d[i][j];} }return ans;
}
void add_id(int x,int y,ll v){v*=x;
//	printf("\nx=%d y=%d v=%lld\n",x,y,v);for(int i=x;i<=n;i+=i&-i){for(int j=y;j<=m;j+=j&-j) id[i][j]+=v;//printf("x=%d y=%d\n",i,j);}
}
ll ask_id(int x,int y){ll ans=0;
//	printf("ask:x=%d y=%d\n------------\n",x,y);for(int i=x;i;i-=i&-i){for(int j=y;j;j-=j&-j){ans+=id[i][j];
//			printf("x=%d y=%d\n",i,j);} }return ans;
}
void add_jd(int x,int y,ll v){v*=y;for(int i=x;i<=n;i+=i&-i){for(int j=y;j<=m;j+=j&-j) jd[i][j]+=v;}
}
ll ask_jd(int x,int y){ll ans=0;for(int i=x;i;i-=i&-i){for(int j=y;j;j-=j&-j){ans+=jd[i][j];} }return ans;
}
void add_ijd(int x,int y,ll v){v*=x*y;for(int i=x;i<=n;i+=i&-i){for(int j=y;j<=m;j+=j&-j) ijd[i][j]+=v;}
}
ll ask_ijd(int x,int y){ll ans=0;for(int i=x;i;i-=i&-i){for(int j=y;j;j-=j&-j){ans+=ijd[i][j];} }return ans;
}void add(int x1,int y1,int x2,int y2,ll v){add_d(x1,y1,v);add_d(x2+1,y2+1,v);add_d(x1,y2+1,-v);add_d(x2+1,y1,-v);add_id(x1,y1,v);add_id(x2+1,y2+1,v);add_id(x1,y2+1,-v);add_id(x2+1,y1,-v);add_jd(x1,y1,v);add_jd(x2+1,y2+1,v);add_jd(x1,y2+1,-v);add_jd(x2+1,y1,-v);add_ijd(x1,y1,v);add_ijd(x2+1,y2+1,v);add_ijd(x1,y2+1,-v);add_ijd(x2+1,y1,-v);
}
ll ask(int x,int y){
//	printf("ok%lld\n",d[2][1]);
//	n--;m--;ll res=(x*y+x+y+1)*ask_d(x,y)+ask_ijd(x,y)-(y+1)*ask_id(x,y)-(x+1)*ask_jd(x,y);
//	printf("x=%d y=%d d=%lld id=%lld jd=%lld ijd=%lld ans=%lld\n",x,y,
//	ask_d(x,y),ask_id(x,y),ask_jd(x,y),ask_ijd(x,y),res);
//	cout<<res<<endl;
//	n++;m++;return res;
}
int main(){
//	printf("%d",sizeof(d)/1024/1024);char flag;scanf("%c%d%d",&flag,&n,&m);
//	add_id(1,1,2);add_id(2,2,1);
//	printf("%lld\n",ask_id(3,3));int a,b,c,e,w;while(scanf(" %c",&flag)!=EOF){if(flag=='L'){scanf("%d%d%d%d%d",&a,&b,&c,&e,&w);add(a,b,c,e,w);
//			printf("okk:%lld\n",ask_d(1,1));}else{scanf("%d%d%d%d",&a,&b,&c,&e);
//			printf("okk:%lld\n",ask_d(1,1));printf("%lld\n",ask(c,e)+ask(a-1,b-1)-ask(a-1,e)-ask(c,b-1));}
//		printf("------------opp:%lld\n",ask_id(3,3));}	return 0;
}
/**/

thanks for reading！