BZOJ 1935: [Shoi2007]Tree 园丁的烦恼
大概 很闲....树状数组+离散化的题 嫌树状数组要操作四个方位 然后写了主席树.....反正sb 跑出来的效率也是....很惨 主席树大概是树状数组的2倍多....内存更是大得多 当做主席树的练习
1935: [Shoi2007]Tree 园丁的烦恼
Time Limit: 15 Sec Memory Limit: 357 MBSubmit: 1555 Solved: 701
[Submit][Status][Discuss]
Description
很久很久以前，在遥远的大陆上有一个美丽的国家。统治着这个美丽国家的国王是一个园艺爱好者，在他的皇家花园里种植着各种奇花异草。有一天国王漫步在花园里，若有所思，他问一个园丁道： “最近我在思索一个问题，如果我们把花坛摆成六个六角形，那么……” “那么本质上它是一个深度优先搜索，陛下”，园丁深深地向国王鞠了一躬。 “嗯……我听说有一种怪物叫九头蛇，它非常贪吃苹果树……” “是的，显然这是一道经典的动态规划题，早在N元4002年我们就已经发现了其中的奥秘了，陛下”。 “该死的，你究竟是什么来头？” “陛下息怒，干我们的这行经常莫名其妙地被问到和OI有关的题目，我也是为了预防万一啊！” 王者的尊严受到了伤害，这是不可容忍的。看来一般的难题是难不倒这位园丁的，国王最后打算用车轮战来消耗他的实力： “年轻人，在我的花园里的每一棵树可以用一个整数坐标来表示，一会儿，我的骑士们会来轮番询问你某一个矩阵内有多少树，如果你不能立即答对，你就准备走人吧！”说完，国王气呼呼地先走了。 这下轮到园丁傻眼了，他没有准备过这样的问题。所幸的是，作为“全国园丁保护联盟”的会长——你，可以成为他的最后一根救命稻草。
Input
第一行有两个整数n，m（0≤n≤500000，1≤m≤500000）。n代表皇家花园的树木的总数，m代表骑士们询问的次数。 文件接下来的n行，每行都有两个整数xi，yi，代表第i棵树的坐标（0≤xi，yi≤10000000）。 文件的最后m行，每行都有四个整数aj，bj，cj，dj，表示第j次询问，其中所问的矩形以（aj，bj）为左下坐标，以（cj，dj）为右上坐标。
Output
共输出m行，每行一个整数，即回答国王以（aj，bj）和（cj，dj）为界的矩形里有多少棵树。
Sample Input
3 1
0 0
0 1
1 0
0 0 1 1
0 0
0 1
1 0
0 0 1 1
Sample Output
3
#include <bits/stdc++.h> #define N 500005 using namespace std; int read(){ int x=0,f=1;char ch=getchar(); while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();} while(isdigit(ch))x=x*10+ch-'0',ch=getchar(); return f*x; } vector<int>x_; vector<int>y_; int rt[N]; typedef struct node{ int x;int y; friend bool operator<(node aa,node bb){ return aa.x<bb.x; } }node; node a[N]; typedef struct List{ int x1,y1,x2,y2;int biao; friend bool operator<(List aa,List bb){ return aa.x2<bb.x2; } }List; List b[N]; typedef struct nod{ int l,r,sum; }nod; nod d[50*N];int cnt; int ans[N]; void update(int &x,int y,int l,int r,int t){ cnt++;x=cnt;d[x]=d[y];d[x].sum++; if(l==r) return ; int mid=(l+r)>>1; if(t<=mid) update(d[x].l,d[y].l,l,mid,t); else update(d[x].r,d[y].r,mid+1,r,t); } int ans1; void querty(int x,int y,int l,int r,int l1,int r1){ if(l1<=l&&r<=r1){ ans1+=(d[y].sum-d[x].sum); return ; } int mid=(l+r)>>1; if(l1<=mid) querty(d[x].l,d[y].l,l,mid,l1,r1); if(r1>mid) querty(d[x].r,d[y].r,mid+1,r,l1,r1); } int main(){ ios::sync_with_stdio(false); int n,m;scanf("%d%d",&n,&m); cnt=0; for(int i=1;i<=n;i++){ scanf("%d%d",&a[i].x,&a[i].y); x_.push_back(a[i].x); y_.push_back(a[i].y); } for(int i=1;i<=m;i++){ scanf("%d%d%d%d",&b[i].x1,&b[i].y1,&b[i].x2,&b[i].y2); x_.push_back(b[i].x1);y_.push_back(b[i].y1); x_.push_back(b[i].x2);y_.push_back(b[i].y2); b[i].biao=i; } sort(x_.begin(),x_.end()); int t1=unique(x_.begin(),x_.end())-x_.begin(); sort(y_.begin(),y_.end()); int t2=unique(y_.begin(),y_.end())-y_.begin(); for(int i=1;i<=n;i++){ a[i].x=lower_bound(x_.begin(),x_.begin()+t1,a[i].x)-x_.begin()+1; a[i].y=lower_bound(y_.begin(),y_.begin()+t2,a[i].y)-y_.begin()+1; } sort(a+1,a+n+1); for(int i=1;i<=m;i++){ b[i].x1=lower_bound(x_.begin(),x_.begin()+t1,b[i].x1)-x_.begin()+1; b[i].y1=lower_bound(y_.begin(),y_.begin()+t2,b[i].y1)-y_.begin()+1; b[i].x2=lower_bound(x_.begin(),x_.begin()+t1,b[i].x2)-x_.begin()+1; b[i].y2=lower_bound(y_.begin(),y_.begin()+t2,b[i].y2)-y_.begin()+1; } sort(b+1,b+m+1);int u=1;int uu=0;int cnt1=0; for(int i=1;i<=m;i++){ for(int j=u;j<=n;j++){ if(a[j].x<=b[i].x2){ cnt1++; if(!rt[a[j].x]) { update(rt[a[j].x],rt[uu],1,t2,a[j].y);uu=a[j].x;} else{ update(rt[a[j].x],rt[a[j].x],1,t2,a[j].y); } } else{ u=j;break; } } if(cnt1==n) u=n+1; ans1=0; int l=1;int r=n;int ans2=0; int tt=b[i].x1-1; while(l<=r){ int mid=(l+r)>>1; if(a[mid].x<=tt){ ans2=mid;l=mid+1; } else r=mid-1; } l=1;r=n;int ans3=0;tt=b[i].x2; while(l<=r){ int mid=(l+r)>>1; if(a[mid].x<=tt){ ans3=mid;l=mid+1; } else r=mid-1; } querty(rt[a[ans2].x],rt[a[ans3].x],1,t2,b[i].y1,b[i].y2); ans[b[i].biao]=ans1; } for(int i=1;i<=m;i++) printf("%d\n",ans[i]); return 0; }
