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;
}

最新文章

  1. mysql在ubuntu下的安装
  2. [NOIP2010初赛]烽火传递+单调队列详细整理
  3. hdu 1556:Color the ball(线段树,区间更新,经典题)
  4. HTML5新特性之WebSocket
  5. 基础学习day08---多态、简单工厂、Object类equals和toString
  6. 2015寒假ACM训练计划
  7. 常见的XMLHttpRequest.status状态码
  8. spring data jpa入门学习
  9. ide远程调试
  10. HDU_1230——火星A+B,加法进制问题
  11. 学习笔记:GLSL Core Tutorial – Vertex Shader(内置变量说明)
  12. Retrofit原理
  13. java导出excel 浏览器直接下载或者或以文件形式导出
  14. ubuntu/debian 多版本gcc共存
  15. 项目thymeleaf
  16. UOJ33 [UR #2] 树上GCD 【点分治】【容斥原理】【分块】
  17. (Alpha)Let&#39;s-技术文档(技术规格说明书)
  18. NOI前训练日记
  19. 【LOJ】#2057. 「TJOI / HEOI2016」游戏
  20. loadrunner日志信息

热门文章

  1. python 3.x 学习笔记9 (面向对象)
  2. OpenGL编程(六)通过三角形绘画出3D模型
  3. servlet中Session的用法
  4. 洛谷P2045 方格取数加强版 最小费用流
  5. Vue过渡与动画
  6. 高并发web系统设计
  7. 安装了python之后Windows的cmd中cd指令无法转换路径怎么办?
  8. jquery.validate动态更改校验规则 【转】
  9. Qt之自定义布局管理器(QBorderLayout)
  10. 将shell脚本运行情况写入Rsyslog日志server