动态规划专题 - 解题报告 - F

首先我们要找如何让状态转移，前提就是需要知道那些状态能干掉j
用压缩后的状态来表示当前的武器，如果i的状态（初始为weapon[0]）可以杀死敌人j，那么就可以升级武器weapon[i]，weapon[i] |= change[j]，所以武器i的下一轮的武器k的状态就是更新后的weapo[i]。
dp就是如果当前武器i无法杀死敌人j，但是武器升级后weapon[i]可以直接杀死，因为i到weapon[i]到杀死j只有一条选择，或者叫在当前状态i想杀敌人j之间还差一个敌人没杀死过，所以dp[i|(1<<j)] += dp[i]，就是我们的状态转移方程。
为什么不用weapon[i]直指j呢？因为会有多种情况导致出现weapon[i]这样的状态，i只是其中之一，我们统计weapon[i]也需要累加dp[i]，所以其他状态等循环到了再加就行

#include<bits/stdc++.h>
#define FOR(a, b, c) for(int a=b; a<=c; a++)
#define maxn 6005
#define maxm 55
#define hrdg 1000000007
#define inf 2147483647
#define llinf 9223372036854775807
#define ll long long
#define pi acos(-1.0)
#define ls p<<1
#define rs p<<1|1
using namespace std;

int n, T, tot;
char start[18], s[18];
int weapon[1<<17], change[1<<17];
ll dp[1<<17];

inline ll read(){
    char c=getchar();long long x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}
void init(){
    memset(weapon, 0, sizeof(weapon));
    memset(change, 0, sizeof(change));
    memset(s, 0, sizeof(s));
    memset(dp, 0, sizeof(dp));
    n = read();
    scanf("%s", start);
    FOR(i, 0, n-1)
        weapon[0] += (start[i] - '0') << i;
    for (int i = 0; i < n; i++)
    {
        scanf("%s", s);
        for (int j = 0; j < n; j++)
            change[i] |= ((s[j] - '0') << j);     //以上是读入完毕
    }
    //FOR(i, 1, n) printf("%d\n", change[i]);
    tot = (1 << n) - 1;
    for (int i = 1; i <= tot; i++)
    {
        weapon[i] = weapon[0];      //一开始先让所有的状态都能杀的人是初始状态
        for (int j = 0; j < n; j++)
            if (i & (1 << j))           //如果能杀j
                weapon[i] |= change[j];      //就能继承它的武器
    }
    //FOR(i, 1, tot) printf("%d\n", weapon[i]);
}

int main()
{
    T = read();
    FOR(cnt, 1, T)
    {
        init();
        dp[0] = 1;          //谁都不杀只有一种情况
        for (int i = 0; i < tot; i++)
            for (int j = 0; j < n; j++)
                if(!(i & (1 << j)) && (weapon[i] & (1 << j)))        //如果当前砍不死这个闸种，但是起了狙就可以
                    dp[i | (1 << j)] += dp[i];					//状态转移
        //FOR(i, 1, tot) printf("dp[%d] = %lld\n", i, dp[i]);
        printf("Case %d: %lld\n", cnt, dp[tot]);
    }
    return 0;
}