Cod sursa(job #2048996)

Utilizator rangal3Tudor Anastasiei rangal3 Data 26 octombrie 2017 19:14:56
Problema Ciclu hamiltonian de cost minim Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.53 kb
#include <fstream>
#include <vector>
#include <cstring>

using namespace std;

const int MAXN = 20;
const int MAXX = (1<<18) + 6;
const int INF = 100000000;

int N, M, Sol;
vector <int> A[MAXN];
int Cost[MAXN][MAXN];
int C[MAXX][MAXN];

int calc(int start, int conf, int last)
{
    if (C[conf][last] == -1)
    {
        C[conf][last] = INF;

        for (size_t i = 0; i < A[last].size(); ++i)
            if (conf & (1<<A[last][i]))
            {
                if (A[last][i] == start && ((conf ^ (1<<last)) != (1<<start))) continue;

                C[conf][last] = min(C[conf][last], calc(start, conf ^ (1<<last), A[last][i]) + Cost[A[last][i]][last]);
            }
    }

    return C[conf][last];
}

int main()
{
    ifstream fin("hamilton.in");
    ofstream fout("hamilton.out");

    fin >> N >> M;

    for (int i = 0; i < N; ++i)
        for (int j = 0; j < N; ++j) Cost[i][j] = INF;

    for (int i = 1; i <= M; ++i)
    {
        int x, y, c;

        fin >> x >> y >> c;

        A[y].push_back(x); // retin muchia invers
        Cost[x][y] = c;
    }

    Sol = INF;

    for (int i = 0; i < N; ++i)
    {
        int nod = i;
        memset(C, -1, sizeof(C));

        C[1<<i][i] = 0;

        for (size_t j = 0; j < A[i].size(); ++j)
        {
            int nodant = A[i][j];
            Sol = min(Sol, calc(nod, (1<<N)-1, nodant) + Cost[nodant][nod]);
        }
    }

    if (Sol != INF) fout << Sol << endl;
    else fout << "Nu exista solutie" << endl;

    return 0;
}