刷题巩固-----DAY19（小国王）

状态机还剩两道题，但是难度比较大，涉及到ac自动机模型，我做这个专栏是为了查漏补缺，难度大的题将放到第二轮刷题进阶计划中

题目链接

活动 - AcWing本课程系统讲解常用算法与数据结构的应用方式与技巧。https://www.acwing.com/problem/content/1066/今天开始状压dp的更新，大概会刷三道题，因为一开始学的时候，学的不是很明白，现在多刷两道题补一下。

状压dp有一个基础知识，就是位运算操作，我们先复习一下位运算

这道题首先枚举在一行当中所有合法的状态，再计算任意一行的合法状态，对应多少个合法状态，

最后再枚举行数，国王数，和合法状态的三重循环计算方案数

ac代码

#include<bits/stdc++.h>
using namespace std;

typedef long long LL;

const int N = 12, M = 1 << 10, K = 110;

int n, m;
vector<int> state;
int cnt[M];
vector<int> head[M];
LL f[N][K][M];

int count(int state)
{
    int res = 0;
    for (int i = 0; i < n; i ++ ) res += state >> i & 1;
    return res;
}

int main()
{
    cin >> n >> m;

    for (int i = 0; i < 1 << n; i ++ )
        if (!(i&i>>1))
        {
            state.push_back(i);
            cnt[i] = count(i);
        }

    for (int i = 0; i < state.size(); i ++ )
        for (int j = 0; j < state.size(); j ++ )
        {
            int a = state[i], b = state[j];
            if ((a & b) == 0 && check(a | b))
                head[i].push_back(j);
        }

    f[0][0][0] = 1;
    for (int i = 1; i <= n + 1; i ++ )
        for (int j = 0; j <= m; j ++ )
            for (int a = 0; a < state.size(); a ++ )
                for (int b : head[a])
                {
                    int c = cnt[state[a]];
                    if (j >= c)
                        f[i][j][a] += f[i - 1][j - c][b];
                }

    cout << f[n + 1][m][0] << endl;

    return 0;
}