[AtCoder Educational DP Contest] J - Sushi(期望dp)

problem

luogu

现有N(1≤N≤300)N(1 ≤ N ≤ 300)N(1N300) 个盘子,编号为1,2,3,…,N1,2,3,…,N1,2,3,,N

iii个盘中放有 ai(1≤ai≤3)a_i(1≤a_i ≤3)ai(1ai3)个寿司。

接下来每次执行以下操作,直至吃完所有的寿司。

从第 1,2,3,…,N1,2,3,…,N1,2,3,,N 个盘子中任选一个盘子,吃掉其中的一个寿司。若没有寿司则不吃。

若将所有寿司吃完,请问此时操作次数的数学期望是多少?

solution

最直接地设 f(a1,a2,a3,...,an):f(a_1,a_2,a_3,...,a_n):f(a1,a2,a3,...,an):iii 盘还剩 aia_iai 个寿司的期望次数。

那么枚举随机到的盘子,有方程:f(a1,a2,...,an)=1+∑i=1n1nf(a1,a2,...,max⁡(ai−1,0),...,an)f(a_1,a_2,...,a_n)=1+\sum_{i=1}^n\frac{1}{n}f(a_1,a_2,...,\max(a_i-1,0),...,a_n)f(a1,a2,...,an)=1+i=1nn1f(a1,a2,...,max(ai1,0),...,an)

显然这个等式不能构成转移方程,因为存在原地打转(第 iii 盘寿司为空时,就变成了自身转移到自身,状态不变)的情况。

由于随机均匀分布,盘子的位置是不重要的,事实上我们真正关注的只有盘子中剩余寿司的数量。

aia_iai 寿司数量又只有四种取值 0/1/2/30/1/2/30/1/2/3

不妨重新设 f(o,i,j,k):f(o,i,j,k):f(o,i,j,k): 当前还剩下 o/i/j/ko/i/j/ko/i/j/k 个盘子中有 0/1/2/30/1/2/30/1/2/3 个寿司。

则有转移:
f(o,i,j,k)=1+onf(o,i,j,k)+inf(o+1,i−1,j,k)+jnf(o,i+1,j−1,k)+knf(o,i,j+1,k−1)f(o,i,j,k)=1+\frac{o}{n}f(o,i,j,k)+\frac{i}{n}f(o+1,i-1,j,k)+\frac{j}{n}f(o,i+1,j-1,k)+\frac{k}{n}f(o,i,j+1,k-1) f(o,i,j,k)=1+nof(o,i,j,k)+nif(o+1,i1,j,k)+njf(o,i+1,j1,k)+nkf(o,i,j+1,k1)

n−onf(o,i,j,k)=1+inf(o+1,i−1,j,k)+jnf(o,i+1,j−1,k)+knf(o,i,j+1,k−1)\frac{n-o}{n}f(o,i,j,k)=1+\frac{i}{n}f(o+1,i-1,j,k)+\frac{j}{n}f(o,i+1,j-1,k)+\frac{k}{n}f(o,i,j+1,k-1) nnof(o,i,j,k)=1+nif(o+1,i1,j,k)+njf(o,i+1,j1,k)+nkf(o,i,j+1,k1)

f(o,i,j,k)=ni+j+k+ii+j+kf(o+1,i−1,j,k)+ji+j+kf(o,i+1,j−1,k)+ki+j+kf(o,i,j+1,k−1)f(o,i,j,k)=\frac{n}{i+j+k}+\frac{i}{i+j+k}f(o+1,i-1,j,k)+\frac{j}{i+j+k}f(o,i+1,j-1,k)+\frac{k}{i+j+k}f(o,i,j+1,k-1) f(o,i,j,k)=i+j+kn+i+j+kif(o+1,i1,j,k)+i+j+kjf(o,i+1,j1,k)+i+j+kkf(o,i,j+1,k1)

这样就不存在状态相同的死循环转移了。

但是现在是 O(n4)O(n^4)O(n4) 的,需要进一步优化。

不难发现,盘子数量是固定不变的,即 o+i+j+k=no+i+j+k=no+i+j+k=n,所以当我们知道了其中任意三个数,就能推出剩下一个数。

f(i,j,k):f(i,j,k):f(i,j,k): 当前还剩下 i/j/ki/j/ki/j/k 个盘子中有 1/2/31/2/31/2/3 个寿司。
f(i,j,k)=ni+j+k+ii+j+kf(i−1,j,k)+ji+j+kf(i+1,j−1,k)+ki+j+kf(i,j+1,k−1)f(i,j,k)=\frac{n}{i+j+k}+\frac{i}{i+j+k}f(i-1,j,k)+\frac{j}{i+j+k}f(i+1,j-1,k)+\frac{k}{i+j+k}f(i,j+1,k-1) f(i,j,k)=i+j+kn+i+j+kif(i1,j,k)+i+j+kjf(i+1,j1,k)+i+j+kkf(i,j+1,k1)
最后还要注意循环枚举的细节:

  • kkk 只用了 k−1k-1k1,当在 (i,j)(i,j)(i,j) 时会问到 j+1j+1j+1,所以 kkk 要在 jjj 循环的外层。
  • 同理 jjj 会用到同 kkk 下的 j+1j+1j+1,但此时是 i−1i-1i1,所以 jjj 循环要在 iii 的外层。
  • 综上我们确定唯一的循环顺序是 k,j,ik,j,ik,j,i

时间复杂度: O(n3)O(n^3)O(n3)

code

#include <bits/stdc++.h>
using namespace std;
#define maxn 305
double f[maxn][maxn][maxn];
int a[5];
int n;
int main() {scanf( "%d", &n );for( int i = 1, x;i <= n;i ++ ) scanf( "%d", &x ), a[x] ++;for( int k = 0;k <= n;k ++ )for( int j = 0;j <= n;j ++ )for( int i = 0;i <= n;i ++ )if( i or j or k ) {if( i ) f[i][j][k] += f[i - 1][j][k] * i / (i + j + k);if( j ) f[i][j][k] += f[i + 1][j - 1][k] * j / (i + j + k);if( k ) f[i][j][k] += f[i][j + 1][k - 1] * k / (i + j + k);f[i][j][k] += n * 1.0 / (i + j + k);}	printf( "%.10f\n", f[a[1]][a[2]][a[3]] );return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/316321.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

cf570 D. Tree Requests

cf570 D. Tree Requests 题意&#xff1a; 给定一个以 1 为根的 n 个结点的树&#xff0c;每个点上有一个字母&#xff08;a-z&#xff09;。每次询问 a, b 查询以 a 为根的子树内深度为 b 的结点上的字母重新排列之后是否能构成回文串。 题解&#xff1a; 回文串形成条件&…

浅谈C#在网络波动时防重复提交

前几天&#xff0c;公司数据库出现了两条相同的数据&#xff0c;而且时间相同&#xff08;毫秒也相同&#xff09;。排查原因&#xff0c;发现是网络波动造成了重复提交。由于网络波动而重复提交的例子也比较多&#xff1a;网络上&#xff0c;防重复提交的方法也很多&#xff0…

[AtCoder Educational DP Contest] V - Subtree(树形dp + 前缀积/后缀积)

problem luogu 给一棵树&#xff0c;对每一个节点染成黑色或白色。 对于每一个节点&#xff0c;求强制把这个节点染成黑色的情况下&#xff0c;所有的黑色节点组成一个联通块的染色方案数&#xff0c;答案对 MMM 取模。 1≤n≤1e5,2≤M≤1e91\le n\le 1e5,2\le M\le 1e91≤n…

P2634 [国家集训队]聪聪可可(点分治做法)

P2634 [国家集训队]聪聪可可 题意&#xff1a; 一颗n个点的树&#xff0c;问其中两点之间的边上数的和加起来是3的倍数的点对有多少个&#xff1f; 输出这样的点对所占比例 题解&#xff1a; 因为是求三的倍数&#xff0c;我们num来记录%30&#xff0c;1&#xff0c;2的数量…

C#并行编程(4):基于任务的并行

C#中的任务Task在C#编程中&#xff0c;实现并行可以直接使用线程&#xff0c;但使用起来很繁琐&#xff1b;也可以使用线程池&#xff0c;线程池很大程度上简化了线程的使用&#xff0c;但是也有着一些局限&#xff0c;比如我们不知道作业什么时候完成&#xff0c;也取不到作业…

[AtCoder Educational DP Contest] W - Intervals(线段树优化dp)

problem 给定 mmm 条规则形如 (li,ri,ai)(l_i,r_i,a_i)(li​,ri​,ai​)&#xff0c;对于一个 01 串&#xff0c;其分数的定义是&#xff1a;对于第 iii 条规则&#xff0c;若该串在 [li,ri][l_i,r_i][li​,ri​] 中至少有一个 1&#xff0c;则该串的分数增加 aia_iai​。 你…

P2634 [国家集训队]聪聪可可(树形dp)

题意&#xff1a; 一颗n个点的树&#xff0c;问其中两点之间的边上数的和加起来是3的倍数的点对有多少个&#xff1f; 输出这样的点对所占比例 题解&#xff1a; 树形dp求解 因为是求长度为3的倍数&#xff0c;模3的结果只有0/1/2&#xff0c;我们可以单独考虑作为一维。 设…

.net core 注入中的三种模式:Singleton、Scoped 和 Transient

从上篇内容不如题的文章《.net core 并发下的线程安全问题》扩展认识.net core注入中的三种模式&#xff1a;Singleton、Scoped 和 Transient我们都知道在 Startup 的 ConfigureServices 可以注入我们想要的服务&#xff0c;那么在注入的时候有三种模式可以选择&#xff0c;那么…

[AtCoder Educational DP Contest] X - Tower(贪心 + dp)

problem luogu 你有 nnn 个箱子&#xff0c;编号从 1 到 n&#xff0c;每个箱子有三个属性&#xff0c;以第 iii 个箱子为例&#xff0c;分别是重量 wiw_iwi​&#xff0c;承重能力 sis_isi​&#xff0c;价值 viv_ivi​。 你想建一座塔&#xff0c;因此需要将一些箱子堆叠起…

P2634 [国家集训队]聪聪可可(树上启发式合并)

P2634 [国家集训队]聪聪可可&#xff08;树上启发式合并&#xff09; 题意&#xff1a; 一颗n个点的树&#xff0c;问其中两点之间的边上数的和加起来是3的倍数的点对有多少个&#xff1f; 输出这样的点对所占比例 题解&#xff1a; 没有修改&#xff0c;统计边长为3的倍数…

【西安活动】 | 4月20日「拥抱开源,又见.NET:云时代 • 新契机」

云计算日渐兴起&#xff0c;成为提升企业效率和生产力的最终解决方案&#xff0c;而云时代也为软件开发模式带来了翻天覆地的变化。可以说 .NET Core就是这个时代催生的产物。自2016年 .NET Core 1.0 发布以来&#xff0c;其强大的生命力让越来越多技术爱好者对她的未来满怀憧憬…

[AtCoder Educational DP Contest] Y - Grid 2(容斥 + dp)

problem luogu 给一个 HWH\times WHW 的网格&#xff0c;每一步只能向右或向下走&#xff0c;给出一些坐标&#xff0c;这些坐标对应的位置不能经过&#xff0c;求从左上角 (1,1)(1,1)(1,1) 走到右下角 (H,W)(H,W)(H,W) 的方案数&#xff0c;答案对 109710^971097 取模。 1≤…

C#并行编程(5):需要知道的异步

异步与并行的联系大家知道“并行”是利用CPU的多个核心或者多个CPU同时执行不同的任务&#xff0c;我们不关心这些任务之间的依赖关系。但是在我们实际的业务中&#xff0c;很多任务之间是相互影响的&#xff0c;比如统计车间全年产量的运算要依赖于各月产量的统计结果。假如你…

[CodeForces1603D] Artistic Partition(四边形不等式 + 决策单调性优化dp + 分治 + 线性筛 + 数论分块)

problem codeforces 对于给定的正整数 l≤l\leql≤&#xff0c;定义 c(l,r)c(l,r)c(l,r) 为满足下列条件的正整数对 (i,j)(i,j)(i,j) 的数量&#xff1a; l≤i≤j≤rl\leq i\leq j\leq rl≤i≤j≤r&#xff1b;gcd⁡(i,j)≥l\gcd(i,j)\geq lgcd(i,j)≥l。 给定正整数 k≤nk\…

P3714 [BJOI2017]树的难题

P3714 [BJOI2017]树的难题 题意&#xff1a; 给你一棵 n 个点的无根树。 树上的每条边具有颜色。一共有 m 种颜色&#xff0c;编号为 1 到 m&#xff0c;第 i 种颜色的权值为 ci。 对于一条树上的简单路径&#xff0c;路径上经过的所有边按顺序组成一个颜色序列&#xff0c;…

从壹开始 [ Id4 ] 之一║ 授权服务器 IdentityServer4 开篇讲计划书

哈喽大家周四好&#xff01;时间过的很快&#xff0c;现在已经是三月份了&#xff0c;我的 IdentityServer4 教程也拖了一定的时间了&#xff0c;正好最近有精力学新东西了&#xff0c;主要中间被小伙伴要求写一个管理后台&#xff0c;目前1.0已经上线&#xff08;《权限后台系…

P2664 树上游戏

P2664 树上游戏 题意&#xff1a; 给一个长度为 n 的颜色序列&#xff0c;定义 s(i,j) 为 i 到 j 的颜色数量。以及 sumi∑j1ns(i,j)sum_{i}\sum_{j1}^{n}s(i,j)sumi​∑j1n​s(i,j) 现在求所有的sumisum_{i}sumi​ 题解&#xff1a; 待补 代码&#xff1a;

[SDOI2019] 热闹的聚会与尴尬的聚会

problem luogu-P5361 他的联系薄上有 nnn 位好友&#xff0c;他们两两之间或者互相认识&#xff0c;或者互相不认识。 小 Q 希望在周六办一个热闹的聚会&#xff0c;再在周日办一个尴尬的聚会。 一场热闹度为 ppp 的聚会请来了任意多位好友&#xff0c;对于每一位到场的好友…

.NET 泛型,详细介绍

今天的文章是因为再给一个朋友讲这个的时候随手记录下整理出来的。说白了就是把前辈们曾经给我吹过的我又吹了出去。泛型&#xff1a;是C# FrameWork 2.0 时代 加入进来的&#xff0c;可以说对与Net开发人员来说泛型是无处不再的&#xff0c;喜欢看源码的同学&#xff0c;可能会…

P2257 YY的GCD

P2257 YY的GCD 题意&#xff1a; 求 1≤x≤N,1≤y≤M1 \leq x \leq N,1 \leq y \leq M1≤x≤N,1≤y≤M 且gcd(x, y) 为质数的 (x,y) 有多少对。 题解&#xff1a; 莫比乌斯反演 代码&#xff1a; #include <bits/stdc.h> #include <unordered_map> #define…