hdu6568.Math——（概率dp)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6568

题目大意

Avin要带着机器人从0走到L，走路速度为每秒一格并且只会在整数格停下，他走路按照以下规则：

如果Avin没有弄丢机器人，则有p的概率弄丢机器人
如果Avin已经弄丢了机器人，则有q的概率发现，特别的，当他到达终点却没有机器人时，会立刻发现
如果Avin发现他弄丢了机器人，则会一直往左走直到遇见机器人，否则他会向右走一格

问Avin从0走到L的期望花费距离是多少？

思路：

设表示从i到i+1的期望花费。

首先，如果没有弄丢机器人，

其次，如果弄丢了机器人，那么就有两种情况

第一：Avin从i到L都没有发现丢弃了机器人，那么到L会立即发现，然后就会从L一直走回来到i，所以这样的概率为

，一共走的距离为（因为有去有回）,所以这里的期望为

第二：Avin在i到L-1的路上发现了机器人已经失踪了，我们枚举发现的位置x，那么发现的概率为

，期望为，因为是枚举x，所有x都满足情况，所以这一部分最后期望为

最后，因为是图上随机游走过程，当你回到i点时，又会重新计算一遍，所以最后整理公式如下：

然后你将g(i)移过去，预处理的前缀和，最后答案就是

#pragma GCC optimize(2)
#pragma GCC optimize(3)
#pragma GCC optimize(4)
#include<bits/stdc++.h>
using namespace std;
#define y1 y11
#define fi first
#define se second
#define pi acos(-1.0)
#define LL long long
//#define mp make_pair
#define pb emplace_back
#define ls rt<<1, l, m
#define rs rt<<1|1, m+1, r
#define ULL unsigned LL
#define pll pair<LL, LL>
#define pli pair<LL, int>
#define pii pair<int, int>
#define piii pair<pii, int>
#define pdd pair<double, double>
#define mem(a, b) memset(a, b, sizeof(a))
#define debug(x) cerr << #x << " = " << x << "\n";
#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
//head

const int N = 1e5 + 10;
int n;
double p, q, _q[N], sum[N];
int main() {
    while(~scanf("%d %lf %lf", &n, &p, &q)) {
        sum[0] = 0;
        _q[0] = 1;
        for (int i = 1; i <= n; ++i) {
            _q[i] = _q[i-1]*(1-q);
            sum[i] = sum[i-1] + _q[i]*i;
        }
        double ans = 0;
        for (int i = 0; i < n; ++i) {
            ans += 1+(2*p*(_q[n-i]*(n-i)+q*sum[n-i-1]))/(1-p);
        }
        printf("%.10f\n", ans);
    }
    return 0;
}