C++(g++ 7.5.0) 解法, 执行用时: 38ms, 内存消耗: 32056K, 提交时间: 2023-04-20 12:00:23

#include <vector>
#include <iomanip>
#include <iostream>
#include <algorithm>
using namespace std;

#define endl '\n'
using LL=long long;
using LD=long double;
constexpr int N=1e3+10;
LD dp[N][N][2];

void solve() {
    int w,b;
    cin>>w>>b;

    for(int i=0;i<=w;i++) {
        for(int j=0;j<=b;j++) {
            LD pw=LD(i)/(i+j);
            LD pb=LD(j)/(i+j);
            
            if(i==0) dp[i][j][0]=0;
            else dp[i][j][0]=pw+(1-dp[i][j-1][1])*pb;

            if(i==0||i+j<=2) dp[i][j][1]=1;
            else {
                if(i) dp[i][j][1]+=pw;
                if(j>=2) dp[i][j][1]+=pb*LD(j-1)/(i+j-1)*(1-dp[i][j-2][0]);
                if(i&&j) dp[i][j][1]+=pb*LD(i)/(i+j-1)*(1-dp[i-1][j-1][0]);
            }
        }
    }
    cout<<fixed<<setprecision(12)<<dp[w][b][0];
}

int main() {
    ios::sync_with_stdio(0);
    cin.tie(nullptr);
    solve();
    return 0;
}

C++ 解法, 执行用时: 18ms, 内存消耗: 8232K, 提交时间: 2022-07-21 00:11:14

#include<bits/stdc++.h>
using namespace std;
const int N=1e3+5;
double f[N][N];
int w,b;
void solve(){
    cin>>w>>b;
    for(int i=1;i<=w;i++) f[i][0]=1;
    for(int i=1;i<=w;i++){
        for(int j=1;j<=b;j++){
            f[i][j]=1.0*i/(i+j);
            if(j>=3) f[i][j]+=1.0*j/(i+j)*(j-1)/(i+j-1)*(j-2)/(j+i-2)*f[i][j-3];
            if(j>=2) f[i][j]+=1.0*j/(i+j)*(j-1)/(i+j-1)*(i)/(j+i-2)*f[i-1][j-2];
        }
    }
    printf("%.10lf\n",f[w][b]);
}
int main(){
    solve();
}