[APIO2009]会议中心（贪心）

P3626 [APIO2009]会议中心

题目描述

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)。 你的任务是帮助销售主管确定应该将会堂租借给哪些公司。

输入输出格式

输入格式：

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

输出格式：

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

输入输出样例

输入样例#1： 复制

4

4 9

9 11

13 19

10 17

输出样例#1： 复制

2

1 3

说明

对于 50%的输入，N≤3000。

在所有输入中，N≤200000。

---

题解

先让我吐槽一波。

今天教练给我们考倍增NOIP（Professional）模拟的考试题目。

第一题就是这一题，就第二题LCA模板？？？？

3，4都是NOI原题。肝了一下午。真棒。

接下来讲题解了。

考场上我想的是把一个线段的左端点排序。我不会告诉你我是从 \(1\) 跳到 \(n\)的

从 \(n\) 遍历到 \(1\) 然后二分最靠近它右端点的另一个线段的左端点。然后如果能跳到 \(i-1\) 的个数小于 \(i\)，那么就把 \(i-1\) 的信息都更新为 \(i\) 的，不然则把 \(i\) 的总值同等于 \(i-1\)。

这样的话我们就能在\(O(nlogn)\)的时间复杂度内搞出第一问了。

但是第二问怎么搞啊。

很明显的一个贪心：当我可以从多个线段跳过来的时候，我肯定是选择那个字典序最小的，那么，我就只需要在发现 \(i\) 的总值和 \(i-1\) 的总值相等时加一个暴力维护一下就好了，不相等就暴力更新。

时间复杂度是玄学的，可以被卡。但是X谷上还是跑的挺快的。

---

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=200001;
const int fuck=99999999;
struct node{
    int l,r,id;
}ch[N];
int n,m,ans,cnt;
int sum[N],pre[N],now=99999999,line[N],minn[N],nn;
int read()
{
    int x=0,w=1;char ch=getchar();
    while(ch>'9'||ch<'0'){if(ch=='-')w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*w;
}

bool cmp(node a,node b){
    return a.l<b.l;
}

bool cmp2(node a,node b){
    return a.id<b.id;
}

int find (int x){
    int l=1,r=n,ans=0;
    while(r>l){
        int mid=(l+r)>>1;
        if(ch[mid].l>x)r=mid,ans=mid;
        else l=mid+1;
    }
    return r;
}

bool judge(int x,int y){
    int a=fuck,b=fuck;
    while(x!=y){
        a=min(ch[x].id,a);b=min(ch[y].id,b);
        x=pre[x];y=pre[y];
    }
    return a<b;
}

void solve()
{
    for(int i=n;i>=1;i--){
        int next=find(ch[i].r);
        sum[i]=1;minn[i]=i;pre[i]=fuck;
        if(ch[next].l>ch[i].r)sum[i]=sum[next]+1,pre[i]=minn[next];
        if(i!=n){
            if(sum[i+1]>sum[i])
            sum[i]=sum[i+1],minn[i]=minn[i+1],pre[i]=pre[i+1];
            if(sum[i+1]==sum[i]&&judge(minn[i+1],minn[i]))
            sum[i]=sum[i+1],minn[i]=minn[i+1],pre[i]=pre[i+1];
        }
    ans=max(ans,sum[i]);
}
    cout<<ans<<endl;
    int x=minn[1];
    while(x!=fuck){
        line[++cnt]=ch[x].id;
        x=pre[x];
    }
    sort(line+1,line+cnt+1);
    for(int i=1;i<=cnt;i++)
    cout<<line[i]<<' ';
}

int main()
{
//	freopen("interval.in","r",stdin);
//	freopen("interval.out","w",stdout);
    memset(minn,0x7f,sizeof(minn));
    n=read();
    for(int i=1;i<=n;i++)pre[i]=fuck;
    for(int i=1;i<=n;i++){
        ch[i].l=read();ch[i].r=read();ch[i].id=i;
    }
    sort(ch+1,ch+n+1,cmp);
    solve();
    return 0;
}