[BZOJ 2109/BZOJ 2535][NOI 2010]航空管制(贪心)

题目链接

http://www.lydsy.com/JudgeOnline/problem.php?id=2109

思路

很容易发现其实就是给了一个有向图，在这个有向图上做拓扑序，并限制某些点$i$必须在某个区间$[1,K_i]$中，第一问就是求一个合法的拓扑序，第二问则是求对于每个点$i$，问$i$在拓扑序中最早在什么位置出现。

网上的做法很多，但是有些做法能过BZOJ 2535却过不了2109?(这些做法的贪心策略应该是错误的)，反正各种做法我都看不懂，毕竟我太弱。。。

引用vfk大爷的做法(http://vfleaking.blog.163.com/blog/static/17480763420134219155120/)：可以把整个有向图反向，在反向图中，某些点$i$的拓扑序被限制最早只能在位置$n-K_i+1$出现。那么在放点$i$之前，就可以随意放很多点了，不用考虑后效性。

那么对于第二问，我们要达到的效果就是在加入一个点之前，加入尽量多的点。对于求点$t$最早什么位置出现，我们就可以看成是在反向图的拓扑序里，这个点最晚什么时候出现。做拓扑排序时无视掉点$t$。将所有点按照$n-K_t$升序排序，然后像做一般的拓扑排序一样，假设当前要找拓扑序中位置$i$的点，那么每次加入队列的入度为0的点必须满足$n-K_t+1\leq i$，表示当前的这个点是可以在位置$i$或更晚的地方出现的。由于所有的点已经排序过了，因此每个时刻队列中的点均能满足从队首到队尾，$n-K_t+1$是升序的，越是能早点出现的点，在拓扑序里也是越早出现。这样生成的拓扑序里，那些不能在点$t$前出现的点均不会出现。而在这个拓扑序里第一个没有填入的点就是$t$，这样做显然可以发现，加入点 t <script type="math/tex" id="MathJax-Element-420">t</script>前，所有能加入的点都被加入进来了，这样第二问就能做出来了。
而第一问的处理也比较简单，用第二问的做法，不无视任何一个点就能得到一个合法的拓扑序了。

代码

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>

#define MAXE 21000
#define MAXV 2100

using namespace std;

struct edge
{
    int u,v,next;
}edges[MAXE];

int head[MAXV],nCount=0;
int inDegree[MAXV],tmpdeg[MAXV];

void AddEdge(int U,int V)
{
    edges[++nCount].u=U;
    edges[nCount].v=V;
    edges[nCount].next=head[U];
    head[U]=nCount;
}

struct Plane
{
    int id,limit;
}plane[MAXV];

int n,m;

bool cmp(Plane a,Plane b)
{
    return a.limit<b.limit;
}

int q[MAXE*2],h,t;
int limit[MAXV];

void work(int x) //在拓扑排序时不填入飞机x
{
    memcpy(tmpdeg,inDegree,sizeof(inDegree));
    h=t=0;
    for(int i=1,p=1;i<=n;i++) //在位置i填入一个飞机
    {
        for(;p<=n&&plane[p].limit<i;p++)
            if(!tmpdeg[plane[p].id]&&plane[p].id!=x)
                q[t++]=plane[p].id;
        if(h<t) //队列不为空
        {
            int u=q[h++]; //将u填入第i个位置
            for(int p1=head[u];p1!=-1;p1=edges[p1].next)
            {
                int v=edges[p1].v;
                tmpdeg[v]--;
                if(!tmpdeg[v]&&v!=x&&limit[v]<i)
                    q[t++]=v;
            }
        }
        else return;
    }
}

int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&plane[i].limit);
        plane[i].limit=n-plane[i].limit;
        plane[i].id=i;
        limit[i]=plane[i].limit;
    }
    sort(plane+1,plane+n+1,cmp);
    for(int i=1;i<=m;i++)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        inDegree[a]++;
        AddEdge(b,a);
    }
    work(0);
    for(int i=t-1;i>=0;i--)
        printf("%d%c",q[i],i==0?'\n':' ');
    for(int i=1;i<=n;i++)
    {
        work(i);
        printf("%d%c",n-t,i==n?'\n':' ');
    }
    return 0;
}