BSOJ：3799 训练场--贪心+倍增

3799 -- 【模拟试题】训练场
Description
　　某中学有很多学生社团，其中电竞社是最受欢迎的一个。该社团中总共有N只游戏战队，但是该中学只有一个游戏训练场馆，每次只能容纳一只战队训练。 
　　每只战队对训练时间都有一定的要求，比如甲战队想要在a到b这段时间训练,乙战队想要在c到d这段时间训练，...... 
　　作为训练场管理员的你总是收到形如(x,y)的询问，意思是查询在x到y这段时间内，最多能满足多少个只战队训练。现在有M个询问摆在你面前，请你快速做出回答！
Input
　　第一行，两个整数N和M。 
　　接下来N行，每行两个整数a和b，表示一只战队训练的起止时间点。 
　　接下来M行，每行两个整数x和y，表示一个询问的起止点。
Output
　　M行，每行一个整数，表示一次询问的答案。
Sample Input
3 2
1 2
2 3
1 3
1 2
1 3
Sample Output
1
2
Hint
【数据范围】 
x < y同时a < b。 0 <= x,y,a,b <= 1,000,000,000 
对于30%的数据,有0 < N,M <= 2000 
对于50%的数据,有0 < N,M <= 50000 

对于100%的数据，有0 < N,M <= 100000

作为当天考试压轴题这道题还是充满了震撼..毕竟当时too young,看到这道题就先去打了个线段树，然后发现根本没法维护..于是30分弃坑了。

前几天看了看覆盖牛棚（内网链接）的题解，发现这道题简直神似，都使用了同样的算法：

倍增。

倍增主要拿来覆盖区间的速度。

现在要二分一个解，对于解的验证，其速度相对较慢，举一个例子：

有线段a(st,ed)，在这条线段后面，找出另一条线段b，使得a与b有交集，且b的右端点靠的最后。

朴素的算法是一条一条线段地枚举，速度为O(n).

O(n)太慢了，考虑优化，将线段排序，获得单调之后考虑二分，每次二分一个mid，指在a后mid个，且与a重合的线段。我们发现这个mid还是要O(mid)走过去，还是较慢。

着重优化O(mid)，记录f[i][j]为i根栏杆后面连接的第2^j根栏杆，再按照朴素的ST推走即可。对于一个mid，只需O(log(mid))就能找到这个元素。

- -感觉讲的很烂，毕竟语文捉鸡。还是回到这个题里面。

对于这个问题，我们有一个贪心策略：被包含的区间肯定比包含其的大区间更优，直接舍弃大区间。证明略。

这样问题就转化成了一个统计问题：在给定区间(x,y)中，统计有多少个互不相交的区间（当然要保证最优）。

数据范围使得O(n)算法不能通过，现在就要用上倍增了。

设f[i][j]为i往右经过2^j个区间后第2^j区间（右端点坐标最小（贪心保证结果最优））的区间编号。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cmath>
using namespace std;
inline int read()
{
	int bj=1;
	char ch=getchar();
	while(ch<'0'||ch>'9')
	{
		if(ch=='-')bj=-1;
		ch=getchar();
	}
	int ret=0;
	while(ch>='0'&&ch<='9')ret=ret*10+ch-'0',ch=getchar();
	return ret*bj;
}
struct node
{
	int st,ed;
}w[100005],a[100005];
int n,m,f[100005][21]={0},top=0;
bool cmp(node x,node y)
{
	return x.st<y.st||(x.st==y.st&&x.ed<y.ed);
}
void init()
{
	n=read();m=read();
	for(int i=1;i<=n;i++)
	{
		w[i].st=read();
		w[i].ed=read();
	}
	sort(w+1,w+n+1,cmp);
	for(int i=1;i<=n;i++)
	{
		while(a[top].ed>w[i].ed&&top)top--;
		a[++top]=w[i];
	}
	n=top;
}

void ST()
{
	int j=1;
	a[n+1]=(node){0x7fffffff,0x7fffffff};
	f[n+1][0]=n+1;
	for(int i=1;i<=n;i++)
	{
		while(j<=n&&a[j].st<a[i].ed)j++;
		f[i][0]=j;
	}
	for(int j=1;j<=log2(n)+1;j++)
	{
		for(int i=1;i<=n+1;i++)f[i][j]=f[f[i][j-1]][j-1];
	}
}
void work()
{
	int x=read(),y=read();
	int l=1,r=n+1,ans;
	while(l<=r)
	{
		int mid=(l+r)>>1;
		if(a[mid].st>=x)
		{
			ans=mid;
			r=mid-1;
		}
		else l=mid+1;
	}
	int k=log2(n)+1,ret=0;
	for(int i=k;i>=0;i--)
	{
		if(a[f[ans][i]].ed<=y)
		{
			ret+=(1<<i);
			ans=f[ans][i];
		}
	}
	printf("%d\n",ret+(a[ans].ed<=y));
}
int main()
{
	init();
	ST();
	while(m--)work();
return 0;
}

注意n+1这个哨兵，不然死得惨。
- -感觉这是我写得最烂的博客了