【bzoj1226】【[SDOI2009]学校食堂Dining】状压dp

（上不了p站我要死了，侵权度娘背锅）

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。

状态没定义清楚，转移没想清楚。最后的代码简直是复刻了黄学长的。。。

因为bi小于7，所以考虑状压。设dp[i][j][s]，表示前i个人都吃了、最后一个吃的人与i的相对位置为j，状态s表示自己及后面7个同学是否吃饭。如果该同学还没有吃，则后面的同学如果要先吃就要满足一个条件（各个bi都要满足）。如果该同学吃了，则不必再讨论了，直接转移到下一个同学。

注意j可以是负数，表示前面的。

初值dp[1][-1][0]=0

答案min(dp[n+1][-8~-1][0])，因为我们不清楚结束的到底是哪一个，所以从n+1来取答案会好处理一些。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define f(a,b,c) (dp[a][b+8][c])
using namespace std;

template <typename T>inline void read(T &res){
    T k=1,x=0;char ch=0;
    while(ch<'0'||ch>'9'){if(ch=='-')k=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    res=x*k;
}

const int N=1000+5;
const int oo=0x3f3f3f3f;

int n,t[N],b[N],dp[N][20][1<<8];
int pow[9];

int cal(int x,int y){
    if(x==0) return 0;
    return t[x]^t[y];
}
void solve(){
    read(n);
    for(int i=1;i<=n;i++) read(t[i]),read(b[i]);
    memset(dp,0x3f,sizeof(dp));
    f(1,-1,0)=0;
    for(int i=1;i<=n;i++){
        for(int s=0;s<pow[8];s++){
            for(int j=-8;j<=7;j++){
                if(f(i,j,s)<oo){
                    if(s&1) f(i+1,j-1,s>>1)=min(f(i+1,j-1,s>>1),f(i,j,s));
                    else{
                        int lim=oo;
                        for(int k=0;k<=7;k++){
                            if((s&pow[k])==0){
                                if(i+k>lim) break;
                                lim=min(lim,i+k+b[i+k]);
                                f(i,k,s|pow[k])=min(f(i,k,s|pow[k]),f(i,j,s)+cal(i+j,i+k));
                            }
                        }
                    }

                }
            }
        }
    }
    int ans=oo;
    for(int i=-8;i<=-1;i++)
        ans=min(ans,f(n+1,i,0));
    printf("%d\n",ans);
}
int main(){
    for(int i=0;i<=8;i++) pow[i]=(1<<i);
    int c;
    read(c);
    while(c--) solve();
    return 0;
}