好题啊!

数树

  • \(\text{opt = 0, 6 pts.}\)

显然答案为 \(y^{n-|E_1∩E_2|}\) 。

  • \(\text{opt = 1, 47 pts.}\)
\[\sum_{E_2}y^{n-|E_1∩E_2|}
\]

考虑该容斥式:\(f(S)=\sum\limits_{T\subset S}\sum\limits_{P\subset T}(-1)^{|T|-|P|}f(P)\) ,则:

\[\sum_{E_2}\sum_{T\subset E_1∩E_2}\sum_{P\subset T}(-1)^{|T|-|P|}y^{n-|P|}
\]

令 \(g(T)\) 表示包含 \(T\) 的 \(E_2\) 个数。

\[\sum_{T\subset E_1}g(T)\sum_{P\subset T}(-1)^{|T|-|P|}y^{n-|P|}
\]
\[\sum_{T\subset E_1}g(T)y^{n-|T|}\sum_{P\subset T}(-y)^{|T|-|P|}
\]
\[\sum_{T\subset E_1}g(T)y^{n-|T|}\sum_{k=0}^{|T|}\binom{|T|}{k}1^k(-y)^{|T|-k}
\]
\[\sum_{T\subset E_1}g(T)y^{n-|T|}(1-y)^{|T|}
\]

考虑一个定理:给定 \(n\) 个点的森林,有 \(k\) 个连通分量,第 \(i\) 个连通分量的大小为 \(a_i\) ,则由 \(\text{prufer}\) 序列得:生成树的个数为 \(n^{k-2}\prod\limits_{i=1}^{k}a_i\) 。

代回原式,记 \(k=n-|T|\) ,表示连通分量个数:

\[\sum_{T\subset E_1}y^k(1-y)^{n-k}n^{k-2}\prod_{i=1}^{k}a_i
\]
\[\frac{(1-y)^n}{n^2}\sum_{T\subset E_1}\prod_{i=1}^{k}\frac{ny}{1-y}a_i
\]

发现一个大小为 \(a_i\) 的连通分量对答案的贡献为 \(Ka_i\) ,其中 \(K=\frac{ny}{1-y}\) 为一个常量。

考虑贡献的组合意义,大小为 \(a_i\) 的连通分量贡献相当于在这个连通分量重选取一个点,产生 \(K\) 的乘积贡献。

据此,可进行 \(dp\) :\(\mathrm{f[u][0/1]}\) 表示在 \(u\) 的子树中,当前连通分量是否已经做出贡献的答案。这里为了方便转移,当前连通分量若已做出贡献,就计入答案。

int f[N][2], K;
void dfs(int u, int fa) {
f[u][0] = 1, f[u][1] = K;
for (auto v: adj[u]) {
if (v == fa) continue;
dfs(v, u);
f[u][1] = (1ll * f[u][0] * f[v][1] + 1ll * f[u][1] * f[v][0] + 1ll * f[u][1] * f[v][1]) % mod;
f[u][0] = (1ll * f[u][0] * f[v][0] + 1ll * f[u][0] * f[v][1]) % mod;
}
}
  • \(\text{opt = 2, 47 pts.}\)
\[\sum_{E_1}\sum_{E_2}y^{n-|E_1∩E_2|}
\]
\[\sum_{E_1}\sum_{E_2}\sum_{T\subset E_1∩E_2}\sum_{P\subset T}(-1)^{|T|-|P|}y^{n-|P|}
\]

这一段类似 \(\text{opt = 2}\) ,直接快进:

\[\sum_{T}g^2(T)y^{n-|T|}(1-y)^{|T|}
\]
\[\sum_{T}y^k(1-y)^{n-k}n^{2(k-2)}\prod_{i=1
}^{k}a_i^2\]
\[\frac{(1-y)^n}{n^4}\sum_{T}\prod_{i=1}^{k}\frac{n^2y}{1-y}a_i^2
\]

换个角度考虑,每个连通分量是无序的,相当于有 \(n\) 个有标号小球扔进 \(k\) 个无标号盒子,盒子不为空。

对于一个装有 \(a\) 个球的盒子,一个生成树对其贡献为 \(\frac{n^2y}{1-y}a^2\) 作为乘积,有 \(a^{a-2}\) 个生成树,因此总乘积贡献为 \(\frac{n^2y}{1-y}a^a\) 。

显然,总方案的 \(EGF\) 等于 \(exp\) 单个盒子的 \(EGF\) ,即:

\[F=\sum\limits_{i=1}^{\infty}\frac{n^2y}{1-y}i^i\frac{x^i}{i!}
\]
\[G=exp(F)
\]

答案即为:

\[\frac{(1-y)^n}{n^4}n![x^n]G
\]

套多项式 \(exp\) 即可。

时间复杂度 \(O(nlogn)\) 。

WC2019的更多相关文章

  1. WC2019滚粗记

    什么?你问WC2019滚粗记在哪里? 抱歉,这篇文章鸽了. 原因? 引用神仙\(yyb\)的话. 恩,想了想还是更一点吧. Day 0 签到海星,我写了个大大的\(Cgod\)有没有人看见啊,然后被广 ...

  2. WC2019游记 && 课件

    WC2019 游记 课件 wc2019.zip_免费高速下载|百度网盘-分享无限制 提取码: un6z day 0 打飞机去广州... 在飞机上刷了爱乐(le)之城, 相当好看... 广二好大! 哈三 ...

  3. WC2019 题目集

    最近写的一些 WC2019 上讲的一些题.还是怕忘了,写点东西记录一下. LOJ2983 「WC2019」数树 题意 本题包含三个问题: 问题 0:已知两棵 \(n\) 个节点的树的形态(两棵树的节点 ...

  4. WC2019 tree

    WC2019唯一一道正常的题,考场上没什么想法,也只拿到了暴力分.搞了一天终于做完了. 前置知识:purfer序,多项式exp或分治FTT. 对于\(type=0\)的,随便维护下,算下联通块即可. ...

  5. WC2019 划水记

    写在前面: 本篇是擅长咕咕咕的\(\text{BLUESKY007}\)同学难得不咕写的游记,将会记录\(WC2019(2019.1.24(Day\ 0)\sim2019.1.30(Day\ 6))\ ...

  6. [WC2019] 数树

    [WC2019] 数树 Zhang_RQ题解(本篇仅概述) 前言 有进步,只做了半天.... 一道具有极强综合性的数数好题! 强大的多合一题目 精确地数学推导和耐心. 有套路又不失心意. 融合了: 算 ...

  7. WC2019 游记

    Day 0 早上奇迹般的六点半起床平常这时候我还没睡呢 早餐在武汉站吃了一碗28的番茄牛肉米线,结果上菜后我把所有非米线的固体(包括番茄和牛肉)全挑出去了 高二大佬:一个愿宰一个愿挨 在高铁上待了四个 ...

  8. WC2019 20天训练

    Day -1 2019.1.2 初步计划: 0x60 图论 std 洛谷提高剩余练习 NOIP2018遗留题解 洛谷省选基础练习 数学: 1.数论 2.组合数学(练习:莫比乌斯反演) 3.概率(练习: ...

  9. PKUWC2019游记&&WC2019游记

    今天好颓,不想写代码了,写写游记 PKUWC2019游记&&WC2019游记 PKUWC2019游记 提前两天就来了中山纪中,考了两天模拟,第一天比较正常,但是可做题只有T3,第二天非 ...

  10. 洛谷 P5206: bzoj 5475: LOJ 2983: [WC2019] 数树

    一道技巧性非常强的计数题,历年WC出得最好(同时可能是比较简单)的题目之一. 题目传送门:洛谷P5206. 题意简述: 给定 \(n, y\). 一张图有 \(|V| = n\) 个点.对于两棵树 \ ...

随机推荐

  1. [FireshellCTF2020]ScreenShooter 1

    此题关键在于理清逻辑,本地将url发送给服务器,服务器请求sereenshooter以后将结果返回 所以应该在服务器查看日志. 发现了PhantomJS 引擎一下 <!DOCTYPE html& ...

  2. .Net Core:Docker无法拉取mcr.microsoft.com相关镜像解决办法

    今天在教同事Docker简单部署Asp.Net Core项目,pull镜像时突然出现下图中的错误: 因为微软在 2018 年五月之后,只会将相关镜像打包发布到 MCR 上.但是 MCR 对国内用户不太 ...

  3. javaweb图书管理系统之账号密码验证登录

    验证账号与密码是否正确功能 一.注册功能 首先,在验证账号与密码是否正确的前提下的,需要先注册一个账号,如果没有账号,就会进不去,也无法验证. 其实,注册功能就是一个添加的功能,仿照我的第一篇文章,往 ...

  4. MySQL中MyISAM和InnoDB引擎的区别

    区别: 1. InnoDB支持事务,MyISAM不支持,对于InnoDB每一条SQL语言都默认封装成事务,自动提交,这样会影响速度,所以最好把多条SQL语言放在begin和commit之间,组成一个事 ...

  5. DRF-认证权限频率

    目录 DRF-认证权限频率 认证 登录接口 认证 权限 作用 使用 频率 作用 使用 认证权限频率+五个接口 模型 视图 序列化器 认证权限频率类 配置文件 路由 DRF-认证权限频率 前后端混合开发 ...

  6. SpringMVC的数据响应-回写数据

    1.直接返回字符串 其他具体代码请访问chilianjie @RequestMapping("/report5") public String save5(HttpServletR ...

  7. Sql递归查询,Sqlserver、Oracle、PG、Mysql

    递归分两种:一种由父项向下级递归,另一种是由子项向上级递归.下面就这两种情况做个简单的处理. 假设有一个表treeview,包含字段 id,parentid,text 分别代表id,上级id,描述字段 ...

  8. 内存之旅——如何提升CMA利用率?

    ​(以下内容来自开发者分享,不代表 OpenHarmony 项目群工作委员会观点)​ 宋远征 李佳伟 OpenAtom OpenHarmony(以下简称"OpenHarmony") ...

  9. 甲骨文严查Java授权,换openJDK要避坑

    背景 外媒The Register报道,甲骨文稽查企业用户,近期开始将把过去看管较松散的Java授权加入. 甲骨文针对标准版Java(Java SE)有2种商业授权.2019年4月甲骨文宣布Java ...

  10. jstl操作session

    1.jstl操作session(添加.删除session中的值)