7-15 关键活动 (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

源码：

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

#define INFINITY 99999999
#define MaxVertexNum 1010 /* maximum number of vertices */
typedef int Vertex;      /* vertices are numbered from 0 to MaxVertexNum-1 */
typedef int WeightType;

typedef struct GNode *PtrToGNode;
struct GNode{
    int Nv;
    int Ne;
    WeightType G[MaxVertexNum][MaxVertexNum];
};
typedef PtrToGNode MGraph;
//边集，用于逆序输出
struct Evode{
    int Head;
    int Tail;
    int lowcost;
}Edge[MaxVertexNum*(MaxVertexNum-1)/2+1];

int outdegree[MaxVertexNum];
MGraph ReadG(void)
{
    MGraph G=(MGraph)malloc(sizeof(struct GNode));
    scanf("%d%d",&G->Nv,&G->Ne);
    int i,j;
    for(i=1;i<=G->Nv;i++)
    {
        for(j=1;j<=G->Nv;j++)
        {
            G->G[i][j]=INFINITY;
            if(i==j)
                G->G[i][j]=0;
        }
    }
    Edge[0].Head=Edge[0].Tail=Edge[0].lowcost;
    for(i=1;i<=G->Ne;i++)
    {
        int a,b,x;
        scanf("%d%d%d",&a,&b,&x);
        G->G[a][b]=x;
        Edge[i].Head=a;
        Edge[i].Tail=b;
        Edge[i].lowcost=x;
    }

    return G;
}
int topo(MGraph G,int *result)
{
    int m=0;

    //初始化入度数组
    int pen[MaxVertexNum];

    int i,j;
    int stack[MaxVertexNum];
    int top=-1;
    for(i=1;i<=G->Nv;i++)
    {
        pen[i]=0;
        outdegree[i]=0;
        for(j=1;j<=G->Nv;j++)
        {
            if(G->G[j][i]!=INFINITY&&i!=j)
                pen[i]++;
            if(G->G[i][j]!=INFINITY&&i!=j)
                outdegree[i]++;
        }
        if(pen[i]==0)stack[++top]=i;
        //把入度为0的点压栈
    }
    while(top!=-1)
    {
        int t;
        result[m++]=t=stack[top--];

        for(int i=1;i<=G->Nv;i++)
        {
            if(G->G[t][i]!=INFINITY&&t!=i)
            {
                pen[i]--;
                if(pen[i]==0)
                    stack[++top]=i;
            }
        }
    }
    if(m<G->Nv)
        return 0;//有环路，无法生成topo序列
    else
        return 1;//成功产生topo序列
}
int least[MaxVertexNum];//交接点最早发生时间
int last[MaxVertexNum];//交接点最迟发生时间
int thing_least[MaxVertexNum*5];//事件最早发生时间
int thing_last[MaxVertexNum*5];//事件最迟发生时间
void key_path(MGraph G,int *result)
{
    int i;
    for(i=1;i<=G->Nv;i++)
    {
        least[i]=0;
        last[i]=INFINITY;
    }
    //按照topo序，生成最早发生时间
    for(i=1;i<=G->Nv;i++)
    {
        int t=result[i-1];
        for(int j=1;j<=G->Nv;j++)
        {
            if(G->G[t][j]!=INFINITY&&t!=j&&G->G[t][j]+least[t]>least[j])
                least[j]=G->G[t][j]+least[t];
        }
    }
    last[i-1]=least[i-1];
    least[0]=-1;
    int ans=least[0];
    for(i=1;i<=G->Nv;i++)
    {
        if(least[i]>ans)ans=least[i];
    }

    printf("%d\n",ans);
    //按照逆topo序，生成最迟发生时间
    int queue[MaxVertexNum];//建立一个简单队列
    int first = -1,rear = -1;

    for(int i=1; i<=G->Nv; i++)//将入度为0的节点压进去
    {
        if(outdegree[i] == 0)
        {
            queue[++rear] = i; //入队
            last[i] =ans;
        }
    }


    while(first < rear)
    {
        int v = queue[++first];//出队
        for( int i=G->Nv; i>=1; i--) //对V的每个邻接点处理
        {
            if(G->G[i][v]!=INFINITY&&i!=v)

            {
                outdegree[i]--;//出度减一
                last[i] = (last[i]<last[v]-G->G[i][v])?last[i]:last[v]-G->G[i][v];
                if(outdegree[i] == 0)
                {
                    queue[++rear] = i;
                }
            }

        }
    }
    int m=1;
    //生成每个事件的最早、最迟发生时间,并将关键事件输出
    for(i=1;i<=G->Nv;i++)
    {
       
        for(int j=G->Ne;j>=1;j--)
        {
            if(Edge[j].Head==i)
            {
                int r=Edge[j].Tail;
          
                if(least[i]==last[r]-G->G[i][r])
                printf("%d->%d\n",i,r);
            }
        }

    }
}
int main()
{
    MGraph G=ReadG();
    int result[MaxVertexNum];
    int f=topo(G,result);

    if(!f)printf("%d\n",f);
    else
        key_path(G,result);

    return 0;
}