bzoj 1226: [SDOI2009]学校食堂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。

Source

状压dp真的是个大坑啊!!!

首先我们设计状态,我们需要知道当前位置i,i附近的没有吃饭的人,上一个吃饭的人;

因为p<=7,那么我们可以状压i附近的人,我们设状态为dp[i][j][k],为i个人之前全部吃了饭,j为i和i后面7个人是否吃饭的状态,k为上一个吃饭的人;

由于p的限制,那么abs(i-k)<=8,那么我们又可以把k换为与i的差值(可以为负),于是空间就开得下了;

然后我们考虑转移:
如果j&1即i吃了饭,那么dp[i][j][k]->dp[i+1][j>>1][k-1];

否则就枚举还没有吃饭的人,dp[i][j][k]+calc(i+k,i+p)->dp[i][j+bin[p]][p];注意放人的时候要避免后面的不合法,所以需要记一个lim,表示最多能放到多远;

//MADE BY QT666
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=100050;
const int xu=8;
int dp[1050][256][20],a[N],b[N],bin[N],T,n;
int calc(int x,int y){
    if(x==0) return 0;
    return (a[x]|a[y])-(a[x]&a[y]);
}
int main(){
    scanf("%d",&T);
    while(T--){
	scanf("%d",&n);
	for(int i=1;i<=n;i++) scanf("%d%d",&a[i],&b[i]);
	bin[0]=1;for(int i=1;i<=8;i++) bin[i]=bin[i-1]<<1;
	memset(dp,127/3,sizeof(dp));
	dp[1][0][-1+xu]=0;
	for(int i=1;i<=n;i++){
	    for(int j=0;j<bin[8];j++)
		for(int k=-8;k<=7;k++){
		    if(dp[i][j][k+xu]==dp[0][0][0]) continue;
		    if(j&1) dp[i+1][j>>1][k-1+xu]=min(dp[i+1][j>>1][k-1+xu],dp[i][j][k+xu]);
		    else{
			int lim=2*n;
			for(int p=0;p<=7;p++){
			    if(!(bin[p]&j)){
				if(i+p>lim) break;
				lim=min(lim,i+p+b[i+p]);
				dp[i][j+bin[p]][p+xu]=min(dp[i][j+bin[p]][p+xu],dp[i][j][k+xu]+calc(i+k,i+p));
			    }
			}
		    }
		}
	}
	int ans=2147483647;
	for(int i=-8;i<=-1;i++) ans=min(ans,dp[n+1][0][i+xu]);
	printf("%d
",ans);
    }
    return 0;
}
原文地址:https://www.cnblogs.com/qt666/p/7682778.html