题目大意:
小A的楼房外有一大片施工工地,工地上有N栋待建的楼房。每天,这片工地上的房子拆了又建、建了又拆。他经常无聊地看着窗外发呆,数自己能够看到多少栋房子。
为了简化问题,我们考虑这些事件发生在一个二维平面上。小A在平面上(0,0)点的位置,第i栋楼房可以用一条连接(i,0)和(i,Hi)的线段表示,其中Hi为第i栋楼房的高度。如果这栋楼房上任何一个高度大于0的点与(0,0)的连线没有与之前的线段相交,那么这栋楼房就被认为是可见的。
施工队的建造总共进行了M天。初始时,所有楼房都还没有开始建造,它们的高度均为0。在第i天,建筑队将会将横坐标为Xi的房屋的高度变为Yi(高度可以比原来大—修建,也可以比原来小—拆除,甚至可以保持不变—建筑队这天什么事也没做)。请你帮小A数数每天在建筑队完工之后,他能看到多少栋楼房?
思路:
这好像也可以说是线段树维护区间最长连续上升的模板题。。。我们用斜率来表示高度,用一个线段树来维护某个区间内的高度最大值和最长连续上升。修改高度的时候只需要普通的线段树更新就好了,但是我们需要更新最长的长度,考虑一个区间[l,r],它的lc和rc的信息都是已知的,那么就是求rc这个区间内大于Max[lc]的连续上升的数的个数。
现在我们开始递归修改,记h为Max[lc]。
1、若右区间的左子区间最大值大于等于h,那么它的右子区间的个数就可以直接求出来,继续递归lc的左子区间即可。
2、若右区间的左子区间最大值小于h,那么他的左子区间的贡献必定为0,所以继续递归右子区间即可。
总的时间复杂度
O(n∗log2n)
O
(
n
∗
log
2
n
)
(比我的分块不知道快了多少。。。)
// luogu-judger-enable-o2
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdio>
using namespace std;
void File(){
freopen("bzoj2957.in","r",stdin);
freopen("bzoj2957_Segment_tree.out","w",stdout);
}
template<typename T>
void read(T &x){
T _=0;char __=getchar();
while(!isdigit(__))__=getchar();
while(isdigit(__))_=(_<<1)+(_<<3)+(__^'0'),__=getchar();
x=_;
}
template<typename T>
void write(T x,char c){
if(x==0){
putchar('0');
putchar(c);
return;
}
T y=10;
int len=1;
while(y<=x)y=(y<<1)+(y<<3),++len;
while(len--){
y=y/10;
putchar(x/y+48);
x=x%y;
}
putchar(c);
}
template<typename T>void chckmax(T &_,T __){_=_>__ ? _ : __;}
template<typename T>void chckmin(T &_,T __){_=_<__ ? _ : __;}
#define REP(i,a,b) for(register int i=a;i<=b;++i)
#define DREP(i,a,b) for(register int i=a;i>=b;--i)
#define MREP(i,x) for(register int i=beg[x];i;i=E[i].last)
#define ll long long
const int maxn=1e5+10;
int n,m;
struct Segment_Tree{
#define mid ((l+r)>>1)
#define lson rt<<1,l,mid
#define rson (rt<<1)|1,mid+1,r
double Max[maxn<<2];
int cnt[maxn<<2];
int cal(int rt,int l,int r,double h){
if(l==r)return Max[rt]>h ? 1 : 0;
if(Max[rt<<1]<=h)return cal(rson,h);
return cal(lson,h)+cnt[rt]-cnt[rt<<1];
}
void update(int rt,int l,int r,int pos,double x){
if(l==r){
Max[rt]=x;
cnt[rt]=1;
return;
}
else if(pos<=mid)update(lson,pos,x);
else update(rson,pos,x);
Max[rt]=Max[rt<<1]>Max[(rt<<1)|1] ? Max[rt<<1] : Max[(rt<<1)|1];
cnt[rt]=cnt[rt<<1]+cal(rson,Max[rt<<1]);
}
}T;
int main(){
//File();
scanf("%d%d",&n,&m);
REP(i,1,m){
int x,y;
scanf("%d%d",&x,&y);
T.update(1,1,n,x,((double)y)/x);
write(T.cnt[1],'\n');
}
return 0;
}