【bzoj4010】[HNOI2015]菜肴制作

题目链接

Description

知名美食家小 A被邀请至ATM 大酒店，为其品评菜肴。

ATM 酒店为小 A 准备了 N 道菜肴，酒店按照为菜肴预估的质量从高到低给予1到N的顺序编号，预估质量最高的菜肴编号为1。由于菜肴之间口味搭配的问题，某些菜肴必须在另一些菜肴之前制作，具体的，一共有 M 条形如“i 号菜肴‘必须’先于 j 号菜肴制作”的限制，我们将这样的限制简写为 < i,j > 。现在，酒店希望能求出一个最优的菜肴的制作顺序，使得小 A能尽量先吃到质量高的菜肴：也就是说，(1)在满足所有限制的前提下，1 号菜肴“尽量”优先制作；(2)在满足所有限制，1号菜肴“尽量”优先制作的前提下，2号菜肴“尽量”优先制作；(3)在满足所有限
制，1号和2号菜肴“尽量”优先的前提下，3号菜肴“尽量”优先制作；(4)在满足所有限制，1 号和 2 号和 3 号菜肴“尽量”优先的前提下，4 号菜肴“尽量”优先制作；(5)以此类推。
例1：共4 道菜肴，两条限制<3,1>、<4,1>，那么制作顺序是 3,4,1,2。例2：共5道菜肴，两条限制<5,2>、 <4,3>，那么制作顺序是 1,5,2,4,3。例1里，首先考虑 1，因为有限制<3,1>和<4,1>，所以只有制作完 3 和 4 后才能制作 1，而根据(3)，3 号又应“尽量”比 4 号优先，所以当前可确定前三道菜的制作顺序是 3,4,1；接下来考虑2，确定最终的制作顺序是 3,4,1,2。例 2里，首先制作 1是不违背限制的；接下来考虑 2 时有<5,2>的限制，所以接下来先制作 5 再制作 2；接下来考虑 3 时有<4,3>的限制，所以接下来先制作 4再制作 3，从而最终的顺序是 ,5,2,4,3。
现在你需要求出这个最优的菜肴制作顺序。无解输出“Impossible!” （不含引号，首字母大写，其余字母小写）

Input

第一行是一个正整数D，表示数据组数。
接下来是D组数据。
对于每组数据：
第一行两个用空格分开的正整数N和M，分别表示菜肴数目和制作顺序限
制的条目数。
接下来M行，每行两个正整数x,y，表示“x号菜肴必须先于y号菜肴制作”的限制。（注意：M条限制中可能存在完全相同的限制）

Output

输出文件仅包含 D 行，每行 N 个整数，表示最优的菜肴制作顺序，或者”Impossible!”表示无解（不含引号）。

Sample Input

3

5 4

5 4

5 3

4 2

3 2

3 3

1 2

2 3

3 1

5 2

5 2

4 3

Sample Output

1 5 3 4 2
Impossible!
1 5 2 4 3

HINT

【样例解释】

第二组数据同时要求菜肴1先于菜肴2制作，菜肴2先于菜肴3制作，菜肴3先于

菜肴1制作，而这是无论如何也不可能满足的，从而导致无解。

100%的数据满足N,M<=100000,D<=3。

题解

等价于反向输出反图中最大字典序的拓扑序。
证明可以看黄学长博客：http://hzwer.com/8642.html。

#include<bits/stdc++.h>
using namespace std;

inline int read(){
    int x = 0, f = 1; char c = getchar();
    while(!isdigit(c)) { if(c == '-') f = -1; c = getchar(); }
    while(isdigit(c)) { x = x * 10 + c - '0'; c = getchar(); }
    return x * f;
}

const int N = 100000 + 10;

int t, m, n, tot, top;
int to[N], nxt[N], hd[N], d[N], ans[N], vis[N];
priority_queue<int, vector<int> > q;

inline void insert(int u, int v){to[++tot] = v; nxt[tot] = hd[u]; hd[u] = tot;}

void init(){
    memset(hd, 0, sizeof(hd));
    memset(d, 0, sizeof(d));
    memset(ans, 0, sizeof(ans));
    tot = top = 0;
    n = read(); m = read();
    for(int i = 1; i <= m; i++){
        int u = read(), v = read();
        insert(v, u); d[u]++;
    }
}

void solve(int x){
    ans[++top] = x; q.pop();
    for(int i = hd[x]; i; i = nxt[i]){
        d[to[i]]--;
        if(!d[to[i]]) q.push(to[i]);
    }
}

void work(){
    init();
    for(int i = 1; i <= n; i++) if(!d[i]) q.push(i);
    while(!q.empty()) solve(q.top());
    if(top != n) puts("Impossible!");
    else for(int i = top; i; i--)
        printf("%d%s", ans[i], i != 1 ? " " : " \n");
}

int main(){
    t = read();
    while(t--){
        work();
    }
}