蒟蒻君的刷题日记Day11（DP专题T2）：CF53E Dead Ends

蒟蒻一枚

已于 2022-02-10 09:32:25 修改

阅读量487

点赞数 3

分类专栏：每天一道CF，蒟蒻的AKCF计划文章标签：动态规划算法 c++ Codeforces 状态压缩

于 2022-02-08 10:08:18 首次发布

本文链接：https://blog.csdn.net/yueyuedog/article/details/122818824

版权

每天一道CF，蒟蒻的AKCF计划专栏收录该内容

17 篇文章 2 订阅

订阅专栏

这篇博客探讨了一道动态规划题目，涉及n≤10的数据规模。作者指出，直接使用状态压栈DP可能会因为顺序不同导致重复计数，从而引发错误答案。解决方案是加入条件判断，确保叶子节点的最高位为最新加入的。博客提供了代码实现，但未通过洛谷平台的审核，作者尝试用集合语言重写代码。最后，博主分享了优化后的代码，并表示这是一道需要深入理解的题目。

摘要由CSDN通过智能技术生成

解题思路

数据规模 $n\le10$ 的 dp 题，显然为状压 dp。

状态定义： $d p [i] [j]$ 为联通情况为 $i$ ,叶子结点情况为 $j$ 时的方案数。

状态转移方程：和正常的状压差不多，这里不好描述，详见代码。

边界处理： $\in [1,n])$ 。

这道紫题真得有那么简单嘛？WA。

如果直接这样 DP，顺序不同的方案也会被重复计算。如以下情况：

在这里插入图片描述

我们发现是因为不同的根节点导致的，所以我们加入一个判断条件：

$j$ 加入时当且仅当最高位是最新加入的。

22.2.9 Upd

洛谷上没审核通过，兔队让我用集合的语言写/kk
重新写了一下：

定义 $d p [i] [j]$ 表示目前生成树状态为 $i$ 时，其中度数为 $1$ 的点的状态为 $j$ 时的方案数。

对于加边，可以从连到度数为 $1$ 的点和其它点的情况中取最大值，状态转移方程见代码。

代码实现

#include <bits/stdc++.h>
using namespace std;
const int N = 1 << 10;
int s[N], head[15], cnt;
long long dp[N][N];
struct edge {
	int v, nxt;
} e[205];
inline void add(int u, int v) {
    e[++cnt] = {v, head[u]};
    head[u] = cnt;
}
int main() {
    ios :: sync_with_stdio(0);
    int n, m, k;
    cin >> n >> m >> k;
    while (m--) {
        int u, v;
        cin >> u >> v;
        add(u, v);
        add(v, u);
    }
    int S = (1 << n) - 1;
    for (int i = 1; i <= S; ++i) {
        for (int j = 1; j <= n; ++j) {
			if (i & (1 << j - 1)) {
				++s[i];
			}
        }
    }
	for (int i = 1; i <= n; ++i) {
		dp[1 << i - 1][1 << i - 1] = 1;
    }
	for (int i = 1; i <= S; ++i) {
		for (int j = i; j; --j &= i) {
			if (!dp[i][j]) {
				continue;
			}
			for (int p = 1; p <= n; ++p) {
				if (i & (1 << p - 1)) {
					for (int l = head[p]; l; l = e[l].nxt) {
						int v = e[l].v, t;
						if (i & (1 << v - 1)) {
                            continue;
                        }
                        t = s[i] == 1 ? (i | (1 << v - 1)) : (j & ~(1 << p - 1) | (1 << v - 1));
						if (!(t >> v)) {
                            dp[i | (1 << v - 1)][t] += dp[i][j];
                        }
					}
				}
			}
		}
    }
    long long res = 0;
	for (int i = 1; i <= S; ++i) {
		res += (s[i] == k) * dp[S][i];
	}
    cout << res << '\n';
    return 0;
}