ZOJ 3551 Bloodsucker 概率DP

题目大意：

就是现在一个城镇里有n - 1个居民, 1个吸血鬼, 每天晚上都有且只会有两个生物见面, 只有当两个见面的物种不同时, 有p的概率人会变成吸血鬼, 为经过多少个晚上城镇里所有的人都会变成吸血鬼, 求这个天数的期望

大致思路：

应该说是最基本的概率题了...

细节见代码吧

代码如下：

Result  :  Accepted     Memory  :  8084 KB     Time  :  150 ms

/*
 * Author: Gatevin
 * Created Time:  2014/12/24 18:35:37
 * File Name: Sora_Kasugano.cpp
 */
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#include<iomanip>
using namespace std;
const double eps(1e-8);
typedef long long lint;

double dp[1000010];
/*
 * 应该说是最最基础的一个概率DP了吧...
 * 用dp[i]表示有i只吸血鬼时到达目标状态的所需要的期望天数
 * 那么dp[i] = (i*(n - i))/(n*(n - 1)/2)*p*(dp[i + 1] + 1) + (1 - (i*(n - i))/(n*(n - 1)/2)*p)*(dp[i] + 1)
 * 解得dp[i] = dp[i + 1] + n*(n - 1)/(2*(n - i)*i*p)
 * 初始时dp[n] = 0
 * dp[1]即为答案
 */

int main()
{
    int t;
    scanf("%d", &t);
    int n;
    double p;
    while(t--)
    {
        scanf("%d %lf", &n, &p);
        dp[n] = 0;
        for(int i = n - 1; i >= 1; i--)
            dp[i] = dp[i + 1] + 1.*n/i*(n - 1)/2/(n - i)/p;
        printf("%.3f\n", dp[1]);
    }  
    return 0;
}