LightOJ ~ 1074 ~ Extended Traffic （SPFA + DFS判点是否在负环中）

题意：T组测试数据，N个点，编号1~N，然后输入这N个点的一个数值。然后又M条边A~B的权值为（B点数值 - A点数值）^3，然后又Q次询问，问1点到X点的最短路为多少？如果①最短路小于3或②不能到达，输出“？”，否则输出最短路的长度。

思路：建图跑SPFA，如果有负权环，那么就跑DFS把可以通过这个负权环达到的点都标记一下，因为他们的距离可以通过负权环无限缩短。如果是被标记的点或者最短路小于-3或者无法到达，就输出“?”，否则输出最短路的长度。

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 1e5 + 5;
const int INF = 0x3f3f3f3f;
struct Edge
{
    int from, to, dist;       //起点,终点,距离
    Edge(int u, int v, int w):from(u), to(v), dist(w) {}
};
struct SPFA
{
    int n, m;                      //结点数,边数(包括反向弧)
    vector<Edge> edges;            //边表。edges[e]和edges[e^1]互为反向弧
    vector<int> G[MAXN];           //邻接表，G[i][j]表示结点i的第j条边在edges数组中的序号
    bool vis[MAXN];                //是否在队列中
    int d[MAXN];                   //Bellman-Ford
    int p[MAXN];                   //上一条弧
    int cnt[MAXN];                 //进队次数
    int cir[MAXN];                 //是否可以通过负权环无限缩短距离

    void init(int n)
    {
        this->n = n;
        edges.clear();
        for (int i = 0; i <= n; i++) G[i].clear();
    }
    void add_edge(int from, int to, int dist)
    {
        edges.push_back(Edge(from, to, dist));
        m = edges.size();
        G[from].push_back(m - 1);
    }

    void dfs(int u)
    {
        cir[u] = true;
        for (int i = 0; i < G[u].size(); i++)
        {
            Edge& e = edges[G[u][i]];
            if (!cir[e.to]) dfs(e.to);
        }
    }

    bool spfa(int s)
    {
        for (int i = 0; i <= n; i++) d[i] = INF;
        memset(vis, 0, sizeof(vis));
        memset(cnt, 0, sizeof(cnt));
        memset(cir, 0, sizeof(cir));
        d[s] = 0; vis[s] = true;

        queue<int> Q;
        Q.push(s);
        while (!Q.empty())
        {
            int u = Q.front(); Q.pop();
            vis[u] = false;
            for (int i = 0; i < G[u].size(); i++)
            {
                Edge& e = edges[G[u][i]];
                if (d[u] < INF && d[e.to] > d[u] + e.dist)
                {
                    d[e.to] = d[u] + e.dist;
                    p[e.to] = G[u][i];
                    if (!vis[e.to])
                    {
                        Q.push(e.to); vis[e.to] = true;
                        if (++cnt[e.to] > n) { dfs(e.to); return false; } //有负环
                    }
                }
            }
        }
        return true;//没有负环
    }
};
SPFA solve;
int a[MAXN];
int main()
{
    int T, n, m, q, CASE = 1;
    scanf("%d", &T);
    while (T--)
    {
        scanf("%d", &n);
        solve.init(n);
        for (int i = 1; i <= n ; i++) scanf("%d", &a[i]);
        scanf("%d", &m);
        while (m--)
        {
            int u, v;
            scanf("%d%d", &u, &v);
            solve.add_edge(u, v, pow(a[v] - a[u], 3));
        }
        solve.spfa(1);
        scanf("%d", &q);
        printf("Case %d:\n", CASE++);
        while (q--)
        {
            int x;
            scanf("%d", &x);
            int ans = solve.d[x];
            bool flag = solve.cir[x];
            if (ans < 3 || ans == INF || flag) printf("?\n");
            else printf("%d\n", ans);
        }
    }
    return 0;
}
/*
2

5
6 7 8 9 10
6
1 2
2 3
3 4
1 5
5 4
4 5
2
4
5

2
10 10
1
1 2
1
2
*/