Cod sursa(job #2852803)

Utilizator ana_madalina_18Radu Ana Madalina ana_madalina_18 Data 19 februarie 2022 16:23:39
Problema Ciclu hamiltonian de cost minim Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 2.3 kb
#include <fstream>
#include <vector>
#include <iostream>
using namespace std;
ifstream fin("hamilton.in");
ofstream fout("hamilton.out");
int n, m;
struct muchie
{
    int nod_curent, nod_urm, cost;
};

vector <muchie> v[20];
vector <muchie> vecini_externi;
int dp[20][300000];
const int MAAX=1e9;
int main()
{
    fin>>n>>m;
    for(int i=1;i<=m;i++)
    {
        int a, b, c;
        fin>>a>>b>>c;
        muchie nou;
        nou.cost=c;
        nou.nod_curent=a;
        nou.nod_urm=b;
        v[a].push_back(nou);
        if(b==0)
        {
            vecini_externi.push_back(nou);
        }
    }
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<=(1<<n)-1;j++)
        {
            dp[i][j]=MAAX;
        }
    }
    //dp[0][0]=0;
    dp[0][1]=0;
    //cout<<dp[1][27]<<' '<<MAAX<<endl;
    for(int i=1;i<=(1<<n)-1;i+=2)
    {
        for(int j=0;j<n;j++)
        {
            for(int x=0;x<v[j].size();x++)
            {
                muchie vecin=v[j][x];
                //cout<<dp[1][27]<<endl;
                if(j==3 && i==27 && v[j][x].nod_urm==2)
                {
                    //cout<<dp[j][i]<<' '<<vecin.cost<<' '<<dp[vecin.nod_urm][i+(1<<vecin.nod_urm)]<<' ';
                }
                int putere=(1<<vecin.nod_urm);
                if((i & putere)==0)
                {
                    if(dp[vecin.nod_urm][i+putere]==10 && vecin.nod_urm==1 && (i+putere)==27)
                    {
                        //cout<<i<<' '<<j<<' '<<vecin.cost<<' '<<dp[j][i]<<endl;
                    }

                    dp[vecin.nod_urm][i+putere]=min(
                        dp[vecin.nod_urm][i+putere],
                        dp[j][i]+vecin.cost
                    );



                }
            }
        }
    }
    //cout<<dp[3][27]<<' ';
    int minim=MAAX;
    for(int i=0;i<vecini_externi.size();i++)
    {
       muchie vecin=vecini_externi[i];
       //cout<<vecin.nod_curent<<' ';
       //cout<<dp[vecin.nod_curent][(1<<n)-1]<<endl;
       if(dp[vecin.nod_curent][(1<<n)-1]+vecin.cost<minim)
        {
            minim=dp[vecin.nod_curent][(1<<n)-1]+vecin.cost;
        }
    }

    if(minim>=MAAX)
    {
        fout<<"Nu exista solutie";
    }
    else
    {
        fout<<minim;
    }
    return 0;
}