【BZOJ-1178】CONVENTION会议中心倍增 + set （神思路好题！）

1178: [Apio2009]CONVENTION会议中心

Time Limit: 15 Sec Memory Limit: 162 MB
Submit: 812 Solved: 323
[Submit][Status][Discuss]

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

HINT

修复数据bug,并新加数据一组By NanoApe 2016.5.11

修复后数据:JudgeOnline/upload/201605/dd.rar

Source

Solution

思路非常神

首先我们如果只要求数目最多，显然可以贪心，即线段覆盖

字典序最小就很麻烦了

方法还是贪心，按照读入的顺序，判断，如果使当前满足，能够达到最优，那么就选当前

问题是如何加速判断的过程

$f[i][j]$表示，以$i$时间开始，选$2^{j}$条线段后，右端点最近在哪

这个信息显然是可以合并的，那么倍增解决

判断的时候，用数据结构维护一下就可以，比如用set查询前驱后继

Code

#include<iostream>

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<cmath>

#include<set>

using namespace std;

int read()

{

    int x=,f=; char ch=getchar();

    while (ch<'' || ch>'') {if (ch=='-') f=-; ch=getchar();}

    while (ch>='' && ch<='') {x=x*+ch-'';  ch=getchar();}

    return x*f;

}

#define MAXN 200010

#define INF 0x3f3f3f3f

int N;

struct CompanyNode{int s,t,id;}c[MAXN],tmp[MAXN],b[MAXN];

bool cmp1(CompanyNode A,CompanyNode B) {return A.s==B.s? A.t>B.t : A.s<B.s;}

bool cmp2(CompanyNode A,CompanyNode B) {return A.id<B.id;}

struct SetNode

{

    int x,f;

    SetNode (int x,int f) : x(x),f(f) {}

    bool operator < (const SetNode & A) const

        {return x<A.x;}

};

set<SetNode>st;

set<SetNode>::iterator lst,rst;

int f[MAXN<<][];

int Cal(int l,int r)

{

    int re=;

    for (int i=; i>=; i--)

        if (f[l][i]<=r) l=f[l][i]+,re+=<<i;  //一定要+1

    return re;

}

int ls[MAXN<<],tp,top;

int main()

{

    //freopen("data.in","r",stdin); freopen("data.out","w",stdout);

    N=read();

    for (int i=; i<=N; i++) ls[++tp]=c[i].s=read(),ls[++tp]=c[i].t=read(),c[i].id=i;

    sort(ls+,ls+tp+);

    top=unique(ls+,ls+tp+)-ls-;

    for (int i=; i<=N; i++)

        c[i].s=lower_bound(ls+,ls+top+,c[i].s)-ls,c[i].t=lower_bound(ls+,ls+top+,c[i].t)-ls;

    sort(c+,c+N+,cmp1);

    int last=INF; tp=;

    for (int i=N; i>=; i--) if (c[i].t<last) tmp[++tp]=c[i],last=c[i].t;

    reverse(tmp+,tmp+tp+);

//    for (int i=1; i<=tp; i++) printf("%d %d %d\n",tmp[i].s,tmp[i].t,tmp[i].id);

    memset(f,0x3f,sizeof(f));

//    printf("%d %d\n",f[1],INF);

    for (int i=top,j=tp; i>=; i--)

        {

            f[i][]=f[i+][];

            if (tmp[j].s==i) f[i][]=min(f[i][],tmp[j].t);

            for (int k=; k<=; k++)

                if (f[i][k-]!=INF)

                    f[i][k]=f[f[i][k-]+][k-];

            while (tmp[j].s==i) j--;

        }

//    for (int i=1; i<=N; i++) printf("%d %d %d\n",c[i].id,c[i].s,c[i].t);

    printf("%d\n",Cal(,top));

    sort(c+,c+N+,cmp2);

    st.insert(SetNode(,-)); st.insert(SetNode(top+,));

    int cnt=;

    for (int i=; i<=N; i++)

        {

            SetNode ls=SetNode(c[i].s,),rs=SetNode(c[i].t,-);

            lst=st.lower_bound(ls); rst=st.upper_bound(rs);

            if (lst!=rst || (*rst).f==-) continue;

            lst--;

            int L=(*lst).x+,R=(*rst).x-;

            if (Cal(L,c[i].s-)+Cal(c[i].t+,R)+==Cal(L,R))

                {

                    st.insert(ls),st.insert(rs);

                    cnt++; if (cnt==) printf("%d",c[i].id); else printf(" %d",c[i].id);

                }

        }

    puts("");

    return ;

}