HDU4598 Difference（差分约束）

题意：

有一个图，给图上每个顶点都赋一个实数Ai。如果存在一个正整数T满足下面两个条件，这个图就是一个"difference"。

1. |Ai| <= T。

2. 如果点i，j构成的边在图中存在，则 |Ai - Aj| >= T；否则 |Ai - Aj| < T。（"<=>" 代表充要条件）

给出图，问这个图是否是一个"difference"。

 

思路：

结合前两个条件，显然每条边的两个端点值都要一正一负，用dfs给每个顶点打上标记colors（1表示值为正，2表示值为负）。

如果图中存在奇数边的环一定无解，可以在dfs的时候进行剪枝，能够让运行时间快一大半。

然后是遍历每个点对（i，j），有下面四种情况：

1. 边ij存在，且点i为正，则 Ai - Aj >= T

2. 边ij存在，且点i为负，则 Aj - Ai >= T

3. 边ij不存在，且点i为正，则 Ai - Aj <= T - 1

4. 边ij不存在，且点i为负，则 Aj - Ai <= T - 1

仅仅是上面的不等式建的图不一定连通，我们可以增加超级源点0，设A0 = 0，那么就有下面两种情况：

1. 点i为正时，0 <= Ai - A0 <= T - 1

2. 点i为负时，0 <= A0 - Ai <= T - 1

T的值不重要，取超过点数即可。

用上面两组不等式建图，然后spfa判负环即可。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
#include <cmath>
#include <list>
#include <cstdlib>
#include <set>
#include <string>
#include <time.h>

using namespace std;

typedef long long ll;
const int inf = 0x3f3f3f3f;
const int maxn = 1005;
struct edg{
    int v, w, nxt;
}G[90005];
int tot, pre[305];
void add(int u, int v, int w) {
    G[tot].v = v;
    G[tot].w = w;
    G[tot].nxt = pre[u];
    pre[u] = tot++;
}
int g[305][305], flag, n, times[305], colors[305], dis[305];
bool vis[305];
void dfs(int x, int k) {
    if (flag) {
        return;
    }
    colors[x] = k;
    vis[x] = true;
    for (int i = 1; i <= n; ++i) {
        if (g[x][i]) {
            if (!colors[i]) {
                dfs(i, 3 - k);
            } else {
                int t = (k + colors[i]) % 2;
                if (!t) {
                    flag = 1;
                }
            }
        }
        if (flag) {
            return ;
        }
    }
}
bool spfa() {
    for (int i = 1; i <= n; ++i) {
        dis[i] = inf;
        vis[i] = false;
        times[i] = 0;
    }
    times[0] = 1;
    vis[0] = true;
    dis[0] = 0;
    queue<int> que;
    que.push(0);
    while (!que.empty()) {
        int u = que.front();
        que.pop();
        vis[u] = false;
        for (int i = pre[u]; ~i; i = G[i].nxt) {
            int v = G[i].v, w = G[i].w;
            if (dis[u] + w < dis[v]) {
                dis[v] = dis[u] + w;
                if (!vis[v]) {
                    if (++times[v] > n + 1) { // 带上0有n+1个点
                        return false;
                    }
                    vis[v] = true;
                    que.push(v);
                }
            }
        }
    }
    return true;
}
int main(){
    int t, T = 305;
    scanf("%d", &t);
    while (t--) {
        scanf("%d", &n);
        for (int i = 1; i <= n; ++i) {
            for (int j = 1; j <= n; ++j) {
                scanf("%1d", &g[i][j]);
            }
        }
        flag = 0;
        memset(colors, 0, sizeof(colors));
        memset(vis, 0, sizeof(vis));
        for (int i = 1; i <= n; ++i) {
            if (!vis[i]) {
                dfs(i, 1);
            }
            if (flag) break;
        }
        if (flag) {
            puts("No");
            continue;
        }
        tot = 0;
        memset(pre, -1, sizeof(pre));
        for (int i = 1; i <= n; ++i) {
            for (int j = i + 1; j <= n; ++j) {
                if (g[i][j]) {
                    if (colors[i] == 1) {
                        add(i, j, -T);           // Ai - Aj >= T
                    } else {
                        add(j, i, -T);           // Aj - Ai >= T
                    }
                } else {
                    if (colors[i] == 1) {
                        add(j, i, T - 1);       // Ai - Aj <= T - 1
                    } else {
                        add(i, j, T - 1);       // Aj - Ai <= T - 1
                    }
                }
            }
            if (colors[i] == 1) {
                add(0, i, T - 1);               // Ai - A0 <= T - 1
                add(i, 0, 0);                   // Ai - A0 >= 0
            } else {
                add(i, 0, T - 1);               // A0 - Ai <= T - 1
                add(0, i, 0);                   // Ai - A0 <= 0
            }
        }
        puts(spfa() ? "Yes" : "No");
    }
    return 0;
}

这道题只用差分约束，不写判奇环也能过，但时间要慢一些。下面一个是带剪枝的，上面一个是不判奇环的。