[bzoj1178]会议中心

考虑用f[i][j]表示以i为起点（i是时间，所以要离散）选$2^j$条线段（这里不是时间）最小的终点，预处理用倍增的方式来求即可
预处理出这个数组后，就可以很快的求出在$[l,r]$的时间内最多能选多少条线段，类似于树上的倍增（注意判断超过最大时间的情况）
然后从前往后枚举每一段时间，考虑能否加入，当且仅当其分割的那个区间的两部分的答案之和+1=原区间的答案，找区间可以用两个set（分别维护起点和终点）来搞

#include<bits/stdc++.h>
using namespace std;
#define N 200005
struct ji{
    int l,r;
    bool operator < (const ji &k)const{
        return (l<k.l)||(l==k.l)&&(r<k.r);
    }
}a[N],aa[N];
set<int>s1,s2;
int n,m,b[N<<1],f[N<<1][21];
void build(){
    memcpy(aa,a,sizeof(a));
    sort(aa+1,aa+n+1);
    aa[n+1].l=aa[n+1].r=m+1;
    for(int i=1,j=1;i<=m;i++){
        while (i>aa[j].l)j++;
        f[i][0]=aa[j].r;
    }
    for(int i=m-1;i;i--)f[i][0]=min(f[i][0],f[i+1][0]);
    for(int i=0;i<=20;i++)f[m+1][i]=m+1;
    for(int j=1;j<=20;j++)
        for(int i=1;i<=m;i++)f[i][j]=f[min(m,f[i][j-1])+1][j-1];
}
int query(int l,int r){
    int ans=0;
    for(int i=20;i>=0;i--)
        if (f[l][i]<=r){
            l=min(m,f[l][i])+1;
            ans+=(1<<i);
        }
    return ans;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d%d",&a[i].l,&a[i].r);
        b[++b[0]]=a[i].l;
        b[++b[0]]=a[i].r;
    }
    sort(b+1,b+2*n+1);
    m=unique(b+1,b+2*n+1)-b-1;
    for(int i=1;i<=n;i++){
        a[i].l=lower_bound(b+1,b+m+1,a[i].l)-b;
        a[i].r=lower_bound(b+1,b+m+1,a[i].r)-b;
    }
    build();
    printf("%d
",query(1,m));
    bool flag=0;
    s1.insert(1);
    s2.insert(m);
    for(int i=1;i<=n;i++){
        if ((!s1.size())||(*s1.begin()>a[i].l)||(*--s2.end()<a[i].r))continue;
        int x=*--s1.upper_bound(a[i].l),y=*s2.lower_bound(a[i].r);
        if ((x!=*--s1.end())&&(*s1.upper_bound(a[i].l)<=y))continue;
        if (query(x,y)==query(x,a[i].l-1)+1+query(a[i].r+1,y)){
            if (a[i].r==y)s2.erase(y);
            else s1.insert(a[i].r+1);
            if (a[i].l==x)s1.erase(x);
            else s2.insert(a[i].l-1);
            if (flag)printf(" ");
            flag=1;
            printf("%d",i);
        }
    }
}