Cod sursa(job #3271430)

Utilizator Alex_AeleneiAlex Aelenei Ioan Alex_Aelenei Data 26 ianuarie 2025 11:25:49
Problema Algoritmul Bellman-Ford Scor 10
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 3.12 kb
#include <iostream>
#include <fstream>
#include <vector>
#include <limits>

using namespace std;

void bellmanFord(int startNode, vector<vector<pair<int, int>>> &costMatrix, vector<int> &answer, vector<int> &predecessor, vector<int> &cycle)
{
    for (int nodeIndex = 1; nodeIndex < costMatrix.size(); nodeIndex++)
    {
        answer[nodeIndex] = numeric_limits<int>::max();
        predecessor[nodeIndex] = 0;
    }

    answer[startNode] = 0;
    for (int nodeIndex = 1; nodeIndex < costMatrix.size() - 1; nodeIndex++)
        for (int sourceNodeIndex = 1; sourceNodeIndex < costMatrix.size(); sourceNodeIndex++)
            for (auto itDestination = costMatrix[sourceNodeIndex].begin(); itDestination != costMatrix[sourceNodeIndex].end(); itDestination++)
                if (answer[sourceNodeIndex] + itDestination->second < answer[itDestination->first])
                {
                    answer[itDestination->first] = answer[sourceNodeIndex] + itDestination->second;
                    predecessor[itDestination->first] = sourceNodeIndex;
                }

    for (int sourceNodeIndex = 1; sourceNodeIndex < costMatrix.size(); sourceNodeIndex++)
        for (auto itDestination = costMatrix[sourceNodeIndex].begin(); itDestination != costMatrix[sourceNodeIndex].end(); itDestination++)
            if (answer[sourceNodeIndex] + itDestination->second < answer[itDestination->first])
            {
                answer[itDestination->first] = answer[sourceNodeIndex] + itDestination->second;

                vector<bool> visited(costMatrix.size(), false);
                visited[sourceNodeIndex] = true;

                int currentNode = itDestination->first;
                while (!visited[currentNode])
                {
                    visited[currentNode] = true;
                    currentNode = predecessor[currentNode];
                }

                cycle.clear();

                cycle.push_back(currentNode);
                sourceNodeIndex = predecessor[currentNode];
                while (sourceNodeIndex != currentNode)
                {
                    cycle.push_back(sourceNodeIndex);
                    sourceNodeIndex = predecessor[sourceNodeIndex];
                }

                return;
            }
}

int main()
{
    ifstream in("bellmanford.in");
    ofstream out("bellmanford.out");

    int nodeCount, vertexCount, source, destination, cost;
    in >> nodeCount >> vertexCount;

    vector<vector<pair<int, int>>> costMatrix(nodeCount + 1, vector<pair<int, int>>(100));
    vector<int> answer(nodeCount + 1, numeric_limits<int>::max()), predecessor(nodeCount + 1, 0), cycle;

    for (int vertexIndex = 1; vertexIndex <= vertexCount; vertexIndex++)
    {
        in >> source >> destination >> cost;
        costMatrix[source].push_back(pair<int, int>(destination, cost));
    }

    answer[1] = 0;
    bellmanFord(1, costMatrix, answer, predecessor, cycle);

    if (!cycle.empty())
    {
        out << "Ciclu negativ!";
        return 0;
    }

    for (int nodeIndex = 2; nodeIndex < answer.size(); nodeIndex++)
        out << ((answer[nodeIndex] == numeric_limits<int>::max()) ? 0 : answer[nodeIndex]) << ' ';

    return 0;
}