hdu 5521 Meeting 2015ACM/ICPC亚洲区沈阳站-重现赛

开始只想到了从1出发和从N出发各求一次最短路，每个点的相遇时间就是从1出发和从0出发最短距离的最大值。后来发现重点是构图==。直接构图会TLE，正解是对每个集合新增一个vertex E_i，for node x in E_i, (x,E_i)=0,(E_i,x)=t。这样就可以保证 for all x,y in E_i, (x,y)=(y,x)=t。因为dijkstra不会往标记过的点走，所以不用担心(x,x)=0的问题。

#include<iostream>
#include<stdio.h>
#include<cstdio>
#include<string>
#include<cmath>
#include<stdlib.h>
#include<algorithm>
#include<string.h>
#include<cstring>
#include<vector>
#include<queue>
#include<map>

using namespace std;
//hdu 5521
const int maxn=2*1e5+10;
const int INF=0x3f3f3f3f;
int T;
int n;
int m;
vector<pair<int,int> >mp[maxn];
int ds[maxn];
int dt[maxn];
vector<int>ans;
bool vis[maxn];
class qnode
{
public:
    int v;
    int c;
    qnode()
    {
        v=0;
        c=0;
    }
    qnode(int a,int b)
    {
        v=a;
        c=b;
    }
    bool operator <(const qnode &r) const
    {
        return c>r.c;
    }
};
void Dijsktra(int d[],int st)
{
    memset(d,0x3f,sizeof(d));
    memset(vis,false,sizeof(vis));
    priority_queue<qnode>que;
    while(!que.empty()) que.pop();
    d[st]=0;
    que.push(qnode(st,0));
    qnode tmp;
    while(!que.empty())
    {
        tmp=que.top();
        que.pop();
        int now=tmp.v;
        if(vis[now]==true) continue;
        vis[now]=true;
        for(int i=0;i<mp[now].size();i++)
        {
            int next=mp[now][i].first;
            int cost=mp[now][i].second;
            if(vis[next]==false&&d[next]>d[now]+cost)
            {
                d[next]=d[now]+cost;
                que.push(qnode(next,d[next]));
            }
        }
    }
//    for(int i=1;i<=n;i++)
//    {
//        cout<<d[i]<<" ";
//    }
//    cout<<endl;
}
int main()
{
    freopen("input.txt","r",stdin);
    scanf("%d",&T);
    for(int ca=1;ca<=T;ca++)
    {
        scanf("%d %d",&n,&m);
        int t=0;
        int s=0;
        memset(ds,0x3f,sizeof(ds));
        memset(dt,0x3f,sizeof(dt));
        memset(mp,0,sizeof(mp));
        ans.clear();
        for(int i=1;i<=m;i++)
        {
            scanf("%d %d",&t,&s);
            for(int j=0;j<s;j++)
            {
                int x=0;
                scanf("%d",&x);
                mp[n+i].push_back(make_pair(x,t));
                mp[x].push_back(make_pair(n+i,0));
            }
        }
        Dijsktra(ds,1);
        Dijsktra(dt,n);
        int res=INF;
        int time[maxn];
        memset(time,0,sizeof(time));
        //cout<<"inf "<<INF<<endl;
        for(int i=1;i<=n;i++)
        {
            time[i]=max(ds[i],dt[i]);
            if(time[i]<INF)
            {
                res=min(res,time[i]);
            }
        }
        for(int i=1;i<=n;i++)
        {
            if(res<INF&&res==time[i])
            {
                ans.push_back(i);
            }
        }
        sort(ans.begin(),ans.end());
        printf("Case #%d: ",ca);
        if(ans.size()!=0)
        {
            printf("%d\n",res);
            for(int i=0;i<ans.size();i++)
            {
                if(i==0) printf("%d",ans[i]);
                else printf(" %d",ans[i]);
            }
            puts("");
        }
        else
        {
            puts("Evil John");
        }
    }
    return 0;
}