洛谷 P2157 [SDOI2009]学校食堂Dining【状态压缩】

不喜欢数据结构，决心练习水dp；

洛谷 P2157 [SDOI2009]学校食堂Dining
题目描述

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

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

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

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

输入输出样例
输入样例#1：
2
5
5 2
4 1
12 0
3 3
2 2
2
5 0
4 0
输出样例#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。

【蒟蒻一开始完全没想状态压缩】
巧妙的是数据范围0~7早已暴露一切。
具体状态转移方程见代码备注栏。

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
inline int read(){
    char ch='*';int sum=0;
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch>='0'&&ch<='9')
        sum*=10,sum+=ch-'0',ch=getchar();
    return sum;
}
#if 0
Writers: G.S.M. && Goes
Waking up
#endif
const int N=1005;
const int INF=0x3f3f3f3f;
int n,C,T[N],B[N];
int dp[N][(1<<8)][20];
/*dp[i][j][k]
第i个人以前的都打上了饭
从i开始往后7个人打饭的状态是j
当前最后一个打饭的是i+k

情况1：第i个人已经打上饭 => j&1==true
f[i+1][j>>1][k-1]=min(self,f[i][j][k])
【比较容易看出这两个式子所表达的含义是等价的】

情况2：第i个人还没有打上饭 => j&1==false
 f[i][ j|(1<<h) ][h]=min(self, f[i][j][k]+ time(i+k,i+h) )
【枚举 从i开始到i+B[i]个人谁先打饭】
*/
inline void gswork(){
    n=read();memset(dp,INF,sizeof(dp));
    memset(T,0,sizeof(T));
    memset(B,0,sizeof(B));
    for(int i=1;i<=n;i++)
        T[i]=read(),B[i]=read();
    dp[1][0][7]=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+8]!=INF){
        if(j&1) dp[i+1][j>>1][k+7]=min(
            dp[i+1][j>>1][k+7],dp[i][j][k+8]);
        else{
            int lir=INF;
            for(int h=0;h<=7;h++)
            if(!((j>>h)&1)){
                if(i+h>lir)break;
                lir=min(lir,i+h+B[i+h]);
                dp[i][ j|(1<<h) ][h+8]=min(
                dp[i][j|(1<<h)][h+8],dp[i][j][k+8]+
                (i+k?(T[i+k]^T[i+h]):0));
            }
        }
    }int res=INF;
    for(int k=0;k<=8;k++)res=min(res,dp[n+1][0][k]);
    printf("%d\n",res); return ;
}
int main()
{
    C=read();
    while(C--)gswork();
    return 0;
}