BZOJ1178 [Apio2009]CONVENTION会议中心

本文作者：ljh2000
作者博客：http://www.cnblogs.com/ljh2000-jump/
转载请注明出处，侵权必究，保留最终解释权！

Description

Siruseri政府建造了一座新的会议中心。许多公司对租 借会议中心的会堂很感兴趣，他们希望能够在里面举行会议。 对于一个客户而言，仅当在开会时能够独自占用整个会堂，他才会租借会堂。会议中心的销售主管认为：最好的策略应该是将会堂租借给尽可能多的客户。显然，有 可能存在不止一种满足要求的策略。 例如下面的例子。总共有4个公司。他们对租借会堂发出了请求，并提出了他们所需占用会堂的起止日期（如下表所示）。 开始日期 结束日期 公司1 4 9 公司2 9 11 公司3 13 19 公司4 10 17 上例中，最多将会堂租借给两家公司。租借策略分别是租给公司1和公司3，或是公司2和公司3，也可以是公司1和公司4。注意会议中心一天最多租借给一个公 司，所以公司1和公司2不能同时租借会议中心，因为他们在第九天重合了。 销售主管为了公平起见，决定按照如下的程序来确定选择何种租借策略：首先，将租借给客户数量最多的策略作为候选，将所有的公司按照他们发出请求的顺序编 号。对于候选策略，将策略中的每家公司的编号按升序排列。最后，选出其中字典序最小1的候选策略作为最终的策略。 例中，会堂最终将被租借给公司1和公司3：3个候选策略是{(1,3),(2,3),(1,4)}。而在字典序中(1,3)<(1,4)< (2,3)。 你的任务是帮助销售主管确定应该将会堂租借给哪些公司。

Input

输入的第一行有一个整数N，表示发出租借会堂申请的公司的个数。第2到第N+1行每行有2个整数。第i+1行的整数表示第i家公司申请租借的起始和终止日期。对于每个公司的申请，起始日期为不小于1的整数，终止日期为不大于10^9的整数。N≤200000

Output

输出的第一行应有一个整数M，表示最多可以租借给多少家公司。第二行应列出M个数，表示最终将会堂租借给哪些公司。

Sample Input

4
4 9
9 11
13 19
10 17

Sample Output

2
1 3

正解：倍增+贪心

解题报告：　　

　　 题目要求从若干个闭区间中取出尽可能多的区间，并且保证方案的字典序最小。

　　 在只有第一问的情况下，这就是一道贪心的经典模型了。但是我们考虑如何保证我们选出来的方案字典序最小。显然我们可以发现一个问题，如果我们从第一个区间开始考虑，假设把它选了，没有使总数量减少，那么我们肯定可以选择它（通过字典序最小的原则很容易想通）。那么我们的任务就是如何快速的检查是否不会使得总数量减少，所以我们需要做的就是快速查询区间的最大放的数量（通过预处理实现）。

　　假设get_ans函数可以快速查询区间的最大放的数量，那么如果我想放入一个[l0,r0]的区间，那么需要满足get_ans(l,r)=get_ans(l,l0-1)+get_ans(r0+1,r)+1，这个式子很容易想通。l是离l0最近的上次插入的区间的右端点，r是离r0最近的上次插入的左端点。开始我还在想怎么维护选择区间之后，这个区间的数量变化情况，事实上我们并不关心怎么变，我们只知道它不会使得数量变化即可。也就是说，我只要知道我能影响的区间怎么变，外面的区间因为没有修改，所以肯定也不会变化，换而言之当前影响区间不变则整个区间的总数量就不会变。

　　 那么我们需要预处理什么呢？首先我们需要在log n的时间内实现快速查询给定范围最大可以选取的区间数量。那么我们把区间重新排序，并且去掉互相包含的区间中较大的部分（显然可以知道这是可行的），这可以先按左端点排序，再用一个右端点单调递增的栈实现。之后对于每个保留下来的区间i，考虑用f[i][j]表示从第i个区间开始往后连续选取2^j个区间之后落在哪个区间上。这个预处理也很简单，n log n。这其实用的就是一个倍增的思想。所以对于每一次查询[l,r]，我可以找到[l,r]中的第一个区间，然后倍增的找区间，不断地往后跳并累加答案。通过上述操作可以实现查询操作。

　　 另外，我们选取了一个区间之后可以用一个set维护一下我选取的区间的端点坐标，方便我快速查询最近的坐标。

　　 整个题目细节比较多，实现起来稍显复杂，需要一点时间。

 //It is made by ljh2000
 #include <iostream>
 #include <cstdlib>
 #include <cstring>
 #include <cstdio>
 #include <cmath>
 #include <algorithm>
 #include <ctime>
 #include <vector>
 #include <queue>
 #include <map>
 #include <set>
 using namespace std;
 typedef long long LL;
 const int inf = (<<);
 const int MAXN = ;
 int n,cnt,end,ans;
 int L[MAXN],nn;
 bool pd[MAXN];
 int f[MAXN][];//f[i][j]表示第i个区间往后取2^j个可行的区间后的可以取到第几个区间
 struct node{
     int l,r;
 }a[MAXN],b[MAXN];
 set<int>S;
 inline bool cmpl(node q,node qq){ if(q.l==qq.l) return q.r>qq.r; return q.l<qq.l; }
 inline int getint()
 {
     int w=,q=; char c=getchar();
     while((c<'' || c>'') && c!='-') c=getchar(); if(c=='-') q=,c=getchar();
     while (c>='' && c<='') w=w*+c-'', c=getchar(); return q ? -w : w;
 }

 inline void build(){
     for(int i=;i<=n;i++) b[i]=a[i]; sort(b+,b+n+,cmpl);
     cnt=; b[++cnt]=b[];
     for(int i=;i<=n;i++) {
     while(b[i].r<=b[cnt].r && cnt>) cnt--;
     b[++cnt]=b[i];
     }
     nn=; for(int i=;i<=cnt;i++) L[++nn]=b[i].l;
     L[++nn]=inf; int xia; for(int i=;i<=;i++) f[cnt+][i]=cnt+;
     for(int i=;i<=cnt;i++) { xia=upper_bound(L+,L+nn+,b[i].r)-L; f[i][]=xia; }
     for(int j=;j<=;j++) for(int i=;i<=cnt;i++) if(f[i][j-] && f[f[i][j-]][j-]) f[i][j]=f[f[i][j-]][j-];
 }

 inline int get_ans(int l,int r){//得到区间[l,r]的可选区间的个数
     if(l>r) return ; if(l>end) return ;//!!!
     int daan=,head,cc; head=lower_bound(L+,L+nn+,l)-L;
     if(head>cnt) return ;//!!!
     if(b[head].r>r) return ;
     for(int i=;i>=;i--) {
     if(f[head][i]==cnt+) continue; if(f[head][i]==) continue;
     cc=b[f[head][i]].r; if(cc>r) continue;
     head=f[head][i]; daan+=(<<i);
     }
     return daan;
 }

 inline void work(){
     n=getint(); for(int i=;i<=n;i++) a[i].l=getint(),a[i].r=getint(),L[++cnt]=a[i].l,L[++cnt]=a[i].r;
     sort(L+,L+cnt+); end=nn=unique(L+,L+cnt+)-L-;
     for(int i=;i<=n;i++) a[i].l=lower_bound(L+,L+nn+,a[i].l)-L,a[i].r=lower_bound(L+,L+nn+,a[i].r)-L;
     build(); int now,nowl,nowr; ans=get_ans(,end); printf("%d\n",ans);
     S.insert(-inf); S.insert(inf);  bool flag=false;
     for(int i=;i<=n;i++) {
     now=*S.lower_bound(a[i].l);
     if(now==inf || !pd[now]) {
         if(now<=a[i].r) continue;
         nowl=*--S.lower_bound(a[i].l); nowl++;
         nowr=*S.lower_bound(a[i].r); nowr--;
         now=get_ans(nowl,a[i].l-)+get_ans(a[i].r+,nowr)+;
         if(now<get_ans(nowl,nowr)) continue;
         S.insert(a[i].l); S.insert(a[i].r); if(a[i].l!=a[i].r) pd[a[i].r]=;//!!!
         if(flag) printf(" ");  printf("%d",i); flag=true;
     }
     }
 }

 int main()
 {
     work();
     return ;
 }