AOE网上的关键路径
Time Limit: 1000MS Memory limit: 65536K
题目描述
一个无环的有向图称为无环图(Directed Acyclic Graph),简称DAG图。
AOE(Activity On Edge)网:顾名思义,用边表示活动的网,当然它也是DAG。与AOV不同,活动都表示在了边上,如下图所示:
如上所示,共有11项活动(11条边),9个事件(9个顶点)。整个工程只有一个开始点和一个完成点。即只有一个入度为零的点(源点)和只有一个出度为零的点(汇点)。
关键路径:是从开始点到完成点的最长路径的长度。路径的长度是边上活动耗费的时间。如上图所示,1 到2 到 5到7到9是关键路径(关键路径不止一条,请输出字典序最小的),权值的和为18。
输入
输出
示例输入
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 8 9 4 7 9 2
示例输出
18 1 2 2 5 5 7 7 9
关键路径的算法:
(1). 输入e条弧<j, k>,建立AOE-网的存储结构。
(2). 从源点v0出发,令ve[0]=0,按拓扑有序求其余各顶点的最早发生时间ve[i](1<=i<=n-1)。如果得到的拓扑有序序列中顶点个数小于网中顶点数n,则说明网中存在环,不能求关键路径,算法终止;否则执行步骤(3)。
(3). 从汇点vn出发,令vl[n-1]=ve[n-1],按逆拓扑有序求其余各顶点的最迟发生时间vl[i](2<=i<=n-2)。
(4). 根据各顶点的ve和vl值,求每条弧s的最早开始时间e(s)和最迟开始时间l(s)。若某条弧满足关系e(s)=l(s),则该弧为关键活动。
PS:(1). 求ve——从ve(0)=0开始向前递推
(2). 求vl——从vl(n-1)=ve(n-1)起向后递推
#include <stdio.h>
#include <stdlib.h>
#include <stack>
#include <string.h>
using namespace std;
typedef struct arcnode//表结点
{
int adj; //存储节点
struct arcnode *next;
int info;//存储权值
} arcnode;
typedef struct vnode//头结点
{
int data;
arcnode *first;
} adjlist[10010];
typedef struct
{
adjlist a;
int vn, an;
} ALG;
int n, m;
int indegree[10010];//记录每个点的入度
int ve[10010]; //最早开始时间
int vl[10010]; //最迟开始时间
void create(ALG &g)
{
int i, j, v1, v2, w;
arcnode *p;
for(i=1; i<=n; i++)
{
g.a[i].first = NULL; //头结点清空
}
for(i=1; i<=m; i++)
{
scanf("%d %d %d", &v1, &v2, &w);
p = new arcnode;
p->adj = v2;
p->info = w;
indegree[v2]++;
p->next = g.a[v1].first;
g.a[v1].first = p;
}
}
stack <int> t; //用于逆拓扑有序时求vl
int topo1(ALG &g)//求 ve(最早开始时间即最大路径长度)
{
int i, j, k;
arcnode *p;
stack <int> s;
for(i=1; i<=n; i++)//入度为0的结点入栈
if(!indegree[i])
s.push(i);
memset(ve, 0, sizeof(ve));
int count = 0;//记录出栈顶点个数
while(!s.empty())
{
j = s.top();
s.pop();
t.push(j);
count++;
for(p=g.a[j].first; p; p=p->next) //删去所有已j为起点的出边
{
k = p->adj;
if(!(--indegree[k]))//将新入度为0的结点入栈
s.push(k);
if(ve[j]+p->info>ve[k]) //更新最大路径长度
ve[k] = ve[j]+p->info;
}
}
if(count<n)
return 0;
else
return 1;
}
int topo2(ALG &g)//求vl(即不影响施工进度的情况下最晚开始的时间)
{
int i, j, k;
arcnode *p;
if(!topo1(g)) return 0;
printf("%d\n", ve[n]);
for(i=1; i<=n; i++) //初始化最迟发生时间
vl[i] = ve[n];
while(!t.empty()) //按拓扑逆序求vl
{
j = t.top();
t.pop();
for(p = g.a[j].first; p; p=p->next)
{
k = p->adj;
if(vl[k] - p->info < vl[j]) //更新最晚开始时间
vl[j] = vl[k] - p->info;
}
}
int flag = 0, a, b;
for(j=1; j<=n; j++) //求关键活动
{
for(p = g.a[j].first; p; p=p->next)
{
k = p->adj;
int e = ve[j];
int l = vl[k] - p->info;
if(e == l)
{
if(flag == 0)//对输出的特殊处理,要求字典序最小
{
a = j;
b = k;
flag = 1;
}
else if(a==j&&b>k)
b = k;
else if(b==j)
{
printf("%d %d\n", a, b);
a = j;
b = k;
}
}
}
}
printf("%d %d\n", a, b);
}
int main()
{
ALG g;
while(~scanf("%d %d", &n, &m))
{
memset(indegree, 0, sizeof(indegree));
create(g);
//topo1(g);
topo2(g);
}
return 0;
}