1. 程式人生 > >lp4390 [BOI2007]Mokia 摩基亞

lp4390 [BOI2007]Mokia 摩基亞

題目概述:
  帶修改的二維數點。

題目思路:
  把一個詢問拆成4個,那麼我們很容易想到cdq分治,直接做即可。
  用一些小技巧可以優化程式碼。老久沒寫cdq了,果然掛了,i=j=1是怎麼回事。

題目程式碼:

#include<cstdio>
#include<algorithm>

const int MaQ1=10000,MaQ2=160000;
#define opt1 (true)
#define opt2 (false)

struct ques{
    int t,x,y,num;
    bool opt;
    friend bool operator <(ques a,ques b){
        
if(a.t!=b.t) return a.t<b.t; return a.x+a.y>b.x+b.y; } }q[MaQ1+MaQ2*4+5],ls[MaQ1+MaQ2*4+5]; int n,m,k; inline int lowbit(int a){return a&-a;} int c[2000005]; void add(int val,int r){ for(;val<=n;val+=lowbit(val)) c[val]+=r; return ; } int que(int val){ int
r; for(r=0;val;val-=lowbit(val)) r+=c[val]; return r; } int st[1000005],con; void clear(int val){ for(;val<=n;val+=lowbit(val)) c[val]=0; return ; } void prout(){ /*printf("\n"); for(int i=1;i<=k;i++){ printf("%d:%d %d %d\n",i,q[i].t,q[i].x,q[i].y); } printf("\n");
*/ return ; } void cdq(int l,int r){ if(l==r) return ; #define mid (l+r>>1) cdq(l,mid);cdq(mid+1,r); //BUG A //int i=1,j=1,k=mid+1; int i=l,j=l,k=mid+1; while(j<=mid&&k<=r){ if(q[j].x<=q[k].x){ if(q[j].opt){ add(q[j].y,q[j].num); st[con++]=q[j].y; } ls[i++]=q[j++]; } else{ if(!q[k].opt) q[k].num+=que(q[k].y); ls[i++]=q[k++]; } } while(j<=mid){ if(q[j].opt){ add(q[j].y,q[j].num); st[con++]=q[j].y; } ls[i++]=q[j++]; } while(k<=r){ if(!q[k].opt) q[k].num+=que(q[k].y); ls[i++]=q[k++]; } #undef mid for(int i=0;i<con;i++) clear(st[i]); con=0; for(int i=l;i<=r;i++) q[i]=ls[i]; prout(); return ; } int main(){ scanf("%*d%d",&n); int opt; while(true){ scanf("%d",&opt); if(opt==3) break; ++m; if(opt==1){ int x,y,z; scanf("%d%d%d",&x,&y,&z); q[++k]= (ques) {m,x,y,z,1}; } else{ int x1,y1,x2,y2; scanf("%d%d%d%d",&x1,&y1,&x2,&y2); q[++k]=(ques) { m,x1-1,y1-1,0,opt2 }; q[++k]=(ques) { m,x1-1,y2,0,opt2 }; q[++k]=(ques) { m,x2,y1-1,0,opt2 }; q[++k]=(ques) { m,x2,y2,0,opt2 }; } } prout(); cdq(1,k); std::sort(q+1,q+k+1); prout(); int now=1; for(int i=1;i<=m;i++){ if(q[now].opt==opt1){++now;continue;} if(q[now].opt==opt2) printf("%d\n",q[now].num+q[now+3].num-q[now+1].num-q[now+2].num),now+=4; } return 0; }