poj3744(概率dp+快速幂)

题意:YYF从坐标1出发,每次走一步的概率是p,走两步的概率是1-p,然后有n个点上有地雷。问YYF走过雷区的概率是多少。


解法:ans[i]表示走到i的概率。概率公式ans[i]=p*ans[i-1]+(1-p)ans[i-2]

转移矩阵:
  ans[i]    | p ,1-p  |    ans[i-1]
            =|            |*
ans[i-1]   | 1 , 0    |   ans[i-2]

根据每个地雷的位置分段计算概率;

代码:

/******************************************************
* author:xiefubao
*******************************************************/
#pragma comment(linker, "/STACK:102400000,102400000")
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <vector>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <string.h>
//freopen ("in.txt" , "r" , stdin);
using namespace std;

#define eps 1e-8
const double pi=acos(-1.0);
typedef long long LL;
const int Max=15;
int m=2;
int n;
double p;
int mine[Max];
struct matrix
{
    double num[3][3];
     matrix(){
         memset(num,0,sizeof num);
     }
};
matrix operator*(const matrix& a,const matrix& b)
{
    matrix ans;
    for(int i=0;i<m;i++)
    for(int j=0;j<m;j++)
    for(int k=0;k<m;k++)
    ans.num[j][k]+=a.num[j][i]*b.num[i][k];
    return ans;
}
matrix pow(matrix mat,int t)
{
    matrix ans;
    ans.num[0][0]=1;
    ans.num[1][1]=1;
    while(t)
    {
        if(t&1)
        ans=ans*mat;
        mat=mat*mat;
        t>>=1;
    }
    return ans;
}
double getans(int t)
{
    if(t<=1)
    return 0;
    if(t==2)
    return 1;
    if(t==3)
    return p;
    t-=3;
    matrix tool;
    tool.num[0][0]=p;
    tool.num[0][1]=1-p;
    tool.num[1][0]=1;
    tool.num[1][1]=0;
    tool=pow(tool,t);
    //cout<<tool.num[0][0]<<" "<<tool.num[0][1]<<" "<<p<<endl;
    return tool.num[0][0]*p+tool.num[0][1];
}
int main()
{
    while(scanf("%d%lf",&n,&p)==2)
    {
        mine[0]=0;
        for(int i=1;i<=n;i++)
        scanf("%d",mine+i);
        sort(mine+1,mine+n+1);
        double out=1;
        for(int i=1;i<=n;i++)
        {
            out*=getans(mine[i]-mine[i-1]);
            out*=(1.0-p);
        }
        printf("%.7lf\n",out);
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值