当前位置: 首页 > news >正文

题解:P3343 [ZJOI2015] 地震后的幻想乡

题意:给出一个图,无重边自环,边权为 \([0,1]\) 内的随机数,问最小生成树最大边权的期望。

做法:

注意到题目中有一个 hint:\(m\) 个随机变量的 \(k\) 小值期望是 \(\frac{k}{m+1}\),考虑怎么使用。

考虑暴力,因为边权是连续而不是离散的,不能直接枚举,所以一个更好的想法是我们枚举边的排名,这样同样是可以做 kruskal 的,然后再利用这个 hint 直接算出答案。

考虑到答案其实是对 \(\frac{k}{m+1}\) 乘上 \(k\) 这条边连通的概率,所以转化一下,变成对 \(\ge k,k\in [0,m-1]\) 的时候仍然不连通的概率。

那么很自然地枚举 \(1\) 所在的集合 \(S\),设补集为 \(S'\),那么要求 \(S\) 内部连通,\(S'\) 随意。但是我们不好计算直接连通的方案数。

所以我们考虑 dp,\(f_{s,i},g_{s,i}\) 分别代表 \(s\) 这个集合用了 \(i\) 条边,不连通/连通的方案数。考虑如何转移。

既然我们要求 \(g\),那么肯定考虑 \(g\) 的转移,但是没啥好的方法。我们注意到 \(f_{s,i}+g_{s,i}=\binom{|E(s)|}{i}\),所以可以转为求 \(f\)

那么 \(f\) 就类似于我们求答案的时候的思路,直接找 lowbit 所在的集合 \(S\),设补集为 \(S'\),枚举我分给 \(S\)\(x\) 条边,那么贡献是:

\[g_{S,x}\binom{|E(S')|}{i-x} \]

解释一下,前者是 \(S\) 连通的方案数,后者是因为 \(S'\) 不能和 \(S\) 中有连边,所以只用管 \(S'\) 内部的,而 \(S'\) 内部不要求连通,随意选就可以。

按上述柿子计算即可,复杂度 \(O(3^nm^2)\)

代码:

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define lowbit(x) (x & (-x))
const int maxn = 15, N = (1 << 11) + 5;
int n, m, d[N];
double f[N][maxn * maxn], g[N][maxn * maxn], C[maxn * maxn][maxn * maxn];
signed main() {cin >> n >> m;for (int i = 1; i <= m; i++) {int x, y; cin >> x >> y;int nw = ((1 << x - 1) | (1 << y - 1));for (int s = 0; s < (1 << n); s++)if((s & nw) == nw)d[s]++;}C[0][0] = 1;for (int i = 1; i <= m; i++) {C[i][0] = 1;for (int j = 1; j <= m; j++)C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]);}for (int s = 0; s < (1 << n); s++) {for (int j = 0; j <= d[s]; j++) {int lb = lowbit(s);for (int t = (s - 1) & s; t; t = (t - 1) & s) {if(!(t & lb))continue;for (int k = 0; k <= min(j, d[t]); k++)f[s][j] = (f[s][j] + g[t][k] * C[d[s ^ t]][j - k]);}g[s][j] = C[d[s]][j] - f[s][j];//	cout << s << " " << j << " " << g[s][j] << endl;}}double ans = 0;for (int i = 0; i <= m; i++)ans += f[(1 << n) - 1][i] / C[m][i];cout << fixed << setprecision(6) << ans / (m + 1) << endl;return 0;
}
http://www.jsqmd.com/news/21154/

相关文章:

  • 暂存:P14214 [COI 2010] 圆圈 / KOLO
  • 树论大封装(直径+重心+中心)
  • QMPlayer2解析
  • 2025年10月广州单位办公室搬家公司全景解析报告,基于专业测评的技术、性能及市场优势深度分析
  • 书评-谋杀黄昏
  • 徐州信息技术服务管理体系认证渠道口碑榜:聚焦机构资质、服务案例及合规性评估
  • 完整教程:【汽车篇】AI深度学习在汽车零部件外观检测——铝铸件中的应用
  • 附加数据文件失败:操作系统错误 5:“5(拒绝访问。)”。 CREATE DATABASE 失败。无法创建列出的某些文件名
  • 20251024- 使用shell脚本分库定时备份MySQL数据
  • 权威调研榜单:东莞工厂装修公司OP3榜单好评深度解析
  • 【Linux】倒计时和进度条完成
  • 2025年口碑好的FPC离型纸,环氧胶离型纸推荐TOP生产厂家
  • 2025年口碑好的数字地磅,电子汽车衡地磅厂家推荐及选择建议
  • 权威调研榜单:四氟换热器生产厂家TOP3榜单好评深度解析
  • 2025年口碑好的食品级贴体盒,榴莲贴体盒实力源头
  • 2025年热门的魔方智能柜,黑金刚智能柜厂家推荐及选择指南
  • 2025年10月进口艺术漆厂家全景解析报告,基于专业测评的技术、性能及市场优势深度分析
  • 2025 年漆包线厂家最新推荐榜,技术实力与市场口碑深度解析,筛选优质品牌助力采购决策
  • 2025年诚信的液压水渠成型机,全自动水渠成型机厂家最新权威推荐榜
  • 2025 年优质销轴厂家最新推荐榜,技术实力与市场口碑深度解析,聚焦高品质连接解决方案发黑 / 异型 / 非标 / 农机销轴公司推荐
  • 最大公约数 gcd
  • 2025年10月扬州公考面试机构全景解析报告,基于专业测评的技术、性能及市场优势深度分析
  • 2025年耐用的陶瓷纤维异性件,硅酸铝纤维陶瓷纤维实力源头加工
  • View root,dirs,files
  • 2025年正规的按动中性笔,多功能中性笔厂家推荐及采购指南
  • 2025年口碑好的空气能地暖管,德国品牌地暖管厂家最新TOP推荐榜
  • 2025 年企业邮箱供应商品牌最新推荐榜,聚焦技术实力与市场口碑深度解析
  • 2025年上海注册公司/上海代理记账公司最新推荐榜,聚焦企业服务品质与特色业务竞争力深度剖析
  • 2025 年接触角测量仪厂家最新推荐榜:聚焦企业专利技术、品质管控及知名客户合作案例的权威解析
  • 2025 年国内无缝钢管厂家最新推荐榜:20#/Q355 系列 / 合金管优质品牌权威测评及选购指南