【BZOJ 2957】 2957: 楼房重建 （线段树）

2957: 楼房重建

Time Limit: 10 Sec  Memory Limit: 256 MB
Submit: 1753  Solved: 841

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

HINT

Source

中国国家队清华集训 2012-2013 第一天

【分析】

　　其实这个很明显是线段树维护的，但是维护方法和求值方法其实都不是很传统，所以我想不到啊。。。

　　分析题目知道先把每栋楼跟原点的三角形的斜率求出来，一个楼有贡献当且仅当他的斜率比他前面的楼的斜率都大，即大于前面的楼的斜率的最大值。

　　线段树维护一个ans和mx，mx表示这个区间的斜率的最大值，ans表示只考虑这个区间的时候的答案。

　　更新ans的时候要带一个外区间的mx，表示前面区间的mx

　　若这个mx>=左孩子的mx，则左边都不能选，只考虑右边，变成子问题。

　　若这个mx<左孩子的mx，则为整个区间的ans-左区间的ans+左区间带mx时的ans（这个也很容易意会吧）

　　【说明我的线段树还是有漏洞啊TAT

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 using namespace std;
 #define Maxn 100010

 double mymax(double x,double y) {return x>y?x:y;}

 struct node
 {
     int l,r,lc,rc;
     double mx;
     int ans;
 }tr[Maxn*];
 int len;

 int build(int l,int r)
 {
     int x=++len;
     tr[x].l=l;tr[x].r=r;tr[x].mx=;tr[x].ans=;
     if(l!=r)
     {
         int mid=(l+r)>>;
         tr[x].lc=build(l,mid);
         tr[x].rc=build(mid+,r);
     }
     else tr[x].lc=tr[x].rc=;
     return x;
 }

 int query(int x,double mx)
 {
     if(tr[x].l==tr[x].r) return tr[x].mx>mx;
     int mid=(tr[x].l+tr[x].r)>>;
     if(tr[tr[x].lc].mx<mx) return query(tr[x].rc,mx);
     return tr[x].ans-tr[tr[x].lc].ans+query(tr[x].lc,mx);
 }

 void change(int x,int y,double z)
 {
     if(tr[x].l==tr[x].r)
     {
         tr[x].mx=z;
         tr[x].ans=;
         return;
     }
     int mid=(tr[x].l+tr[x].r)>>;
     if(y<=mid) change(tr[x].lc,y,z);
     else change(tr[x].rc,y,z);
     tr[x].mx=mymax(tr[tr[x].lc].mx,tr[tr[x].rc].mx);
     tr[x].ans=tr[tr[x].lc].ans+query(tr[x].rc,tr[tr[x].lc].mx);
 }

 int main()
 {
     int n,m;
     scanf("%d%d",&n,&m);
     build(,n);
     for(int i=;i<=m;i++)
     {
         int x,y;
         scanf("%d%d",&x,&y);
         change(,x,y*1.0/x);
         printf("%d\n",tr[].ans);
     }
     return ;
 }

2017-03-25 09:50:54