\(\mathcal{Description}\)

  Link.

  给定一棵 \(n\) 个点的树,其中 \(2|n\),你需要把这些点两两配对,并把每对点间的路径染色。求使得所有边被染色的方案数,对 \(10^9+7\) 取模。

  \(n\le5000\)。

\(\mathcal{Solution}\)

  容斥,令 \(f(S)\) 表示钦定边集 \(S\) 全部为被覆盖的方案数。显然答案为:

\[\sum_{S\subseteq E}(-1)^{|S|}f(S)
\]

  \(S\) 的断边相当于把原树切分为联通块,且块内可以任意配对。令 \(g(n)\) 表示大小为 \(n\) 的块任意配对的方案数,则:

\[g(n)=\begin{cases}
[n=2]&n\le2\\
(n-1)g(n-2)&\text{otherwise}
\end{cases}
\]

  理解上,考虑在 \(g(n-2)\) 的基础上新加两个点。则可以拆开原来的一对点和这两个点配对,方案数 \(\frac{n-2}2\times2\);这两个点亦可直接配对,方案数 \(1\)。

  利用树上 DP 求解,令 \(h(u,i)\) 表示 \(u\) 子树内有 \(i\) 个点与 \(u\) 连通的方案数。树上背包转移:

\[h(u,i)=\sum h(v,j)h(w,i-j)
\]

  特别地,\(h(u,0)\) 表示 \(u\) 与父亲被切断,即 \(u\) 成为独立连通块。这是考虑容斥系数进行转移:

\[h(u,0)=-\sum h(u,i)g(i)
\]

  答案即为 \(-g(root,0)\)。

\(\mathcal{Code}\)

  1. #include <cstdio>
  2. const int MAXN = 5000, MOD = 1e9 + 7;
  3. int n, ecnt, head[MAXN + 5], siz[MAXN + 5];
  4. int g[MAXN + 5], f[MAXN + 5][MAXN + 5];
  5. struct Edge { int to, nxt; } graph[MAXN * 2 + 5];
  6. inline void addeq ( int& a, const int b ) { if ( ( a += b ) >= MOD ) a -= MOD; }
  7. inline void link ( const int s, const int t ) {
  8. graph[++ ecnt] = { t, head[s] };
  9. head[s] = ecnt;
  10. }
  11. inline void solve ( const int u, const int fa ) {
  12. static int tmp[MAXN + 5];
  13. f[u][1] = siz[u] = 1;
  14. for ( int i = head[u], v; i; i = graph[i].nxt ) {
  15. if ( ( v = graph[i].to ) ^ fa ) {
  16. solve ( v, u );
  17. for ( int j = 1; j <= siz[u] + siz[v]; ++ j ) tmp[j] = 0;
  18. for ( int j = 0; j <= siz[v]; ++ j ) {
  19. for ( int k = 1; k <= siz[u]; ++ k ) {
  20. addeq ( tmp[j + k], 1ll * f[v][j] * f[u][k] % MOD );
  21. }
  22. }
  23. for ( int j = 1; j <= siz[u] + siz[v]; ++ j ) f[u][j] = tmp[j];
  24. siz[u] += siz[v];
  25. }
  26. }
  27. for ( int i = 2; i <= siz[u]; i += 2 ) addeq ( f[u][0], 1ll * f[u][i] * g[i] % MOD );
  28. f[u][0] = ( MOD - f[u][0] ) % MOD;
  29. }
  30. int main () {
  31. scanf ( "%d", &n );
  32. for ( int i = 1, u, v; i < n; ++ i ) {
  33. scanf ( "%d %d", &u, &v );
  34. link ( u, v ), link ( v, u );
  35. }
  36. g[0] = 1;
  37. for ( int i = 2; i <= n; i += 2 ) g[i] = ( i - 1ll ) * g[i - 2] % MOD;
  38. solve ( 1, 0 );
  39. printf ( "%d\n", ( MOD - f[1][0] ) % MOD );
  40. return 0;
  41. }

Solution -「ARC 101E」「AT 4352」Ribbons on Tree的更多相关文章

  1. Solution Set -「ARC 107」

    「ARC 107A」Simple Math   Link.   答案为: \[\frac{a(a+1)\cdot b(b+1)\cdot c(c+1)}{8} \] 「ARC 107B」Quadrup ...

  2. 「ARC 139F」Many Xor Optimization Problems【线性做法,踩标】

    「ARC 139F」Many Xor Optimization Problems 对于一个长为 \(n\) 的序列 \(a\),我们记 \(f(a)\) 表示从 \(a\) 中选取若干数,可以得到的最 ...

  3. Solution -「CTS 2019」「洛谷 P5404」氪金手游

    \(\mathcal{Description}\)   Link.   有 \(n\) 张卡牌,第 \(i\) 张的权值 \(w_i\in\{1,2,3\}\),且取值为 \(k\) 的概率正比于 \ ...

  4. 【翻译】西川善司的「实验做出的游戏图形」「GUILTY GEAR Xrd -SIGN-」中实现的「纯卡通动画的实时3D图形」的秘密,后篇

    http://www.4gamer.net/games/216/G021678/20140714079/     连载第2回的本回,  Arc System Works开发的格斗游戏「GUILTY G ...

  5. 「题解」「美团 CodeM 资格赛」跳格子

    目录 「题解」「美团 CodeM 资格赛」跳格子 题目描述 考场思路 思路分析及正解代码 「题解」「美团 CodeM 资格赛」跳格子 今天真的考自闭了... \(T1\) 花了 \(2h\) 都没有搞 ...

  6. Android内存管理(4)*官方教程 含「高效内存的16条策略」 Managing Your App's Memory

    Managing Your App's Memory In this document How Android Manages Memory Sharing Memory Allocating and ...

  7. SSH连接时出现「WARNING: REMOTE HOST IDENTIFICATION HAS CHANGED!」解决办法

    用ssh來操控github,沒想到連線時,出現「WARNING: REMOTE HOST IDENTIFICATION HAS CHANGED!」,後面還有一大串英文,這時當然要向Google大神求助 ...

  8. 「Windows MFC 」「Edit Control」 控件

    「Windows MFC 」「Edit Control」 控件

  9. 「ZJOI2019」&「十二省联考 2019」题解索引

    「ZJOI2019」&「十二省联考 2019」题解索引 「ZJOI2019」 「ZJOI2019」线段树 「ZJOI2019」Minimax 搜索 「十二省联考 2019」 「十二省联考 20 ...

随机推荐

  1. spring boot 使用 AOP 的正确姿势 --- 心得

    1.前言 向spring boot转型,所有的配置基本上是用注解完成 ,以前使用spring MVC 需要写一大堆xml文件来配置. 基本上没什么变化,但是有些地方需要注意: 环绕通知不要使用异常捕获 ...

  2. console.log(a)和console.log(window.a)的区别?

    console.log(window.l); //undefined console.log(l); //Uncaught ReferenceError: l is not defined js对于未 ...

  3. 使用.NET 6开发TodoList应用(27)——实现API的Swagger文档化

    系列导航及源代码 使用.NET 6开发TodoList应用文章索引 需求 在日常开发中,我们需要给前端提供文档化的API接口定义,甚至需要模拟架设一个fake服务用来调试接口字段.或者对于后端开发人员 ...

  4. 【hexo指南】hexo配置ER图流程图时序图插件

    偏技术的文章有时会用到各种图形,一般来说可以做好图然后截图放到文章中就好了,虽然但图片本身也很小,但存一大堆图片占用空间总觉得不是很好. mermaid mermaid官方网站 mermaid支持很多 ...

  5. C# 季节判断

    编写一个控制台应用程序,可根据输入的月份判断所在季节. 代码如下 using System; using System.Collections.Generic; using System.Linq; ...

  6. 【刷题-LeetCode】207. Course Schedule

    Course Schedule There are a total of numCourses courses you have to take, labeled from 0 to numCours ...

  7. thinkpad s5 电源功率不足提示

    相关答案 作者:路灯瓜 链接:https://www.zhihu.com/question/47551448/answer/122578101 来源:知乎 著作权归作者所有.商业转载请联系作者获得授权 ...

  8. java接口概述及特点

    1 package face_09; 2 3 4 5 6 abstract class AbsDemo{ 7 abstract void show1(); 8 abstract void show2( ...

  9. MySQL存储引擎(最全面的概括)

    目录 一:MySQL存储引擎 1.什么是存储引擎? 2.查看存储引擎信息 二:MySQL支持的存储引擎 1.存储引擎 三:innoDB存储引擎 1.特性 2.存储结构 3.优缺点.适用场景 四:MyI ...

  10. pyhon笔记入门

    人生苦短,我用Python 博客园精华区01-15 23:46 (一)认识Python Python背景介绍 Python的格言: Life is short,use python. (人生苦短,我用 ...