好题啊!

数树

  • \(\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. [翻译]Service workers:PWA背后的英雄

    原文地址:https://medium.freecodecamp.org/service-workers-the-little-heroes-behind-progressive-web-apps-4 ...

  2. Java/C++实现命令模式---多次撤销和撤回

    某系统需要提供一个命令集合(注:可以使用链表,栈等集合对象实现),用于存储一系列命令对象,并通过该命令集合实现多次undo()和redo()操作,可以使用加法运算来模拟实现.\ 类图: Java代码: ...

  3. Java中JSONArray转换成int[]的办法

    今天写项目的时候要做一个MyBatis的带IN子句的删除,于是用一个整型数组来保存待删除数据的ID 从前端将JSON字符串搞过来之后如何将JSONArray转换成int类型数组就成了个问题 下面是我的 ...

  4. maven安装和配置阿里云镜像(各种详细配置)

    maven安装和详细配置 提示:下面是maven3.6.3版本百度云链接,记住maven3.3以上版本必须安装jdk1.7及以上版本,否则会出错. 链接:https://pan.baidu.com/s ...

  5. 关于#pragma 和 _pragma

    首先要明确 #pragma 和_Pragma 是什么 这两个都是出自于c/c++ 的 ,其中#pragma 是预处理指令(preProcess directive ) ,#pragma是用来向编译器传 ...

  6. mpvue使用scss

    安装scss 安装命令如下,不带版本号可能会导致报错 npm i sass-loader@7.3.1 -D npm i node-sass@4.14.1 -D 然后修改 build 文件夹下的 web ...

  7. vivo 短视频推荐去重服务的设计实践

    一.概述 1.1 业务背景 vivo短视频在视频推荐时需要对用户已经看过的视频进行过滤去重,避免给用户重复推荐同一个视频影响体验.在一次推荐请求处理流程中,会基于用户兴趣进行视频召回,大约召回2000 ...

  8. SpringMVC-设置编码过滤器

    1.接上文->springmvc获取请求参数链接 2.在web.xml配置编码过滤器 <!-- 配置编码过滤器--> <filter> <filter-name&g ...

  9. 【直播回顾】OpenHarmony知识赋能第四期第三课——I2C驱动开发

    3月24日晚上19点,​知识赋能第四期直播的第三节课<OpenHarmony标准系统HDF框架之I2C驱动开发>​,在OpenHarmony开发者成长计划社群内成功举办. 本期课程,由​拓 ...

  10. IP协议/地址(IPv4&IPv6)概要

    IP协议/地址(IPv4&IPv6)概要 IP协议 什么是IP协议 IP是Internet Protocol(网际互连协议)的缩写,是TCP/IP体系中的网络层协议. [1] 协议的特征 无连 ...