Pagini recente » Cod sursa (job #2037781) | Cod sursa (job #297859) | Cod sursa (job #1443725) | Cod sursa (job #2030219) | Cod sursa (job #2136711)
#include<cstdio>
#include<vector>
#define MAX_N 2000
#define MAX_K 15
#define oo 0x3f3f3f3f
using namespace std;
vector<pair<int,int> >g[MAX_N + 1];
int dp[MAX_N + 1][MAX_N + 1], dist[MAX_N + 1], h[MAX_N + 1], pos[MAX_N + 1], v[MAX_K + 1], st[MAX_K + 1];
int n, m, k, heapSize, minCost;
bool used[MAX_N + 1];
FILE *fout = fopen("ubuntzei.out","w");
void readGraph() {
int i, x, y, cost;
FILE *fin = fopen("ubuntzei.in","r");
fscanf(fin,"%d%d%d",&n,&m,&k);
for(i = 0; i < k; i++)
fscanf(fin,"%d",&v[i]);
for(i = 0; i < m; i++) {
fscanf(fin,"%d%d%d",&x,&y,&cost);
g[x].push_back(make_pair(y,cost));
g[y].push_back(make_pair(x,cost));
}
fclose(fin);
}
inline void Swap(int i, int j) {
int aux = h[i];
h[i] = h[j];
h[j] = aux;
aux = pos[h[i]];
pos[h[i]] = pos[h[j]];
pos[h[j]] = aux;
}
inline void heapUp(int i) {
if(dist[h[i/2]] <= dist[h[i]]) return;
Swap(i,i/2);
heapUp(i/2);
}
inline void heapDown(int i) {
int dleftSon, drightSon;
if(2 * i > heapSize) return;
dleftSon = dist[h[2*i]];
if(2 * i + 1 <= heapSize)
drightSon = dist[h[2*i+1]];
else drightSon = dleftSon + 1;
if(dleftSon <= drightSon) {
if(dist[h[i]] <= dleftSon) return;
Swap(i,2*i);
heapDown(2*i);
} else {
if(dist[h[i]] <= drightSon) return;
Swap(i,2*i+1);
heapDown(2*i+1);
}
}
void Dijkstra(int source) {
vector<pair<int,int> >::iterator it;
int i, node;
for(i = 1; i <= n; i++) {
dist[i] = oo;
h[i] = pos[i] = i;
}
dist[source] = 0;
Swap(1,source);
heapSize = n;
for(i = 1; i <= n; i++) {
node = h[1];
Swap(1,heapSize);
heapSize--;
heapDown(1);
for(it = g[node].begin(); it != g[node].end(); it++) {
if(dist[(*it).first] > dist[node] + (*it).second) {
dist[(*it).first] = dist[node] + (*it).second;
heapUp(pos[(*it).first]);
}
}
}
for(i = 1; i <= n; i++)
dp[source][i] = dp[i][source] = dist[i];
}
void precalcDistances() {
Dijkstra(1);
Dijkstra(n);
for(int i = 0; i < k; i++)
Dijkstra(v[i]);
}
void GenPerm(int K) {
if(K == k) {
int cost = 0;
for(int i = 0; i < k - 1; i++)
cost += dp[v[st[i]]][v[st[i+1]]];
cost += dp[1][v[st[0]]] + dp[v[st[K-1]]][n];
if(cost < minCost)
minCost = cost;
} else {
for(int i = 0; i < k; i++)
if(!used[v[i]]) {
st[K] = i;
used[v[i]] = true;
GenPerm(K+1);
used[v[i]] = false;
}
}
}
int main() {
readGraph();
if(k == 0) {
Dijkstra(1);
fprintf(fout,"%d\n",dist[n]);
} else {
precalcDistances();
minCost = oo;
GenPerm(0);
fprintf(fout,"%d\n",minCost);
}
return 0;
}