码迷,mamicode.com
首页 > Windows程序 > 详细

BZOJ1178 [Apio2009]CONVENTION会议中心 贪心 set

时间:2017-08-18 20:04:16      阅读:169      评论:0      收藏:0      [点我收藏+]

标签:scan   处理   线段   区间   print   cst   namespace   注意   std   

欢迎访问~原文出处——博客园-zhouzhendong

去博客园看该题解


题目传送门 - BZOJ1178


题意概括

  一堆线段,现在取出最多条数,并输出字典序最小的方案。


 

题解

  这题好坑。

  首先,注意一点,最后不能有多余的空格。

  第一问就是基础的线段覆盖。

  关键在于第二问。

  我们要准备一个函数—— Get_Ans(L,R),用来求解L~R这个区间内,最多可以取多少线段。

  这个可以O(log n)解决。前提是倍增预处理。

  然后就是关键部分。

  我们按照字典序从小到大,能选择就选择。

  怎么判断是否可以选择?只需要在之前选择的线段空隙中可以放这个线段,并且满足Get_Ans(这个线段到它左边的第一条线段之前)+Get_Ans(它到它右边的第一个线段之前)+1 = Get_Ans(他左右两端的),那么就可以放进去。


代码

 

#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <set>
using namespace std;
const int N=200000+5,Inf=1e9;
struct seg{
    int L,R;
    bool operator < (const seg x) const {
        if (L==x.L)
            return R>x.R;
        return L<x.L;
    }
}s[N],s2[N];
int n,rn,Hash[N*2],hs,p[N][20],L[N];
bool alive[N],vis[N*2];
set <int> S;
void HASH(){
    sort(Hash+1,Hash+hs+1);
    int hs_=1;
    for (int i=2;i<=hs;i++)
        if (Hash[i]!=Hash[i-1])
            Hash[++hs_]=Hash[i];
    hs=hs_;
}
void Delete(){
    for (int i=1;i<=n;i++)
        s2[i]=s[i];
    sort(s+1,s+n+1);
    memset(alive,true,sizeof alive);
    int MinR=Inf,n_=0;
    for (int i=n;i>=1;i--)
        if (MinR>s[i].R)
            MinR=s[i].R;
        else
            alive[i]=0;
    for (int i=1;i<=n;i++)
        if (alive[i])
            s[++n_]=s[i];
    n=n_;
}
void Get_Pos(){
    memset(p,0,sizeof p);
    for (int i=1;i<=n;i++)
        L[i]=s[i].L;
    L[n+1]=Inf;
    for (int i=1;i<=n;i++){
        p[i][0]=upper_bound(L+1,L+n+2,s[i].R)-L;
        if (p[i][0]==n+1)
            p[i][0]=0;
    }
    for (int i=1;i<=18;i++)
        for (int j=1;j<=n;j++)
            p[j][i]=p[p[j][i-1]][i-1];
}
int Get_Ans(int Le,int Ri){
    if (Le>Ri||Le>hs)
        return 0;
    int ans=1,pos=lower_bound(L+1,L+n+2,Le)-L;
    if (pos>n||Ri<s[pos].R)
        return 0;
    for (int i=18;i>=0;i--){
        if (!p[pos][i]||s[p[pos][i]].R>Ri)
            continue;
        pos=p[pos][i];
        ans+=1<<i;
    }
    return ans;
}
int main(){
    scanf("%d",&n),rn=n;
    for (int i=1;i<=n;i++){
        scanf("%d%d",&s[i].L,&s[i].R);
        Hash[++hs]=s[i].L;
        Hash[++hs]=s[i].R;
    }
    HASH();
    for (int i=1;i<=n;i++){
        s[i].L=lower_bound(Hash+1,Hash+hs+1,s[i].L)-Hash;
        s[i].R=lower_bound(Hash+1,Hash+hs+1,s[i].R)-Hash;
    }
    Delete();
    Get_Pos();
    printf("%d\n",Get_Ans(1,hs));
    S.clear();
    S.insert(Inf),S.insert(-Inf);
    bool empty_block=0;
    for (int i=1;i<=rn;i++){
        int now=*S.lower_bound(s2[i].L);
        if (now==Inf||!vis[now]){
            if (now<=s2[i].R)
                continue;
            int left=*--S.lower_bound(s2[i].L);
            int right=*S.lower_bound(s2[i].R);
            int val=Get_Ans(left+1,s2[i].L-1)+1+Get_Ans(s2[i].R+1,right-1);
            if (val<Get_Ans(left+1,right-1))
                continue;
            S.insert(s2[i].L);
            S.insert(s2[i].R);
            if (s2[i].L!=s2[i].R)
                vis[s2[i].R]=1;
            if (empty_block)
                putchar(‘ ‘);
            printf("%d",i);
            empty_block=1;
        }
    }
    return 0;
}

 

  

 

BZOJ1178 [Apio2009]CONVENTION会议中心 贪心 set

标签:scan   处理   线段   区间   print   cst   namespace   注意   std   

原文地址:http://www.cnblogs.com/zhouzhendong/p/7391164.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!