Поэтому я пытался реализовать алгоритм Дейкстры для кратчайшего пути в ориентированном графе, используя списки смежности, но я не знаю, по какой причине, он не выводит результаты (выводит минимальное расстояние как 0 для всех узлов) ,
Код, который я написал:
#include <fstream>
#include <functional>
#include <climits>
#include <vector>
#include <queue>
#include <list>
using namespace std;
struct node {
int vertex;
int weight;
node(int v, int w) : vertex(v), weight(w) { };
node() { }
};
class CompareGreater {
public:
bool const operator()(node &nodeX, node &nodeY) {
return (nodeX.weight > nodeY.weight) ;
}
};
vector< list<node> > adj;
vector<int> weights;
priority_queue<node, vector<node>, CompareGreater> Q;
int nrVertices, nrEdges;
void readData();
void Dijkstra(node);
void writeData();
int main(int argc, char *argv[]) {
readData();
Dijkstra(node(1, 0));
writeData();
return 0;
}
void readData() {
fstream in("dijkstra.in", ios::in);
int nodeX, nodeY, weight;
in >> nrVertices >> nrEdges;
adj.resize(nrVertices+1);
weights.resize(nrVertices+1);
for (int i = 1; i <= nrVertices; ++i) {
weights.push_back(INT_MAX);
}
for (int i = 1; i <= nrEdges; ++i) {
in >> nodeX >> nodeY >> weight;
adj[nodeX].push_back(node(nodeY, weight));
}
in.close();
}
void Dijkstra(node startNode) {
node currentNode;
weights[startNode.vertex] = 0;
Q.push(startNode);
while (!Q.empty()) {
currentNode = Q.top();
Q.pop();
if (currentNode.weight <= weights[currentNode.vertex]) {
for (list<node>::iterator it = adj[currentNode.vertex].begin(); it != adj[currentNode.vertex].end(); ++it) {
if (weights[it->vertex] > weights[currentNode.vertex] + it->weight) {
weights[it->vertex] = weights[currentNode.vertex] + it->weight;
Q.push(node((it->vertex), weights[it->vertex]));
}
}
}
}
}
void writeData() {
fstream out("dijkstra.out", ios::out);
weights.resize(nrVertices+1);
for (vector<int>::iterator it = weights.begin()+1; it != weights.end(); ++it) {
out << (*it) << " ";
}
out.close();
}
Входные данные были:
5 7
1 2 10
1 3 2
1 5 100
2 4 3
3 2 5
4 3 15
4 5 5
Это означает, что есть 5 узлов, 7 дуг (направленных ребер), и дуги существуют от 1 до 2 со стоимостью 10, от 1 до 3 со стоимостью 2 и так далее.
Однако вывод неправильный. Я понятия не имею, где может произойти сбой программы. Я взял основную идею отсюда: http://community.topcoder.com/tc?module=Static&d1 = учебники&d2 = standardTemplateLibrary2 # dijkstra1
(В конце он дает идею алгоритма Дейкстры с использованием priority_queue).
Заранее спасибо.
Raul
Проблема в линии
weights.resize(nrVertices+1);
в readData()
, Это устанавливает вектор с nrVertices+1
элементы значения 0. Позже вы добавляете фактические значения, которые вы хотите к этому вектору, используя weights.push_back(INT_MAX);
,
В фактическом алгоритме Дейкстры все интересное weights
Таким образом, 0, а не INT_MAX
ты хочешь.
Заменить строку на
weights.resize(1);
(просто чтобы убедиться, что индекс 1 действительно ссылается на первый элемент — кажется, вы используете 1 в качестве первого индекса вместо 0), и это может сработать.
Других решений пока нет …