08-图9 关键活动 (30分)

题目描述

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

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

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

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

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

输入格式:

输入第1行给出两个正整数N(≤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

解题思路

这道题是上一道题的升级版，上一题只需求出关键路径长度，而这一题还需要打印出所有的关键路径。注意，关键路径可能不只一条，输出时需要打印所有可构成关键路径的边。由于需要输出具体的关键路径，这道题在上道题的最早完成时间数组的基础上还需要添加一个最晚完成时间数组，为了构造这个数组，因此还需要一个反向的邻接表。
那么稍微总结一下这道题主要用到的数据结构：正向邻接表graph1，反向邻接表graph2，最早完成时间数组eTime[]，最晚完成时间数组lTime[]。最后当通过拓扑排序形成最终的eTime[]和lTime[]时，只需要检测是否满足关键路径条件lTime[e->target] - e->weight == eTime[i]然后打印路径即可。

代码

#include <stdio.h>
#include <stdlib.h>

#define MAXSIZE 101

struct ENode {
    int target;
    int weight;
    struct ENode *next;
};
typedef struct ENode *Edge;

typedef struct {
    Edge firstEdge;
} List[MAXSIZE];

struct Graph {
    int vertexCount;
    int edgeCount;
    List header;
};
typedef struct Graph *LGraph;

LGraph graph1, graph2;  //graph1是正向邻接表，graph2是反向邻接表

void createGraph();
void topSort();

int main() {
    createGraph();
    topSort();
    return 0;
}

void createGraph() {
    int N, M;
    scanf("%d %d", &N, &M);

    graph1 = (LGraph) malloc(sizeof(struct Graph));
    graph2 = (LGraph) malloc(sizeof(struct Graph));
    graph1->vertexCount = N; graph2->vertexCount = N;
    graph1->edgeCount = M; graph2->edgeCount = M;
   

    for (int i = 1; i <= N; i++) {
        graph1->header[i].firstEdge = NULL;
        graph2->header[i].firstEdge = NULL;
    }
    
    for (int i = 0; i < M; i++) {
        int x, y, weight;
        scanf("%d %d %d", &x, &y, &weight);

        Edge newEdge1 = (Edge) malloc(sizeof(struct ENode));
        newEdge1->target = y;
        newEdge1->weight = weight;
        newEdge1->next = graph1->header[x].firstEdge;
        graph1->header[x].firstEdge = newEdge1;

        Edge newEdge2 = (Edge) malloc(sizeof(struct ENode));
        newEdge2->target = x;
        newEdge2->weight = weight;
        newEdge2->next = graph2->header[y].firstEdge;
        graph2->header[y].firstEdge = newEdge2;
    }
}

void topSort() {
    int counter = 0;                //判断是否有回路
    int indegree[MAXSIZE] = {0};    //入度
    int outdegree[MAXSIZE] = {0};   //出度
    int eTime[MAXSIZE] = {0};       //最早完成时间
    int lTime[MAXSIZE] = {0};       //最晚完成时间
    int queue[MAXSIZE], front = 0, rear = 0;    //队列

    for (int i = 1; i <= graph1->vertexCount; i++)
        for (Edge e = graph1->header[i].firstEdge; e; e = e->next)
            indegree[e->target]++;  //初始化入度

    for (int i = 1; i <= graph1->vertexCount; i++)
        if (indegree[i] == 0)
            queue[rear++] = i;      //初始化队列

    while (front != rear) {
        int vertex = queue[front++];
        counter++;
        for (Edge e = graph1->header[vertex].firstEdge; e; e = e->next) {
            if (eTime[vertex] + e->weight > eTime[e->target])
                eTime[e->target] = eTime[vertex] + e->weight;
            if (--indegree[e->target] == 0) queue[rear++] = e->target;
        } 
    }
    //如果没有环路
    if (counter == graph1->vertexCount) {
        int max = 0;
        for (int i = 1; i <= counter; i++)
            if (eTime[i] > max)
                max = eTime[i];
        
        front = rear = 0;               //重置队列
        for (int i = 1; i <= graph2->vertexCount; i++)
            lTime[i] = max;             //初始化最晚完成时间

        for (int i = 1; i <= graph2->vertexCount; i++)
            for (Edge e = graph2->header[i].firstEdge; e; e = e->next)
                outdegree[e->target]++; //初始化出度
        
        for (int i = 1; i <= graph2->vertexCount; i++)
            if (outdegree[i] == 0)
                queue[rear++] = i;      //初始化队列

        while (front != rear) {
            int vertex = queue[front++];
            for (Edge e = graph2->header[vertex].firstEdge; e; e = e->next) {
                if (lTime[e->target] > lTime[vertex] - e->weight)
                    lTime[e->target] = lTime[vertex] - e->weight;
                if (--outdegree[e->target] == 0) queue[rear++] = e->target;
            }
        }

        printf("%d
", max);
        for (int i = 1; i <= counter; i++)
            for (Edge e = graph1->header[i].firstEdge; e; e = e->next)
                if (lTime[e->target] - e->weight == eTime[i])
                    printf("%d->%d
", i, e->target);
    } else {
        printf("0
");
    }
}