ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

【概率期望 && 动态规划】 Card Collector

2019-08-10 19:08:33  阅读:265  来源: 互联网

标签:double sum lt Card && pi Collector dp 式子


题目传送门

题目描述:

在这里插入图片描述


Solution

首先知道这道题的数据范围只有20,很显然的想到了状态压缩。

  • dp[S]S10dp[S]表示|S|这个集合到卡牌全部购买的概率期望(1表示买,0表示不买)dp[S]表示∣S∣这个集合到卡牌全部购买的概率期望(1表示买,0表示不买)
    我们能够得到以下很显然的式子:
    dp[S]=iSp[i]dp[S{i}]dp[S] = \sum_{i\in|S| } p[i]*dp[S∪\{i\}]dp[S]=i∈∣S∣∑​p[i]∗dp[S∪{i}]
  • 但是,这个式子很显然会出锅,因为如果i已经出现过,那么就会重复算,我们需要将1和0之间对答案的关系分开考虑
  • 将分开考虑后的式子进行移项,能够得到以下等式:
    (1iSpi)dps=1+iSp[i]dp[i(S1&lt;&lt;i)](1 - \sum_{i\in |S|} p_i)*dp_s = 1+\sum_{i∉|S|}p[i]*dp[i|(S|1&lt;&lt;i)](1−i∈∣S∣∑​pi​)∗dps​=1+i∈/​∣S∣∑​p[i]∗dp[i∣(S∣1<<i)]
  • 这个式子的左边的意思就是如果当前i在集合中,就不用抽,直接乘上抽不到的概率(应该还有一个1-抽倒得总概率),而右边表示当前状态可以由第i张牌抽到的状态转移过来,由于是由重点向起点转移,所以这样显然是对的,这个等式可以用移项推到。
  • 我们只需要在做的过程中处理出(1iSpi)(1 - \sum_{i\in |S|} p_i)(1−∑i∈∣S∣​pi​)以及iSp[i]dp[i(S1&lt;&lt;i)]\sum_{i∉|S|}p[i]*dp[i|(S|1&lt;&lt;i)]∑i∈/​∣S∣​p[i]∗dp[i∣(S∣1<<i)],就可以把左边的除到右边去得到答案即可。
    //ps:E=1pE=\frac1pE=p1​

Code

#include<bits/stdc++.h>
using namespace std;
int n;
double p[100100];
double unp = 0;
double dp[1<<21];
void work(){
	unp = 0.000;
	for (int i=0;i<n;i++)
	  scanf("%lf",&p[i]),unp+=p[i];
	unp = 1-unp;
	dp[(1<<n)-1] = 0;
	for (int i=(1<<n)-2;i>=0;i--){
		double isp = 0,s = 1;
		for (int j=0;j<n;j++)
		  if (i&(1<<j)) isp+=p[j];
		  else s+=p[j]*dp[i|(1<<j)];
		dp[i] = (s)/(1-unp-isp);
	}
	printf("%.4lf\n",dp[0]);
}

int main(){
	while (scanf("%d",&n)!=EOF) work();
	return 0;
}

标签:double,sum,lt,Card,&&,pi,Collector,dp,式子
来源: https://blog.csdn.net/huang_ke_hai/article/details/99083817

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有