POJ 3682 King Arthur's Birthday Celebration 概率DP

最新推荐文章于 2020-01-21 10:22:22 发布

_NMfloat_

最新推荐文章于 2020-01-21 10:22:22 发布

阅读量1.2k

点赞数

文章标签： poj dp

本文链接：https://blog.csdn.net/NMfloat/article/details/50650489

版权

题意：亚瑟王掷一枚硬币，概率p正面向上，概率1-p反面朝上，现在亚瑟王要掷 $k$ 次正面朝上，第 $i$ 次掷硬币时花费 $2*i-1$ 。问：期望要掷多少枚硬币才能达到 $k$ 次正面朝上，以及达到 $k$ 次正面朝上时的花费。

思路：

$E i 表示已经掷到 i 次正面朝上了，还要再掷多少次才能到达 k 次正面朝上。$ $E_i表示已经掷到i次正面朝上了，还要再掷多少次才能到达k次正面朝上。$

$E i = p E i + 1 + (1 - p) E i + 1$ $E_i = p E_{i+1}+(1-p)E_i+1$

$很容易解出期望掷 k p 次达到目的。$ $很容易解出期望掷\frac{k}{p}次达到目的。$

$F i 表示已经掷到 i 次正面朝上了，还要再花费多少才能到达 k 次正面朝上。$ $F_i表示已经掷到i次正面朝上了，还要再花费多少才能到达k次正面朝上。$

$F i = p F i + 1 + (1 - p) F i + 2 * E i - 1$ $F_i = pF_{i+1}+(1-p)F_i+2*E_i-1$

$+ 2 * E i - 1 是因为从 i 次正面朝上的时候掷下一次的时候是在第 E i 次掷的。$ $+2*E_i-1是因为从i次正面朝上的时候掷下一次的时候是在第E_i次掷的。$

$解出花费为 n 2 + n p 2 - n p$ $解出花费为\frac{n^2+n}{p^2}-\frac{n}{p}$

http://poj.org/problem?id=3682

/*********************************************
    Problem : POJ 3682
    Author  : NMfloat
    InkTime (c) NM . All Rights Reserved .
********************************************/

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <cmath>
#include <ctime>
#include <cstdio>
#include <vector>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>

#define rep(i,a,b)  for(int i = (a) ; i <= (b) ; i ++) //遍历
#define rrep(i,a,b) for(int i = (b) ; i >= (a) ; i --) //反向遍历
#define repS(it,p) for(auto it = p.begin() ; it != p.end() ; it ++) //遍历一个STL容器
#define repE(p,u) for(Edge * p = G[u].first ; p ; p = p -> next) //遍历u所连接的点
#define cls(a,x)   memset(a,x,sizeof(a))
#define eps 1e-8

using namespace std;

const int MOD = 1e9+7;
const int INF = 0x3f3f3f3f;
const int MAXN = 1e5+5;
const int MAXE = 2e5+5;

typedef long long LL;
typedef unsigned long long ULL;

int T,n,m;

int fx[] = {0,1,-1,0,0};
int fy[] = {0,0,0,-1,1};

double p;

void input() {
    scanf("%lf",&p);
}

void solve() {
    printf("%.3f %.3f\n",n/p,((n*n+n)/p-n)/p);
}

int main(void) {
    //freopen("a.in","r",stdin);
    while(scanf("%d",&n),n) {
        input();
        solve();
    }
    return 0;
}