Cod sursa(job #2244398)

Utilizator RaduXD1Nicolae Radu RaduXD1 Data 22 septembrie 2018 17:58:01
Problema Ubuntzei Scor 20
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.73 kb
#include <fstream>
#include <queue>
#include <algorithm>
using namespace std;
ifstream fin("ubuntzei.in");
ofstream fout("ubuntzei.out");
vector<pair <int, int> > v[2001];
vector<pair <int, int> >::iterator it;
priority_queue<pair<int,int> ,  vector<pair<int,int> > , greater<pair<int,int> > > h;
int n,m,i,a,b,c,cost,nod,vecin,f[131072],nr,d[131072],mini;
int ma[18][18],fa[2001],k,l[18],r,j;



void solve(int poz)
{
    for(i=1;i<=n;i++)
    {
        d[i]=2000000000;
        f[i]=0;
    }
    d[poz]=0;
    h.push(make_pair(d[poz], poz));
    while(!h.empty())
    {
        nod=h.top().second;
        h.pop();
        while(f[nod]==1)
        {
            if(h.empty())
                break;
            nod=h.top().second;
            h.pop();
        }
        f[nod]=1;
        nr--;
        for(it=v[nod].begin();it!=v[nod].end();it++)
        {
            vecin=it->first;
            cost=it->second;
            if(f[vecin]==0&&d[vecin]>d[nod]+cost)
            {
                d[vecin]=d[nod]+cost;
                h.push(make_pair(d[vecin], vecin));
            }
        }
    }
    for(i=1;i<=n;i++)
        if(i!=poz&&fa[i]!=0)
        {
            ma[fa[poz]][fa[i]]=min(d[i],ma[fa[poz]][fa[i]]);
            ma[fa[i]][fa[poz]]=min(d[i],ma[fa[i]][fa[poz]]);
        }
}

int main()
{
    fin>>n>>m;
    fin>>k;
    for(i=1;i<=k;i++)
        fin>>l[i];
    sort(l+1, l+k+1);
    if(l[k]!=n)
        l[++k]=n;
    if(l[1]!=1)
        l[++k]=1;
    sort(l+1, l+k+1);
    for(i=1;i<=k;i++) for(j=1;j<=k;j++) ma[i][j]=2000000000;
    for(i=1;i<=k;i++)
        fa[l[i]]=i;
    for(i=1;i<=m;i++)
    {
        fin>>a>>b>>c;
        v[a].push_back(make_pair(b, c));
        v[b].push_back(make_pair(a, c));
    }
    for(j=1;j<=k;j++)
        solve(l[j]);
    for(i=1;i<k;i++)
        r+=(1<<i);
    for(i=0;i<=r;i++)
    {
        d[i]=2000000000;
        f[i]=0;
    }
    d[r]=0;
    h.push(make_pair(d[r], r));
    while(f[0]==0)
    {
        nod=h.top().second;
        h.pop();
        while(f[nod]==1)
        {
            nod=h.top().second;
            h.pop();
        }
        f[nod]=1;
        for(i=0;i<k;i++)
            if((nod&(1<<i))==0)
                for(j=0;j<k;j++)
                    if((nod&(1<<j))!=0)
                    {
                        vecin=nod-(1<<j);
                        cost=ma[i+1][j+1];
                        if(f[vecin]==0&&d[vecin]>d[nod]+cost)
                        {
                            d[vecin]=d[nod]+cost;
                            h.push(make_pair(d[vecin], vecin));
                        }
                    }
    }
    fout<<d[0];
    fin.close();
    fout.close();
    return 0;
}