1. 基本概念
把用顶点表示事件,用弧表示活动,弧的权值表示活动所需要的时间的有向无环图称为边表示活动的网,简称AOE-网。
在AOE-网中存在唯一的、入度为0的顶点,称为源点;存在唯一的、出度为0的顶点,称为汇点。从源点到汇点的最长路径长度即为完成整个工程任务所需的时间,该路径称为关键路径。关键路径上的活动称为关键活动。这些活动中的任意一项活动未能按期完成,则整个工程的完成时间就要推迟;相反如果能加快关键活动的进度,则整个工程可以提前完成。
例如,上图所示的AOE-网中,共有v0,v1,v2,…,v8这9个事件。v0为源点,v8为汇点。v0到v8的关键路径(最长路径)有两条:v0,v1,v4,v6,v8或v0,v1,v4,v7,v8,长度均为18。关键活动为a1,a4,a7,a10或a1,a4,a8,a11。关键活动a1计划6天完成,若a1提前两天完成,则整个工程也可以提前两天完成。
几个基本定义:
拓扑排序见:https://blog.csdn.net/weixin_51450101/article/details/123018532?spm=1001.2014.3001.5501
① 事件vi的最早发生时间ve(i):从源点到顶点vi的最长路径长度。
求ve(i)可以从源点开始,按拓扑顺序向汇点逆推:
ve(0)=0;
ve(i)=Max{ve(k)+dut(<k,i>)},dut(<k,i>)表示与弧<k,i>对应活动的持续时间
② 时间vi的最晚发生时间vl(i):在保证汇点按其最早发生时间发生的这一前提下,求事件vi的最晚发生时间。即在求出ve(i)的基础上,从汇点开始,按逆拓扑顺序向源点递推,求出vl(i):
vl(n-1)=ve(n-1);
vl(i)=Min{vl(k)-dut(<i,k>)},dut(<i,k>)表示与弧<i,k>对应活动的持续时间
③ 活动ai的最早开始时间e(i):如果活动ai对应的弧为<j,k>,则e(i)等于从源点到顶点 j 的最长路径长度,即e(i)=ve(j)。
④ 活动ai的最晚开始时间l(i):如果活动ai对应的弧为<j,k>,其持续时间是dut(<j,k>),则有l(i)=vl(k)-dut(<j,k>)。
⑤ 活动ai的松弛时间(时间余量):ai的最晚开始时间与最早开始时间之差,即l(i)-e(i)。显然松弛时间为0的活动为关键活动。
2. 求关键路径的基本过程及计算步骤
基本过程:
① 对图中顶点进行拓扑排序,在排序过程中按逆拓扑序列求出每个事件的最早发生时间ve(i);
② 按逆拓扑序列求每个事件的最晚发生时间vl(i);
③ 求出每个活动ai的最早开始时间ei(i)和最晚开始时间l(i);
④ 找出e(i)=l(i)的活动ai即为关键活动。
计算步骤:
以上图所示AOE-网为例:
① 计算各顶点的最早开始时间
ve(0) = 0
ve(1) = Max {ve(0) + dut(<0,1>)} = 6
ve(2) = Max {ve(0) + dut(<0,2>)} = 4
ve(3) = Max {ve(0) + dut(<0,3>)} = 5
ve(4) = Max {ve(1) + dut(<1,4>),ve(2) + dut(<2,4>)} = 7
ve(5) = Max {ve(3) + dut(< 3,5>)} = 7
ve(6) = Max {ve(4) + dut(<4,6>)} = 16
ve(7) = Max {ve(4) + dut(<4,7>),ve(5) + dut(<5,7>)} = 14
ve(8) = Max {ve(6) + dut(<6,8>),ve(7) + dut(<7,8>)} = 18
② 计算各顶点的最迟开始时间
vl(8) = ve(8) = 18
vl(7) = Min {vl(8) - dut<7,8>} = 14
vl(6) = Min {vl(8) - dut<6,8>} = 16
vl(5) = Min {vl(7) - dut<5,7>} = 10
vl(4) = Min {vl(7) - dut<4,7>,vl(6) - dut(<4,6>)} = 7
vl(3) = Min {vl(5) - dut<3,5>} = 8
vl(2) = Min {vl(4) - dut<2,4>} = 6
vl(1) = Min {vl(4) - dut<1,4>} = 6
vl(0) = Min {vl(3) - dut(<0,3>),vl(2) - dut<0,2>,vl(1) - dut(<0,1>)} = 0
③ 计算各活动的最早开始时间
e(a1) = ve(0) = 0
e(a2) = ve(0) = 0
e(a3) = ve(0) = 0
e(a4) = ve(1) = 6
e(a5) = ve(2) = 4
e(a6) = ve(3) = 5
e(a7) = ve(4) = 7
e(a8) = ve(4) = 7
e(a9) = ve(5) = 7
e(a10) = ve(6) = 16
e(a11) = ve(7) = 14
④ 计算各活动的最迟开始时间
l(a11) = vl(8) - dut(<7,8>) = 14
l(a10) = vl(8) - dut(<6,8>) = 16
l(a9) = vl(7) - dut(<5,7>) = 10
l(a8) = vl(7) - dut(<4,7>) = 7
l(a7) =vl(6) - dut(<4,6>) = 7
l(a6) =vl(5) - dut(< 3,5>) = 8
l(a5) =vl(4) - dut(<2,4>) = 6
l(a4) =vl(4) - dut(<1,4>) = 6
l(a3) =vl(3) - dut(<0,3>) = 3
l(a2) =vl(2) - dut(<0,2>) = 2
l(a1) =vl(1) - dut(<0,1>) = 0
⑤ 计算活动ai的松弛时间
l(a1) - e(a1) = 0
l(a2) - e(a2) = 2
l(a3) - e(a3) = 3
l(a4) - e(a4) = 0
l(a5) - e(a5) = 2
l(a6) - e(a6) = 3
l(a7) - e(a7) = 0
l(a8) - e(a8) = 0
l(a9) - e(a9) = 3
l(a10) - e(a10) = 0
l(a11) - e(a11) = 0
松弛时间为0的活动即为关键活动,由计算结果得关键活动为:a1,a4,a7,a8,a10,a11
关键路径为:v0,v1,v4,v6,v8 或 v0,v1,v4,v7,v8
3. 实现代码及运行结果
求逆拓扑序列及事件最早发生时间ve(i)
int indegree[MAX_VERTEX_NUM]; //存放各顶点入度数
int ve[MAX_VERTEX_NUM]; //各顶点的最早发生时间
void FindID(AdjList G); //求各顶点入度算法
/*拓扑排序,求逆拓扑序列及事件最早发生时间ve(i)*/
int TopoSort(AdjList G, SeqStack* T) {
//栈T用于生成逆拓扑序列
int i, k, count = 0;
SeqStack S; //栈S用于存放入度为0的顶点
ArcNode* p;
FindID(G); //求各顶点入度
InitStack(T); //初始化栈T
InitStack(&S); //初始化栈S
for (i = 0; i < G.vexnum; i++) {
if (indegree[i] == 0)
Push(&S, i); //将入度为0的顶点入栈S
}
for (i = 0; i < G.vexnum; i++)
ve[i] = 0; //初始化最早发生时间
while (!IsEmpty(&S)) {
i = S.elem[S.top];
Pop(&S);
count++;
Push(T, i); //按拓扑顺序进入栈T
p = G.vertex[i].firstarc;
while (p != NULL) {
k = p->adjvex;
indegree[k]--; //i号顶点的每个邻接点的入度减1
if (indegree[k] == 0)
Push(&S, k); //若入度减为0则入栈
if (ve[i] + p->weight > ve[k]) //按拓扑顺序计算事件最早发生时间
ve[k] = ve[i] + p->weight;
p = p->nextarc;
}
}
if (count < G.vexnum) {
printf("该图存在回路!\n");
return FALSE;
}
else
return TRUE;
}
关键路径算法
/*关键路径算法*/
void CriticalPath(AdjList G, SeqStack* T) {
ArcNode* p;
int i, j, k, dut, ei, li;
char tag;
int vl[MAX_VERTEX_NUM]; //每个顶点的最迟发生时间
TopoSort(G, T); //① 求事件最早发生时间和逆拓扑序列栈T
for (i = 0; i < G.vexnum; i++)
vl[i] = ve[G.vexnum - 1]; //将各事件的最晚发生时间初始化为汇点的最早发生时间
while (!IsEmpty(T)) {
//② 按逆拓扑顺序求各顶点的最晚发生时间vl值
j = T->elem[T->top];
Pop(T);
p = G.vertex[j].firstarc;
while (p != NULL) {
k = p->adjvex;
dut = p->weight;
if (vl[k] - dut < vl[j])
vl[j] = vl[k] - dut;
p = p->nextarc;
}
}
for (j = 0; j < G.vexnum; j++) {
//③ 求各活动的最早开始时间ei和最晚开始时间li
p = G.vertex[j].firstarc;
while (p != NULL) {
k = p->adjvex;
dut = p->weight;
ei = ve[j];
li = vl[k] - dut;
if (ei == li) //④ 输出关键活动及对应的路径
printf("a%d v%c->v%c\n", p->activity, G.vertex[j].data, G.vertex[k].data);
p = p->nextarc;
}
}
}
完整实现代码
/*AOE-网的关键路径*/
# include<stdio.h>
# include<malloc.h>
# define MAX_VERTEX_NUM 20
# define TRUE 1
# define FALSE 0
/*AOE-网的邻接表表示法*/
typedef char VertexData;
//弧结点结构
typedef struct ArcNode {
int adjvex; //该弧指向顶点的位置
struct ArcNode* nextarc; //指向下一条弧的指针
int activity; //弧表示的活动
int weight; //权值
}ArcNode;
//表头结点结构
typedef struct VertexNode {
VertexData data; //顶点数据
ArcNode* firstarc; //指向该顶点的第一条弧的指针
}VertexNode;
//邻接表结构
typedef struct {
VertexNode vertex[MAX_VERTEX_NUM];
int vexnum, arcnum; //图的顶点数和弧数
}AdjList;
/*求顶点位置*/
int LocateVertex(AdjList* G, VertexData v) {
int k;
for (k = 0; k < G->vexnum; k++) {
if (G->vertex[k].data == v)
break;
}
return k;
}
/*创建AOE-网的邻接表*/
int CreateAdjList(AdjList* G) {
int i, j, k, activity, weight;
VertexData v1, v2;
ArcNode* p;
printf("输入图的顶点数和弧数:"); //输入图的顶点数和弧数
scanf("%d%d", &G->vexnum, &G->arcnum);
printf("输入图的顶点:");
for (i = 0; i < G->vexnum; i++) {
//输入图的顶点,初始化顶点结点
scanf(" %c", &(G->vertex[i].data));
G->vertex[i].firstarc = NULL;
}
for (k = 0; k < G->arcnum; k++) {
printf("输入第%d条弧的两个顶点、弧表示的活动及权值:", k + 1);
scanf(" %c %c %d %d", &v1, &v2, &activity, &weight); //输入一条弧的两个顶点、弧表示的活动及权值
i = LocateVertex(G, v1);
j = LocateVertex(G, v2);
p = (ArcNode*)malloc(sizeof(ArcNode)); //申请新弧结点
p->activity = activity;
p->weight = weight;
p->adjvex = j;
p->nextarc = G->vertex[i].firstarc;
G->vertex[i].firstarc = p;
}
}
/*顺序栈的存储结构*/
typedef struct {
int elem[MAX_VERTEX_NUM]; //用于存放栈中元素的一维数组
int top; //存放栈顶元素的下标,top为-1表示空栈
}SeqStack;
/*初始化顺序栈*/
void InitStack(SeqStack* S) {
S->top = -1;
}
/*判空*/
int IsEmpty(SeqStack* S) {
if (S->top == -1) //栈为空
return TRUE;
else
return FALSE;
}
/*顺序栈进栈*/
int Push(SeqStack* S, int x) {
if (S->top == MAX_VERTEX_NUM - 1) //栈已满
return FALSE;
S->top++;
S->elem[S->top] = x; //x进栈
return TRUE;
}
/*顺序栈出栈*/
int Pop(SeqStack* S) {
if (S->top == -1) //栈为空
return FALSE;
S->top--;
return TRUE;
}
int indegree[MAX_VERTEX_NUM]; //存放各顶点入度数
int ve[MAX_VERTEX_NUM]; //各顶点的最早发生时间
/*求各顶点入度算法*/
void FindID(AdjList G) {
int i;
ArcNode* p;
for (i = 0; i < G.vexnum; i++)
indegree[i] = 0;
for (i = 0; i < G.vexnum; i++) {
p = G.vertex[i].firstarc;
while (p != NULL) {
indegree[p->adjvex]++;
p = p->nextarc;
}
}
}
/*拓扑排序,求逆拓扑序列及事件最早发生时间ve(i)*/
int TopoSort(AdjList G, SeqStack* T) {
//栈T用于生成逆拓扑序列
int i, k, count = 0;
SeqStack S; //栈S用于存放入度为0的顶点
ArcNode* p;
FindID(G); //求各顶点入度
InitStack(T); //初始化栈T
InitStack(&S); //初始化栈S
for (i = 0; i < G.vexnum; i++) {
if (indegree[i] == 0)
Push(&S, i); //将入度为0的顶点入栈S
}
for (i = 0; i < G.vexnum; i++)
ve[i] = 0; //初始化最早发生时间
while (!IsEmpty(&S)) {
i = S.elem[S.top];
Pop(&S);
count++;
Push(T, i); //按拓扑顺序进入栈T
p = G.vertex[i].firstarc;
while (p != NULL) {
k = p->adjvex;
indegree[k]--; //i号顶点的每个邻接点的入度减1
if (indegree[k] == 0)
Push(&S, k); //若入度减为0则入栈
if (ve[i] + p->weight > ve[k]) //按拓扑顺序计算事件最早发生时间
ve[k] = ve[i] + p->weight;
p = p->nextarc;
}
}
printf("\n事件最早发生时间为:"); //输出事件的最早发生时间
for (i = 0; i < G.vexnum; i++)
printf("%d ", ve[i]);
if (count < G.vexnum) {
printf("该图存在回路!\n");
return FALSE;
}
else
return TRUE;
}
/*关键路径算法*/
void CriticalPath(AdjList G, SeqStack* T) {
ArcNode* p;
int i, j, k, dut, ei, li;
char tag;
int vl[MAX_VERTEX_NUM]; //每个顶点的最迟发生时间
TopoSort(G, T); //① 求事件最早发生时间和逆拓扑序列栈T
for (i = 0; i < G.vexnum; i++)
vl[i] = ve[G.vexnum - 1]; //将各事件的最晚发生时间初始化为汇点的最早发生时间
while (!IsEmpty(T)) {
//② 按逆拓扑顺序求各顶点的最晚发生时间vl值
j = T->elem[T->top];
Pop(T);
p = G.vertex[j].firstarc;
while (p != NULL) {
k = p->adjvex;
dut = p->weight;
if (vl[k] - dut < vl[j])
vl[j] = vl[k] - dut;
p = p->nextarc;
}
}
printf("\n事件最晚发生时间为:"); //输出事件的最晚发生时间
for (i = 0; i < G.vexnum; i++)
printf("%d ", vl[i]);
printf("\n关键活动及对应的弧为:\n");
for (j = 0; j < G.vexnum; j++) {
//③ 求各活动的最早开始时间ei和最晚开始时间li
p = G.vertex[j].firstarc;
while (p != NULL) {
k = p->adjvex;
dut = p->weight;
ei = ve[j];
li = vl[k] - dut;
if (ei == li) //④ 输出关键活动及对应的路径
printf("a%d v%c->v%c\n", p->activity, G.vertex[j].data, G.vertex[k].data);
p = p->nextarc;
}
}
}
int main() {
AdjList G;
SeqStack T;
CreateAdjList(&G);
CriticalPath(G, &T);
return 0;
}
运行结果
参考:耿国华《数据结构——用C语言描述(第二版)》
更多数据结构内容关注我的《数据结构》专栏:https://blog.csdn.net/weixin_51450101/category_11514538.html?spm=1001.2014.3001.5482