【题解】P3349 [ZJOI2016]小星星 - 子集dp - 容斥
P3349 [ZJOI2016]小星星
声明:本博客所有题解都参照了网络资料或其他博客,仅为博主想加深理解而写,如有疑问欢迎与博主讨论✧。٩(ˊᗜˋ)و✧*。
题目描述
小 \(Y\) 是一个心灵手巧的女孩子,她喜欢手工制作一些小饰品。她有 \(n\) 颗小星星,用 \(m\) 条彩色的细线串了起来,每条细线连着两颗小星星。
有一天她发现,她的饰品被破坏了,很多细线都被拆掉了。这个饰品只剩下了 \(n-1\) 条细线,但通过这些细线,这颗小星星还是被串在一起,也就是这些小星星通过这些细线形成了树。小 \(Y\) 找到了这个饰品的设计图纸,她想知道现在饰品中的小星星对应着原来图纸上的哪些小星星。如果现在饰品中两颗小星星有细线相连,那么要求对应的小星星原来的图纸上也有细线相连。小 \(Y\) 想知道有多少种可能的对应方式。
只有你告诉了她正确的答案,她才会把小饰品做为礼物送给你呢。
输入格式
第一行包含个 \(2\) 正整数 \(n,m\),表示原来的饰品中小星星的个数和细线的条数。接下来 \(m\) 行,每行包含 \(2\) 个正整数 \(u,v\),表示原来的饰品中小星星 \(u\) 和 \(v\) 通过细线连了起来。这里的小星星从 \(1\) 开始标号。保证 \(u≠v\),且每对小星星之间最多只有一条细线相连。接下来 \(n-1\) 行,每行包含个 \(2\) 正整数 \(u,v\),表示现在的饰品中小星星 \(u\) 和 \(v\) 通过细线连了起来。保证这些小星星通过细线可以串在一起。\(n<=17,m<=n*(n-1)/2\)
输出格式
输出共 \(1\) 行,包含一个整数表示可能的对应方式的数量。如果不存在可行的对应方式则输出 \(0\)。
\(Solution\)
从一个暴力开始考虑,设 \(f[i][j][S]\) 表示,以 \(i\) 为树根的子树,当 \(i\) 映射到 \(j\) 上,映射的子集为 \(S\) 时的方案数。
不断合并 \(i\) 的子树来转移,复杂度貌似是 \(O(3^n n)\) 的,但是不知道为什么 \(emm\)
考虑如何优化
如果我们去除 \(S\) 这一维,那么答案中会有一些不合法的情况,但这些不合法的情况仅为多个点映射到同一个点,不会出现两个点之间没有边相连的情况
于是我们可以做容斥,去除掉这些不合法情况
首先枚举删除一个点,在剩下的图上做一遍 \(f[i][j]\) 的 \(dp\),此时找到的状态至少有两个点映射在了同一个点上(因为图上只有 \(n - 1\) 个点,而树上有 \(n\) 个)
我们把这些状态都删除,但是会发现,会重复删除掉三个点,四个点等映射在同一个点上的情况,于是再加上它们
然后就可以做容斥了
\(Code\)
代码其实挺简单的,但是要注意在做 \(dp\) 的时候,把儿子循环放外面,先做儿子的 \(dp\) 再枚举 \(j\),不然会超时(不过这个逻辑应该挺简单,我真是降智了)
#include<bits/stdc++.h>
#define ll long long
#define F(i, x, y) for(int i = x; i <= y; ++i)
using namespace std;
int read();
const int N = 400;
int n, m;
int mp[N][N];
int head[N], cnt, ver[N], next[N];
ll ban[N], f[N][N];
ll ans;
void add(int x, int y)
{
ver[++ cnt] = y, next[cnt] = head[x], head[x] = cnt;
}
void dfs(int x, int fa)
{
F(j, 1, n) f[x][j] = 1;
for(int i = head[x]; i; i = next[i])
if(ver[i] != fa)
{
dfs(ver[i], x);
F(j, 1, n)
if(! ban[j])
{
ll res = 0;
F(v, 1, n) if(! ban[v]) res += f[ver[i]][v] * mp[v][j];
f[x][j] *= res;
}
}
}
int main()
{
n = read(), m = read();
for(int i = 1, u, v; i <= m; ++ i) u = read(), v = read(), mp[u][v] = mp[v][u] = 1;
for(int i = 1, u, v; i <= n - 1; ++ i) u = read(), v = read(), add(u, v), add(v, u);
F(k, 0, (1 << n) - 1) // 可以直接枚举状态,不用按顺序,这样简单很多
{
F(i, 1, n) ban[i] = 0;
F(i, 1, n) F(j, 1, n) f[i][j] = 0;
int size = 0; ll res = 0;
for(int i = 1, s = k; s; ++ i, s >>= 1) if(s & 1) ban[i] = 1, ++ size;
dfs(1, 0);
F(i, 1, n) res += f[1][i];
if(size & 1) ans -= res;
else ans += res;
}
printf("%lld\n", ans);
return 0;
}
int read()
{
int x = 0, f = 1;
char c = getchar();
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;
}
最新文章
- MYSQL 5.0 USING BTREE 错误
- JS 获取客户端操作系统
- 在ssh中利用Solr服务建立的界面化站内搜索---solr2
- java基础--相等
- 图片过大导致OOM
- 使用SQL Server存储ASP.NET Session变量
- Android ViewPager的简单实现
- 解决cocoapods在64位iOS7系统以下的警告问题
- JavaScript 错误 - Throw、Try 和 Catch
- pci 相关资料
- Sencha Extjs4.2 皮肤制作
- JDBC连接SQL Server 2005步骤详解
- 测试对bug如何分析和定位
- 我发起了一个 .Net 开源 跨平台 GUI (界面开发框架)项目 HtmlCore
- 剑指offer编程题Java实现——面试题14调整数组顺序使奇数位于偶数之前
- verilog语法实例学习(8)
- Material Design系列第二篇——Getting Started
- python 实现判断一个用户输入字符串是否是小数的小程序
- 析构方法 deinit
- Azure 和 Linux
热门文章
- coding++:解决Not allowed to load local resource错误-SpringBoot配置虚拟路径
- RabbitMQ AMQP 事务机制
- js原生模拟new 关键字
- docker 本地镜像导入导出 compose安装
- AngularJS插件使用---angular-cookies.js
- 爬虫简介和requests模块
- Aggressive cows(二分法)
- 错误:Attempt to resolve method: [XXX方法] on undefined variable or class name: [XXX类]的解决(IDEA打包jar问题)
- LeetCode 题解 | 面试题57 - II. 和为s的连续正数序列
- MTK Android修改System分区