【BZOJ1226】[SDOI2009] 学校食堂
题目描述
小F 的学校在城市的一个偏僻角落，所有学生都只好在学校吃饭。学校有一个食堂，虽然简陋，但食堂大厨总能做出让同学们满意的菜肴。当然，不同的人口味也不一定相同，但每个人的口味都可以用一个非负整数表示。 由于人手不够，食堂每次只能为一个人做菜。做每道菜所需的时间是和前一道菜有关的，若前一道菜的对应的口味是a，这一道为b，则做这道菜所需的时间为（a or b）-（a and b），而做第一道菜是不需要计算时间的。其中，or 和and 表示整数逐位或运算及逐位与运算，C语言中对应的运算符为“|”和“&”。
学生数目相对于这个学校还是比较多的，吃饭做菜往往就会花去不少时间。因此，学校食堂偶尔会不按照大家的排队顺序做菜，以缩短总的进餐时间。
虽然同学们能够理解学校食堂的这种做法，不过每个同学还是有一定容忍度的。也就是说，队伍中的第i 个同学，最多允许紧跟他身后的Bi 个人先拿到饭菜。一旦在此之后的任意同学比当前同学先拿到饭，当前同学将会十分愤怒。因此，食堂做菜还得照顾到同学们的情绪。 现在，小F 想知道在满足所有人的容忍度这一前提下，自己的学校食堂做完这些菜最少需要多少时间。
输入输出格式
输入格式：
第一行包含一个正整数C，表示测试点的数据组数。 每组数据的第一行包含一个正整数N，表示同学数。 每组数据的第二行起共N行，每行包含两个用空格分隔的非负整数Ti和Bi，表示按队伍顺序从前往后的每个同学所需的菜的口味和这个同学的忍受度。 每组数据之间没有多余空行。
输出格式：
包含C行，每行一个整数，表示对应数据中食堂完成所有菜所需的最少时间。
输入输出样例
说明
对于第一组数据：
同学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呢，一开始真的一直想不到状态转移的方法。
（别告诉我你完全没有想到状态压缩，B[i]才7诶！）
一个人既和前面的人有关系，又和后面的人有关系，真是没法转移了。。。
但是总归是有办法的
我们开一个DP数组
f[i][st][k]
表示此时考虑到第 i 个人，从第i个人开始（包括他自己）后面七个人的吃饭状态是st（1表示吃了，0表示没吃），前一个吃饭的人与i的距离为k（因而k的范围就是 [-8,7]对吧^.^）所用时间的最小值。
所以此时的第i个人不一定吃了哦。
思考下面两种情况：
1. 此时枚举到的i，已经吃过了（st里i对应的位置就为1），那么后面的人不管按什么顺序吃，i都不会生气。
所以这时可以直接更新i+1的状态即 f[i+1][st>>1][k-1]=min(f[i+1][st>>1][k-1],f[i][st][k])
（向后移一个人，st要右移，前一个吃完的人离i+1的距离变成k-1）
（i同学：“你们爱咋咋，反正我已经吃完了，i+1同学，我可以直接更新你的状态了。”）
2.此时枚举到的i，还没有吃饭 ，那么他最多只能忍受后面B[i]个人比他先吃到饭，所以这时候就可以枚举从i开始后面的人谁先吃饭了。
所以此时用lx来枚举后面的人（0<=lx<=7），f[i][ st|tw[lx] ] [lx]=min(f[i][st][k]+Time(i+k,i+lx),f[i][ st|tw[lx] ][lx])
不过还要注意后面的人的B[i+lx]的值是不一样的，我们要开一个变量rx来记录最右能枚举到的lx，并随时更新rx（不然会有很多人生气的！），当lx>rx时， 适时地break就可以了。
（i同学：“我还没吃饭！我后面B[i]个人后的人谁敢先吃饭，我就让你WA”）
最后的答案当然就存在 f[n][1][k]里啦 （-8<=k<=0）
B[i]最多只会有7，而n也才1000，复杂度是很优秀的呢。
1 #include<iostream> 2 #include<cstdio> 3 #include<cstdlib> 4 #include<cstring> 5 #include<algorithm> 6 7 #define For(i,a,b) for(register int i=a;i<=b;++i) 8 #define Re register 9 #define inf 0x7f7f7f 10 #define Nm(i) i+8 11 #define Pn putchar('\n') 12 13 using namespace std; 14 const int N=1e3+10,ST=(1<<8)+10; 15 int f[N][ST][20]; 16 int n,m,T[N],B[N],tw[20],nf; 17 18 inline void read(int &v){ 19 v=0; 20 char c=getchar(); 21 while(c<'0'||c>'9')c=getchar(); 22 while(c>='0'&&c<='9')v=v*10+c-'0',c=getchar(); 23 } 24 void write(int x){ 25 if(x>9)write(x/10); 26 int xx=x%10; 27 putchar(xx+'0'); 28 } 29 int TM(int x,int y){ 30 if(x==0)return 0; 31 return T[x]^T[y]; 32 } 33 34 int main(){ 35 freopen("dining.in","r",stdin); 36 freopen("dining.out","w",stdout); 37 int Tx; read(Tx); 38 For(i,0,10)tw[i]=1<<i; 39 while(Tx--){ 40 read(n); 41 For(i,1,n) read(T[i]),read(B[i]); 42 43 memset(f,inf,sizeof(f)); 44 nf=f[0][0][0]; 45 T[0]=B[0]=0; 46 47 f[1][0][Nm(-1)]=0; 48 For(i,1,n) For(st,0,tw[8]-1) For(k,-8,7){ 49 int Fx=f[i][st][Nm(k)]; 50 if(Fx==nf)continue; 51 //cout<<Fx<<endl; 52 if(st&1){ 53 f[i+1][st>>1][Nm(k-1)]=min(f[i+1][st>>1][Nm(k-1)],Fx); 54 }else{ 55 int mxR=nf; 56 For(lx,0,7){ 57 if(st&tw[lx])continue; 58 if(i+lx>mxR)break; 59 mxR=min(mxR,i+lx+B[i+lx]); 60 61 f[i][ st|tw[lx] ][ Nm(lx) ] =min( f[i][ st|tw[lx] ][ Nm(lx) ],Fx+TM(i+k,i+lx) ); 62 63 } 64 } 65 } 66 int fn=nf; 67 For(i,-8,0)fn=min(fn,f[n][1][Nm(i)]); 68 write(fn); Pn; 69 } 70 fclose(stdin); fclose(stdout); 71 return 0; 72 73 }
求资瓷～
