Cod sursa(job #2470002)

Utilizator AndreiD31Dragan Andrei AndreiD31 Data 8 octombrie 2019 16:32:35
Problema Ciclu hamiltonian de cost minim Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.99 kb
#include <bits/stdc++.h>

using namespace std;

ifstream f("hamilton.in");
ofstream g("hamilton.out");

int n,m,i,x,y,sol,cost,j,v,k,vecin[20][20],ct[20][20],dp[2][262150][20];
int main()
{
    f>>n>>m;

    for(i=1; i<=m; i++)
    {
        f>>x>>y>>cost;
        ct[x][y]=cost;
        vecin[x][y]=1;
    }

    for(i=0;i<(1<<n);i++)
    {
        for(j=0;j<n;j++)
            dp[0][i][j]=INT_MAX/2; // ca sa fac minim, initializez matricea dp
    }

    // Initializez nodul 0
    dp[0][(1<<0)][0]=0;

    //Matricea "dp" poate fi privita ca o matrice bidimensionala.
    // Trebuie sa incep neaparat cu nodurile de la 0, pentru ca altfel se strica configuratia binara cand fac pentru j

    // Cum e ciclu, pornesc doar de la un nod i(fie el 0) si schimb doar j si k
    for(j=0;j<(1<<n);j++) // pentru fiecare configuratie de j. (daca j in baza 2 este de exemplu 0110, inseamna ca de la i la k iau nodurile 1 si 2
    {
        for(k=0;k<n;k++) // pentru fiecare nod, ca sa fac de la i la k
        {
            if((j&(1<<k))) // daca nodul k se afla in configuratia lui j, adica daca il iau si pe el
            {
                for(v=0;v<n;v++) // verific vecinii intermediari
                {
                    if(vecin[v][k]==1 && ((j&(1<<v)))) // iau vecinul doar daca se afla in configuratia lui j, daca se afla pe drumul ales de la nodul i la nodul k
                    {
                        dp[0][j][k]=min(dp[0][j][k], dp[0][(j^(1<<k))][v]+ct[v][k]); // aici e dinamica
                    }
                }
            }
        }
    }


    sol=INT_MAX;
    for(i=0;i<n;i++)
    {
        if(vecin[i][0]==1 && dp[0][(1<<n)-1][i]!=INT_MAX/2) // daca este vecin si daca are o valoare (daca inca este INT_MAX/2 atunci nu ma intereseaza pentru ca inseamna ca nu exista solutie in acel caz)
        {
            sol=min(sol,dp[0][(1<<n)-1][i]+ct[i][0]);
        }
    }
    if(sol==INT_MAX)g<<"Nu exista solutie"<<'\n';
    else g<<sol;

    return 0;
}