码迷,mamicode.com
首页 > 其他好文 > 详细

codeforces-148D-Bag of mice-概率DP

时间:2014-07-22 22:36:54      阅读:205      评论:0      收藏:0      [点我收藏+]

标签:blog   os   io   2014   art   for   

dp[x][y]:如今有x个白老鼠,y个黑老鼠,公主赢的概率。

那么:

假设公主直接拿到白老鼠,概率为x/(x+y),公主赢。

假设公主拿到黑老鼠,概率为y/(x+y),那么公主假设想赢,龙必须拿到黑老鼠,概率为(y-1)/(x+y-1);

那么逃跑的老鼠为黑老鼠的概率为(y-2)/(x+y-2),为白老鼠的概率为(x)/(x+y-2);

那么dp[x][y]=x/(x+y)+y/(x+y) * (y-1)/(x+y-1) * ( (y-2)/(x+y-2)  * dp[x][y-3]   +  (x)/(x+y-2)  *  dp[x-1][y-2] );

记忆化深搜也行,直接递推DP也行。

#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stdlib.h>
#include<vector>
#include<cmath>
#include<queue>
#include<set>
#include<stack>
using namespace std;
#define N 100010
#define LL long long
#define INF 0xfffffff
#define maxn 1100
double dp[maxn][maxn];
double dos(int x,int y)
{
    if(x<0||y<0)return 0;
    if(x==0)return 0;
    if(y==0)return 1;
    if(dp[x][y]>-0.5)return dp[x][y];
    double a,b;
    a=1.0*x;
    b=1.0*y;
    dp[x][y]=a/(a+b);
    if(x+y>=3)
    {
        dp[x][y]+=(b*(b-1)/((a+b)*(a+b-1)))*
        (dos(x-1,y-2)*a/(a+b-2)+dos(x,y-3)*(b-2)/(a+b-2));
    }
  //  cout<<x<<" "<<y<<" "<<dp[x][y]<<endl;
    return dp[x][y];
}
int main()
{
    int a,b;
    while(~scanf("%d%d",&a,&b))
    {
        memset(dp,-1,sizeof(dp));
        printf("%.10lf\n",dos(a,b));
    }
    return 0;
}











codeforces-148D-Bag of mice-概率DP,布布扣,bubuko.com

codeforces-148D-Bag of mice-概率DP

标签:blog   os   io   2014   art   for   

原文地址:http://www.cnblogs.com/mengfanrong/p/3860743.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!