HDU-5781 ATM Mechine(概率DP)

题目大意:某个未知整数x等概率的分布在[0,k]中。每次你都可以从这个整数中减去一个任意整数y,如果x>=y,那么x=x-y,操作次数累计加1;否则,将会受到一次错误提示。当错误提示超过w次,将会对你的人生产生影响。现在,你的任务是将x逐步变为0,求最少操作次数的期望值。
题目分析:概率DP求期望。定义状态dp(k,w)表示整数分布在[0,k],错误提示次数上限为w时的最少操作次数的期望。
则dp(k,w)=min(p1*dp(k-y,w)+p2*(y-1,w-1))+1,其中p1、p2分别为k>=y、k<y的概率,p1=(k-y+1)/(k+1)、p2=y/(k+1)。因为你非常聪明,所以你每次都会二分的选择要减掉的整数y。根据题目的数据规模,你最多会操作log2(k)+1次,所以你被错误提示的次数最多log2(k)次。这样,便大大减少了状态数目,使得上述方程能够得以实现。

 

参考代码:

# include<bits/stdc++.h>
using namespace std;

const int N=2000;
const double inf=1e9;

double dp[N+5][13];

void init()
{
	for(int i=0;i<=12;++i)
		dp[0][i]=0;
	for(int i=1;i<=N;++i)
		for(int j=0;j<=12;++j) dp[i][j]=inf;
	for(int i=1;i<=N;++i) for(int j=1;j<=12;++j)
		for(int k=1;k<=i;++k)
			dp[i][j]=min(dp[i][j],1.0*(i-k+1)/(i+1)*dp[i-k][j]+1.0*k/(i+1)*dp[k-1][j-1]+1.0);
}

int main()
{
	init();
	int k,w;
	while(~scanf("%d%d",&k,&w))
	{
		w=min(w,12);
		printf("%.6lf\n",dp[k][w]);
	}
	return 0;
}

  

  

  

posted @ 2016-08-06 09:27  20143605  阅读(221)  评论(0编辑  收藏  举报