7-4 关键活动 （拓扑序列+关键路径）代码知识点详解

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

比如完成一个专业的所有课程学习和毕业设计可以看成一个本科生要完成的一项工程，各门课程可以看成是子任务。有些课程可以同时开设，比如英语和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

---

拓扑序列和关键路径的难点主要是拓扑序列的顺序特殊性，它与我们往常认识的在邻接矩阵中从0-N的顺序不同。

所以要明白拓扑序列和固定行序和列序的关系。也正因此，我们写了GetIndex()函数找到他们之间的关系。

完整代码：

#include <iostream>
#include <stdlib.h>
#include <stack>
using namespace std;

#define INFINITY 32767

typedef struct 
{
    int *vex;
    int **arc;
    int numVertexes,numEdge;
}MGraph;

MGraph* InitGraph(int numVer,int numEdge)   //初始化图
{
    MGraph *G = (MGraph*)malloc(sizeof(MGraph));
    G->vex = (int*)malloc(sizeof(int)*numVer);
    G->arc = (int**)malloc(sizeof(int*)*numVer);
    for(int i=0;i<numVer;i++)
    {
        G->vex[i] = i+1;
        G->arc[i] = (int*)malloc(sizeof(int)*numVer);
        for(int j=0;j<numVer;j++)
        {
            G->arc[i][j]=INFINITY;
        }
    }
    G->numVertexes=numVer;
    G->numEdge=numEdge;

    return G;
}

int* FindInDgree(MGraph *G)  //计算各个顶点的入度，为拓扑排序做准备
{
    int *Dgree = (int*)malloc(sizeof(int)*G->numVertexes);
    for(int i=0;i<G->numVertexes;i++)
    {
        Dgree[i]=0;
    }

    for(int i=0;i<G->numVertexes;i++)
    {
        for(int j=0;j<G->numVertexes;j++)
        {
            if(G->arc[i][j]!=INFINITY)
            {
                Dgree[j]++;
            }
        }
    }
    return Dgree;
}

int* ToplogicalSort(MGraph *G)  //拓扑排序
{
    int index=0;
    stack<int> S;
    int *TopSort = (int*)malloc(sizeof(int)*G->numVertexes);
    int *inDgree = FindInDgree(G);
    for(int i=0;i<G->numVertexes;i++)
    {
        if(inDgree[i]==0)
            S.push(i);
    }

    while(S.empty()==0)
    {
        int vex = S.top();
        TopSort[index++]=vex;
        S.pop();
        for(int i=0;i<G->numVertexes;i++)
        {
            if(G->arc[vex][i]!=INFINITY)
            {
                inDgree[i]--;
                if(inDgree[i]==0)
                    S.push(i);
            } 
        }
    }

    if(index<G->numVertexes)
    {
        cout<<0;
        exit(0);
    }

   return TopSort;
}

int GetIndex(MGraph *G,int* top,int i)
{
    int j;
    for(j=0;j<G->numVertexes;j++)
    {
        if(top[j]==i)
            break;
    }

    return j;
}

void CriticalPath(MGraph *G)  //关键路径中，求最早开始时间
{
    int *TopSort = ToplogicalSort(G);
    int *early = (int*)malloc(sizeof(int)*G->numVertexes);
    int *late = (int*)malloc(sizeof(int)*G->numVertexes);
    for(int i=0;i<G->numVertexes;i++)
    {
        early[i] = 0;
        late[i] = 0;
    }

    //计算最早开始时间
    for(int i=0;i<G->numVertexes;i++)
    {
        int max=0;
        for(int j=0;j<G->numVertexes;j++)
        {
            if(G->arc[j][TopSort[i]]!=INFINITY)
            {
                int index = GetIndex(G,TopSort,j);
                if(early[index] + G->arc[j][TopSort[i]] > max)
                    max = early[index] + G->arc[j][TopSort[i]];
            }
        }
        early[i] = max;
    }
    cout<<early[(G->numVertexes)-1]<<endl;
    //计算最晚开始时间
    late[(G->numVertexes)-1] = early[(G->numVertexes)-1];
    for(int i = (G->numVertexes) - 2;i>=0;i--)
    {
        int min = INFINITY;
        for(int j=0;j<G->numVertexes;j++)
        {
            if(G->arc[TopSort[i]][j]!=INFINITY)
            {
                int index = GetIndex(G,TopSort,j);
                if(late[index] - G->arc[TopSort[i]][j] < min)
                    min = late[index] - G->arc[TopSort[i]][j];
            }

        }

        late[i] = min;
    }

    //开始求关键路径,遍历整个图的边
    for(int i=0;i<G->numVertexes;i++)
    {
        for(int j=0;j<G->numVertexes;j++)
        {
            if(G->arc[i][j]!=INFINITY)
            {
                int start = GetIndex(G,TopSort,i);
                int end = GetIndex(G,TopSort,j);
                if((late[end] - G->arc[i][j]) - early[start] == 0) //满足条件的边即是关键路径
                {
                    //cout<<i+1<<"->"<<j+1<<endl;
                    cout<<G->vex[i]<<"->"<<G->vex[j]<<endl;


                }
            }
        }
    }

}


int main()
{
    int N,M;
    MGraph *G;

    cin>>N>>M;
    G = InitGraph(N,M);

    for(int i=0;i<M;i++)
    {
        int x,y,time;
        cin>>x>>y>>time;
        G->arc[x-1][y-1] = time;
    }

    //ToplogicalSort(G);

    CriticalPath(G);

    return 0;
}