AcWing 393. 雇佣收银员（差分约束）

Description

一家超市要每天 24 小时营业，为了满足营业需求，需要雇佣一大批收银员。

已知不同时间段需要的收银员数量不同，为了能够雇佣尽可能少的人员，从而减少成本，这家超市的经理请你来帮忙出谋划策。

经理为你提供了一个各个时间段收银员最小需求数量的清单 R(0),R(1),R(2),…,R(23)。

R(0) 表示午夜00:00 到凌晨 01:00 的最小需求数量，R(1) 表示凌晨 01:00 到凌晨 02:00 的最小需求数量，以此类推。

一共有 N 个合格的申请人申请岗位，第 i 个申请人可以从 ti 时刻开始连续工作 8 小时。

收银员之间不存在替换，一定会完整地工作 8 小时，收银台的数量一定足够。

现在给定你收银员的需求清单，请你计算最少需要雇佣多少名收银员。

Input

第一行包含一个不超过 20 的整数，表示测试数据的组数。

对于每组测试数据，第一行包含 24 个整数，分别表示 R(0),R(1),R(2),…,R(23)(0≤R(i)≤1000)。

第二行包含整数 N(0≤N≤1000)。

接下来 N 行，每行包含一个整数 ti(0≤ti≤23)。

Output

每组数据输出一个结果，每个结果占一行。

如果没有满足需求的安排，输出 No Solution。

Sample Input

1
1 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
5
0
23
22
1
10

Sample Output

1

思路：利用差分约束求最小值：答案的最小值等于答案的最大上界，因此我们需要找出题目中带有‘ ≥ ’符号的不等式。

利用前缀和S[i]表示前i小时即[00:00, i - 1:00]雇佣的收银员人数，因此可以找出S[i]满足以下不等式：
1. 每个时间段雇佣收银员的人数 ≥ 0：S[i] - S[i - 1] ≥ 0 → S[i] ≥ S[i - 1] + 0 → add(i - 1, i, 0);
2. 每个时间段雇佣收银员的人数 ≤ 该时间段收银员的申请人数：
    S[i] - S[i - 1] ≤ r[i]  →  S[i - 1] ≥ S[i] - r[i]  →  add(i, i - 1, -r[i]);
3. 需要满足题目要求，使得需求清单的时间段上有收银员工作：
  对于7 ≤ i ≤ 23 的时间段需要满足：
  S[i + 1] - S[i + 1 - 8] ≥ need[i]  →  S[i + 1] ≥ S[i + 1 - 8] +need[i]  →  add(i + 1 - 8, i + 1, need[i]);
  对于0 ≤ i ≤ 6 的时间段需要满足：
  S[i + 1] + S[24] - S[16 + i + 1] ≥ need[i]  →  S[i + 1] ≥ S[16 + i + 1] -S[24] + need[i] → add(16 + i + 1, i + 1, -S[24] + need[i]);

注意：这里我们并不知道S[24]的取值，但我们可以知道0≤N≤1000，所以0≤S[24]≤1000。因此我们只需要从小到大遍历S[24]的取值，如果出现一个取值满足以上不等式且不发生矛盾，则该值为S[24]的最小取值，即为答案[00:00, 24 - 1:00]时间段雇佣的收银员人数。

4.在遍历时，我们还需要添加S[24] = i 这条不等式，我们可将其转化为S[24] ≥ i，S[24] ≤ i，又因为S[0]没有被用到，其值为0，因此S[24] ≥ S[0] + i，S[24] ≤ S[0] + i  →  add(0, 24, i); add(24, 0, -i);
 

---

注意：我们空出来的S[0]刚好可以作为一个超级源点，因此求spfa时从0出发即可
证明：因为add(i-1, i, 0)，所以有0→1→2→…→50000，因此0可以到达每一个点，满足超级源点的性质

---

代码如下：

#include<bits/stdc++.h>
using namespace std;

#define N 30
#define M 1000

int h[N], e[M], w[M], ne[M], idx;
void add(int a, int b, int c) {
	e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx++;
}

int need[N], r[N];
void build(int s24) {
	memset(h, -1, sizeof h);
	idx = 0;

	for (int i = 1; i <= 24; i++)add(i - 1, i, 0);
	for (int i = 1; i <= 24; i++)add(i, i - 1, -r[i]);
	for (int i = 8; i <= 24; i++)add(i - 8, i, need[i]);
	for (int i = 1; i <= 7; i++)add(16 + i, i, -s24 + need[i]);
	add(0, 24, s24), add(24, 0, -s24);
}

int dist[N], st[N], cnt[N];
bool spfa(int s24) {
	build(s24);
	memset(dist, -0x3f, sizeof dist);
	memset(st, 0, sizeof st);
	memset(cnt, 0, sizeof cnt);

	queue<int>q;
	q.push(0);
	dist[0] = 0;

	while (q.size()) {
		int t = q.front();
		q.pop();
		st[t] = 0;

		for (int i = h[t]; i != -1; i = ne[i]) {
			int j = e[i];
			if (dist[j] < dist[t] + w[i]) {
				dist[j] = dist[t] + w[i];
				cnt[j] = cnt[t] + 1;
				if (cnt[j] == 25)return false;
				if (!st[j]) {
					q.push(j);
					st[j] = 1;
				}
			}
		}
	}
	return true;
}

int main() {

	int t, n;
	scanf("%d", &t);
	while (t--) {
		
		for (int i = 1; i <= 24; i++)scanf("%d", &need[i]);
	
		memset(r, 0, sizeof r);
		scanf("%d", &n);
		while (n--) {
			int a;
			scanf("%d", &a);
			r[a + 1]++;
		}

		bool flag = false;
		for (int i = 0; i < 1001; i++) 
			if (spfa(i)) {
				printf("%d\n", i);
				flag = true;
				break;
			}
		if (!flag)printf("No Solution\n");
	}

	return 0;
}