LibreOJ #2007.「SCOI2015」国旗计划 贪心

题意

有一个长度为m的环和n个区间，对每个区间分别求，若必须要选这个区间，则至少需要选多少个区间才能覆盖完环上的每一条边。
$n\le200000,m\le10^9$

分析

先离散化，然后破环成链，然后把链倍长。
对于每个点$i$，求出$F_i$表示用左端点不超过$i$的区间能到达的最右边的点是多少，这么可以$O(n)$求出。
然后对于每一个区间，从该区间的右端点$r_i$开始跳，直到跳到某个不小于$r_i+m$的点，步数即为答案。
如果我们把$F_i$看成是$i$的父亲的话，这实际上是一棵树。
注意到每个点的答案最多为最优值+1，那么设$T=1的答案-1$，然后从根开始dfs整棵树，记录一下深搜栈，对每个点从$T$开始枚举答案，然后就可以在栈中$O(1)$判断了。
时间复杂度$O(n)$。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>

const int N=800005;

int n,m,cnt,last[N],fa[N],ans[N],lef[N],rig[N],stack[N],top,T,w[N];
struct edge{int to,next;}e[N];
std::vector<int> vec[N];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void addedge(int u,int v)
{
    e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
}

int query(int x)
{
    int ans=0;
    while (x<m+1) x=fa[x],ans++;
    return ans;
}

void dfs(int x)
{
    stack[++top]=x;
    if (x<=m)
    {
        for (int i=T;;i++)
            if (stack[top-i]>=x+m) {ans[x]=i;break;}
    }
    for (int i=last[x];i;i=e[i].next) dfs(e[i].to);
    top--;
}

int main()
{
    n=read();m=read();m=0;
    for (int i=1;i<=n;i++)
    {
        lef[i]=read();w[++m]=lef[i];
        rig[i]=read();w[++m]=rig[i];
    }
    std::sort(w+1,w+m+1);
    m=std::unique(w+1,w+m+1)-w-1;
    for (int i=1;i<=n;i++)
    {
        lef[i]=std::lower_bound(w+1,w+m+1,lef[i])-w;
        rig[i]=std::lower_bound(w+1,w+m+1,rig[i])-w;
        if (lef[i]<rig[i])
        {
            vec[lef[i]].push_back(rig[i]);
            vec[lef[i]+m].push_back(rig[i]+m);
        }
        else
        {
            vec[lef[i]].push_back(rig[i]+m);
            vec[1].push_back(rig[i]);
            vec[lef[i]+m].push_back(m*2);
        }
    }
    int mx=0;
    for (int i=1;i<=m*2;i++)
    {
        for (int j=0;j<vec[i].size();j++) mx=std::max(mx,vec[i][j]);
        if (mx!=i) fa[i]=mx,addedge(mx,i);
    }
    T=query(1)-1;
    dfs(m*2);
    for (int i=1;i<=n;i++) printf("%d ",ans[rig[i]]);
    return 0;
}