bzoj 1226 [SDOI2009]学校食堂Dining（状压DP）

 

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

       设f[i][j][k]表示前i-1个人都吃了饭，i..i+7号人吃饭状况为j，最后一个吃饭的与i相距k

       如果j&1==1，则第i个人已经吃饭可以转移到f[i+1][j>>1][k-1]

       否则，枚举j内还没有吃饭的第l个人，有转移f[i][j^(1<<i)][l]=f[i][j][k]+(a[i+l]|a[i+k])-(a[i+l]&a[i+k])

       另外x|y-x&y=x^y，可以用集合理解理解，并集-交集=两个集合所特有的之并，与xor的定义不谋而合 :)。

【代码】

#include<cstdio>
#include<cstring>
#include<iostream>
#define FOR(a,b,c) for(int a=(b);a<=(c);a++)
#define f(a,b,c) (g[a][b][c+8])
using  namespace std;

const int N = 1e3+10;
const int INF = 1e9;

int g[N][1<<8][16];
int a[N],b[N],bin[20],T,n;

int calc(int i,int j) {
    if(!i) return 0;
    else return (a[i]|a[j])-(a[i]&a[j]);
}

int main() {
    bin[0]=1;
    FOR(i,1,10) bin[i]=bin[i-1]<<1;
    scanf("%d",&T);
    while(T--) {
        scanf("%d",&n);
        FOR(i,1,n) scanf("%d%d",&a[i],&b[i]);
        FOR(i,1,n+1) FOR(j,0,bin[8]-1) FOR(k,-8,7) f(i,j,k)=INF;
        f(1,0,-1)=0;
        FOR(i,1,n) FOR(j,0,bin[8]-1) FOR(k,-8,7)
            if(f(i,j,k)<INF) {
                if(j&1) f(i+1,j>>1,k-1)=min(f(i+1,j>>1,k-1),f(i,j,k));
                else {
                    int r=INF;
                    FOR(l,0,7) if((j&bin[l])==0) {
                        if(i+l>r) break;
                        r=min(r,i+l+b[i+l]);
                        f(i,j^bin[l],l)=min(f(i,j^bin[l],l),f(i,j,k)+calc(i+k,i+l));
                    }
                }
            }
        int ans=INF;
        FOR(k,-8,-1) ans=min(ans,f(n+1,0,k));
        printf("%d
",ans);
    }
    return 0;
}