【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

HINT

Source

中国国家队清华集训 2012-2013 第一天

【解析】

线段树

先求出每栋楼到原点的斜率，线段树记录区间最大值c和能看见的楼数d

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
struct node
{
	int l,r,lc,rc,d;
	double c;
}tr[410000];int len;
void bt(int l,int r)
{
	len++;int now=len;
	tr[now].l=l;tr[now].r=r;
	tr[now].lc=tr[now].rc=-1;
	tr[now].c=-1.0;tr[now].d=0;
	if(l==r) return ;
	else
	{
		int mid=(l+r)/2;
		tr[now].lc=len+1;bt(l,mid);
		tr[now].rc=len+1;bt(mid+1,r);
	}
}
int findmax(int now,double k)
{
	if(tr[now].l==tr[now].r)
	{
		if(k<tr[now].c) return 1;
		else return 0;
	}
	int lc=tr[now].lc,rc=tr[now].rc;
	if(k>tr[lc].c) return findmax(rc,k);
	return tr[now].d-tr[lc].d+findmax(lc,k);
}
void change(int now,int x,double k)
{
	if(tr[now].l==tr[now].r){tr[now].c=k;tr[now].d=1;return ;}
	int lc=tr[now].lc,rc=tr[now].rc;
	int mid=(tr[now].l+tr[now].r)/2;
	if(x<=mid) change(lc,x,k);
	else if(x>mid) change(rc,x,k);
	if(tr[lc].c>tr[rc].c){tr[now].c=tr[lc].c;tr[now].d=tr[lc].d;}
	else{tr[now].c=tr[rc].c;tr[now].d=tr[lc].d+findmax(rc,tr[lc].c);}
}
int main()
{
	int n,m,x,y;
	scanf("%d%d",&n,&m);
	len=0;bt(1,n);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&x,&y);
		change(1,x,(double)y/x);
		printf("%d\n",tr[1].d);
	}
	return 0;
}

线段树

先求出每栋楼到原点的斜率，线段树记录区间最大值c和能看见的楼数d