搬运大佬的博客:http://blog.csdn.net/calabash_boy/article/details/77487605
题解:
回想最短路算法,首先排除掉N^3的的那个,然后剩下SPFA和Dijkstra跑多次的复杂度比较能接受,这两个其实是差不多的东西,由于边权都是正的,就上Dijkstra吧。
基础版的Dijkstra是单源多汇的,但是本题是多源多汇,但是Dij他是单源的……等等。。。Dij也可以多源呀,只要开一个超源0,用长度为0的边连接到各个起点,在把每个终点用长度为0的边连接到超汇n+1,这样0 - n+1的最短路就是从所有的起点到所有的终点路径中最短的。那么我们要想办法把真正最短答案的起点 分到起点集合中,把真正的最优终点放到终点集合,其他的随便放哪里都行。emmmm随机化算法随机分组。。。期望做4次可以得到正确答案。。。。
官解:
按照点的标号的每个二进制位分组,最多分20次(准确的说是17次)。每次会把某一位不同的点分开到起点和终点集,然后再起点终点互换,再做一次。
正确性在于:
对于任意两个点u和v,u和v是不同的点,必然有至少一个位不同,因此至少有一次他们被分到了各自正确的集合中,得到了正确答案,其他的答案都比他要大。
#include<bits/stdc++.h>
using namespace std;
#define LL long long
const LL INF = 0x3f3f3f3f3f3f3f3fLL;
const int maxn = 100005;
LL dis[maxn];
int head[maxn], a[maxn], cnt, n, m;
bool vis[maxn];
struct edge
{
int u, v, nxt;
LL w;
}Edge[maxn*2], edg[maxn];
struct node
{
int point;LL distance;
node(){}
node(int _point, LL _distance) {point = _point; distance = _distance;}
bool operator < (const node &other)const
{
return distance > other.distance;
}
};
void add_Edge(int u, int v, int w)
{
Edge[cnt].v = v;
Edge[cnt].w = w;
Edge[cnt].nxt = head[u];
head[u] = cnt++;
}
void Dijkstra(int s)
{
memset(dis, 0x3f, sizeof(dis));
memset(vis, false, sizeof(vis));
dis[s] = 0;
priority_queue<node> q;
q.push(node(s, dis[s]));
while(!q.empty())
{
node now = q.top();
q.pop();
if(vis[now.point])
continue;
vis[now.point] = true;
for(int i = head[now.point]; i != -1; i = Edge[i].nxt)
{
int to = Edge[i].v;
if(dis[to] > dis[now.point] + Edge[i].w)
{
dis[to] = dis[now.point] + Edge[i].w;
q.push(node(to, dis[to]));
}
}
}
}
int main()
{
int T, k;scanf("%d", &T);
for(int kase=1; kase<=T; kase++)
{
scanf("%d%d", &n, &m);
for(int i=1; i<=m; i++) scanf("%d%d%I64d", &edg[i].u, &edg[i].v, &edg[i].w);
scanf("%d", &k);
for(int i=1; i <=k; i++) scanf("%d", &a[i]);
LL ans=INF;
for(int t=0; t<20; t++)
{
cnt=0;
memset(head, -1, sizeof(head));
for(int i=1; i<=m; i++) add_Edge(edg[i].u, edg[i].v, edg[i].w);
for(int i=1; i<=k; i++)
{
if((a[i]>>t)&1) add_Edge(0, a[i], 0);
else add_Edge(a[i], n+1, 0);
}
Dijkstra(0);
ans=min(ans, dis[n+1]);
cnt=0;
memset(head, -1, sizeof(head));
for(int i=1; i<=m; i++) add_Edge(edg[i].u, edg[i].v, edg[i].w);
for(int i=1; i<=k; i++)
{
if(!((a[i]>>t)&1)) add_Edge(0, a[i], 0);
else add_Edge(a[i], n+1, 0);
}
Dijkstra(0);
ans=min(ans, dis[n+1]);
}
printf("Case #%d: %I64d\n", kase, ans);
}
return 0;
}