タイトル効果:長方形N * N所与は、Qが操作を有し、各操作は、単一の点の重みとサブ矩形を変更またはクエリ矩形であってもよいです。
問題のソリューション:独立して、オフライン操作とサポートのトピックの間で、互いの変更の処理するために適しCDQパーティション。
明らかに、ツリーの木のカバーを行うことができます上記の条件を満たしているが、コードの量が、私はQAQを書くことができませんでした、離陸します。。
私たちは、この問題を解決するためにCDQ用分割統治を学びました。それはすなわち、動作間の時間ディメンションがあることを見出し:現在のクエリ操作への貢献を生成する機能は、お問い合わせの前に操作を変更するための時間でなければなりません。体重維持サブ矩形二次元プレフィックスと、すなわち、処理することができる:メンテナンス(1,1)及び重みは(M、N)の間とすることができます。三次元の半順序の問題:このように、すなわち点と、右に対応する点のX1 <X2及びY1 <Y2全て(X1、Y1)を満足するように求めて、所定の時間ディメンション前提で問題に形質転換しました。
コードは以下の通りです
#include <cstdio>
using namespace std;
const int maxn=2e6+10;
const int maxq=2e5+10;
int n,s,bit[maxn],ans[maxq];
struct node{int x,y,val,opt,aid;}q[maxq],tmp[maxq];// opt: 0 -> modify 1 -> query
int cnt,tot;
void modify(int pos,int val){
for(int i=pos;i<=n;i+=i&-i)bit[i]+=val;
}
int query(int pos){
int ret=0;
for(int i=pos;i;i-=i&-i)ret+=bit[i];
return ret;
}
void cdq(int l,int r){
if(l==r)return;
int mid=l+r>>1;
cdq(l,mid),cdq(mid+1,r);
int x=l,y=mid+1;
for(int i=l;i<=r;i++){
if(y>r||(x<=mid&&q[x].x<=q[y].x)){
tmp[i]=q[x];
if(q[x].opt==0)modify(q[x].y,q[x].val);
++x;
}else{
tmp[i]=q[y];
if(q[y].opt==1)ans[q[y].aid]+=query(q[y].y)*q[y].val;
++y;
}
}
for(int i=l;i<=mid;i++)if(q[i].opt==0)modify(q[i].y,-q[i].val);
for(int i=l;i<=r;i++)q[i]=tmp[i];
}
void add(int x,int y,int val,int opt){
++tot,q[tot].x=x,q[tot].y=y,q[tot].val=val,q[tot].opt=opt;
if(q[tot].opt==1)q[tot].aid=cnt;
}
void read_and_parse(){
scanf("%d%d",&s,&n);
while(1){
int opt;scanf("%d",&opt);
if(opt==3)break;
if(opt==1){
int x,y,a;
scanf("%d%d%d",&x,&y,&a);
add(x,y,a,0);
}
if(opt==2){
int x1,x2,y1,y2;
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
++cnt,add(x2,y2,1,1),add(x2,y1-1,-1,1),add(x1-1,y2,-1,1),add(x1-1,y1-1,1,1);
}
}
}
void solve(){
cdq(1,tot);
for(int i=1;i<=cnt;i++)printf("%d\n",ans[i]);
}
int main(){
read_and_parse();
solve();
return 0;
}