Advertisement
atmoPunk

Untitled

Jun 22nd, 2018
85
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 2.12 KB | None | 0 0
  1. #include <algorithm>
  2. #include <iostream>
  3. #include <vector>
  4.  
  5. const int INF = 1000000000;
  6.  
  7. struct edge {
  8.   int from;
  9.   int to;
  10.   int cost;
  11. };
  12.  
  13. int n; // Кол-во вершин
  14. int e; // Кол-во рёбер
  15. std::vector<edge> gr(e);
  16. /* Ввод графа */
  17.  
  18. void FordBellman(int start) {
  19.   std::vector<int> d(n, INF); // Цены пути
  20.   std::vector<int> p(n, -1);  // Предки вершин
  21.   d[start] = 0;
  22.   int x;
  23.   for (int i = 0; i < n;
  24.        ++i) { // Выполняем n раз, чтобы найти отрицательный цикл
  25.     for (int j = 0; j < e; ++j) {
  26.       if (d[gr[j].from] < INF) {
  27.         if (d[gr[j].to] > d[gr[j].from] + gr[j].cost) {
  28.           d[gr[j].to] = d[gr[j].from] + gr[j].cost;
  29.           // если есть отрицательный цикл, чтобы не ушли за границы инта
  30.           d[gr[j].to] = std::max(d[gr[j].to], -INF);
  31.           p[gr[j].to] = gr[j].from;
  32.           x = gr[j].to; // Запоминаем последнее сделанное изменение
  33.         }
  34.       }
  35.     }
  36.     if (x == -1) { // Не было сделано ни одного изменения - дальнейшие итерации
  37.                    // бессмысленны
  38.       break;
  39.     }
  40.   }
  41.  
  42.   if (x == -1) {
  43.     // На последней фазе не было изменений - нет отрицательных циклов
  44.     std::cout << "No negative cycles" << std::endl;
  45.   } else {
  46.     std::vector<int> path;
  47.     int y = x; // x - либо лежит на отрицательном цикле, либо достижима из него
  48.     for (int i = 0; i < n; ++i) {
  49.       y = p[y];
  50.     } // y - точно лежит на цикле
  51.     std::vector<int> negcycle;
  52.     for (int cur = y;; cur = p[cur]) {
  53.       path.push_back(cur);
  54.       if (cur == y && path.size() > 1) {
  55.         break;
  56.       }
  57.     }
  58.     std::reverse(path.begin(), path.end());
  59.     std::cout << "Negative cycle" << std::endl;
  60.     for (int i = 0; i < path.size(); ++i) {
  61.       std::cout << path[i] << " ";
  62.     }
  63.     std::cout << std::endl;
  64.   }
  65. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement