[SDOI2009]学校食堂（状态压缩）

题目描述

小F 的学校在城市的一个偏僻角落，所有学生都只好在学校吃饭。学校有一个食堂，虽然简陋，但食堂大厨总能做出让同学们满意的菜肴。当然，不同的人口味也不一定相同，但每个人的口味都可以用一个非负整数表示。由于人手不够，食堂每次只能为一个人做菜。做每道菜所需的时间是和前一道菜有关的，若前一道菜的对应的口味是a，这一道为b，则做这道菜所需的时间为（a or b）-（a and b），而做第一道菜是不需要计算时间的。其中，or 和and 表示整数逐位或运算及逐位与运算，C语言中对应的运算符为“|”和“&”。

学生数目相对于这个学校还是比较多的，吃饭做菜往往就会花去不少时间。因此，学校食堂偶尔会不按照大家的排队顺序做菜，以缩短总的进餐时间。

虽然同学们能够理解学校食堂的这种做法，不过每个同学还是有一定容忍度的。也就是说，队伍中的第i 个同学，最多允许紧跟他身后的Bi 个人先拿到饭菜。一旦在此之后的任意同学比当前同学先拿到饭，当前同学将会十分愤怒。因此，食堂做菜还得照顾到同学们的情绪。现在，小F 想知道在满足所有人的容忍度这一前提下，自己的学校食堂做完这些菜最少需要多少时间。

输入输出格式

输入格式：

第一行包含一个正整数C，表示测试点的数据组数。每组数据的第一行包含一个正整数N，表示同学数。每组数据的第二行起共N行，每行包含两个用空格分隔的非负整数Ti和Bi，表示按队伍顺序从前往后的每个同学所需的菜的口味和这个同学的忍受度。每组数据之间没有多余空行。

输出格式：

包含C行，每行一个整数，表示对应数据中食堂完成所有菜所需的最少时间。

输入输出样例

输入样例#1：复制

输出样例#1：复制

16
1

说明

对于第一组数据：

同学1允许同学2或同学3在他之前拿到菜；同学2允许同学3在他之前拿到菜；同学3比较小气，他必须比他后面的同学先拿菜……

一种最优的方案是按同学3、同学2、同学1、同学4、同学5做菜，每道菜所需的时间分别是0、8、1、6及1。

【数据规模和约定】

对于30%的数据，满足1 ≤ N ≤ 20。

对于100%的数据，满足1 ≤ N ≤ 1,000，0 ≤ Ti ≤ 1,000，0 ≤ Bi ≤ 7，1 ≤ C ≤ 5。

存在30%的数据，满足0 ≤ Bi ≤ 1。

存在65%的数据，满足0 ≤ Bi ≤ 5。

存在45%的数据，满足0 ≤ Ti ≤ 130。】

状态：

　　f[i][j][k] 表示：前面i-1个人全都打完饭了，第i个人以及这个人后面的七个人是否打过饭的状态（什么？可能不到七个人？后面会处理的！），最后一名打过饭的同学的相对编号是k（-8 ≤ k ≤ 7），即当前最后一名打过饭的同学是i+k。

状态转移方程（我为人人）：

　　如果当前第i个人打过饭了if(j & 1)：

　　f[i+1][j>>1][k-1] = min(f[i][j][k]);（这里不用增加时间，是因为已经打过饭了）

　　如果当前第i个人没打过饭else

　　枚举h = 0...7，只要满足之前同学的容忍度（用limit维护，如下），且第i+h个人没打过饭，就可以更新：

　　f[i][j|(1<<h)[h] = min(f[i][j][k] + t(i+h, i+k));（最后一个打饭的人更新为i+h）

这个dp有点假，是从前往后推的，更像是递推。

转自https://www.cnblogs.com/Lubixiaosi-Zhaocao/p/9877531.html

 1 #include <bits/stdc++.h>
 2 
 3 using namespace std;
 4 const int INF = 0x3f3f3f3f;
 5 const int MAX_N = 1000 + 5;
 6 const int MAX_B = 8;
 7 
 8 int C, N;
 9 int T[MAX_N], B[MAX_N];
10 int f[MAX_N][1<<MAX_B][20];
11 
12 void input()
13 {
14     cin >> N;
15     for (int i = 1; i <= N; i++)
16         scanf("%d%d", T+i, B+i);
17 }
18 
19 void tomin(int& a, int b)
20 {
21     a = min(a, b);
22 }
23 void tomax(int& a, int b)
24 {
25     a = max(a, b);
26 }
27 
28 inline int cal(int inda, int indb)
29 {
30     if (inda == 0 || indb == 0)
31         return 0;
32     return T[inda] ^ T[indb];
33 }
34 
35 void dp()
36 {
37     memset(f, 0x3f, sizeof f);
38     f[1][0][7] = 0;
39     for (int i = 1; i <= N; i++)
40         for (int j = 0; j < (1<<MAX_B); j++)
41             for (int k = -8; k <= 7; k++) if (f[i][j][k+8] != INF) {
42                 if (j & 1) {
43                     tomin(f[i+1][j>>1][k+7], f[i][j][k+8]);
44                 }
45                 else {
46                     int limit = INF;
47                     for (int h = 0; h <= 7; h++) if (!((j>>h)&1)) {
48                         int nxt = i+h;
49                         if (nxt > limit)
50                             break;
51                         tomin(limit, nxt+B[nxt]);
52                         tomin(f[i][j|(1<<h)][h+8], f[i][j][k+8] + cal(nxt, i+k));
53                     }
54                 }
55             }
56 
57     int ans = INF;
58     for (int k = 0; k < 8; k++)
59         tomin(ans, f[N+1][0][k]);
60     cout << ans << endl;
61 }
62 
63 int main()
64 {
65     cin >> C;
66     while (C--) {
67         input();
68         dp();
69     }
70     return 0;
71 }