HDU 6166 && 2017 多校训练：Senior Pan（最短路）

题意：

有一张n个点m条边的有向图，还有一个包含k个点的点集，求出这个点集中任意两点间最短路的最小值

官方题解看不懂。。

如果一条边的两个端点都在这个集合中，就将这条边直接删掉（中间记录下最小值）以后不会再用

之后就可以愉快的SPFA了，len[i]表示到i点的最短路

随便设集合中的一个点x1为起点，令len[x1]=0然后求最短路，跑完之后记录所有的len[xi] (xi>1 && xi∈k)看哪个更小，之后再将集合中的第二个点x2的len[]置为0并加入队列，继续spfa，然后记录所有的len[xi] (xi>2 && xi∈k)看哪个更小，如此操作直到集合中所有点的len[]全为0

不过这样只处理了ai到aj (i<j)的最小值，所以还要初始化len[]倒过来再跑一次spfa，操作和上面一样

例如样例k=3，集合中的点为1,3,5

①删掉所有1,3,5之间的边，ans = min(ans, len(u, v))  (e(u, v)∈k)

②求出1号点到所有点的单源最短路,  然后ans = min(ans, len[3], len[5])

③令len[3]=0，将3号点加入队列继续spfa，ans = min(ans, len[5])

④没必要再加最后5号点了，初始化len数组

⑤求出5号点到所有点的单源最短路,  然后ans = min(ans, len[3], len[1])

⑥令len[3]=0，将3号点加入队列继续spfa，ans = min(ans, len[1])

结束，输出ans

看上去貌似是跑了2*k次spfa其实只跑了2次，因为你并没有清空len[]数组，所以只是相当于手动放缩了k次而已

复杂度O((n+m)log(n+m))

#include<stdio.h>
#include<vector>
#include<string.h>
#include<queue>
#include<algorithm>
using namespace std;
#define LL long long
typedef struct
{
LL v;
LL len;
}Road;
Road now;
vector<Road> G[100005];
queue<LL> q;
LL God[100005], vis[100005], sum[100005];
int main(void)
{
LL T, i, n, m, j, x, y, len, k, bet, cas = 1;
scanf("%lld", &T);
while(T--)
{
scanf("%lld%lld", &n, &m);
for(i=1;i<=n;i++)
G[i].clear();
for(i=1;i<=m;i++)
{
scanf("%lld%lld%lld", &x, &y, &len);
now.len = len, now.v = y;
G[x].push_back(now);
}
memset(God, 0, sizeof(God));
scanf("%lld", &k);
for(i=1;i<=k;i++)
{
scanf("%lld", &x);
God[x] = 1;
}
bet = 100000ll*100000+5;
for(i=1;i<=n;i++)
{
if(!God[i])  continue;
for(j=0;j<G[i].size();j++)
{
y = G[i][j].v;
if(God[y])
{
bet = min(bet, G[i][j].len);
G[i].erase(G[i].begin()+j);
j--;
}
}
}
memset(vis, 0, sizeof(vis));
memset(sum, 62, sizeof(sum));
for(i=1;i<=n;i++)
{
if(God[i])
{
vis[i] = 1;
q.push(i);
bet = min(bet, sum[i]);
sum[i] = 0;
while(q.empty()==0)
{
x = q.front();
q.pop();
vis[x] = 0;
for(j=0;j<G[x].size();j++)
{
now = G[x][j];
if(sum[x]+now.len<sum[now.v])
{
sum[now.v] = sum[x]+now.len;
if(vis[now.v]==0)
{
vis[now.v] = 1;
q.push(now.v);
}
}
}
}
}
}
memset(sum, 62, sizeof(sum));
for(i=n;i>=1;i--)
{
if(God[i])
{
vis[i] = 1;
q.push(i);
bet = min(bet, sum[i]);
sum[i] = 0;
while(q.empty()==0)
{
x = q.front();
q.pop();
vis[x] = 0;
for(j=0;j<G[x].size();j++)
{
now = G[x][j];
if(sum[x]+now.len<sum[now.v])
{
sum[now.v] = sum[x]+now.len;
if(vis[now.v]==0)
{
vis[now.v] = 1;
q.push(now.v);
}
}
}
}
}
}
printf("Case #%lld: %lld\n", cas++, bet);
}
return 0;
}