[BZOJ]1935:[Shoi2007]Tree 园丁的烦恼

题解: 离线排序对y坐标离散化后差分用树状数组统计贡献

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <stack>
#include <queue>
#include <cmath>
#include <set>
#include <map>
#define mp make_pair
#define pb push_back
#define pii pair<int,int>
#define link(x) for(edge *j=h[x];j;j=j->next)
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
const int MAXN=5e5+10;
const double eps=1e-8;
#define ll long long
using namespace std;
struct edge{int t,v;edge*next;}e[MAXN<<1],*h[MAXN],*o=e;
void add(int x,int y,int vul){o->t=y;o->v=vul;o->next=h[x];h[x]=o++;}
ll read(){
    ll 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 x*f;
}



int ans[MAXN];
typedef struct node{
    int x,y;
    friend bool operator<(node aa,node bb){return aa.x<bb.x;}
}node;
node d[MAXN];
vector<int>vec;
int base;

typedef struct Node{
    int x,y1,y2,id;
    friend bool operator<(Node aa,Node bb){return aa.x<bb.x;}
}Node;
Node que[MAXN<<1];

int sum[MAXN*3];
int get_id(int x){return x&(-x);}
void update(int x){
    for(int i=x;i<=base;i+=get_id(i))sum[i]++;
}
int Sum(int x){
    int ans1=0;
    for(int i=x;i>0;i-=get_id(i))ans1+=sum[i];
    return ans1;
}

int main(){
    int n=read();int m=read();
    inc(i,1,n)d[i].x=read(),d[i].y=read(),vec.pb(d[i].y);
    int x1,y1,x2,y2;
    int cnt=0;
    inc(i,1,m){
	x1=read();y1=read();x2=read();y2=read();
	que[++cnt]=(Node){x2,y1,y2,i};
	que[++cnt]=(Node){x1-1,y1,y2,-i};
	vec.pb(y1);vec.pb(y2);
    }
    sort(vec.begin(),vec.end());
    base=unique(vec.begin(),vec.end())-vec.begin();
    sort(d+1,d+n+1);
    sort(que+1,que+cnt+1);
    int tot=1;
    inc(i,1,cnt){
	while(tot<=n&&d[tot].x<=que[i].x){
	    int t1=lower_bound(vec.begin(),vec.begin()+base,d[tot].y)-vec.begin()+1;
	    update(t1);
	    tot++;
	}
	y1=lower_bound(vec.begin(),vec.begin()+base,que[i].y1)-vec.begin()+1;
	y2=lower_bound(vec.begin(),vec.begin()+base,que[i].y2)-vec.begin()+1;
	if(que[i].id<0)ans[-que[i].id]-=(Sum(y2)-Sum(y1-1));
	else ans[que[i].id]+=(Sum(y2)-Sum(y1-1));
    }
    inc(i,1,m)printf("%d
",ans[i]);
}

1935: [Shoi2007]Tree 园丁的烦恼

Time Limit: 15 Sec Memory Limit: 357 MB
Submit: 2077 Solved: 953
[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