BZOJ 2957 楼房重建

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

　　原以为不过是道水题，是道线段树裸题。不就是个长度为n的序列，m个操作，每个操作单点修改（delta=y*1.0/x，斜率）然后询问root罢了。

　　不用pushdown，但是update确实让人很迷。我最开始以为自己的update是O（lg n）的，但最后居然成了O（n）的（逻辑不严谨）。啊啊啊啊啊！

　　注意到题目要求，nd->sum=nd->ls->sum+find(nd->rs,nd->ls->max);原来居然还存了一段线段的min值，结果发现并不需要。find(nd,val)指一段线段中比val大的能看见的总数。递归分解，若nd->ls->max<=nd，那么左边就注定为0了，再去find（nd->rs,val）。否则右边没有影响，则return nd->num-nd->ls->num+find（nd->ls,val），这样就可以保证每次范围减半，就是O（lg n）的。

　　代码如下：

 1 /**************************************************************
 2     Problem: 2957
 3     User: Doggu
 4     Language: C++
 5     Result: Accepted
 6     Time:1424 ms
 7     Memory:4732 kb
 8 ****************************************************************/
 9  
10 #include <cstdio>
11 #include <cstring>
12 #include <algorithm>
13  
14 template<class T>inline void readin(T &res) {
15     static char ch;T flag=1;
16     while((ch=getchar())<‘0‘||ch>‘9‘)if(ch==‘-‘)flag=-1;
17     res=ch-48;while((ch=getchar())>=‘0‘&&ch<=‘9‘)res=(res<<1)+(res<<3)+ch-48;res*=flag;
18 }
19  
20 const int N = 100100;
21 struct Node {
22     double max;
23     int num;Node *ls, *rs;
24 }pool[N<<1], *tail=pool, *root;
25  
26 Node *build(int lf,int rg) {
27     Node *nd=tail++;
28     if(lf==rg) {
29         nd->max=nd->num=0;
30     } else {
31         int mid=(lf+rg)>>1;
32         nd->ls=build(lf,mid);
33         nd->rs=build(mid+1,rg);
34         nd->max=nd->num=0;
35     }
36     return nd;
37 }
38 int find(Node *nd,int lf,int rg,double delta) {
39     if(lf==rg) {
40         return nd->max>delta;
41     }
42     int mid=(lf+rg)>>1, ans=0;
43     if(nd->ls->max<=delta) return find(nd->rs,mid+1,rg,delta);
44     return nd->num-nd->ls->num+find(nd->ls,lf,mid,delta);
45 }
46 void modify(Node *nd,int lf,int rg,int pos,double delta) {
47     if(lf==rg) {
48         nd->max=delta;
49         nd->num=1;return ;
50     }
51     int mid=(lf+rg)>>1;
52     if(pos<=mid) modify(nd->ls,lf,mid,pos,delta);
53     else modify(nd->rs,mid+1,rg,pos,delta);
54     nd->num=nd->ls->num+find(nd->rs,mid+1,rg,nd->ls->max), nd->max=std::max(nd->ls->max,nd->rs->max);
55 }
56  
57 int main() {
58     int n, m, x, y;
59     readin(n);readin(m);
60     root=build(1,n);
61     for( int i = 1; i <= m; i++ ) {
62         readin(x);readin(y);
63         modify(root,1,n,x,y*1.0/x);
64         printf("%d\n",root->num);
65     }
66     return 0;
67 }

线段树

　　不过，天算不如人算。没有想到，居然可以直接上暴力！当然，前提是要分块。每一次把修楼点所属的块暴力重构。这是有策略的，扫一遍保证单调性，是O（sqrt（n））的，之前做O（sqrt（n）log n）便会TLE。之后枚举每个块，upper_bound之前的最大值，加上并更新最值，而这是O（sqrt（n）log n）的。

　　我用vector来存块，因为不熟悉，所以卡了很久。调试时间比线段树久，但确实是不可多得的绝佳经历。

　　代码如下（似乎比线段树短）：

 1 /**************************************************************
 2     Problem: 2957
 3     User: Doggu
 4     Language: C++
 5     Result: Accepted
 6     Time:4052 ms
 7     Memory:3976 kb
 8 ****************************************************************/
 9  
10 #include <cstdio>
11 #include <cmath>
12 #include <algorithm>
13 #include <vector>
14 template<class T>inline void readin(T &res) {
15     static char ch;T flag=1;
16     while((ch=getchar())<‘0‘||ch>‘9‘)if(ch==‘-‘)flag=-1;
17     res=ch-48;while((ch=getchar())>=‘0‘&&ch<=‘9‘)res=(res<<1)+(res<<3)+ch-48;res*=flag;
18 }
19 const int N = 100100;
20 int belong[N], pos, bound, num[N], tnum;
21 double aa[N], tmp;
22 std::vector<double> bb[N];
23 int main() {
24     int n, m, S, x, y;
25     readin(n);readin(m);
26     S=(int)floor(sqrt(n*log(n)/log(2)/2));
27     pos=1-S;for( int i = 1; i <= n; i++ ) {
28         belong[i]=pos;
29         if(i-pos>=S) pos=i, num[i]=++tnum, belong[i]=pos;
30     }
31     for( int i = 1; i <= m; i++ ) {
32         readin(x);readin(y);
33         aa[x]=y*1.0/x;
34         bound=std::min(belong[x]+S-1,n);
35         tmp=0;bb[num[belong[x]]].clear();for( int i = belong[x]; i <= bound; i++ ) if(aa[i]>tmp) bb[num[belong[x]]].push_back(aa[i]),tmp=aa[i];
36         tmp=0;int ans=0;
37         for( int i = 1; i <= tnum; i++ ) {
38             bound=bb[i].size()-1;
39             pos=std::upper_bound(bb[i].begin(),bb[i].end(),tmp)-bb[i].begin();
40             ans+=bound+1-pos;
41             if(bound>=0) tmp=std::max(tmp,bb[i][bound]);
42         }
43         printf("%d\n",ans);
44     }
45     return 0;
46 }

暴力分块+vector+upper_bound

　　恩。最后讨论一下速度。似乎暴力只慢了3倍左右。　　

线段树	4732 kb	1424 ms	1523 B
暴力分块+vector+upper_bound	3976 kb	4052 ms	1232 B

BZOJ 2957 楼房重建

2957: 楼房重建

Description

Input

Output

Sample Input

Sample Output

Source