***HDU 1529 - Cashier Employment(差分约束)

题目：

http://acm.hdu.edu.cn/showproblem.php?pid=1529

题意：

给出一天中24个小时，每个小时至少需要的工人数，有m个工人，每个工人有其开始时间，每个工人持续工作8小时。求出完成一天的工作最少需要多少人。

思路：

num[i]为从i时刻开始工作的人数

x[i] 为 i时刻实际工作的人数

r[i]为 i 时刻至少需要工作的人数

Sn表示1~n这段时间开始工作的员工数

每个时刻实际工作人数 x[i] >= 0 , 得 Si - S(i-1) >= 0

每个时刻实际工作人数 x[i] <= num[i] , 得 Si - S(i-1) <= num[i]   S(i-1) - Si >= -num[i]

在 i(24>=i>=8) 时刻可以工作的人数 为 x[i-7]+x[i-6]+x[i-5]+x[i-4]+x[i-3]+x[i-2]+x[i-1]+x[i] >= r[i] , 得 Si - S(i-8) >= r[i]

在i(1<=i<8)时刻可以工作的人数为 x[i-7]+x[i-6]+x[i-5]+x[i-4]+x[i-3]+x[i-2]+x[i-1]+x[i] >= r[i] ，得S(24) + Si - S(i+16) >=r[i] ;

x[1]+x[2]+...+x[24] >= 0 得到  S(24) - S0 >= 0 ;

然后用spfa求最长路

AC.

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>

using namespace std;
const int maxn = 25;
const int inf = 0x3f3f3f3f;
int val[maxn], num[maxn];

int tol, head[maxn];
struct Edge{
    int to, w, next;
}edge[maxn*5];
void addedge(int u, int v, int w)
{
    edge[tol].to = v;
    edge[tol].w = w;
    edge[tol].next = head[u];
    head[u] = tol++;
}

int vis[maxn], dis[maxn], cnt[maxn];
bool spfa(int s)
{
    for(int i = 0; i < 25; ++i) {
        vis[i] = 0;
        cnt[i] = 0;
        dis[i] = -inf;
    }
    queue<int> que;
    vis[s] = 1;
    dis[s] = 0;
    que.push(s);

    while(!que.empty()) {
        int u = que.front(); que.pop();
        if(++cnt[u] > 23) return false;
        vis[u] = 0;
        for(int i = head[u]; ~i; i = edge[i].next) {
            int v = edge[i].to, w = edge[i].w;
            if(dis[v] < dis[u] + w) {
                dis[v] = dis[u] + w;
                if(!vis[v]) {
                    vis[v] = 1;
                    que.push(v);
                }
            }
        }
    }
    return true;
}
void init()
{
    tol = 0;
    memset(head, -1, sizeof(head));
}
int main()
{
    //freopen("in", "r", stdin);
    int T;
    scanf("%d", &T);
    while(T--) {
        for(int i = 1; i <= 24; ++i) {
            scanf("%d", &val[i]);
        }
        int m;
        scanf("%d", &m);
        memset(num, 0, sizeof(num));
        for(int i = 0; i < m; ++i) {
            int c;
            scanf("%d", &c);
            num[c+1]++;
        }

        int ans = inf;
        for(int j = 0; j <= m; ++j) {
            init();
            for(int i = 1; i <= 24; ++i) {
                if(i > 0) {
                    addedge(i-1, i, 0);
                    addedge(i, i-1, -num[i]);
                }
                if(i >= 8) {
                    addedge(i-8, i, val[i]);
                }
                else {
                    addedge(i+16, i, val[i]-j);
                }
            }
            addedge(0, 24, j);
            if(spfa(0)) {
                ans = j;
                break;
            }
        }

        if(ans == inf) printf("No Solution\n");
        else
            printf("%d\n", ans);
    }
    return 0;
}