数据结构与算法——7-11 关键活动 (30分)

7-11 关键活动 (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为该任务开始和完成涉及的交接点编号。关键活动输出的顺序规则是：任务开始的交接点编号小者优先，起点编号相同时，与输入时任务的顺序相反。

思路：这道题就是拓补排序 关键路径问题 关键就是判断是否是关键路径，也就是当最早开始时间和最晚开始时间相等时 就意味着没有空闲，该活动是关键活动，活动间的路径就是关键路径
所以就要求出最早发生时间 earliest time 和最晚发生时间 latest time
求最早发生时间就是从头到尾找拓补序列的过程
earliest[k]=0 当 k=0时
earliest[k]=max{earliest[i]+len<vi,vk>},当k不等于0且<vi,vk>属于P[k]时
P[k]表示所有到达顶点vk的弧的集合 len<vi,vk>表示弧<vi,vk>上的权值

求最晚发生时间就是把拓补排序倒过来进行的
latest[k]=0 当 k=n-1时
latest[k]=max{latest[j]-len<vk,vj>},当k<n-1 且<vk,vj>属于S[k]时
S[k]表示所有从顶点vk出发的弧的集合 len<vk,vj>表示弧<vk,vj>上的权值

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

#define MAXVER 105
#define INFINITY 65535

int G[MAXVER][MAXVER];  //图
int earliest[MAXVER];      //最早发生时间
int latest[MAXVER];      //最迟发生时间
int in[MAXVER];        //入度
int out[MAXVER];       //出度
int num_vertex,num_edge;            //顶点数目 ，边数目

void CreatGraph();
int EarliestTime();
void LatestTime(int Scost);
int FindMax( int a,int b);
int FindMin( int a,int b);

int main()
{
    int flag;
    int i,j;
    scanf("%d %d",&num_vertex,&num_edge);

    CreatGraph();//建图
    flag = EarliestTime();
    if( flag==-1)
    {
        printf("0\n");
    }
    else
    {
        printf("%d\n",flag);
        LatestTime( flag );
        for( i=1; i<=num_vertex; i++)
        {
            if(earliest[i] != latest[i])
                continue;
            for( j=num_vertex; j>=1 ; j--)
            {
                if( G[i][j]>=0 && earliest[j]==latest[j] &&latest[j]-G[i][j]==earliest[i])
                {
                    //i,j均在关键路径上且相邻
                    printf("%d->%d\n",i,j);
                }
            }
        }

    }
    return 0;
}

void CreatGraph()
{
    int i,j;
    int s,d,weight;

    for( i=1; i<=num_vertex; i++)
    {
        for( j=1; j<=num_vertex; j++)
        {
            G[i][j] = -1;
        }
        earliest[i] = 0;
        latest[i] = INFINITY;
        in[i] = 0;
        out[i] = 0;
    }
    for( i=0; i<num_edge; i++)
    {
        scanf("%d %d %d",&s,&d,&weight);
        G[s][d] = weight;   //有向边
        in[d] ++;
        out[s]++;
    }

}

int EarliestTime()
{
    int queue[num_vertex];
    int first =-1,rear = -1;//定义队列
    int count=0;
    int i;
    int temp,ret=0;

    for( i=1; i<=num_vertex; i++)
    {
        if( in[i]==0)
        {
            //如果入度为0则入队
            queue[++rear] = i;
        }
    }

    while( first<rear)    //判断队是否为空
    {
        temp = queue[++first];   //出队
        count++;
        for( i=1; i<=num_vertex; i++)
        {
            if( G[temp][i]>=0 )
            {
                in[i]--;
                earliest[i] = FindMax( earliest[i],earliest[temp]+G[temp][i]);
                if( in[i]==0)
                {
                    queue[++rear] = i;//如果入度减少过后为0，则入队
                }
            }
        }
    }
    if( count!=num_vertex)
    {
        ret = -1;
    }
    else
    {
        ret = earliest[1];
        for( i=2; i<=num_vertex; i++)
        {
            if(earliest[i] > ret)
            {
                //找出最大的earliest[i]
                ret = earliest[i];
            }
        }
    }

    return ret;
}

void LatestTime(int Sweight)
{
    int i;
    int queue[MAXVER];
    int first=-1,rear=-1;
    int temp;

    for( i=1; i<=num_vertex; i++)
    {
        if( out[i]==0)
        {
            queue[++rear] = i;
            latest[i] = Sweight;
        }
    }

    while( first<rear )
    {
        temp = queue[++first];
        for( i=num_vertex; i>=1; i--)
        {
            if( G[i][temp]>=0)
            {
                latest[i] = FindMin( latest[i],latest[temp]-G[i][temp]);
                out[i]--;
                if(out[i]==0)
                {
                    queue[++rear] = i;
                }
            }
        }

    }

}

int FindMax( int a,int b)
{
    return (a>b)?a:b;
}
int FindMin( int a,int b)
{
    return (a>b)? b:a;
}