【题解】【BZOJ2957楼房重建】（线段树）

Description

小A的楼房外有一大片施工工地，工地上有N栋待建的楼房。每天，这片工地上的房子拆了又建、建了又拆。他经常无聊地看着窗外发呆，数自己能够看到多少栋房子。
　　为了简化问题，我们考虑这些事件发生在一个二维平面上。小A在平面上(0,0)点的位置，第i栋楼房可以用一条连接(i,0)和(i,Hi)的线段表示，其中Hi为第i栋楼房的高度。如果这栋楼房上任何一个高度大于0的点与(0,0)的连线没有与之前的线段相交，那么这栋楼房就被认为是可见的。
　　施工队的建造总共进行了M天。初始时，所有楼房都还没有开始建造，它们的高度均为0。在第i天，建筑队将会将横坐标为Xi的房屋的高度变为Yi(高度可以比原来大—修建，也可以比原来小—拆除，甚至可以保持不变—建筑队这天什么事也没做)。请你帮小A数数每天在建筑队完工之后，他能看到多少栋楼房？
Input
　　第一行两个正整数N,M
　　接下来M行，每行两个正整数Xi,Yi
Output
　　M行，第i行一个整数表示第i天过后小A能看到的楼房有多少栋
Sample Input
3 4
2 4
3 6
1 1000000000
1 1
Sample Output
1
1
1
2
数据约定
　　对于所有的数据1<=Xi<=N，1<=Yi<=10^9
N,M<=100000

这道题一看就是线段树（当然分块也是可以的），单点修改，到是没有什么问题，最主要的还是要想出如何统计答案；

我们在每个节点统计两个信息：ans[t],mx[t]，分别表示只看该区间的答案，和这个区间斜率的最大值
mx数组还是很好转移的：mx[t]=max(mx[lson[t]],mx[rson[t]]);
然后就是ans数组的区间合并了：
显然左边的区间肯定是会被保留的，且它最后一个楼房就是斜率最大的那个；
然后我们就基于这个进行分类讨论：
若右区间的斜率最大都没有左区间的大，那么答案就是左区间，即:ans[t]=ans[lsont]];
若右区间的最大斜率比左区间的最大斜率大，我们就可以往下递归左儿子，知道右区间的左儿子的最大斜率比左区间小，就递归右儿子，统计答案即可。
那么这样做的正确性怎么保证的呢？？
这其实就是一种类似于二分查找的方法，找到第一个斜率大于左区间最大斜率的节点，要明确：二分查找的前提是数列有序，而统计在ans中的节点的斜率是单调递增的。所以正确性显然。

就看一下代码吧：

#include<bits/stdc++.h>

using namespace std;

const int maxn = 101000; 
int n,m,lson[maxn*4],rson[maxn*4],ans[maxn*4],ndnum;

double mx[maxn*4];

int cut(int t,int l,int r,double bz)
{
	if(l==r) return mx[t]>bz;
	int mid=(l+r)/2;
	if(mx[lson[t]]<=bz) return cut(rson[t],mid+1,r,bz);
	return cut(lson[t],l,mid,bz)+ans[t]-ans[lson[t]];
}
inline void update(int t,int l,int r)
{
	mx[t]=max(mx[lson[t]],mx[rson[t]]);
	if(mx[rson[t]]<mx[lson[t]]) ans[t]=ans[lson[t]];
	else 
		ans[t]=ans[lson[t]]+cut(rson[t],((l+r)/2)+1,r,mx[lson[t]]);
}


void build(int &t,int l,int r)
{
	t=++ndnum;
	if(l==r) return;
	int mid=(l+r)/2;
	build(lson[t],l,mid);
	build(rson[t],mid+1,r);
}

void modify(int t,int l,int r,int pos,double val)
{
	if(l==r)
	{
		ans[t]=1;
		mx[t]=val;
		return; 
	}
	int mid=(l+r)/2;
	if(pos<=mid) modify(lson[t],l,mid,pos,val);
	else modify(rson[t],mid+1,r,pos,val);
	update(t,l,r);
}
int main()
{
	scanf("%d%d",&n,&m);
	int root;
	build(root,1,n);
	for(int i=1;i<=m;++i)
	{
		int x,y;
		scanf("%d%d",&x,&y);
		double yy=y*1.0/x;
		modify(root,1,n,x,yy);
		printf("%d\n",ans[1]);
	}
	return 0;
}
/*
3 4
2 4
3 6
1 1000000000
1 1
*/