[BZOJ1226][SDOI2009]学校食堂(状压DP)

状压DP
f[i][j][k]表示第1个人到第i1个人已经打完饭,第i个人以及后面7个人是否打饭的状态为j,当前最后一个打饭的人的编号为i+kk的范围为87,所以用数组存时要加上8),那么转移为:
如果j&1为真,就表示第i个人已经打完饭,i之后的7个人中,还没打饭的人就再也不会插入到第i个人前面了。所以这时候可以转移到f[i+1][j>>1][k1],即f[i+1][j>>1][k1]=min(f[i+1][j>>1][k1],f[i][j][k]),不需要累积时间(因为在j&1为真的情况下,f[i][j][k]f[i+1][j>>1][k1]的意义是一样的)。
而为什么意义是一样的呢?因为可以看出,最后一个打饭的人的编号为(i+1)+(k1)=i+k,和f[i][j][k]表示的一样。而第i个人也打完了饭,所以满足「第1个人到第i个人已经打完饭」这个条件。而j>>1就是说i之后的第1个人就是i+1之后的第0个人(就是i+1本人),i之后的第2个人就是i+1之后的第1个人,i之后的第3个人就是i+1之后的第2个人,…。这样就可以看出意义一样了。
j&1为假时,是没办法转移到f[i+1]的(因为i+1之前的人还有i没有打完饭)。但是这时候可以把i以及i之后的7个人中选出一个人打饭,也就是枚举h07f[i][j|(1<<h)][h]=min(f[i][j|(1<<h)][h],f[i][j][k]+time(i+k,i+h)),其中time(i,j)表示如果上一个人编号为i,当前的人编号为j,那么做编号为j的人的菜需要的时间。 当然,这个转移需要考虑到忍耐度的问题。这样,在ii之后的7个人,不是每一个还未打饭的人都可以先打饭的。因为编号在他之前的所有未打饭的人的忍耐度必须能忍受这个人在他们之前打饭。所以,在这里用了一个变量r来统计了一下,表示到目前为止的未打饭的人的忍受范围(注意,不是忍耐度,忍受范围是指能忍受在其之前打饭的最大位置)的最小值,对于任何一个人,如果i+h>r,就表示他无法满足编号在他之前的所有人的忍受范围,就不要考虑这个人了。代码实现如下:

lir = INF;
for (h = 0; h <= 7; h++) if (!((j >> h) & 1)) {
    if (i + h > lir) break;
    chkmin(lir, i + h + B[i + h]);
    chkmin(f[i][j | (1 << h)][h + 8], f[i][j][k + 8] +
    (i + k ? (T[i + k] ^ T[i + h]) : 0));
}

其中lir为上面提到的统计变量rchkmin(a,b)为如果b<a则把a赋为b
最后答案即为min(f[n+1][0][k],8<=k<=0)
完整代码:

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
inline int read() {
    int res = 0; bool bo = 0; char c;
    while (((c = getchar()) < '0' || c > '9') && c != '-');
    if (c == '-') bo = 1; else res = c - 48;
    while ((c = getchar()) >= '0' && c <= '9')
        res = (res << 3) + (res << 1) + (c - 48);
    return bo ? ~res + 1 : res;
}
const int N = 1005, INF = 0x3f3f3f3f;
int n, T[N], B[N], f[N][1 << 8][20];
void chkmin(int &a, int b) {a = min(a, b);}
void work() {
    int i, j, k, h, lir;
    n = read(); for (i = 1; i <= n; i++)
        T[i] = read(), B[i] = read();
    memset(f, INF, sizeof(f)); f[1][0][7] = 0;
    for (i = 1; i <= n; i++) for (j = 0; j < (1 << 8); j++)
    for (k = -8; k <= 7; k++) if (f[i][j][k + 8] != INF) {
        if (j & 1) chkmin(f[i + 1][j >> 1][k + 7], f[i][j][k + 8]);
        else {
            lir = INF;
            for (h = 0; h <= 7; h++) if (!((j >> h) & 1)) {
                if (i + h > lir) break;
                chkmin(lir, i + h + B[i + h]);
                chkmin(f[i][j | (1 << h)][h + 8], f[i][j][k + 8] +
                (i + k ? (T[i + k] ^ T[i + h]) : 0));
            }
        }
    }
    int res = INF; for (k = 0; k <= 8; k++)
        res = min(res, f[n + 1][0][k]);
    printf("%d\n", res);
}
int main() {
    int T = read();
    while (T--) work();
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值