[Codeforces 666B] World Tour

[题目链接]

https://codeforces.com/contest/666/problem/B

[算法]

首先 ， 用BFS求出任意两点的最短路径

然后 ， 我们用f[i][0-2]表示从i出发到达的最远三点 ， g[i][0-2]表示到i距离最远的三个点

枚举b和c ， 然后在枚举3 * 3个点对 ， 从中选出最优的a和d即可

时间复杂度 ： O(N^2)

[代码]

#include<bits/stdc++.h>
using namespace std;
#define MAXN 3010
#define MAXM 5010
const int inf = 2e9;

int tot , n , m;
int head[MAXN];
int f[MAXN][],g[MAXN][],dist[MAXN][MAXN];

struct info
{
        int a , b , c , d;
} res;
struct edge
{
        int to , nxt;
} e[MAXM];

template <typename T> inline void chkmax(T &x,T y) { x = max(x,y); }
template <typename T> inline void chkmin(T &x,T y) { x = min(x,y); }
template <typename T> inline void read(T &x)
{
    T f = ; x = ;
    char c = getchar();
    for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
    for (; isdigit(c); c = getchar()) x = (x << ) + (x << ) + c - '';
    x *= f;
}
inline void addedge(int u,int v)
{
        tot++;
        e[tot] = (edge){v,head[u]};
        head[u] = tot;
}

int main()
{

        read(n); read(m);
        for (int i = ; i <= m; i++)
        {
                int u , v;
                read(u); read(v);
                addedge(u,v);
        }
        for (int i = ; i <= n; i++)
        {
                for (int j = ; j <= n; j++) dist[i][j] = inf;
                queue< int > q;
                q.push(i);
                dist[i][i] = ;
                while (!q.empty())
                {
                        int cur = q.front();
                        q.pop();
                        for (int j = head[cur]; j; j = e[j].nxt)
                        {
                                int v = e[j].to;
                                if (dist[i][cur] +  < dist[i][v])
                                {
                                        dist[i][v] = dist[i][cur] + ;
                                        q.push(v);
                                }
                        }
                }
        }
        for (int i = ; i <= n; i++)
        {
                for (int j = ; j <= n; j++)
                {
                        if (i == j || dist[i][j] == inf) continue;
                        int pos = j;
                        for (int k = ; k < ; k++)
                        {
                                if (f[i][k] ==  || dist[i][pos] > dist[i][f[i][k]])
                                        swap(pos,f[i][k]);
                        }
                }
                for (int j = ; j <= n; j++)
                {
                        if (i == j || dist[j][i] == inf) continue;
                        int pos = j;
                        for (int k = ; k < ; k++)
                        {
                                if (g[i][k] ==  || dist[pos][i] > dist[g[i][k]][i])
                                        swap(pos,g[i][k]);
                        }
                }
        }
        int ans = ;
        for (int b = ; b <= n; b++)
        {
                for (int c = ; c <= n; c++)
                {
                        if (b == c || dist[b][c] == inf) continue;
                        for (int t1 = ; t1 < ; t1++)
                        {
                                int a = g[b][t1];
                                if (a ==  || a == c) continue;
                                for (int t2 = ; t2 < ; t2++)
                                {
                                        int d = f[c][t2];
                                        if (d ==  || d == b || d == a) continue;
                                        if (dist[a][b] + dist[b][c] + dist[c][d] > ans)
                                        {
                                                ans = dist[a][b] + dist[b][c] + dist[c][d];
                                                res = (info){a,b,c,d};
                                        }
                                }
                        }
                }
        }
        printf("%d %d %d %d\n",res.a,res.b,res.c,res.d);

        return ;

}