Fast Matrix Operations

题目：https://cn.vjudge.net/problem/UVA-11992

给定一个r*c(r<=20,r*c<=1e6)的矩阵，其元素都是0，现在对其子矩阵进行操作。

1 x1 y1 x2 y2 val 表示将(x1,y1,x2,y2)（x1<=x2,y1<=y2）子矩阵中的所有元素加val；

2 x1 y1 x2 y2 val 表示将(x1,y1,x2,y2)（x1<=x2,y1<=y2）子矩阵中的所有元素变为val；

3 x1 y1 x2 y2 val 表示输出(x1,y1,x2,y2)（x1<=x2,y1<=y2）子矩阵中的所有元素的和，最小值和最大值。

共有m次查询(1<=m<=20000)

 

好吧 二维线段树维护

为了优化这里的置成v，可以看成是*0+v，加上v可以看成*1+v；

打20棵线段树维护就可以解决

其实我不会打代码

只好。。。

//一维做法见http://blog.csdn.net/a601025382s/article/details/3835653（好像很麻烦）

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
long long read()
{long long x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;}
struct twtree{
int l,r;
int add,mul;
int sum,minn,maxn;
void init(){add=0;mul=1;}
inline int mid(){return (l+r)>>1;}
inline int len(){return r-l+1;}
void tch(int m){sum*=m;minn*=m;maxn*=m;mul*=m;add*=m;}
void tadd(int a){minn+=a;maxn+=a;sum+=a*len();add+=a;}
}T[21][500005];//e..
void pushdown(int num,int rt){
T[num][rt<<1].tch(T[num][rt].mul);
T[num][rt<<1|1].tch(T[num][rt].mul);
T[num][rt<<1].tadd(T[num][rt].add);
T[num][rt<<1|1].tadd(T[num][rt].add);
T[num][rt].init();
}
void pushup(int num,int rt){
T[num][rt].sum=T[num][rt<<1].sum+T[num][rt<<1|1].sum;
T[num][rt].maxn=max(T[num][rt<<1].maxn,T[num][rt<<1|1].maxn);
T[num][rt].minn=min(T[num][rt<<1].minn,T[num][rt<<1|1].minn);
}
void build(int num,int rt,int l,int r){
T[num][rt].l=l;
T[num][rt].r=r;
T[num][rt].minn=T[num][rt].maxn=T[num][rt].sum=0;
T[num][rt].init();
if(l==r){
return;
}
int mid=T[num][rt].mid();
build(num,rt<<1,l,mid);
build(num,rt<<1|1,mid+1,r);
pushup(num,rt);
}
void cheng(int num,int rt,int l,int r,int mul,int add){
if(T[num][rt].l==l&&T[num][rt].r==r){
T[num][rt].tch(mul);
T[num][rt].tadd(add);
return;
}
pushdown(num,rt);
int mid=T[num][rt].mid();
if(mid>=r)cheng(num,rt<<1,l,r,mul,add);
else if(mid<l)cheng(num,rt<<1|1,l,r,mul,add);
else {
cheng(num,rt<<1,l,mid,mul,add);
cheng(num,rt<<1|1,mid+1,r,mul,add);
}
pushup(num,rt);
}
twtree query(int num,int rt,int l,int r){
if(T[num][rt].l==l&&T[num][rt].r==r)
return T[num][rt];
pushdown(num,rt);
int mid=T[num][rt].mid();
if(mid>=r)return query(num,rt<<1,l,r);
else if(mid<l)return query(num,rt<<1|1,l,r);
else{
twtree t1=query(num,rt<<1,l,mid);
twtree t2=query(num,rt<<1|1,mid+1,r);
twtree t;
t.maxn=max(t1.maxn,t2.maxn);
t.minn=min(t1.minn,t2.minn);
t.sum=t1.sum+t2.sum;
return t;
}
}
int main(){
int n,m,r;
while(~scanf("%d",&r)){
n=read(),m=read();
for(int i=1;i<=r;i++)
build(i,1,1,n);
for(int i=1;i<=m;i++){
int t=read(),x1=read(),y1=read(),x2=read(),y2=read();;
if(t==1){
int v=read();
for(int i=x1;i<=x2;i++)
cheng(i,1,y1,y2,1,v);
}
else if(t==2){
int v=read();
for(int i=x1;i<=x2;i++)
cheng(i,1,y1,y2,0,v);
}
else {
int sum=0,minn=999999999,maxn=0;
twtree ans;
for(int i=x1;i<=x2;i++){
ans=query(i,1,y1,y2);
sum+=ans.sum;
minn=min(minn,ans.minn);
maxn=max(maxn,ans.maxn);
}
cout<<sum<<" "<<minn<<" "<<maxn<<endl;
}
}
}
return 0;
}
//二维线段树模板

转载于:https://www.cnblogs.com/0degreeofsail/p/7146070.html