Алгоритм Дейкстры - в с ++? - PullRequest
12 голосов
/ 10 августа 2010

за последние четыре дня я пытаюсь понять алгоритм Дейкстры.Но я не могу.У меня есть вектор точек.Из этого я создал матрицу затрат.Но я не знаю, как сделать алгоритм Дейкстры.Источники доступны в сети, но я не из области компьютерных наук, поэтому я не могу их понять.Я пытаюсь сделать такую ​​функцию

vector<int> dijkstra(costMatrix[][])
{
  ....
  ....
  return vector<int>pathPointindex
}

main()
{
    vector<Point> availablePoints;
    costMatrix[][]=createCostMatrix();
    vector<int> indexes=dijikstra(costMatrix)
    for(int i=0;i<indexes.size();i++)
       cout << "path points are " << availablePoints[indexes[i]] << endl;
}

Если кто-нибудь, не могли бы вы опубликовать код.Я не ленивый.Но мой проект уже превысил крайний срок один день назад.Теперь я потерял надежду понять логику.Теперь просто я хочу функцию.«Человек в беде - это действительно ангел».

РЕДАКТИРОВАТЬ: Особая благодарность "Loki astari" за отличный ответ

Ответы [ 6 ]

40 голосов
/ 10 августа 2010

алгоритм Дейкстры

По-английски:

Это алгоритм поиска кратчайшего маршрута из точки А в точку Б.
В вычислительном отношении мы упрощаем маршрут к графу, состоящему из узлов и дуг. Каждый узел представляет промежуточную точку, в то время как каждая дуга соединяет два узла и имеет (неотрицательный) вес, представляющий стоимость прохождения между двумя узлами.

Для реализации алгоритма вам понадобятся два списка:

  • закончено: набор (узел, стоимость), в котором вы вычислили минимальную стоимость, которую нужно достичь.
  • работает: отсортированный список (узел, стоимость), которые были проверены.

Алгоритм:

working.addNode(Start, 0); // No cost to get to start.

for( (node, cost) = working.popHead(); node != End; (node,cost) = working.popHead())
{
    // If we have already processed this node ignore it.
    if (finished.find(node))
    {    continue;
    }

    // We have just removed a node from working.
    // Because it is the top of the list it is guaranteed to be the shortest route to
    // the node. If there is another route to the node it must go through one of the
    // other nodes in the working list which means the cost to reach it will be higher
    // (because working is sorted). Thus we have found the shortest route to the node.

    // As we have found the shortest route to the node save it in finished.
    finished.addNode(node,cost);

    // For each arc leading from this node we need to find where we can get to.
    foreach(arc in node.arcs())
    {
        dest = arc.dest();
        if (NOT (finished.find(dest)))
        {
            // If the node is already in finished then we don't need to worry about it
            // as that will be the shortest route other wise calculate the cost and add
            // this new node to the working list.
            destCost = arc.cost() + cost;
            working.addNode(dest,destCost); // Note. Working is sorted list
        }
    }
} 

Так что, если вы думаете об этом алгоритме. Скажем, вы едете из Лондона в Манчестер.

finished = {} // empty.
working  = { (London,0) }

Использование следующей матрицы затрат:

                  L    S    O    B    N    M    W
(L) ondon         -    50   60   100  130  -    -
(S) outhampton    50   -    70   -    -    -    -
(O) xford         60   70   -    50   -    200  -
(B) irmingham     100  -    50   -    -    80   110
(N) orwich        130  -    -    -    -    -    -
(M) anchester     -    -    200  80   -    -    80
Ne(W) castle      -    -    -    110  -    80   -

Теперь вы вынимаете Лондон из рабочего списка (как он есть во главе) и помещаете его в готовый список. Затем добавьте в рабочий список все города, непосредственно связанные с Лондоном.

finished = { (London,0) }
working  = { (Southampton, 50), (Oxford, 60), (Birmingham, 100), (Norwich,130) }

Рассмотрим города в рабочем наборе, внешний край пузыря которого расширился от Лондона. Задача алгоритма Дейкстры состоит в том, чтобы продолжать расширять пузырь, пока мы не достигнем Манчестера (без повторения каких-либо шагов, которые мы уже предприняли). Таким образом, пузырь всегда расширяется наружу, и мы всегда расширяем самую маленькую часть пузыря.

Итак, следующий шаг - взять заголовок списка и повторить. Из Саутгемптона есть только два пункта назначения. Вернуться в Лондон (который мы отбрасываем, как он есть в готовом списке) и Оксфорд Стоимость проезда в Оксфорд составляет 50 + стоимость от Саутгемптона до Оксфорда (поэтому обратите внимание, что он дважды находится в рабочем списке, но не волнуйтесь, позже мы откажемся от него, так как это не оптимальный маршрут).

finished = { (London,0), (Southampton,50) }
working  = { (Oxford, 60), (Birmingham, 100), (Oxford, 120), (Norwich,130) }

Так повторите цикл. Главой списка является Оксфорд. Из Оксфорда мы можем поехать в Манчестер (200), Бирмингем (50) или обратно в Лондон (60) или Саутгемптон (помните, что нам нужно добавить стоимость проезда в Оксфорд к каждой из этих затрат выше. отправился в Саутгемптон, но мы уже нашли кратчайший путь в Саутгемптон, поэтому обработка не требуется). Это оставит нас с:

finished = { (London,0), (Southampton,50), (Oxford, 60) }
working  = { (Birmingham, 100), (Birmingham,110), (Oxford, 120), (Norwich,130), (Manchester,200)}

Обратите внимание, что у нас есть Манчестер в рабочем списке (это наш пункт назначения). Но нам нужно продолжать работать, поскольку мы можем найти более короткий маршрут. Так что теперь мы расширяемся из Бирмингема. Оттуда мы можем пойти в Оксфорд (50), Манчестер (80), Лондон (100), Ньюкасл (110). Если добавить стоимость поездки в Бирмингем, это дает нам:

finished = { (London,0), (Southampton,50), (Oxford, 60), (Birmingham, 100) }
working  = { (Birmingham,110), (Oxford, 120), (Norwich,130), {Manchester, 180), (Manchester,200), (Newcastle, 210)}

Следующие два узла. Оксфорд и Бирмингем уже в готовом списке, поэтому мы можем их игнорировать. Таким образом, если нет маршрута из Норвича в Манчестер, который составляет менее 50 миль, мы доберемся до Манчестера в итерации после этого.

13 голосов
/ 10 августа 2010

Я советую вам взглянуть на TopCoder учебник, который очень прагматичен.Вам нужно будет выяснить, как работает очередь приоритетов STL, и убедиться, что на вашем графике нет negative edge weights.

Достаточно полная реализация может быть найдена здесь .Вам нужно будет добавить к нему вектор Path и реализовать метод RecoverPath, чтобы получить узлы на пути от источника к приемнику.Чтобы использовать это решение, вам также необходимо преобразовать adjacency matrix в adjacency list следующим образом:

for (int i=0;i<nNodes;i++)
    for (int j=0;j<nNodes;j++){
        if (costMatrix[i][j] != NO_EDGE_VALUE){
            G[i].pb(make_pair(j,costMatrix[i],[j]));
        }
    }

РЕДАКТИРОВАТЬ: Если ваш график плотный, я бы посоветовал вам использовать Ford Bellman *Алгоритм 1015 * намного проще и не должен быть намного медленнее.

EDIT2: чтобы вычислить путь, вы должны добавить в заголовок

int P[MAX]; /*array with links to parents*/
for(i=0; i<=nodes; i++) P[i] = -1; /*magic unset value*/

// dijkstra
while(!Q.empty()) {
    ....
    if(!F[v] && D[u]+w < D[v]) {
        D[v] = D[u] + w;
        /*By setting P[v] value we will remember what is the 
          previous node on path from source */
        P[v] = u; // <-- added line
        Q.push(pii(v, D[v]));
    }
    ...
}

Затем вы должны добавить метод RecoverPath (он работаеттолько если путь существует)

vector<int> RecoverPath(int src, int dest){
    vector<int> path;
    int v = dest;
    while (v != src) {
        path.push_back(v);
        v = P[v];
    }
    path.push_back(src);
    std::reverse(path.begin(),path.end());
    return path;
}
5 голосов
/ 21 марта 2014
#include <iostream>
#include <vector>
#include <string>
#include <list>
#include <limits> 
#include <set>
#include <utility>
#include <algorithm> 
#include <iterator>

using namespace std;


typedef int vertex_t;
typedef double weight_t;

const weight_t max_weight = numeric_limits<double>::infinity();

struct neighbor {
    vertex_t target;
    weight_t weight;
    neighbor(vertex_t arg_target, weight_t arg_weight)
        : target(arg_target), weight(arg_weight) { }
};

typedef vector<vector<neighbor> > adjacency_list_t;

// Computing the shortest pathway

void
DijkstraComputePaths(vertex_t source,
                     const adjacency_list_t &adjacency_list,
                     vector<weight_t> &min_distance,
                     vector<vertex_t> &previous)
{
    int n = adjacency_list.size();
    min_distance.clear();
    min_distance.resize(n, max_weight);
    min_distance[source] = 0;
    previous.clear();
    previous.resize(n, -1);
    set<pair<weight_t, vertex_t> > vertex_queue;
    vertex_queue.insert(make_pair(min_distance[source], source));

    while (!vertex_queue.empty())
    {
        weight_t dist = vertex_queue.begin()->first;
        vertex_t u = vertex_queue.begin()->second;
        vertex_queue.erase(vertex_queue.begin());

        // Visit each edge exiting u
        const vector<neighbor> &neighbors = adjacency_list[u];
        for (vector<neighbor>::const_iterator neighbor_iter = neighbors.begin();
            neighbor_iter != neighbors.end();
            neighbor_iter++)
        {
            vertex_t v = neighbor_iter->target;
            weight_t weight = neighbor_iter->weight;
            weight_t distance_through_u = dist + weight;
            if (distance_through_u < min_distance[v]) {
                vertex_queue.erase(make_pair(min_distance[v], v));

                min_distance[v] = distance_through_u;
                previous[v] = u;
                vertex_queue.insert(make_pair(min_distance[v], v));

            }
        }
    } // while
}
2 голосов
/ 06 августа 2012

Основная идея алгоритма Дейкстры довольно проста: предположим, у вас есть набор точек с известными кратчайшими путями к заданной точке A. Затем давайте предположим, что мы хотим добавить новую точку C в набор.Позволяет определить, какие точки из набора связаны с точкой, которую мы хотим добавить.Пусть это будет точка B (i). Поэтому для всех точек B (i) мы будем новыми, чтобы найти сумму расстояний между A, B (i) и B (i) и C. Наименьшее из этих расстояний будет минимальным междуА и С.

1 голос
/ 10 августа 2010

Реализация на с ++

#include <cstdio>
#include <cstring>
#include <set>
#include <vector>
using namespace std;

#define pb push_back
#define mp make_pair
#define MAXN 50100
#define INF 1000000000

int N, M, d[MAXN]; vector<int> G[MAXN], C[MAXN];
set< pair<int, int> > T;

void solve(void)
{
    int i, j, k, val, x;

    for(i = 2; i <= N; i++) d[i] = INF;
    T.insert( mp(0, 1) );

    while( T.size() > 0 )
    {
        val = (*T.begin()).first, x = (*T.begin()).second;
        T.erase(*T.begin());
        for(i = 0; i < G[x].size(); i++)
         if(d[ G[x][i] ] > val + C[x][i] )
            d[ G[x][i] ] = val + C[x][i], T.insert(mp(d[G[x][i]],G[x][i]));
    }
}

int main(void)
{
    freopen("dijkstra.in", "rt", stdin);
    freopen("dijkstra.out", "wt", stdout);

    int i, a, b, c;

    scanf("%d %d\n", &N, &M);

    for(i = 1; i <= M; i++)
        scanf("%d %d %d\n", &a, &b, &c), G[a].pb(b), C[a].pb(c);

    solve();

    for(i = 2; i <= N; i++)
        printf("%d ", d[i] == INF ? 0 : d[i]);

    return 0;
}
0 голосов
/ 08 января 2017
#include<iostream>
#include<vector>
#include<algorithm>
#include<map>
#include<queue>


using namespace std;

const size_t INT_MAX = 0xFFFFFFFF; // or any other value representing infinite distance.

сначала создайте ребро структуры, содержащее индекс исходного узла, индекс целевого узла и ребро "weight" (длина).

struct edge { size_t from; size_t to; size_t length; };

определяет узел класса, содержащий ребра для соседних соседей.

class Node 
{ 
public:
    void AddNeighborEdge( edge _NeighborEdge ) { m_neighborsEdges.push_back( _NeighborEdge ); } 
    vector<edge>::iterator FirstNeighborEdge() { return  m_neighborsEdges.begin(); }
    vector<edge>::iterator LastNeighborEdge() { return  m_neighborsEdges.end(); }

private: 
     vector<edge>  m_neighborsEdges; 
};

класс NeighborsDistanceUpdator будет использоваться в качестве «функтора» алгоритмом for_each для итеративного прохождения и обновления минимального расстояния от текущего узла в графе до соседних соседей.

class NeighborsDistanceUpdator
{
public:
    NeighborsDistanceUpdator( vector<size_t>& _min_distance_from_source, queue< size_t >& _nodes_to_visit ) : m_min_distance_from_source( _min_distance_from_source ),                                                              m_nodes_to_visit( _nodes_to_visit ) 
                                                            {} 
    void operator()( edge& _edge )
    {   
        size_t from = _edge.from;
        size_t to = _edge.to;

        if ( m_min_distance_from_source[ to ] > m_min_distance_from_source[ from ] + _edge.length ) 
        {
            m_min_distance_from_source[ to ] = m_min_distance_from_source[ from ] + _edge.length;
            m_nodes_to_visit.push( to );
        }
    }    

private:
    vector<size_t>& m_min_distance_from_source;
    queue< size_t >& m_nodes_to_visit;
};

Что касается алгоритма Дейкстры, просто запустите все узлы на графике и для каждого узла обновите минимальное расстояние от источника (если меньше), сохранив соседние узлы для посещения.

size_t dijkstra( map< size_t, Node  >& _graph, size_t _sourceIndex, size_t _targetIndex ) 
{
    vector<size_t> min_distance_from_source( _graph.size(), INT_MAX );
    min_distance_from_source[ _sourceIndex ] = 0;
    queue< size_t > nodes_to_visit;
    nodes_to_visit.push( _sourceIndex );
    NeighborsDistanceUpdator neighborsDistanceUpdator( min_distance_from_source, nodes_to_visit );

    while ( ! nodes_to_visit.empty() ) 
    {

        size_t currNodeIndex = nodes_to_visit.front();

        if ( currNodeIndex ==  _targetIndex ) return min_distance_from_source[ currNodeIndex ];

        nodes_to_visit.pop();

        vector<edge>::iterator firstNeighborEdge= _graph[ currNodeIndex ].FirstNeighborEdge();
        vector<edge>::iterator lastNeighborEdge= _graph[ currNodeIndex ].LastNeighborEdge();

        for_each( firstNeighborEdge, lastNeighborEdge, neighborsDistanceUpdator );
    }
    return INT_MAX;
}

тест ...

int main()
{
    Node node1;
    Node node2;
    Node node3;
    Node node4;

    map< size_t, Node > graph;
    edge ed;

    ed.from = 0;
    ed.to = 1;
    ed.length = 1;
    node1.AddNeighborEdge( ed );

    cout << "node: " << 0 << " to: " << ed.to ;
    cout << " lenth: " << ed.length << endl << endl;

    ed.from = 0;        
    ed.to = 2;
    ed.length = 4;
    node1.AddNeighborEdge( ed );
    graph.insert( make_pair( 0, node1 ) );

    cout << "node: " << 0 << " to: " << ed.to ;
    cout << " lenth: " << ed.length << endl << endl;

    ed.from = 1;
    ed.to = 2;
    ed.length = 1;
    node2.AddNeighborEdge( ed );

    cout << "node: " << 1 << " to: " << ed.to ;
    cout << " lenth: " << ed.length << endl << endl;

    ed.from = 1;
    ed.to = 3;
    ed.length = 3;
    node2.AddNeighborEdge( ed );
    graph.insert( make_pair( 1, node2 ) );

    cout << "node: " << 1 << " to: " << ed.to ;
    cout << " lenth: " << ed.length << endl << endl;

    ed.from = 2;
    ed.to = 3;
    ed.length = 1;
    node3.AddNeighborEdge( ed );
    graph.insert( make_pair( 2, node3 ) );

    cout << "node: " << 2 << " to: " << ed.to ;
    cout << " lenth: " << ed.length << endl << endl;

    ed.from = 3;
    ed.to = INT_MAX;
    ed.length = INT_MAX;
    node3.AddNeighborEdge( ed );
    graph.insert( make_pair( 3, node4 ) );

    cout << "node: " << 2 << " to: " << ed.to ;
    cout << " lenth: " << ed.length << endl << endl;


    cout << "min length from: 1 to 4 = " << dijkstra( graph, 0,3  ) << endl;
}
...