[BZOJ1226][SDOI2009]学校食堂Dining-状压DP

学校食堂Dining

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是一种神奇的东西。
它最大的特点是：难想、难写、难调……

---

思路:
首先我们把目光投向数据范围。
恩，Bi <= 7，看来钦定状压DP了。
那么设f[i][j][k]表示i之前的所有人已打完饭，i及其后边7个人打饭的状况为j，上一个打饭的人距离i的距离为k。

考虑如何转移:

若第i个人在j中显示打过了饭，那么f[i][j][k]可以转移到f[i+1][j>>1][k-1]。
j>>1是因为新加入j的人显然没有打过饭，k-1是因为向后移了一个位置，离上一个打饭的人(显然这个人会在在i以及i的后面)更近了一步。

若第i个人在j中显示没有打饭，那么可以在j所表示的8个人中枚举下一个打饭的人，设为第i+l个人，则可以转移到f[i][j|(1 << l)][l]。
最后对所有的f[n+1][0][i]取min即可~

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>

using namespace std;

inline int read()
{
    int x=0;char ch=getchar();
    while(ch<'0' || '9'<ch)ch=getchar();
    while('0'<=ch && ch<='9')x=x*10+(ch^48),ch=getchar();
    return x;
}

inline int minn(int a,int b){if(a<b)return a;return b;}

const int N=1009;
const int Inf=1e9;

int b[N],t[N],bit[19];
int n,f[N][1<<8][19];

int main()
{
    bit[0]=1;
    for(int i=1;i<=8;i++)
        bit[i]=bit[i-1]<<1;

    int c=read();
    while(c--)
    {
        n=read();
        for(int i=1;i<=n;i++)
        {
            t[i]=read();
            b[i]=read();
        }

        for(int i=0,e=1<<8;i<=n+1;i++)
            for(int j=0;j<e;j++)
                for(int k=0;k<=15;k++)
                    f[i][j][k]=Inf;

        f[1][0][7]=0;

        for(int i=1,e=1<<8;i<=n;i++)
            for(int j=0;j<e;j++)
                for(int k=0;k<=15;k++)
                    if(f[i][j][k]!=Inf)
                    {
                        if(j&1)
                            f[i+1][j>>1][k-1]=minn(f[i+1][j>>1][k-1],f[i][j][k]);
                        else
                        {
                            int mval=Inf;
                            for(int l=0;l<=7;l++)
                                if(!(j&(1<<l)))
                                {
                                    if(i+l>mval)
                                        break;
                                    mval=minn(mval,i+l+b[i+l]);
                                    f[i][j|(1<<l)][l+8]=minn(f[i][j|(1<<l)][l+8],f[i][j][k]+(i+k-8==0?0:t[i+k-8]^t[i+l]));
                                }
                        }
                    }
        int ans=Inf;

        for(int i=0;i<=8;i++)
            ans=minn(ans,f[n+1][0][i]);
        printf("%d\n",ans);
    }

    return 0;
}