Light OJ 1050 - Marbles（概率DP）

题目大意：

一个包裹里有蓝色和红色的弹珠，在这个包裹里有奇数个弹珠，你先走， 你先从背包里随机的拿走一个弹珠，拿走每个弹珠的可能性是一样的。然后Jim从背包里拿走一个蓝色的弹珠，如果没有蓝色的弹珠让Jim拿走，那么Jim赢，如果最终从包裹里移走的是蓝的弹珠，那么你赢，否则Jim赢。给你蓝色和红的背包的里弹珠的个数，求赢的概率。

 

dp[R][B] 代表 R个红球， B个白球的时候赢得概率。当遇到 (R+B)%2 == 0 得时候直接 dp[R][B] = dp[R][B-1]就OK了。

这里采用的是记忆化搜索，所以用起来比较方便。

 

 

 ==============================================================================================================

 

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<map>
using namespace std;
typedef long long LL;
const int INF = 1e9+7;
const int MAXN = 555;
double dp[MAXN][MAXN];

double DFS(int R,int B)
{
    if(dp[R][B] != -1)
        return dp[R][B];
    if(R >= B || B == 0)
        return dp[R][B] = 0;
    if(R == 0)
        return dp[R][B] = 1;
    double k = 1.0/(R + B);
    if( (R+B)%2 )
        dp[R][B] =  1.0*R*k*DFS(R-1,B) + 1.0*B*k*DFS(R, B-1);
    else
        dp[R][B] = DFS(R, B-1);

    return dp[R][B];
}

int main()
{
    int T, cas = 1, n, R, B;
    for(int i=0; i<=500; i++)
    for(int j=0; j<=500; j++)
    dp[i][j] = dp[i][j] = -1;

    scanf("%d", &T);
    while(T --)
    {
        scanf("%d %d", &R, &B);

        printf("Case %d: %.9lf\n",cas ++, DFS(R, B) );
    }

    return 0;
}