在C ++中计算DAG的关键路径

时间:2011-05-21 08:25:32

标签: c++ path graph-theory dynamic-programming directed-acyclic-graphs

我正在根据this algorithm对另一篇文章计算图像DAG的关键路径。我的老师要求实现aarray,我简化了家庭作业陈述,通过以下方式实现了一个简单的图表阵列。

enter image description here

这是我的代码,其中我有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个顶点的距离。

1 个答案:

答案 0 :(得分:1)

您的队列包含顶点。您的数组uvd由边数编号。 所以你不能写

first = Q.front();
... u[first] ...

因为first是一个顶点。

更一般地说,如果使用有意义的变量名,您的代码将更容易阅读(并且错误将更加明显)。 first不是很明确(首先是什么?),而uvd也非常神秘。

写下像

这样的东西
cur_vertex = todo.front()
distance[dest[cur_vertex]] = std::max(distance[dest[cur_vertex]], 
    distance[source[cur_vertex]]+ weight[cur_vertex]);

会立即提出一个问题:顶点的来源,那是什么? (这里我们使用变量名作为正确类型检查的替代.ADA程序员会声明两种不同的整数类型,以避免顶点和边数之间的混淆。)

另一个问题:first的后继者的循环在哪里?它在伪代码中,但不在您的源代码中。