2016多校第五场 1001 HDU 5781 DP

HDU 5781 DP

题意：
从银行取钱，钱数为随机数[0,w]。如果当前取钱数大于真实的钱数，就会报错。报错次数有限制，求取钱次数的期望值。
思路：
首先用二分试，发现第三个样例过不了。
赛中的时候不是主要负责的题，给了一个dp[i][j]的公式后就没管了。
然而这样的dp过程是10^9的。
实际上，本题是有trick的。假设每次都按二分取，发现最多15次不合法的取法，就可以把所有状态走到。那么大于15的不合法数对答案是没有有效贡献的。于是复杂度编程降低，就可以做了。
勤劳的队友赛中的时候打表发现有极值，然而用三分莫名炸了所以最后也没过。
代码：

    #include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <string>
#include <vector>
#include <map>
using namespace std;
const int MAXN = 2000 + 5;
#define LL long long
#define inf (1e12)
double dp[MAXN][18];
void init()
{
    for(int i = 0 ; i < 18 ; i++) dp[0][i] = 0;
    for(int i = 1 ; i < MAXN ; i++) dp[i][0] = inf;
    for(int j = 1 ; j < 18 ; j++) {
        for(int i = 1 ; i < MAXN ; i++) {
            dp[i][j] = -1;
            for(int k = 1 ; k <= i ; k++) {
                double temp = dp[k - 1][j - 1] * k / (i + 1) + dp[i - k][j] * (i + 1 - k) / (i + 1) + 1;
                if(dp[i][j] == -1 || dp[i][j] > temp) dp[i][j] = temp;
            }
        }
    }
}
int main()
{
    init();
    int n, m;
    while(scanf("%d%d", &n, &m) != EOF) {
        m = min(m, 17);
        printf("%.6f\n", dp[n][m]);
    }
    return 0;
}