根据this algorithm在另一篇文章中,我正在计算图像DAG的关键路径。我的老师要求实现数组,简化家庭作业陈述,这是通过数组实现的简单图形。
这是我的代码,其中有3个数组v,u和d,分别表示边的起点节点,边的终点节点和每对顶点之间的距离,如上图所示。在图像的图形中,项目的持续时间等于25,对应于与关键路径的距离之和。
我的代码无法根据this link的伪代码来计算距离
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <queue>
#include <iostream>
#include <algorithm>
using namespace std;
int main (){
int numberVertex=6; //number of vertex
int numberActivities=9;//number of edges
int i, j;
/*vertices are of the form (v, u) */
//indegree of each vertex (that is, count the number of edges entering them)
int indegree [6]={0,0,0,0,0,0};
int v[9]={0,0,1,1,2,4,3,3,3};//array v represent the starting vertex of de edge
int u[9]={1,2,2,3,4,5,4,5,5};//array u represent the coming vertex of de edge
int d[9]={5,6,3,8,2,12,0,1,4};//array d represent the time of de activity (v,u)
int project_duration=0;//project duration
/*# Compute the indegree for each vertex v from the graph:
for each neighbor u of v: indegree[u] += 1*/
for (j=0; j<numberActivities; j++){
indegree[u[j]]++;
}
for (j=0;j<numberVertex; j++)
printf ("indegree %d= %d\n",j,indegree[j] );
queue<int> Q; //queue Q = empty queue
int distance [numberVertex];
memset(distance, 0, sizeof(int) * numberVertex);//distance = array filled with zeroes
//for each vertex v:
//if indegree[v] = 0:
//insert v on Q
for (j=0; j<numberVertex; j++)
{
if (indegree[j]==0)
Q.push(v[j]);
}
int first;
//printf ("first in the queue=%d\n", Q.front());
/*for each neighbor u of v:
d istance[u] = max(distance[u], distance[v] + time(v, u))
indegree[u] -= 1
if indegree[u] = 0:
insert u on Q
*/
while (!Q.empty()){ //while Q is not empty:
first= Q.front (); //v = get front element from Q
Q.pop(); //delete de first from queue
distance[u[first]]=std::max(distance[u[first]],
distance[v[first]]+ d[first]);
indegree[u[first]]-=1;
if (indegree[u[first]]==0){
Q.push(u[first]);
}
}
for (j=0; j<numberVertex; j++)
{
printf ("dist [%d]= %d\n", j, distance[j]);
}
/*Now, select the vertex x with the largest distance.
This is the minimum total project_duration.*/
printf ("Total Project Duration %d\n", project_duration);
return (0);
}
我在做错什么或如何解决代码以告诉我项目的持续时间是多少(对应于与关键路径的距离之和)?只能计算到前三个顶点的距离。
最佳答案
您的队列包含顶点。您的数组u
,v
,d
由边号索引。
所以你不能写
first = Q.front();
... u[first] ...
因为
first
是一个顶点。更一般而言,如果您使用有意义的变量名,您的代码将更容易阅读(错误将更加明显)。
first
不是很明确(首先是什么?),而且u
,v
和d
也是很隐秘的。写像
cur_vertex = todo.front()
distance[dest[cur_vertex]] = std::max(distance[dest[cur_vertex]],
distance[source[cur_vertex]]+ weight[cur_vertex]);
马上会提出一个问题:顶点的来源是什么?
(在这里,我们使用变量名来代替适当的类型检查。ADA程序员将声明两个不同的整数类型,以避免顶点和边号之间的混淆。)
另一个问题:
first
后继者的循环在哪里进行?它在伪代码中,但不在您的源代码中。