HDU 4546简单的概率 dp

/*
题意:杭州邀请赛的一道题;有标记着1~n的环,初始时有个机器人在1的位置,有m个操作,
每个操作是使机器人随机的顺时针或逆时针的走w步,求所有操作过后机器人停在区间[l,r]内的概率;
思路:挺明显的概率dp,首先可以很直观的定义出状态dp[i][j],表示第j次操作后停留在第i位置的概率,
转移方程是 dp[(i + w) % n][j] = dp[i[j-1] * 0.5;
dp[i-w][j]=dp[i][j-1]*0.5;其中i-w不断的加n使它大于等于0;
但是题目操作数比较大,这样会超内存,观察发现每次转移只和上一次操作有关,所以可以用个滚动数组
来节省空间;优化空间后转移方程是 dp[(i + w) % n][j&1] += dp[i][!(j&1)] * 0.5;
dp[x][j&1] += dp[i][!(j&1)] * 0.5;x=i-w,x也要叠加到非负;
最后答案就是所有操作完后,位置l到r的概率总和。

*/
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;

#define manx 209
double dp[manx][2];
const double eps=0.0;

int main(){
    int n,m,l,r;
    while(cin>>n>>m>>l>>r){
        if(n==0 && m==0 && l==0 && r==0) break;
        int w,i,j;
        for(i=0;i<n;i++) dp[i][0]=dp[i][1]=0;
        dp[0][0] = 1.0;
        for(j=1;j<=m;j++){
            scanf("%d",&w);
            for(i=0;i<n;i++){
                dp[(i+w)%n][j&1] += dp[i][!(j&1)]*0.5;
                int ans = i - w;
                if(ans<0) ans += n;
                dp[ans][j&1] += dp[i][!(j&1)]*0.5;
            }
            for(i=0;i<=n;i++)
            	dp[i][!(j&1)] = eps;
        }
        double res = 0.0;
        l--,r--;
        for(i=0;i<n;i++)
        	if(i>=l && i<=r) res += dp[i][m&1];
       	printf("%.4lf\n",res);
    }
}

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值