[BZOJ 2957]楼房重建(线段树)

Description

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

为了简化问题,我们考虑这些事件发生在一个二维平面上。小A在平面上(0,0)点的位置,第i栋楼房可以用一条连接(i,0)和(i,Hi)的线段表示,其中Hi为第i栋楼房的高度。如果这栋楼房上任何一个高度大于0的点与(0,0)的连线没有与之前的线段相交,那么这栋楼房就被认为是可见的。
施工队的建造总共进行了M天。初始时,所有楼房都还没有开始建造,它们的高度均为0。在第i天,建筑队将会将横坐标为Xi的房屋的高度变为Yi(高度可以比原来大---修建,也可以比原来小---拆除,甚至可以保持不变---建筑队这天什么事也没做)。请你帮小A数数每天在建筑队完工之后,他能看到多少栋楼房?

Solution

把每栋楼房高度角的tan值放在线段树上维护,然后就变成了一个动态求最长上升子序列的问题

于是线段树每个节点都维护max和ans,这个ans可以由左子树的ans+右子树中大于左边max的最长上升子序列得到,后面的这一部分可以在右子树进行二分

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#define MAXN 100005
using namespace std;
int n,m;
int read()
{
    int x=0,f=1;char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}
struct Node
{
    int l,r,ans;
    double maxn;
}t[MAXN*4];
void build(int idx,int l,int r)
{
    t[idx].l=l,t[idx].r=r,t[idx].maxn=t[idx].ans=0;
    if(l==r)return;
    int mid=(l+r)>>1;
    build(idx<<1,l,mid),build(idx<<1|1,mid+1,r);
}
int update(int idx,double v)
{
    if(t[idx].l==t[idx].r)return t[idx].maxn>v;
    if(t[idx<<1].maxn<=v)return update(idx<<1|1,v);
    return t[idx].ans-t[idx<<1].ans+update(idx<<1,v);
}
void add(int idx,int p,double x)
{
    if(t[idx].l==t[idx].r){t[idx].maxn=x;t[idx].ans=1;return;}
    int mid=(t[idx].l+t[idx].r)>>1;
    if(p<=mid)add(idx<<1,p,x);else add(idx<<1|1,p,x);
    t[idx].maxn=max(t[idx<<1].maxn,t[idx<<1|1].maxn);
    t[idx].ans=t[idx<<1].ans+update(idx<<1|1,t[idx<<1].maxn);
}
int main()
{
    n=read(),m=read();
    build(1,1,n);
    for(int i=1;i<=m;i++)
    {
        int x=read(),y=read();
        add(1,x,(double)y/x);
        printf("%d
",t[1].ans);
    }
    return 0;
}
原文地址:https://www.cnblogs.com/Zars19/p/6964199.html