hdu 6166 Senior Pan

地址：http://acm.split.hdu.edu.cn/showproblem.php?pid=6166

题目：

Senior Pan

Time Limit: 12000/6000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 245    Accepted Submission(s): 71

Problem Description

Senior Pan fails in his discrete math exam again. So he asks Master ZKC to give him graph theory problems everyday.
The task is simple : ZKC will give Pan a directed graph every time, and selects some nodes from that graph, you can calculate the minimum distance of every pair of nodes chosen in these nodes and now ZKC only cares about the minimum among them. That is still too hard for poor Pan, so he asks you for help.

 

Input

The first line contains one integer T, represents the number of Test Cases.1≤T≤5.Then T Test Cases, for each Test Cases, the first line contains two integers n,m representing the number of nodes and the number of edges.1≤n,m≤100000
Then m lines follow. Each line contains three integers xi,yi representing an edge, and vi representing its length.1≤xi,yi≤n,1≤vi≤100000
Then one line contains one integer K, the number of nodes that Master Dong selects out.1≤K≤n
The following line contains K unique integers ai, the nodes that Master Dong selects out.1≤ai≤n,ai!=aj

 

Output

For every Test Case, output one integer: the answer

 

Sample Input

1 5 6 1 2 1 2 3 3 3 1 3 2 5 1 2 4 2 4 3 1 3 1 3 5

 

Sample Output

Case #1: 2

 

Source

2017 Multi-University Training Contest - Team 9

 

思路：

　　官方题解看的有点迷，看了半天才懂。

　　就是把k个关键点分为两个点集st，se，两个点集间的最短距离可以通过两次dijkstra或者spfa求出来。

　　所以现在的关键问题是怎么划分点集，使得任意点对都会被划分到不同点集中。

　　一次划分显然是不可能的，所以可进行有限次划分满足上面的要求。

　　考虑到任意不同数对x，y。x，y在二进制下必有一位不相同。

　　所以可以通过二进制位进行20次划分，这样就可以覆盖所有情况。

　　同时还有一种随机做法，就是进行t次random_shuffle用来划分点集。

#include <bits/stdc++.h>

using namespace std;

#define MP make_pair
#define PB push_back
typedef long long LL;
typedef pair<int,int> PII;
const double eps=1e-8;
const double pi=acos(-1.0);
const int K=1e5+7;
const int mod=1e9+7;

vector<PII >mp[K];
vector<int>key,st,se;
LL dis[K];
int n,m,k;

LL spfa(vector<int> &s,vector<int> &e)
{
    queue<int>q;
    fill(dis+1,dis+n+1,1e12);
    for(auto &x:s)
        q.push(x),dis[x]=0;
    while(q.size())
    {
        int u=q.front();q.pop();
        for(auto &v:mp[u])
        if(dis[v.first]>dis[u]+v.second)
        {
            dis[v.first]=dis[u]+v.second;
            q.push(v.first);
        }
    }
    LL ret=1e15;
    for(auto &x:e)
        ret=min(ret,dis[x]);
    return ret;
}
int main(void)
{
    //freopen("in.acm","r",stdin);
    int t,cs=1;cin>>t;
    while(t--)
    {
        memset(mp,0,sizeof mp);
        key.clear();
        LL ans=1e15;
        scanf("%d%d",&n,&m);
        for(int i=1,x,y,z;i<=m;i++)
            scanf("%d%d%d",&x,&y,&z),mp[x].PB(MP(y,z));
        scanf("%d",&k);
        for(int i=1,x;i<=k;i++)
            scanf("%d",&x),key.PB(x);
        for(int i=0;i<=20;i++)
        {
            st.clear(),se.clear();
            for(auto &x:key)
            if(x&(1<<i))
                st.push_back(x);
            else
                se.push_back(x);
            ans=min(ans,min(spfa(st,se),spfa(se,st)));
        }
        printf("Case #%d: %lld
",cs++,ans);
    }
    return 0;
}