题目链接:bzoj2957
题目大意:
小A的楼房外有一大片施工工地,工地上有N栋待建的楼房。每天,这片工地上的房子拆了又建、建了又拆。他经常无聊地看着窗外发呆,数自己能够看到多少栋房子。
为了简化问题,我们考虑这些事件发生在一个二维平面上。小A在平面上(0,0)点的位置,第i栋楼房可以用一条连接(i,0)和(i,Hi)的线段表示,其中Hi为第i栋楼房的高度。如果这栋楼房上任何一个高度大于0的点与(0,0)的连线没有与之前的线段相交,那么这栋楼房就被认为是可见的。
施工队的建造总共进行了M天。初始时,所有楼房都还没有开始建造,它们的高度均为0。在第i天,建筑队将会将横坐标为Xi的房屋的高度变为Yi(高度可以比原来大—修建,也可以比原来小—拆除,甚至可以保持不变—建筑队这天什么事也没做)。请你帮小A数数每天在建筑队完工之后,他能看到多少栋楼房?
题解:
线段树
mx[i]维护i这个点的区间中的最大斜率,ans[i]维护i这个点的区间中能被看到的楼房数
改的话直接修改。对一段区间now来说,左半区间lc的ans是一定可以对它产生贡献的,那么设左半区间的斜率的最大值为mv,递归进右半区间rc看比mv大的楼房个数。
同样将右区间分成左子区间和右子区间。
如果mv≥mx[rc],即右半区间的最大值没有左半区间的最大值大,那么右半区间的所有楼房都不能被看到了,就return 0.
否则看左子区间的mx与mv的比较:
若mv≥mx[左子区间],那么左子区间就没有贡献,就递归进右子区间看比mv大的个数。
若mx[左子区间]≥mv,那么右子区间对区间now的贡献跟它对右区间的贡献是一样的,因为左区间并没有对它造成影响即ans[rc]-ans[左子区间]
然后就递归进左子区间看比mv大的个数。
都加起来就好了。
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define maxn 400010
double mx[maxn];int ans[maxn];
double mymax(double x,double y){return (x>y)?x:y;}
int calc(int now,double mv,int l,int r)
{
if (l==r) return (mx[now]>mv);
if (mv>=mx[now]) return 0;
int mid=(l+r)>>1,lc=now*2,rc=now*2+1;
if (mv>=mx[lc]) return calc(rc,mv,mid+1,r);
return ans[now]-ans[lc]+calc(lc,mv,l,mid);
}
void change(int now,int l,int r,int x,double k)
{
if (l==r)
{
mx[now]=k;ans[now]=1;
return;
}
int mid=(l+r)>>1,lc=now*2,rc=now*2+1;
if (x<=mid) change(lc,l,mid,x,k);
else change(rc,mid+1,r,x,k);
ans[now]=ans[lc]+calc(rc,mx[lc],mid+1,r);
mx[now]=mymax(mx[lc],mx[rc]);
}
int main()
{
//freopen("a.in","r",stdin);
//freopen("a.out","w",stdout);
int n,m,i,x,y;
scanf("%d%d",&n,&m);
for (i=1;i<=m;i++)
{
scanf("%d%d",&x,&y);
double sp=y*1.0/x;
change(1,1,n,x,sp);
printf("%d\n",ans[1]);
}
return 0;
}