08-图9 关键活动（浙大数据结构）

中国大学MOOC-陈越、何钦铭-数据结构-2022夏

2 天

08-图9 关键活动

分数 30

作者 DS课程组

单位 浙江大学

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

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

// 错误已更正，结构体加一个flag来判断边的先后顺序；
/**
    e[r],l[r]:分别表示活动 ar 的最早开始时间和最迟开始时间；
    当e[r] == l[r] 说明活动 ar 是关键活动；

    ve[i],vl[i] : 表示事件 i 的最早发生时间和最迟发生时间；

    事件 Vi 经过活动 ar 到达 事件 Vj；那么 e[r]=ve[i], l[r]=vl[j]-length[r];
*/

#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>

using namespace std;

struct ENode
{
    int v;
    int cost;
    int flag; //记录边数输入的顺序，输出保证起点相同时，与输入时的任务相反
};

typedef pair<int,int> PII;

const int maxn = 110;

vector<ENode> Adj[maxn];
int ve[maxn],vl[maxn];
int inDegree[maxn];

struct Node
{
    int u,v,flag;
};

vector<Node> res;
stack<int> topOrder;  //用栈存储拓朴排序的结果，方便后面计算vl 数组；

int Nv,Ne;

bool cmp(Node &a,Node &b)
{
    if(a.u==b.u)
        return  a.flag>b.flag;
    else
        return a.u < b.u;
}

void Read(); //读入数据
bool topSort();
int CriticalPath(); //关键路径

int main()
{
    Read();

     int ans = CriticalPath();

    cout << ans << endl;

    if(ans!=0)
    {
        sort(res.begin(),res.end(),cmp);
        for(auto a:res)
            printf("%d->%d\n",a.u,a.v);
    }

    return 0;
}

void Read() //读入数据
{
    int num=1;
    cin >> Nv >> Ne;
    for(int i=0;i<Ne;++i)
    {
        int u,v,cost;
        cin >> u >> v >> cost;
        Adj[u].push_back({v,cost,num++});
        inDegree[v]++;
    }
}

bool topSort()
{
    queue<int> q;
    for(int i=1;i<=Nv;++i)
        if(inDegree[i]==0)
            q.push(i);

    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        topOrder.push(u);

        for(int i=0;i<Adj[u].size();++i)
        {
            int v=Adj[u][i].v;
            inDegree[v]--;
            if(inDegree[v]==0)
                q.push(v);

            if(ve[u]+Adj[u][i].cost > ve[v])   //通过前驱节点 ve[il] ~ ve[ik]的值更新 ve[j]的值：
                ve[v] = ve[u]+Adj[u][i].cost;  //ve[j] = max{ve[ip] + length[rp]};
        }
    }

    if(topOrder.size() == Nv)
        return true;
    else
        return false;
}

int CriticalPath() //关键路径
{
    if(topSort() == false)
        return 0;

    int Maxval=*max_element(ve+1,ve+Nv+1);
    fill(vl,vl+Nv+1,Maxval);

    while(!topOrder.empty())
    {
        int u=topOrder.top();
        topOrder.pop();

        for(int i=0;i<Adj[u].size();++i)    //通过 vi的后继节点 vjl~vjk 来推出vi的值：
        {
            int v=Adj[u][i].v;              //vi = min{vl[jp]-length[rp]};
            int cost = Adj[u][i].cost;
            if(vl[v]-cost < vl[u])
                vl[u] = vl[v]-cost;
        }
    }

    for(int u=1;u<=Nv;++u)
        for(int i=0;i<Adj[u].size();++i)
        {
            int v=Adj[u][i].v;
            int cost = Adj[u][i].cost;
            int e=ve[u],l=vl[v]-cost;

            if(e==l)
                res.push_back({u,v,Adj[u][v].flag});
        }

    return Maxval;
}