【】 BZOJ1176 Bugs

Заголовок эффект: дается прямоугольная N * N, Q есть операции, каждая операция может быть прямоугольной или изменить опрашивать один весы точек и суб-прямоугольник.

Решение проблемы: CDQ раздел подходит для обработки независимо друг от друга модификаций между операционным и поддержкой тем , в автономном режиме.
Встречайте выше условие эксплуатации , очевидно , может сделать древовидные дерева крышки, но объем кода взлетать, я не мог писать QAQ. ,
Мы научились CDQ использовать разделяй и властвуй , чтобы решить эту проблему. Было установлено , что существует измерение времени между операциями, а именно: способность генерировать вклад в текущие операции запроса должно быть время , чтобы изменить операции перед запросом. поддержания веса суб-прямоугольник и двумерная префикс и может быть обработан, а именно: техническое обслуживание (1,1) и весовые коэффициенты могут быть между (п, т). Таким образом, превращается в проблему при заданном времени измерения помещения, стремясь удовлетворить 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;
} 

рекомендация

отwww.cnblogs.com/wzj-xhjbk/p/11109649.html