4444: [Scoi2015]国旗计划

4444: [Scoi2015]国旗计划

Time Limit: 20 Sec  Memory Limit: 256 MB
Submit: 485  Solved: 232

Description

A国正在开展一项伟大的计划——国旗计划。这项计划的内容是边防战士手举国旗环绕边境线奔袭一圈。这

项计划需要多名边防战士以接力的形式共同完成，为此，国土安全局已经挑选了N名优秀的边防战上作为这

项计划的候选人。

A国幅员辽阔，边境线上设有M个边防站，顺时针编号1至M。每名边防战士常驻两个边防站，并且善于

在这两个边防站之间长途奔袭，我们称这两个边防站之间的路程是这个边防战士的奔袭区间。n名边防战士

都是精心挑选的，身体素质极佳，所以每名边防战士的奔袭区间都不会被其他边防战士的奔袭区间所包含。

现在，国十安全局局长希望知道，至少需要多少名边防战士，才能使得他们的奔袭区间覆盖全部的边境线，

从而顺利地完成国旗计划。不仅如此，安全局局长还希望知道更详细的信息：对于每一名边防战士，在他必

须参加国旗计划的前提下，至少需要多少名边防战士才能覆盖全部边境线，从而顺利地完成国旗计划。

 

 

Input

第1行，包含2个正整数N,M，分别表示边防战士数量和边防站数量。

随后n行，每行包含2个正整数。其中第i行包含的两个正整数Ci、Di分别表示i号边防战士常驻的两个边防站编号，

Ci号边防站沿顺时针方向至Di号边防站力他的奔袭区间。数据保证整个边境线都是可被覆盖的。

 

 

Output

输出数据仅1行，需要包含n个正整数。其中，第j个正整数表示j号边防战士必须参加的前提下至少需要

多少名边防战士才能顺利地完成国旗计划

 

Sample Input

4 8
2 5
4 7
6 1
7 3

Sample Output

3 3 4 3

HINT

 

 n≤2×10^5，M< 10^9，1≤Ci，Di≤M

分析

每一个区间都不能被其他区间所包含，也就是如果l_i<l_j那么一定满足r_i<r_j，然后就可以贪心一下。 

对于区间[l,r]，可以在区间找到一个l_i，使得这个l_i的右端点r_i最大，然后继续寻找。 

这样可以按照左端点排序，然后每一个战士要找的下一个战士都是确定的。这样预处理出来。

求解有多少战士的话，可以用倍增，f[i][j]表示从i点走2^j步到达的点，如果f[i][j]所到达的点没有超出m，就可以从i点走2^j步，到达下一个点。 继续寻找。

code

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<set>

using namespace std;

const int MAXN = 400100;
const int INF = 0x7fffffff;

struct Peo{
    int l,r,id;
    bool operator < (const Peo &a) const
    {
        return l < a.l;
    }
}p[MAXN];

int f[MAXN][30],ans[MAXN];//f[i][j]表示从i点走2^j步到达的点 
int tot,n,m;
set< pair<int,int> >s;

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 work(int x)
{
    int d = p[x].id, L = p[x].l+m;
    for (int i=20; i>=0; --i)
        if (f[x][i])
            if (p[f[x][i]].r<L) 
                x = f[x][i],ans[d] += (1<<i);
}

int main()
{
    n = read(),m = read();
    for (int x,y,i=1; i<=n; ++i)
    {
        x = read();y = read();
        if (y<x) y += m;
        p[++tot].l = x; p[tot].r = y; p[tot].id = i;
        p[++tot].l = x+m;p[tot].r = y+m;    
    }
    
    sort(p+1,p+tot+1);
    s.insert(make_pair(p[tot].l,tot));
    for (int i=tot-1; i; --i)
    {
        f[i][0] = (--s.upper_bound(make_pair(p[i].r,INF)))->second; 
        s.insert(make_pair(p[i].l,i));
    }
    
    for (int i=1; i<=20; ++i)
        for (int j=tot; j; --j)
            f[j][i] = f[f[j][i-1]][i-1];
    for (int i=1; i<=tot; ++i)
        if (p[i].id) work(i);
        
    for (int i=1; i<=n; ++i)
         printf("%d ",ans[i]+2);
    return 0;
}