bzoj 1395
很久很久以前，在遥远的大陆上有一个美丽的国家。统治着这个美丽国家的国王是一个园艺爱好者，在他的皇家花园里种植着各种奇花异草。有一天国王漫步在花园里，若有所思，他问一个园丁道： “最近我在思索一个问题，如果我们把花坛摆成六个六角形，那么……” “那么本质上它是一个深度优先搜索，陛下”，园丁深深地向国王鞠了一躬。 “嗯……我听说有一种怪物叫九头蛇，它非常贪吃苹果树……” “是的，显然这是一道经典的动态规划题，早在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 Input3 1
0 0
0 1
1 0
0 0 1 1
0 0
0 1
1 0
0 0 1 1
Sample Output
3
注意，这道题，一看就是可以用二维树状数组来做的，但是这样必定mle，范围太大了，但是由于查询只有5e6次，所以开1维树状数组，然后坐标离散化，是可以做到！
坐标离散化，也就是根据y的坐标顺序排列，把他们分别看作1，2，3，4，5；方法可以看挑战，上面有。
但是，注意我们的离散化数组要开原来数据范围的5倍，为什么，因为输入的点要一次，然后查询时候，因为我们的树状数组是一个“假”的二维数组，所以这时候必定要算分块，一共四个
树状数组这里借鉴了hzwer大佬的模板，bzoj'上的很多题看起来都很悬，有的时候居然卡常数！！！！！今天两个简易树状数组把我坑死了。
1 #include <iostream> 2 #include <cstdio> 3 #include <algorithm> 4 using namespace std; 5 int n,m; 6 int a[500005],b[500005],c[500003],d[500003],y[500003],x[500003]; 7 int total,howmanyq; 8 inline int lowbit(int x){return x&(-x);} 9 struct data{ 10 int x,y,id,foo; 11 } q[2500005]; 12 int disc[1500004];int t[2500005]; 13 int ans[500005][5]; 14 inline bool operator<(data a,data b) 15 { 16 return a.x<b.x||(a.x==b.x&&a.foo<b.foo); 17 } 18 inline void add(int x,int y) 19 { 20 for(int i=x;i<=total;i+=lowbit(i)) 21 t[i]+=y; 22 } 23 inline int query(int x) 24 { 25 int sum=0; 26 for(int i=x;i;i-=lowbit(i)) 27 sum+=t[i]; 28 return sum; 29 } 30 int find(int x) 31 { 32 return lower_bound(disc,disc+total,x)-disc+1; 33 } 34 35 int main() 36 { 37 scanf("%d%d",&n,&m); 38 for(int i=1;i<=n;i++) 39 { 40 scanf("%d%d",&x[i],&y[i]); 41 disc[total]=y[i]; 42 total++; 43 } 44 for(int i=1;i<=m;i++) 45 { 46 scanf("%d%d%d%d",&a[i],&b[i],&c[i],&d[i]); 47 disc[total]=b[i]; 48 total++; 49 disc[total]=d[i]; 50 total++; 51 } 52 sort(disc,disc+total); 53 for(int i=1;i<=n;i++) 54 { 55 y[i]=find(y[i]); 56 howmanyq++; 57 q[howmanyq].x=x[i]; 58 q[howmanyq].y=y[i]; 59 } 60 for(int i=1;i<=m;i++) 61 { 62 b[i]=find(b[i]);d[i]=find(d[i]); 63 q[++howmanyq].x=c[i];q[howmanyq].y=d[i];q[howmanyq].id=i;q[howmanyq].foo=1; 64 q[++howmanyq].x=a[i]-1;q[howmanyq].y=d[i];q[howmanyq].id=i;q[howmanyq].foo=2; 65 q[++howmanyq].x=c[i];q[howmanyq].y=b[i]-1;q[howmanyq].id=i;q[howmanyq].foo=3; 66 q[++howmanyq].x=a[i]-1;q[howmanyq].y=b[i]-1;q[howmanyq].id=i;q[howmanyq].foo=4; 67 } 68 sort(q+1,q+howmanyq+1); 69 for(int i=1;i<=howmanyq;i++) 70 { 71 if(!q[i].foo)add(q[i].y,1); 72 else 73 { 74 int t=query(q[i].y); 75 ans[q[i].id][q[i].foo]=t; 76 } 77 } 78 for(int i=1;i<=m;i++) 79 { 80 int t=ans[i][1]+ans[i][4]-ans[i][2]-ans[i][3]; 81 printf("%d\n",t); 82 } 83 return 0; 84 }
