代码:
#include<bits/stdc++.h>
using namespace std;
double dp[2000][2000];
int n,m;
int main(){
cin>>n>>m;
for(int i = 0;i<=n;i++){
for(int j = 0;j<=m;j++){
if(i==0 && j==0)
dp[i][j] = 0;
else if(i==0)
dp[i][j] = 0;
else if(j==0)
dp[i][j] = 1;
else
dp[i][j] = 1.0*i/(i+j) + 1.0*j/(i+j) * (j-1)/(i+j-1)* ( ( j-2>=0 ? 1.0*i/(i+j-2)*dp[i-1][j-2] : 0) + ( j-3>=0 ? 1.0*(j-2)/(i+j-2)*dp[i][j-3] : 0 ));
/// cout<<"i:"<<i<<" j:"<<j<<" :"<<dp[i][j]<<endl;
}
}
printf("%.4lf\n",dp[n][m]);
}
/*dp[n][m] =
1 中了 n/(m+n);
(A抽了一张不中奖,剩下n张中奖,m-1张不中奖) 然后B也抽了一张不中奖
剩下n张中奖,m-2张不中奖.B可以随便摸
B如果摸到中奖,则剩下n-1和m-2 B如果摸到不中奖,则剩下n和m-3
2 没中 m/(n+m) * (m-1)/(n+m-1)* { dp[n-1][m-2] + dp[n][m-3] }
*/
该题源于360笔试题。
考试时推出公式,但是边界写错了。。。。真想学好dp。。。