POJ 3744 Scout YYF I [概率DP]

原创 2015年11月20日 13:02:04

Description

YYF is a couragous scout. Now he is on a dangerous mission which is to penetrate into the enemy's base. After overcoming a series difficulties, YYF is now at the start of enemy's famous "mine road". This is a very long road, on which there are numbers of mines. At first, YYF is at step one. For each step after that, YYF will walk one step with a probability of p, or jump two step with a probality of 1- p. Here is the task, given the place of each mine, please calculate the probality that YYF can go through the "mine road" safely.

Input

The input contains many test cases ended with EOF
Each test case contains two lines. 
The First line of each test case is N (1 ≤ N ≤ 10) and p (0.25 ≤ p ≤ 0.75) seperated by a single blank, standing for the number of mines and the probability to walk one step. 
The Second line of each test case is N integer standing for the place of N mines. Each integer is in the range of [1, 100000000].

Output

For each test case, output the probabilty in a single line with the precision to 7 digits after the decimal point.

题意:

一条路上每个点标号为1,2,3....,你一开始站在1往正方向出发, 有P的概率走到下一格,(1-P)的概率走到下一格的下一格,现在路上可能会有一些地雷,给出N个地雷的坐标,问你不踩地雷走过这条路的概率。

范围:

地雷数小于等于10,坐标小于10忆

解法:

OXOOOOX

若O表示可以走的,X表示有地雷,那么走过这段路可以分成三个部分

第一部分:坐标1->坐标3 这一阶段达成的概率是1-p

第二部分:坐标3->坐标6

第三部分:坐标6->坐标8 这一阶段达成的概率是1-p

可以发现,只要求出第二部分的概率即可,这一部分的概率可以DP

记F[i]为走过长度为i的没有地雷的路程的概率,那么公式为f[i]=f[i-1]*p+f[i-2]*(1.0-p);

但是因为I可能很大,需要考虑对它进一步分析,最后发现i趋于无穷时,F[I]是有极限的,实际过程中这个极限很快就会达到。

所以处理出大约前1000个F[I],大于1000的取F[1000],然后按阶段处理即可。

代码:

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<iostream>
#include<stdlib.h>
#include<set>
#include<map>
#include<queue>
#include<vector>
#include<bitset>
#pragma comment(linker, "/STACK:1024000000,1024000000")
template <class T>
bool scanff(T &ret){ //Faster Input
    char c; int sgn; T bit=0.1;
    if(c=getchar(),c==EOF) return 0;
    while(c!='-'&&c!='.'&&(c<'0'||c>'9')) c=getchar();
    sgn=(c=='-')?-1:1;
    ret=(c=='-')?0:(c-'0');
    while(c=getchar(),c>='0'&&c<='9') ret=ret*10+(c-'0');
    if(c==' '||c=='\n'){ ret*=sgn; return 1; }
    while(c=getchar(),c>='0'&&c<='9') ret+=(c-'0')*bit,bit/=10;
    ret*=sgn;
    return 1;
}
#define inf 1073741823
#define llinf 4611686018427387903LL
#define PI acos(-1.0)
#define lth (th<<1)
#define rth (th<<1|1)
#define rep(i,a,b) for(int i=int(a);i<=int(b);i++)
#define drep(i,a,b) for(int i=int(a);i>=int(b);i--)
#define gson(i,root) for(int i=ptx[root];~i;i=ed[i].next)
#define tdata int testnum;scanff(testnum);for(int cas=1;cas<=testnum;cas++)
#define mem(x,val) memset(x,val,sizeof(x))
#define mkp(a,b) make_pair(a,b)
#define findx(x) lower_bound(b+1,b+1+bn,x)-b
#define pb(x) push_back(x)
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;

int n,a[1111];
double p,f[1111];
int main(){
    while(scanf("%d%lf",&n,&p)!=EOF){
        mem(f,0);
        f[1]=1.0;
        rep(i,2,1000){
           f[i]=f[i-1]*p+f[i-2]*(1.0-p);
        }
        double ans=1.0;
        rep(i,1,n)scanff(a[i]);
        sort(a+1,a+1+n);
        int an=unique(a+1,a+1+n)-a-1;
        a[0]=0;
        rep(i,1,an){
            int x=a[i];
            int dis=min(a[i]-a[i-1]-1,1000);
            ans=ans*f[dis]*(1.0-p);
        }
        printf("%.7f\n",ans);
    }
    return 0;
}


版权声明:文為博主原創文章,未經博主允許不得轉載。

相关文章推荐

POJ 3744 Scout YYF I(概率DP)

传送门 Scout YYF I Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 7221 Accepted...

POJ 3744-Scout YYF I (概率DP+矩阵快速幂)

Scout YYF I Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 8894   Accepted: 2607...

POJ 3744 Scout YYF I (概率DP+矩阵快速幂)

题目地址:POJ 3744 一个线性概率DP递推式。dp[i]=p*dp[i-1]+(1-p)*dp[i-2]。但是i的值太大。所以可以分成n次,每一次中间过程的纯递推过程用矩阵快速幂来优化。只要想...

poj3744 Scout YYF I 概率dp+矩阵快速幂

题意:在一条路上布置了一些地雷,初始的时候你在起点1处,

poj 3744 Scout YYF I(概率dp+矩阵快速幂)

题意:有个人从1位置开始走,每次走1歩的概率为p,走两步的概率为1-p,在N个位置上有地雷,问这个人不踩到地雷通过的概率。 思路:可以看出安全走到某一位置的概率f(x)=f(x-1)*p+f(x-2...
  • qian99
  • qian99
  • 2014年03月17日 20:48
  • 558

poj 3744 Scout YYF I(矩阵优化概率DP)

Scout YYF I Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 5153   Ac...
  • caduca
  • caduca
  • 2014年11月03日 16:13
  • 537

POJ-3744 Scout YYF I (概率DP&&矩阵快速幂)

设dp[i]表示安全到达点i处的概率,则状态转移方程为:①点i处有雷,dp[i]=0;②点i出无雷,dp[i]=dp[i-1]*p+dp[i-2]*(1-p) 由于雷坐标的范围很大,而雷的数量很少,所...

POJ 3744 Scout YYF I (概率dp+矩阵优化)

Language: Default Scout YYF I Time Limit: 1000MS   Memory Limit: 65536K Total Su...

POJ 3744 Scout YYF I 矩阵快速幂+概率dp

题意:   yyf每次能走一步或二步概率为p和1-p 路上有n个地雷  分布在di位置上 问yyf能安全通过这条路的概率 分析 可以容易的得出 dp【i】=p*dp[i-1]+(1-p)*d...

POJ 3744 Scout YYF I (矩阵优化,分段计算的概率DP)

参考: http://www.cnblogs.com/kuangbin/archive/2012/10/02/2710586.html 此题是一个用矩阵优化的求概率的题目。 主要思想是分...
  • guognib
  • guognib
  • 2013年10月26日 22:29
  • 578
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:POJ 3744 Scout YYF I [概率DP]
举报原因:
原因补充:

(最多只允许输入30个字)