UESTC 2015dp专题 F 邱老师看电影 概率dp

邱老师看电影

Time Limit: 20 Sec  Memory Limit: 256 MB

题目连接

http://acm.uestc.edu.cn/#/contest/show/65

Description

一天邱老师心血来潮想去看电影，但是邱老师的妹子想去逛街，他们谁也没有办法说服对方，于是准备来玩一个游戏来决定听谁的。

邱老师找来w只白鼠和b只黑鼠，邱老师和妹子轮流从袋子里面抓老鼠，谁先抓到白色老鼠谁就赢。

但是有酱神在旁边捣乱，邱老师每抓一只老鼠出来，酱神就偷偷的也从里面抓一只出来，这3个人抓出来的老鼠都是随机的。

如果袋子里没有白老鼠，且之前没有人拿到白老鼠的时候，邱老师胜。

为了体现绅士精神，邱老师让妹子先抓，那么妹子赢的概率是多少呐？

wij 千 克的蛋糕，Gorwin从左上角(1,1)的格子开始走，走到右下角(n,m)的格子。在每一步中，Gorwin可以向右或者向下走，即是：Gorwin 站在(i,j)的时候，她可以走向(i+1,j)或者(i,j+1) (然而她不能走出这个花园)。 当Gorwin到达一个格子的时候，她可以把那个格子里面的蛋糕吃完或者不吃。但是，她不能只吃一部分。她的胃不是那么大，所以她最多只能吃K千克的蛋 糕。现在，Gorwin站在左上角，她在看蛋糕园的地图，想要找出一条路，能够使得她吃到的蛋糕最多的一条路。请你来帮帮忙。

Input

只有两个数字 w和b w<=1000 b<=1000
 

Output

输出妹子赢的概率 保留9位小数

Sample Input

1 3

Sample Output

0.500000000

HINT

题意

题解：

概率dp

代码：

//qscqesze
#include <cstdio>
#include <cmath>
#include <cstring>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <typeinfo>
#include <fstream>
#include <map>
#include <stack>
typedef long long ll;
using namespace std;
//freopen("D.in","r",stdin);
//freopen("D.out","w",stdout);
#define sspeed ios_base::sync_with_stdio(0);cin.tie(0)
#define maxn 200001
#define mod 10007
#define eps 1e-9
int Num;
char CH[20];
//const int inf=0x7fffffff;   //нчоч╢С
const int inf=0x3f3f3f3f;
/*

inline void P(int x)
{
    Num=0;if(!x){putchar('0');puts("");return;}
    while(x>0)CH[++Num]=x%10,x/=10;
    while(Num)putchar(CH[Num--]+48);
    puts("");
}
*/
inline ll read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
inline void P(int x)
{
    Num=0;if(!x){putchar('0');puts("");return;}
    while(x>0)CH[++Num]=x%10,x/=10;
    while(Num)putchar(CH[Num--]+48);
    puts("");
}
//**************************************************************************************

double dp[2000][2000];
double dfs(int w,int b)
{
    if(w==0) return 0;
    if(b==0&&w!=0) return 1;
    if(dp[w][b]>=0) return dp[w][b];
    dp[w][b]=w*1.0/(w+b);
    if(b>=2)
    {
        double tb=b*1.0/(w+b)*(b-1)*1.0/(w+b-1);
        if(b>=3)
            dp[w][b]+=tb*dfs(w,b-3)*(b-2)*1.0/(w+b-2);
        if(w>=1)
            dp[w][b]+=tb*dfs(w-1,b-2)*(w*1.0)/(w+b-2);
    }
    return dp[w][b];
}

int w,b;

int main()
{
    memset(dp,-1,sizeof(dp));
    scanf("%d%d",&w,&b);
    printf("%.9f
",dfs(w,b));
}