●BZOJ 4008 [HNOI2015]亚瑟王

Posted *ZJ

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了●BZOJ 4008 [HNOI2015]亚瑟王相关的知识,希望对你有一定的参考价值。

题链:

http://www.lydsy.com/JudgeOnline/problem.php?id=4008
题解:

概率dp,神仙题
如果我们可以求出每种牌被取到的概率f,那么最后期望造成的伤害也就很好计算了。
定义dp[i][j]表示有j轮游戏在1~i中的某张牌处就结束的概率。
那么此时我们考虑dp[i][j]会怎样对f[i+1]造成贡献:
只剩下了R-j轮游戏进行到了第i+1张牌,怎么计算这种情况下第i+1张牌发动技能的概率g呢?
(令p为其发动技能的概率,并给这R-j轮游戏重新依次编号为1,2,……,R-j)
显然有:g=p+(1-p)*p+(1-p)^2*p+p……+(1-p)^(R-j-1)*p
上式表示重新编号后的在第1轮发动技能的概率+在第2轮发动技能的概率+……+在第R-j轮发动技能的概率。
然而不需要这么麻烦的计算,因为上面的g=1-(1-p)^(R-j),(自己YY为什么是对的吧)
然后把对f[i+1]进行贡献:f[i+1]+=dp[i][j]*g

接下来考虑如何转移dp[i][j]:
1.这R-j轮可以进行到第i+1张牌的机会都没有让其发动技能:
dp[i+1][j]+=dp[i][j]*(1-p)^(R-j)
2.这R-j轮可以进行到第i+1张牌的机会让其发动了一次技能:
dp[i+1][j+1]+=dp[i][j]*(1-(1-p)^(R-j))

然后就是不断转移dp的同时去求出f[]数组。

(真的是神仙题,题解都看了好久,好像第一次遇到这种定义了一个莫名其妙的dp状态去辅助求出另外一个东西从而得出答案的题。。。)


代码:

 

#include<bits/stdc++.h>
using namespace std;
double dp[250][150],p[250],f[250],ans;
int d[250];
int N,R,Case;
double fastpow(double a,int b){
	double ret=1;
	for(;b;a=a*a,b>>=1)
		if(b&1) ret*=a;
	return ret;
}
int main(){
	for(scanf("%d",&Case);Case;Case--){
		scanf("%d%d",&N,&R);
		for(int i=1;i<=N;i++) scanf("%lf%d",&p[i],&d[i]),f[i]=0;
		for(int i=0;i<=N;i++) for(int j=0;j<=R;j++) dp[i][j]=0;
		dp[0][0]=1; ans=0;
		for(int i=0;i<N;i++)
			for(int j=0;j<=R;j++){
				double k=fastpow(1-p[i+1],R-j);
				dp[i+1][j]+=dp[i][j]*k;
				if(j+1<=R){
					dp[i+1][j+1]+=dp[i][j]*(1-k);
					f[i+1]+=dp[i][j]*(1-k);
				}
			}
		for(int i=1;i<=N;i++) ans+=f[i]*d[i];
		printf("%.10lf\n",ans);
	}
	return 0;
}

 

  

 

以上是关于●BZOJ 4008 [HNOI2015]亚瑟王的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ-4008: [HNOI2015]亚瑟王 (概率期望DP)

BZOJ4008 [HNOI2015]亚瑟王

BZOJ 4008 HNOI2015 亚瑟王

bzoj4008 HNOI2015—亚瑟王

BZOJ 4008: [HNOI2015]亚瑟王( dp )

[BZOJ 4008][HNOI2015]亚瑟王(期望Dp)