CF #349 div1 B. World Tour

题目链接：http://codeforces.com/problemset/problem/666/B

大意是给一张有向图，选取四个点，使得走这四个点，任意两个点之间走最短路，总距离最长。

3000个点直接floyd肯定不行，但是注意到任意每条路距离都是1，其实可以枚举所有源点作bfs，求出距离数组。

然后对于每个点求出以这个点为入点和出点距离最大的3个点。再根据这个信息，枚举四个点中的中间两个，再枚举这两个点距离他们最远的那3*3种情况，判断是否有重复，没有重复的话，更新答案。

#include <iostream>
#include <vector>
#include <algorithm>
#include <string>
#include <cstring>
#include <cstdio>
#include <math.h>
#include <queue>
#include <stack>
#include <map>
#include <set>
using namespace std;


const int N=3005;

vector<int>edge[N];
int dis[N][N];
int n;
void bfs(int s) {
    queue< pair<int,int> >que;
    que.push(make_pair(s,0));
    dis[s][s]=0;
    while (!que.empty()) {
        pair<int,int> now=que.front();
        que.pop();
        int u=now.first,d=now.second;
        for (int i=0;i<edge[u].size();i++) {
            int v=edge[u][i];
            if (dis[s][v]>=0) continue;
            dis[s][v]=d+1;
            que.push(make_pair(v,d+1));
        }
    }
}
pair<int,int> in[N][3],out[N][3];
void upd(pair<int,int> x,pair<int,int> *y) {
    for (int i=0;i<3;i++) {
        if (x>y[i]) {
            for (int j=2;j>i;j--) {
                y[j]=y[j-1];
            }
            y[i]=x;
            break;
        }
    }
}
void solve() {
    memset(dis,-1,sizeof dis);
    for (int i=1;i<=n;i++){
        bfs(i);
        for (int j=0;j<3;j++) {
            in[i][j]=make_pair(-1,-1);
            out[i][j]=make_pair(-1,-1);
        }
    }
    for (int i=1;i<=n;i++) {
        for (int j=1;j<=n;j++) {
            if (i==j) continue;
            if (dis[i][j]==-1) continue;
            upd(make_pair(dis[i][j],i),in[j]);
            upd(make_pair(dis[i][j],j),out[i]);
        }
    }
    int best=0;
    int a=-1,b=-1,c=-1,d=-1;
    for (int i=1;i<=n;i++) {
        for (int j=1;j<=n;j++) {
            if (i==j) continue;
            if (dis[i][j]==-1) continue;
            for (int p=0;p<3;p++) {
                if (in[i][p].first==-1) continue;
                int ii=in[i][p].second;
                if (ii==j) continue;
                for (int q=0;q<3;q++) {
                    if (out[j][q].first==-1) continue;
                    int jj=out[j][q].second;
                    if (jj==i||jj==ii) continue;

                    int tot=dis[ii][i]+dis[i][j]+dis[j][jj];
                    if (tot>best) {
                        best=tot;
                        a=ii;b=i;c=j;d=jj;
                    }
                }
            }
        }
    }
    printf("%d %d %d %d
",a,b,c,d);
}
int main () {
    int m;
    scanf("%d %d",&n,&m);
    for (int i=1;i<=m;i++) {
        int u,v;
        scanf("%d %d",&u,&v);
        edge[u].push_back(v);
    }
    solve();
    return 0;
}