原题链接:https://www.luogu.com.cn/problem/P4198
楼房重建
题目描述
小 A 的楼房外有一大片施工工地,工地上有 N N N 栋待建的楼房。每天,这片工地上的房子拆了又建、建了又拆。他经常无聊地看着窗外发呆,数自己能够看到多少栋房子。
为了简化问题,我们考虑这些事件发生在一个二维平面上。小 A 在平面上 ( 0 , 0 ) (0,0) (0,0) 点的位置,第 i i i 栋楼房可以用一条连接 ( i , 0 ) (i,0) (i,0) 和 ( i , H i ) (i,H_i) (i,Hi) 的线段表示,其中 H i H_i Hi 为第 i i i 栋楼房的高度。如果这栋楼房上任何一个高度大于 0 0 0 的点与 ( 0 , 0 ) (0,0) (0,0) 的连线没有与之前的线段相交,那么这栋楼房就被认为是可见的。
施工队的建造总共进行了 M M M 天。初始时,所有楼房都还没有开始建造,它们的高度均为 0 0 0。在第 i i i 天,建筑队将会将横坐标为 X i X_i Xi 的房屋的高度变为 Y i Y_i Yi(高度可以比原来大—修建,也可以比原来小—拆除,甚至可以保持不变—建筑队这天什么事也没做)。请你帮小 A 数数每天在建筑队完工之后,他能看到多少栋楼房?
输入格式
第一行两个正整数 N , M N,M N,M。
接下来 M M M 行,每行两个正整数 X i , Y i X_i,Y_i Xi,Yi 。
输出格式
M M M 行,第 i i i 行一个整数表示第 i i i 天过后小 A 能看到的楼房有多少栋。
输入输出样例
输入 #1
3 4
2 4
3 6
1 1000000000
1 1
输出 #1
1
1
1
2
说明/提示
对于 100 % 100\% 100% 的数据, 1 ≤ X i ≤ N 1 \le X_i \le N 1≤Xi≤N, 1 ≤ Y i ≤ 1 0 9 1 \le Y_i \le 10^9 1≤Yi≤109, 1 ≤ N , M ≤ 1 0 5 1\le N,M \le 10^5 1≤N,M≤105。
题解
首先,要判断一个楼房是否能被看到,肯定不能去判断什么“是否有交点”,稍微一想,就可以想到使用斜率,当某座楼左侧没有斜率比它大的楼时,它就可以被看到。
当把问题“斜率化”之后,问题就变为求一个上升子序列的长度,该子序列中的每个元素大于原序列中位于它左侧的任意元素,需要支持单点修改和全局查询,于是尝试线段树。
因为只有单点修改和全局查询,所以整棵线段树不需要什么打标记、pushdown
什么的,只要能合并子区间信息就行了。
那么现在考虑,已知两个子节点的上升子序列长度,如何得到父节点的最长上升子序列长度。
先看左儿子,因为是从左侧向右侧观测,所以左儿子的上升子序列一定是父节点上升子序列的一部分,而右儿子的上升子序列要拼接到左儿子的序列上,就应该是从大于左儿子上升子序列最右端的值开始接到左儿子的上升子序列上。要找到这个“划分点”,就需要讨论左儿子的最大值
k
k
k和右儿子
v
v
v中值的关系。如果
k
≥
k\ge
k≥ tree[v].mx
,那么右儿子的序列会被全部挡住,直接返回左儿子即可;如果
k
<
k<
k< tree[v].mx
,则进一步讨论
k
k
k和右儿子的左儿子的关系:若
k
≥
k\ge
k≥tree[ls].mx
那么右儿子的左儿子是会被挡住的,直接递归右儿子;反之,“划分点”就在右儿子的左儿子里,就递归左儿子,找到断点后加上右儿子的右儿子中属于右儿子本身的上升子序列的部分(比较拗口,用符号表示就是不能加tree[rs].len
,而是加上tree[v].len-tree[ls].len
)。
由于这个更新区间信息的过程是递归的,所以总复杂度为 O ( n l o g 2 2 n ) O(nlog_2^2n) O(nlog22n)。
代码
#include<bits/stdc++.h>
#define ls v<<1
#define rs v<<1|1
using namespace std;
const int M=1e5+5;
struct node{int le,ri,len;double mx;}tree[M<<2];
int n,m;
void build(int v,int le,int ri)
{
tree[v].le=le,tree[v].ri=ri;
if(le==ri)return;
int mid=le+ri>>1;
build(ls,le,mid);build(rs,mid+1,ri);
}
int ask(int v,double k)
{
if(tree[v].mx<=k)return 0;
if(tree[v].le==tree[v].ri)return tree[v].mx>k;
if(tree[ls].mx<=k)return ask(rs,k);
return ask(ls,k)+tree[v].len-tree[ls].len;
}
void modify(int v,int x,double k)
{
if(tree[v].le==tree[v].ri){tree[v].mx=k,tree[v].len=1;return;}
if(x<=tree[ls].ri)modify(ls,x,k);
else modify(rs,x,k);
tree[v].mx=max(tree[ls].mx,tree[rs].mx);
tree[v].len=tree[ls].len+ask(rs,tree[ls].mx);
}
void in(){scanf("%d%d",&n,&m);}
void ac()
{
build(1,1,n);
for(int i=1,x,y;i<=m;++i)
{
scanf("%d%d",&x,&y);
modify(1,x,1.0*y/x);
printf("%d\n",tree[1].len);
}
}
int main()
{
in(),ac();
system("pause");
}