概率dp

题目链接

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<vector>
#include<queue>
#include<math.h>
using namespace std;
typedef long long LL;
LL p;
LL f[100004];
double dp[1100][1100];
int main()
{
    int w,b;
    scanf("%d %d",&w,&b);
    int i,j;
    for(i=1;i<=w;i++)
        dp[i][0]=1;
    for(i=1;i<b;i++)
        dp[0][i]=0;
    for(i=1;i<=w;i++)
        for(j=1;j<=b;j++)
    {
        dp[i][j]+=(i*1.0)/(i+j);
        if(j>=2)
            dp[i][j]+=(j*1.0)/(i+j)*((j-1)*1.0)/(i+j-1)*(i*1.0)/(i+j-2)*dp[i-1][j-2];
        if(j>=3)
            dp[i][j]+=(j*1.0)/(i+j)*((j-1)*1.0)/(i+j-1)*((j-2)*1.0)/(i+j-2)*dp[i][j-3];
    }
    printf("%.9f\n",dp[w][b]);
    return 0;
}