20150227题解

关键子工程（project）
题目：
在大型工程的施工前，我们把整个工程划分为若干个子工程，并把这些子工程编号为1、2、……、N；这样划分之后，子工程之间就会有一些依赖关系，即一些子工程必须在某些子工程完成之后才能施工。由于子工程之间有相互依赖关系，因此有两个任务需要我们去完成：首先，我们需要计算整个工程最少的完成时间；同时，由于一些不可预测的客观因素会使某些子工程延期，因此我们必须知道哪些子工程的延期会影响整个工程的延期，我们把有这种特征的子工程称为关键子工程，因此第二个任务就是找出所有的关键子工程，以便集中精力管理好这些子工程，尽量避免这些子工程延期，达到用最快的速度完成整个工程。为了便于编程，现在我们假设：
（1）根据预算，每一个子工程都有一个完成时间。
（2）子工程之间的依赖关系是：部分子工程必须在一些子工程完成之后才开工。
（3）只要满足子工程间的依赖关系，在任何时刻可以有任何多个子工程同时在施工，也既同时施工的子工程个数不受限制。
（4）整个工程的完成是指：所有子工程的完成。
例如，有五个子工程的工程规划表：
序号 完成时间 子工程1 子工程2 子工程3 子工程4 子工程5
子工程1 5 0 0 0 0
子工程2 4 0 0 0 0
子工程3 12 0 0 0 0
子工程4 7 1 1 0 0
子工程5 2 1 1 1 1
其中，表格中第I+1行J+2列的值如为0表示“子工程I”可以在“子工程J”没完成前施工，为1表示“子工程I”必须在“子工程J”完成后才能施工。上述工程最快完成时间为14天，其中子工程1、3、4、5为关键子工程。

又例如，有五个子工程的工程规划表：
序号 完成时间 子工程1 子工程2 子工程3 子工程4 子工程5
子工程1 5 0 1 0 0
子工程2 4 0 0 0 0
子工程3 12 0 0 1 0
子工程4 7 1 1 0 0
子工程5 2 1 1 1 1
上述的子工程划分不合理，因为无法安排子工程1，3，4的施工。

输入数据：
第1行为N，N是子工程的总个数，N≤200。
第2行为N个正整数，分别代表子工程1、2、……、N的完成时间。
第3行到N+2行，每行有N-1个0或1。其中的第I+2行的这些0，1，分别表示“子工程I”与子工程1、2、…、I-1、I+1、…N的依赖关系，（I=1、2、……、N）。每行数据之间均用一个空格分开。

输出数据：
如子工程划分不合理，则输出-1；
如子工程划分合理，则用两行输出：第1行为整个工程最少的完成时间。第2行为按由小到大顺序输出所有关键子工程的编号。
样例：
输入文件名：project. in
5
5 4 12 7 2
0 0 0 0
0 0 0 0
0 0 0 0
1 1 0 0
1 1 1 1
输出文件名：project. out
14
1 3 4 5

题解：
本题是一个经典的dp题，题目大意是完成三个任务：
1、给定一个拓扑图，判断是否有环；
2、给定一个无环拓扑图，求总最早完成时间；
3、求无环拓扑图中是否有子工程可以延迟。
N<=200的规模搜索至少是很勉强，但更快的是记忆化搜索。
对第一个任务，可以简单地拓扑排序，如果有环输出“-1”并退出，这样至少可得10分；
第二个任务则是记忆化搜索，（本质上还是dp）若第x个子工程依赖别的子工程i，则递归计算i；若i已被算出，直接返回答案；x的最短完成时间即为它所依赖的所有子工程中完成时间最长的加上x本身的完成时间；
第三个任务依然是记忆化搜索。考虑第x个子工程，它能被延迟当且仅当它的最晚完成时间小于等于全部{【依赖x的所有子工程】在不影响总最早完成时间的情况下的最晚开始时间}的最小值。（dp方向与第二个任务相反）
初始化最晚完成时间为总最早完成时间；（要从后往前推）
若有工程i依赖x，则递归计算i；若i已被算出，直接返回答案；x的最晚完成时间即为所有依赖x的子工程最晚开始时间的最小值。（该值等于min{所有依赖x的子工程i最晚完成时间减去i工程所需要的时间}）
若某工程的最早完成时间与最晚完成时间相等，则不能被延迟，应输出。
CODE

#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
int n, map[201][201], t[201], sumt;          
int imp[201], top;      
int ealy[201], lst[201];
queue <int> Q;
inline void init() {    
    scanf("%d", &n);
    for(int i = 1; i <= n; ++i) scanf("%d", &t[i]);
    for(int i = 1; i <= n; ++i) 
        for(int j = 1; j <= n; ++j)
            if(i == j) map[i][j] = 0;
            else {
                scanf("%d", &map[i][j]);
                if(map[i][j]) ++map[i][0];
            }
}
inline bool topo() {    
    int din[201], sigma = n;