Rating HDU-4870 概率DP

HDU-4870

solution

先发明天再写

code

/*SiberianSquirrel*/
/*CuteKiloFish*/
#include<bits/stdc++.h>
using namespace std;

double p;
double dp[30];

void solve(double res = 0) {
    dp[0] = 1 / p;
    dp[1] = dp[0] / p;
    res = dp[0] + dp[1];
    for(int i = 2; i < 20; ++ i) {
        dp[i] = 1 + (1 - p) / p * (1 + dp[i - 2] + dp[i - 1]);
        res += dp[i];
    }
    cout << fixed << setprecision(6) << res * 2 - dp[19] << endl;
}

int main() {
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(nullptr);
#ifdef ACM_LOCAL
    freopen("input", "r", stdin);
    freopen("output", "w", stdout);
#endif
    while(cin >> p)
        solve();
    return 0;
}