线段树-楼房重建-洛谷-P4198

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/weixin_37517391/article/details/82969725

楼房重建

题目大意

小A的楼房外有一大片施工工地,工地上有N栋待建的楼房。每天,这片工地上的房子拆了又建、建了又拆。他经常无聊地看着窗外发呆,数自己能够看到多少栋房子。

为了简化问题,我们考虑这些事件发生在一个二维平面上。小A在平面上(0,0)点的位置,第i栋楼房可以用一条连接(i,0)和(i,Hi)的线段表示,其中Hi为第i栋楼房的高度。如果这栋楼房上任何一个高度大于0的点与(0,0)的连线没有与之前的线段相交,那么这栋楼房就被认为是可见的。

施工队的建造总共进行了M天。初始时,所有楼房都还没有开始建造,它们的高度均为0。在第i天,建筑队将会将横坐标为Xi的房屋的高度变为Yi(高度可以比原来大—修建,也可以比原来小—拆除,甚至可以保持不变—建筑队这天什么事也没做)。请你帮小A数数每天在建筑队完工之后,他能看到多少栋楼房?

题目解答

给每个楼房表示成从原点到其的斜率,我们要求的就是从原点开始,斜率严格增大的序列长度.

用线段树维护.

维护区间 [ l , r ] [l,r] 中以 a [ l ] a[l] 开始的递增序列的长度.

合并的时候,判断 [ l , m i d ] [l,mid] 的最大值与 a [ m i d + 1 ] a[mid+1] 之间的关系.

如果 m a x ( a [ l , m i d ] ) < a [ m i d + 1 ] max(a[l,mid]) < a[mid+1] ,则 s e g [ o ] = s e g [ l c h ] + s e g [ r c h ] seg[o] = seg[lch] + seg[rch]

否的话,需要在 [ m i d + 1 , r ] [mid+1,r] 中寻找递增序列的第一个数 > m a x ( a [ l , m i d ] ) >max(a[l,mid]) 的递增序列的长度,用 f i n d ( r c h , m i d + 1 , r , m a x ( a [ l , m i d ] ) ) find(rch,mid+1,r,max(a[l,mid])) 表示,那么合并就是
s e g [ o ] = s e g [ l c h ] + f i n d ( r c h , m i d + 1 , r , m a x ( a [ l , m i d ] ) ) seg[o] = seg[lch] + find(rch,mid+1,r,max(a[l,mid]))

下面我们需要实现 f i n d find 操作

边界判断,如果 f i n d find 的区间长度为 1 1 ,那么直接判断范围即可.

我们需要对传给 f i n d ( o , l , r , m x ) find(o,l,r,mx) 的区间 [ l , r ] [l,r] 再分一次,递归查找.

如果 m x m a x ( a [ l , r ] ) mx \le max(a[l,r]) ,那么 r e s = s e g [ o ] s e g [ o . l c h ] + f i n d ( o . l c h , l , m i d , x ) res = seg[o] - seg[o.lch] + find(o.lch,l,mid,x)
如果 m x > m a x ( a [ l , r ] ) mx > max(a[l,r]) ,那么 r e s = f i n d ( o . r c h , m i d + 1 , r , x ) res = find(o.rch,mid+1,r,x)
因为 s e g seg 是从下网上合并的,所以传给 f i n d find s e g [ o ] seg[o] 是计算好了的,所以可以直接使用.

时间复杂度是 O ( n l o g n l o g n ) O(nlognlogn)

代码

#include <iostream>

const int N = 100007;
int n,m;

double T[N<<2],M[N<<2];


inline int find(int o,int l,int r,double v) {
    if(l == r) 
        return v < M[o];
    int mid = (l + r) >> 1;
    if(v < M[o<<1]) 
        return find(o<<1,l,mid,v) + T[o] - T[o<<1];
    else 
        return find(o<<1|1,mid+1,r,v);
}
inline void maintain(int o,int l,int r) {
    M[o] = std::max(M[o<<1],M[o<<1|1]);
    int mid = (l + r) >> 1;
    T[o] = T[o<<1] + find(o<<1|1,mid+1,r,M[o<<1]);

}
inline void change(int o,int l,int r,int pos,double val) {
    if(l == r) {
        M[o] = val;T[o] = 1;
        return ;
    }
    int mid = (l + r) >> 1;
    if(pos <= mid) change(o<<1,l,mid,pos,val);
    else change(o<<1|1,mid+1,r,pos,val);
    maintain(o,l,r);
}
inline void build(int o,int l,int r) {
    if(l == r) {
        M[o] = 0;T[o] = 1;
        return ;	
    }
    int mid = (l + r) >> 1;
    build(o<<1,l,mid);
    build(o<<1|1,mid+1,r);
}


int main() {
    std::ios::sync_with_stdio(false);
    std::cin >> n >> m;
    build(1,1,n);
    while(m--) {
        int pos;double val;
        std::cin >> pos >> val;
        change(1,1,n,pos,val/pos);
        std::cout << find(1,1,n,0.) << std::endl;
    }
    return 0;
}

猜你喜欢

转载自blog.csdn.net/weixin_37517391/article/details/82969725