BZOJ 2957: 楼房重建

Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 3287 Solved: 1549
[Submit][Status][Discuss]
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

解题思路

考试题,当时暴力70分。正解线段树,貌似分块也能做。发现这道题跟斜率有关,如果要被看见就必须斜率大于前面所有的,所以就形成一个单调递增的序列,我们用线段树来维护这个序列。设cnt[x]是x这个区间的最长单调序列,还需要一个mx数组来记录每个区间的斜率最大值,每次进行修改时,首先change操作找到这个点进行修改,因为这个点只影响它右边的,所以直接将左边的cnt加入答案,继续count操作递归数右边的值。右边的又可以看成一个子问题,如果右边这个区间的左边的最大值没有x大,就直接递归右边区间的右边。否则将右边的答案加上递归左边。(有点绕。。。。)

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>

using namespace std;
const int MAXN = 400005;

inline int rd(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=(x<<1)+(x<<3)+ch-48;ch=getchar();}
    return x*f;
}

int n,m,cnt[MAXN];
double mx[MAXN];

inline int count(int l,int r,int x,double k){
    if(l==r) return mx[x]>k;
    int mid=(l+r)>>1;
    if(mx[x<<1]<=k) return count(mid+1,r,x<<1|1,k);
    else return cnt[x]-cnt[x<<1]+count(l,mid,x<<1,k);
    //注意cnt[x]-cnt[x<<1]不能直接写成cnt[x<<1|1],因为cnt数组表示的只是
    //每个区间的递增序列,与其他区间无关,所以左边可能挡住右边。
}

inline void change(int l,int r,int x,int q,double k){
    if(l==r){
        mx[x]=k;
        cnt[x]=1;
        return;
    }
    int mid=(l+r)>>1;
    if(q<=mid) change(l,mid,x<<1,q,k);
    else change(mid+1,r,x<<1|1,q,k);
    mx[x]=mx[x<<1]>mx[x<<1|1]?mx[x<<1]:mx[x<<1|1];
    cnt[x]=cnt[x<<1]+count(mid+1,r,x<<1|1,mx[x<<1]);
}

int main(){
    n=rd();m=rd();
    for(register int i=1;i<=m;i++){
        int x=rd(),y=rd();
        double k=(double)y/(double)x;
        change(1,n,1,x,k);
        printf("%d\n",cnt[1]);
    }
}
posted @ 2018-07-03 17:21  Monster_Qi  阅读(130)  评论(0编辑  收藏  举报