\(\mathcal{Description}\)

  Link.

  给定一棵 \(n\) 个点的树,反复随机选取一条边,合并其两端两点,新点编号在两端两点等概率选取。问每个点留到最后的概率。

  \(n\le50\)。

\(\mathcal{Solution}\)

  推荐 @ywy_c_asm博客 owo。

  所有的操作方案数是 \((n-1)!\),我们可以按删边顺序看做一个长度为 \(n-1\) 的序列。对于每个点分别计算答案,把当前要算的点提为根(记为 \(r\)),我们只需要求出 \(r\) 在所有操作序列中存活的概率和除以 \((n-1)!\) 即可。

  令 \(f(u,i)\) 为 \(r\) 已经走到 \(u\),\(u\) 子树内还剩下 \(i\) 条边没删(没加入删边序列),最终 \(u\)(即 \(r\))存活的概率和。显然答案为 \(f(r,n-1)\),边界 \(f(leaf,0)=1\)。


  第一步,考虑儿子 \(v\) 与 \(u\) 合并。相当于需要考虑边 \((u,v)\) 在操作序列中的位置。粗略来说,若 \(r\) 没走到 \(u\),我们并不关心 \(u\) 号结点的生死;而 \(r\) 到 \(u\) 后,\(u\) 就必须存活。

  定义辅助状态 \(g(u,i)\) 表示 \(u\) 子树内以及 \(u\) 的父边还剩下 \(i\) 条边没删,最终 \(u\) 存活的概率和。现在我们要计算 \(g(u,i)\)。

  第一类,\((u,v)\) 保留到 \(r\) 到达 \(u\) 后再删,那么就涉及到 \(u\) 点存活的概率。于是有转移:

\[g(u,i)=\frac{1}2\sum_{v\in son_u\land j\in[0,i)}f(v,j)
\]

  第二类,\((u,v)\) 在 \(r\) 到达 \(u\) 之前就删,那就很随意啦—— \(v\) 子树中已删除了 \(siz_v-1-i\) 条边,我们把 \((u,v)\) 随便插进一个位置就好,即:

\[g(u,i)=(siz_v-i)f(v,i)
\]

  上两类转移贡献之和即为最终的 \(g(u,i)\)。


  考虑合并,始终记住删除的“序列意义”——保留的边(状态第二维)在删边序列的右端,其它的边在删边序列的左端。合并两个删边序列,仍需要保证这一点,那么分别用组合数合并已删除的左端序列和待删除的右端序列即可。下是 @ywy_c_asm 博客的一张图 owo(红色已删除,蓝色待删除):

  答案呼之欲出啦:

\[f'(u,i+j)=\sum_{i,j}\binom{i+j}{i}\binom{(siz_u-1)+siz_v-i-j}{(siz_u-1)-i}f(u,i)g(v,j)
\]

  两个组合数分别对应分配已删和待删的方案数。

  最终,复杂度 \(\mathcal O(n^4)\) 解决了这道毒瘤 DP qwq。

\(\mathcal{Code}\)

#include <cstdio>
#include <cstring> const int MAXN = 50;
int n, ecnt, head[MAXN + 5], siz[MAXN + 5];
double fac[MAXN + 5];
double f[MAXN + 5][MAXN + 5];
double g[MAXN + 5], h[MAXN + 5]; struct Edge { int to, nxt; } graph[MAXN * 2 + 5]; inline void link ( const int s, const int t ) {
graph[++ ecnt] = { t, head[s] };
head[s] = ecnt;
} inline void init () {
fac[0] = 1;
for ( int i = 1; i <= n; ++ i ) fac[i] = fac[i - 1] * i;
} inline double comb ( const int n, const int m ) {
return n < m ? 0 : fac[n] / fac[m] / fac[n - m];
} inline void solve ( const int u, const int fa ) {
f[u][0] = siz[u] = 1;
for ( int i = head[u], v; i; i = graph[i].nxt ) {
if ( ( v = graph[i].to ) ^ fa ) {
solve ( v, u );
for ( int j = 0; j <= siz[v]; ++ j ) {
g[j] = 0;
for ( int k = 0; k < j; ++ k ) g[j] += 0.5 * f[v][k];
g[j] += ( siz[v] - j ) * f[v][j];
}
for ( int j = 0; j <= siz[v] + siz[u]; ++ j ) h[j] = 0;
for ( int j = 0; j < siz[u]; ++ j ) {
for ( int k = 0; k <= siz[v]; ++ k ) {
h[j + k] += f[u][j] * g[k] * comb ( j + k, j )
* comb ( siz[u] + siz[v] - 1 - j - k, siz[u] - 1 - j );
}
}
siz[u] += siz[v];
for ( int j = 0; j <= siz[u]; ++ j ) f[u][j] = h[j];
}
}
} int main () {
scanf ( "%d", &n ), init ();
for ( int i = 1, u, v; i < n; ++ i ) {
scanf ( "%d %d", &u, &v );
link ( u, v ), link ( v, u );
}
for ( int i = 1; i <= n; ++ i ) {
memset ( f, 0, sizeof f );
solve ( i, 0 );
printf ( "%.12f\n", f[i][n - 1] / fac[n - 1] );
}
return 0;
}

Solution -「CF 1060F」Shrinking Tree的更多相关文章

  1. Solution -「CF 1491H」Yuezheng Ling and Dynamic Tree

    \(\mathcal{Description}\)   Link. 做题原因:题目名.   给定一个长度 \(n-1\) 的序列 \(\{a_2,a_3,\cdots,a_n\}\),其描述了一棵 \ ...

  2. Solution -「CF 1375G」Tree Modification

    \(\mathcal{Description}\)   Link.   给定一棵 \(n\) 个结点的树,每次操作选择三个结点 \(a,b,c\),满足 \((a,b),(b,c)\in E\),并令 ...

  3. Solution -「CF 1342E」Placing Rooks

    \(\mathcal{Description}\)   Link.   在一个 \(n\times n\) 的国际象棋棋盘上摆 \(n\) 个车,求满足: 所有格子都可以被攻击到. 恰好存在 \(k\ ...

  4. Solution -「CF 1237E」Balanced Binary Search Trees

    \(\mathcal{Description}\)   Link.   定义棵点权为 \(1\sim n\) 的二叉搜索树 \(T\) 是 好树,当且仅当: 除去最深的所有叶子后,\(T\) 是满的: ...

  5. Solution -「CF 494C」Helping People

    \(\mathcal{Description}\)   Link.   给定序列 \(\{a_n\}\) 和 \(m\) 个操作,第 \(i\) 个操作有 \(p_i\) 的概率将 \([l_i,r_ ...

  6. Solution -「CF 793G」Oleg and Chess

    \(\mathcal{Description}\)   Link.   给一个 \(n\times n\) 的棋盘,其中 \(q\) 个互不重叠的子矩阵被禁止放棋.问最多能放多少个互不能攻击的车.   ...

  7. Solution -「CF 1622F」Quadratic Set

    \(\mathscr{Description}\)   Link.   求 \(S\subseteq\{1,2,\dots,n\}\),使得 \(\prod_{i\in S}i\) 是完全平方数,并最 ...

  8. Solution -「CF 923F」Public Service

    \(\mathscr{Description}\)   Link.   给定两棵含 \(n\) 个结点的树 \(T_1=(V_1,E_1),T_2=(V_2,E_2)\),求一个双射 \(\varph ...

  9. Solution -「CF 923E」Perpetual Subtraction

    \(\mathcal{Description}\)   Link.   有一个整数 \(x\in[0,n]\),初始时以 \(p_i\) 的概率取值 \(i\).进行 \(m\) 轮变换,每次均匀随机 ...

随机推荐

  1. 强化学习实战 | 自定义gym环境之显示字符串

    如果想用强化学习去实现扫雷.2048这种带有数字提示信息的游戏,自然是希望自定义 gym 环境时能把字符显示出来.上网查了很久,没有找到gym自带的图形工具Viewer可以显示字符串的信息,反而是通过 ...

  2. LINUX学习-Nginx+JDK+Tomcat+Mysql搭建JAVA WEB

    一.安装JDK环境 1.先到  https://www.oracle.com/java/technologies/javase-jdk8-downloads.html  下载linux jdk-8u2 ...

  3. SYCOJ304末尾0的个数

    https://oj.shiyancang.cn/Problem/304.html 首先数据范围不可能算出来的,那么就要看数的性质. 0是怎么来的首先我们知道,有一个0,就必然会有一个5和2. n!在 ...

  4. 腾讯 TKE 厉害了!用 eBPF绕过 conntrack 优化K8s Service,性能提升40%

    Kubernetes Service[1] 用于实现集群中业务之间的互相调用和负载均衡,目前社区的实现主要有userspace,iptables和IPVS三种模式.IPVS模式的性能最好,但依然有优化 ...

  5. POJ 1664 放苹果 (递推思想)

    原题链接:http://poj.org/problem?id=1664 思路:苹果m个,盘子n个.假设 f ( m , n ) 代表 m 个苹果,n个盘子有 f ( m , n ) 种放法. 根据 n ...

  6. 1011day-人口普查系统

    1.Browse.jsp <%@ page language="java" contentType="text/html; charset=UTF-8" ...

  7. 解决windows下因为防火墙无法通过go get 下载gin的问题

    使用: go get -u github.com/gin-gonic/gin 出现以下错误: unrecognized import path "gopkg.in/yaml.v2" ...

  8. C++线程基础笔记(一)

    标准写法: #include<iostream> #include<thread> using namespace std; void MyThread() { cout &l ...

  9. jdk11+安装(win)

    jdk11+安装(win) 官网下载 官网下载地址:https://adoptopenjdk.net/index.html 选择合适的版本 安装 运行下载的 MSI 包 下一步 选择安装位置,下一步 ...

  10. JavaWeb后端工程师技能图