Voici une implémentation de l'algorithme de Dijkstra je l'ai écrit à partir du pseudocode dans le Wikipedia article. Pour un graphe avec environ 40 000 nœuds et 80 000 arêtes, il faut compter 3 ou 4 minutes. Est-ce que c'est quelque chose comme le bon ordre de grandeur? Si non, quel est le problème avec ma mise en œuvre?Performance de l'implémentation de l'algorithme de Dijkstra
struct DijkstraVertex {
int index;
vector<int> adj;
vector<double> weights;
double dist;
int prev;
bool opt;
DijkstraVertex(int vertexIndex, vector<int> adjacentVertices, vector<double> edgeWeights) {
index = vertexIndex;
adj = adjacentVertices;
weights = edgeWeights;
dist = numeric_limits<double>::infinity();
prev = -1; // "undefined" node
opt = false; // unoptimized node
}
};
void dijsktra(vector<DijkstraVertex*> graph, int source, vector<double> &dist, vector<int> &prev) {
vector<DijkstraVertex*> Q(G); // set of unoptimized nodes
G[source]->dist = 0;
while (!Q.empty()) {
sort(Q.begin(), Q.end(), dijkstraDistComp); // sort nodes in Q by dist from source
DijkstraVertex* u = Q.front(); // u = node in Q with lowest dist
u->opt = true;
Q.erase(Q.begin());
if (u->dist == numeric_limits<double>::infinity()) {
break; // all remaining vertices are inaccessible from the source
}
for (int i = 0; i < (signed)u->adj.size(); i++) { // for each neighbour of u not in Q
DijkstraVertex* v = G[u->adj[i]];
if (!v->opt) {
double alt = u->dist + u->weights[i];
if (alt < v->dist) {
v->dist = alt;
v->prev = u->index;
}
}
}
}
for (int i = 0; i < (signed)G.size(); i++) {
assert(G[i] != NULL);
dist.push_back(G[i]->dist); // transfer data to dist for output
prev.push_back(G[i]->prev); // transfer data to prev for output
}
}
Merci pour votre réponse. J'ai l'impression que ma mise en œuvre actuelle n'est pas outrageusement mauvaise, et qu'avec vos suggestions, je pourrais m'attendre à un temps d'exécution de 1 à 3 minutes pour un problème avec 40 000 nœuds. Des exécutions plus proches de 30 secondes ou 1 seconde ne sont pas raisonnables. Est-ce vrai? – zoo