Codeforce148D（概率dp）

题意：

给出w个白球b个黑球，公主和龙没人轮流去球，当然公主要最先取，最先取出白球的赢。求：公主赢的概率。

题解：

状态方程：
dp[i][j]+=1.0*i/(i+j);
dp[i][j]+=dp[i][j-3]*j/(i+j)*(j-1)/(i+j-1)*(j-2)/(i+j-2);  if(j>=3)
dp[i][j]+=dp[i-1][j-2]*j/(i+j)*(j-1)/(i+j-1)*i/(i+j-2);     if(i>=1&&j>=2)

从未知退出已知是常见概率dp推导方法。

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
using namespace std;
//typedef __int64 lld;
#define oo 0x3f3f3f3f
#define Mod 1000000007
#define maxn 1010
#define Exp 1e-9
double dp[maxn][maxn];

int main()
{
    int w,b;
    while(scanf("%d %d",&w,&b)!=EOF)
    {
        memset(dp,0,sizeof dp);
        for(int i=1;i<=w;i++)
            dp[i][0]=1.0;
        for(int j=0;j<=b;j++)
            dp[0][j]=0.0;
        for(int i=1;i<=w;i++)
        {
            for(int j=1;j<=b;j++)
            {
                dp[i][j]+=1.0*i/(i+j);
                if(j>=3)
                    dp[i][j]+=dp[i][j-3]*j/(i+j)*(j-1)/(i+j-1)*(j-2)/(i+j-2);
                if(i>=1&&j>=2)
                    dp[i][j]+=dp[i-1][j-2]*j/(i+j)*(j-1)/(i+j-1)*i/(i+j-2);
            }
        }
        printf("%.9lf\n",dp[w][b]);
    }
    return 0;
}