我正在根据this algorithm对另一篇文章计算图像DAG的关键路径。我的老师要求实现aarray,我简化了家庭作业陈述,通过以下方式实现了一个简单的图表阵列。
这是我的代码,其中我有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);
}
我做错了什么或如何解决代码告诉我项目的持续时间(对应于距离关键路径的距离之和)?只能计算到前3个顶点的距离。
答案 0 :(得分:1)
您的队列包含顶点。您的数组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
的后继者的循环在哪里?它在伪代码中,但不在您的源代码中。