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

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

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

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

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

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

输入格式:

输入第1行给出两个正整数$N$($\le 100$)和$M$，其中$N$是任务交接点（即衔接相互依赖的两个子任务的节点，例如：若任务2要在任务1完成后才开始，则两任务之间必有一个交接点）的数量。交接点按1~$N$编号，$M$是子任务的数量，依次编号为1~$M$。随后$M$行，每行给出了3个正整数，分别是该任务开始和完成涉及的交接点编号以及该任务所需的时间，整数间用空格分隔。

输出格式:

如果任务调度不可行，则输出0；否则第1行输出完成整个工程项目需要的时间，第2行开始输出所有关键活动，每个关键活动占一行，按格式“V->W”输出，其中V和W为该任务开始和完成涉及的交接点编号。关键活动输出的顺序规则是：任务开始的交接点编号小者优先，起点编号相同时，与输入时任务的顺序相反。

序号

	输入	输出
1	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
2	9 11 1 2 6 1 3 4 1 4 5 2 5 1 3 5 1 4 6 2 5 7 9 5 8 7 6 8 4 7 9 2 8 9 4	18 1->2 2->5 5->8 5->7 7->9 8->9
3	11 14 1 2 4 1 3 3 2 4 5 3 4 3 4 5 1 4 6 6 5 7 5 6 7 2 8 3 7 9 3 7 9 10 6 4 10 2 10 6 5 6 11 4	21 3->4 4->10 6->11 8->3 9->3 10->6
4	4 5 1 2 4 2 3 5 3 4 6 4 2 3 4 1 2	0

关键路径求法：点击打开链接

百度文库里讲的很清楚。

坑点：就是有多起点和多终点的情况存在

#include<iostream>
#include<queue>
using namespace std;
#include<stdio.h>
#include<vector>
#include<string.h>
#include<algorithm>
#define MS(a,b) memset(a,b,sizeof(a))
int n,m,dfn[200][200],dis[200][200];
int vl[200],ve[200],in[200],out[200];
int mat[200][200],sum,ans;
vector<int>G[300],G1[300];
struct node
{
    int s,e,DFN;
}p[300];
bool cmp(node a,node b)
{  if(a.s==b.s)
     return a.DFN>b.DFN;
   return a.s<b.s;
}
void bfs_early()
{  int i;
    queue<int>q;
    for(i=1;i<=n;i++)
        if(!in[i])
           q.push(i);
    while(!q.empty())
    {
        int cur=q.front();
        q.pop();
     for(i=0;i<G[cur].size();i++)
     {
         int v=G[cur][i];
         ve[v]=max(ve[v],ve[cur]+dis[cur][v]);
         in[v]--;
        if(!in[v])
        {
           q.push(v);
           ans++;
        }

     }
    }
}
void bfs_later()
{
  int i;
  queue<int>q;
  for(i=1;i<=n;i++)
    if(!out[i])
      q.push(i),vl[i]=sum;
  while(!q.empty())
  {
    int cur=q.front();
    q.pop();
    for(i=0;i<G1[cur].size();i++)
    {
        int v=G1[cur][i];
        vl[v]=min(vl[v],vl[cur]-dis[cur][v]);
        out[v]--;
        if(!out[v])
            q.push(v);
    }
  }
}
int main()
{
   int i,ss,ee,d;
   ans=0 ;
  cin>>n>>m;
  MS(in,0);
  MS(out,0);
  MS(mat,0);
  MS(ve,0);
  for(i=1;i<=m;i++)
  {
    cin>>ss>>ee>>d;
    G[ss].push_back(ee);
    G1[ee].push_back(ss);
    mat[ss][ee]=1;
    dis[ss][ee]=dis[ee][ss]=d;
    in[ee]++;out[ss]++;
    dfn[ss][ee]=i;
  }
  queue<int>q;
  for(i=1;i<=n;i++)
     if(!in[i])
      {
        q.push(i);
        ans++;
      }
   for(i=1;i<=n;i++)
      vl[i]=99999;
    bfs_early();
   if(ans!=n)cout<<"0\n";
   else
   {

         sum=0;
     for(i=1;i<=n;i++)
        if(!out[i])
          sum=max(sum,ve[i]);//多终点时全部活动完成所需要的时间。
       bfs_later();
       cout<<sum<<endl;
       int j,k=0;
       for(i=1;i<=n;i++)
       {
         if(vl[i]==ve[i])
         {
             for(j=1;j<=n;j++)
              if(mat[i][j]&&vl[j]==ve[j]&&vl[j]-dis[i][j]==ve[i])
//一开始vl[j]-dis[i][j]==ve[i]这个条件没写，错了好多次，我以为求出关键结点就可以啦，百度文库的资料最后求关键活动没仔细看。。
                 p[k].s=i,p[k].e=j,p[k].DFN=dfn[i][j],k++;
         }
       }
        sort(p,p+k,cmp);
       for(i=0;i<k;i++)
         cout<<p[i].s<<"->"<<p[i].e<<endl;
   }
    return 0;
}