BZOJ 1226: [SDOI2009]学校食堂Dining 状压DP

[Submit][Status][Discuss]

Description

小F 的学校在城市的一个偏僻角落，所有学生都只好在学校吃饭。学校有一个食堂，虽然简陋，但食堂大厨总能做出让同学们满意的菜肴。当然，不同的人口味也不一定相同，但每个人的口味都可以用一个非负整数表示。由于人手不够，食堂每次只能为一个人做菜。做每道菜所需的时间是和前一道菜有关的，若前一道菜的对应的口味是a，这一道为b，则做这道菜所需的时间为（a or b）-（a and b），而做第一道菜是不需要计算时间的。其中，or 和and 表示整数逐位或运算及逐位与运算，C语言中对应的运算符为“|”和“&”。学生数目相对于这个学校还是比较多的，吃饭做菜往往就会花去不少时间。因此，学校食堂偶尔会不按照大家的排队顺序做菜，以缩短总的进餐时间。虽然同学们能够理解学校食堂的这种做法，不过每个同学还是有一定容忍度的。也就是说，队伍中的第i 个同学，最多允许紧跟他身后的Bi 个人先拿到饭菜。一旦在此之后的任意同学比当前同学先拿到饭，当前同学将会十分愤怒。因此，食堂做菜还得照顾到同学们的情绪。现在，小F 想知道在满足所有人的容忍度这一前提下，自己的学校食堂做完这些菜最少需要多少时间。

Input

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

Output

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

Sample Input

2

5

5 2

4 1

12 0

3 3

2 2

2

5 0

4 0

Sample Output

16

1

HINT

对于第一组数据：同学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。

解法：真心不会，竟然是状压DP，真心跪服。思路是抄hzwer神牛的。f(i, j, k)表示前i-1人都吃过饭，j表示i与i之后7人

的吃饭情况，k表示上一个吃饭的人与i的相对位置。

那么转移大概是这样子的。

如果j&1=1, f(i, j, k) -> f(i+1, j>>1, k-1)

否则枚举j集合内没吃饭的人

f(, j, k) -> f(i, j+bin[l], l)

复杂度是O(n*7*7*2^7)

//BZOJ 1226

#include <bits/stdc++.h>
using namespace std;
const int inf = 0x3f3f3f3f;
#define dp(i,j,k) ge[i][j][k+8]
///dp[i][j][k]表示前i-1人都吃过饭，j表示i与i之后7人的吃饭情况，k表示上一个吃饭的人与i的相对位置

int t[1005], b[1005], ge[1005][257][18];
int T, n, ans;
void init()
{
memset(ge, inf, sizeof(ge));
}
inline int sum(int x, int y)
{
return x==0?0:t[x]^t[y];
}

int main()
{
scanf("%d", &T);
while(T--)
{
init();
scanf("%d", &n);
for(int i=1; i<= n; i++)
{

scanf("%d%d", &t[i],&b[i]);
}
dp(1, 0, -1) = 0;
for(int i = 1; i <= n; i++)
{
for(int j = 0; j < (1<<8); j++)
{
for(int k=-8; k<=7; k++)
{
if(dp(i, j, k) < inf)
{
if(j&1) dp(i+1, j>>1, k-1) = min(dp(i, j, k), dp(i+1, j>>1, k-1));
else
{
int lim = inf;
for(int l = 0; l <= 7; l++)
{
if(!(j&(1<<l)))
{
if(i+l>lim) break;
lim = min(lim, i+l+b[i+l]);
dp(i, j+(1<<l), l) = min(dp(i, j+(1<<l), l), dp(i, j, k)+sum(i+k, i+l));
}
}
}
}
}
}
}
ans = inf;
for(int i = -8; i<= 7; i++)
{
ans = min(ans, dp(n+1, 0, i));
}
printf("%d\n", ans);
}
return 0;
}