7-11 关键活动

假定一个工程项目由一组子任务构成，子任务之间有的可以并行执行，有的必须在完成了其它一些子任务后才能执行。“任务调度”包括一组子任务、以及每个子任务可以执行所依赖的子任务集。

比如完成一个专业的所有课程学习和毕业设计可以看成一个本科生要完成的一项工程，各门课程可以看成是子任务。有些课程可以同时开设，比如英语和C程序设计，它们没有必须先修哪门的约束；有些课程则不可以同时开设，因为它们有先后的依赖关系，比如C程序设计和数据结构两门课，必须先学习前者。

但是需要注意的是，对一组子任务，并不是任意的任务调度都是一个可行的方案。比如方案中存在“子任务A依赖于子任务B，子任务B依赖于子任务C，子任务C又依赖于子任务A”，那么这三个任务哪个都不能先执行，这就是一个不可行的方案。

任务调度问题中，如果还给出了完成每个子任务需要的时间，则我们可以算出完成整个工程需要的最短时间。在这些子任务中，有些任务即使推迟几天完成，也不会影响全局的工期；但是有些任务必须准时完成，否则整个项目的工期就要因此延误，这种任务就叫“关键活动”。

请编写程序判定一个给定的工程项目的任务调度是否可行；如果该调度方案可行，则计算完成整个工程项目需要的最短时间，并输出所有的关键活动。

输入格式:

输入第1行给出两个正整数$N$($\le 100$)和$M$，其中$N$是任务交接点（即衔接相互依赖的两个子任务的节点，例如：若任务2要在任务1完成后才开始，则两任务之间必有一个交接点）的数量。交接点按1~$N$编号，$M$是子任务的数量，依次编号为1~$M$。随后$M$行，每行给出了3个正整数，分别是该任务开始和完成涉及的交接点编号以及该任务所需的时间，整数间用空格分隔。

输出格式:

如果任务调度不可行，则输出0；否则第1行输出完成整个工程项目需要的时间，第2行开始输出所有关键活动，每个关键活动占一行，按格式“V->W”输出，其中V和W为该任务开始和完成涉及的交接点编号。关键活动输出的顺序规则是：任务开始的交接点编号小者优先，起点编号相同时，与输入时任务的顺序相反。

输入样例:

7 8
1 2 4
1 3 3
2 4 5
3 4 3
4 5 1
4 6 6
5 7 5
6 7 2

输出样例:

17
1->2
2->4
4->6
6->7

参考代码：

#include <cstdlib>
#include <queue>
#include <iostream>
#include <stack>
using namespace std;
/*
关键活动——拓扑排序+最短路径
邻接矩阵
*/
#define MaxNum 100
#define INFINITY 65533
struct GNode{
    int Nv;
    int Ne;
    int G[MaxNum][MaxNum];
};
typedef struct GNode *Graph;
/*
记录Top排序后的结果，其中
change[1]表示任务开始顶点
change[G->Nv]表示任务结束顶点
*/
int change[MaxNum] = {0,};

void INDEGREE(Graph G, int Indegree[])
{
    for (int V = 1; V <= G->Nv; V++){
        for (int i = 1; i <= G->Nv; i++){
            if (G->G[i][V] < INFINITY)
                Indegree[V]++;
        }
    }
}

void KeyActivity(int Ecost[], int Lcost[], Graph G)
{
    //以下求最早发生时间
    for (int V = 1; V <= G->Nv; V++){
        for (int i = 1; i <= G->Nv; i++){
            if (G->G[change[V]][i] < INFINITY)
                if (Ecost[change[V]] + G->G[change[V]][i] > Ecost[i]){
                    Ecost[i] = Ecost[change[V]] + G->G[change[V]][i];
                }
        }
    }
    //以下求最晚发生时间
    for (int i = 1; i <= G->Nv; i++)
        Lcost[i] = Ecost[change[G->Nv]];
    printf("%d\n", Ecost[change[G->Nv]]);
    for (int V = G->Nv; V >= 1; V--){
        for (int i = G->Nv; i >= 1; i--){
            if (G->G[change[V]][i] < INFINITY)
                if (Lcost[i] - G->G[change[V]][i] < Lcost[change[V]]){
                    Lcost[change[V]] = Lcost[i] - G->G[change[V]][i];
                }
        }
    }
    for (int i = 1; i <= G->Nv; i++){
        for (int j = 1; j <= G->Nv; j++)
            if (G->G[i][j] < INFINITY){
                int e = Ecost[i];
                int l = Lcost[j] - G->G[i][j];
                if ( e == l)
                    printf("%d->%d\n", i, j);
            }
    }
}

int TopSort(Graph G, int Indegree[])
//用来判断任务的调度是否可行
{
    int cnt = 0;
    int sum = 0;
    queue<int> q;
    for (int i = 1; i <= G->Nv; i++){
    //计算每个顶点的入度
        if (Indegree[i] == 0){
            q.push(i);
            cnt++;
        }
    }
    //int max[MaxNum] = {0,};
    int V;
    while (!q.empty()){
        V = q.front();
        q.pop();
        sum++;
        change[sum] = V;//将拓扑排序结果存入change中
        for (int i = 1; i <= G->Nv; i++)
            if (G->G[V][i] < INFINITY){//对于V的每个邻接点
                Indegree[i]--;
                if (Indegree[i] == 0){
                    q.push(i);
                    cnt++;
                }
            }
    }
    if (cnt < G->Nv)
        return 0;
    else 
        return 1;
}


int main(int argc, char const *argv[])
{
    int N, M, time;
    scanf("%d %d\n", &N, &M);
    int a, b;
    Graph G = (Graph)malloc(sizeof(struct GNode));
    G->Nv = N;
    G->Ne = M;
    for (int i = 1; i <= N; i++)
        for (int j = 1; j <= N; j++)
            G->G[i][j] = INFINITY;
    for (int i = 0; i < M; i++){
        scanf("%d %d %d", &a, &b, &time);
        //有向图
        G->G[a][b] = time; 
    }
    int Indegree[MaxNum] = {0,};
    int Ecost[MaxNum] = {0,};
    int Lcost[MaxNum] = {0,};
    INDEGREE(G, Indegree);
    int check = TopSort(G, Indegree);
    if (check == 0)
        printf("0\n");
    else
        KeyActivity(Ecost, Lcost, G);
    system("pause");
    return 0;
}