Pagini recente » Cod sursa (job #1933927) | Cod sursa (job #1723063) | Cod sursa (job #3004299) | Cod sursa (job #2941815) | Cod sursa (job #2949998)
#include <fstream>
#include <vector>
using namespace std;
ifstream f("ubuntzei.in");
ofstream g("ubuntzei.out");
int n,h[2011],k,r,w[22],dist[2011],m,x,y,z,ver[2011];
int oo=999999999;
int b[21][21],ab[20][(1<<16)];
vector<pair<int, int>>a[2001];
void update(int poz)
{
while(dist[h[poz]]<dist[h[poz/2]])
{
swap(ver[h[poz]],ver[h[poz/2]]);
swap(h[poz],h[poz/2]);
poz/=2;
}
}
void erase_heap()
{
swap(ver[h[1]],ver[h[k]]);
swap(h[1],h[k]);
k--;
int ok=1,cc=2,tt=1;
while(cc<=k&&ok)
{
ok=0;
if(cc!=k)
{
if(dist[h[cc]]>dist[h[cc+1]])
cc++;
}
if(dist[h[cc]]<dist[h[tt]])
{
swap(ver[h[tt]],ver[h[cc]]);
swap(h[tt],h[cc]);
tt=cc;
cc*=2;
ok=1;
}
}
}
void add_heap(int val)
{
h[++k]=val;
int poz=k;
update(poz);
}
void dijkstra(int s,int TTT)
{
for(int i=1; i<=n; i++)
{
ver[i]=-1;
dist[i]=oo;
}
k=0;
dist[s]=0;
ver[s]=1;
h[++k]=s;
while(k)
{
int u=h[1];
erase_heap();
int cost=dist[u];
for(int i=0; i<a[u].size(); i++)
{
int x1=a[u][i].first;
int y1=a[u][i].second;
if(dist[x1]>y1+cost)
{
if(dist[x1]==oo)
{
ver[x1]=k+1;
dist[x1]=y1+cost;
add_heap(x1);
}
else
{
dist[x1]=y1+cost;
update(ver[x1]);
}
}
}
}
for(int i=0;i<=r;i++)
if(w[i]!=s)
b[TTT][i]=b[i][TTT]=dist[w[i]];
}
void solutie()
{
for(int i=0; i<r; i++)
for(int j=0; j<(1<<r); j++)
ab[i][j]=oo;
ab[0][1]=0;
for(int j=1; j<(1<<r); j++)
for(int i=0; i<r; i++)
if(j&(1<<i))
for(int k=0; k<r; k++)
if((1<<k)&j)
ab[i][j]=min(ab[i][j],b[k][i]+ab[k][(j^(1<<i))]);
g<<ab[r-1][(1<<r)-1];
}
int main()
{
f>>n>>m>>r;
for(int i=1;i<=r;i++)
f>>w[i];
w[++r]=n;
w[0]=1;
for(int i=1; i<=m; i++)
{
f>>x>>y>>z;
a[x].push_back(make_pair(y,z));
a[y].push_back(make_pair(x,z));
}
for(int i=0;i<=r;i++)
for(int j=0;j<=r;j++)
b[i][j]=oo;
for(int i=0;i<=r;i++)
dijkstra(w[i],i);
r++;
solutie();
return 0;
}