HDU 4576 Robot （概率 & 期望）

Robot

Time Limit: 8000/4000 MS (Java/Others)    Memory Limit: 102400/102400 K (Java/Others)
Total Submission(s): 779    Accepted Submission(s): 304

Problem Description

Michael has a telecontrol robot. One day he put the robot on a loop with n cells. The cells are numbered from 1 to n clockwise.



At first the robot is in cell 1. Then Michael uses a remote control to send m commands to the robot. A command will make the robot walk some distance. Unfortunately the direction part on the remote control is broken, so for every command the robot will chose a direction(clockwise or anticlockwise) randomly with equal possibility, and then walk w cells forward.
Michael wants to know the possibility of the robot stopping in the cell that cell number >= l and <= r after m commands.

 

 

Input

There are multiple test cases. 
Each test case contains several lines.
The first line contains four integers: above mentioned n(1≤n≤200) ,m(0≤m≤1,000,000),l,r(1≤l≤r≤n).
Then m lines follow, each representing a command. A command is a integer w(1≤w≤100) representing the cell length the robot will walk for this command.  
The input end with n=0,m=0,l=0,r=0. You should not process this test case.

 

 

Output

For each test case in the input, you should output a line with the expected possibility. Output should be round to 4 digits after decimal points.

 

 

Sample Input

3 1 1 2 1 5 2 4 4 1 2 0 0 0 0

 

 

Sample Output

0.5000 0.2500

 

 

Source

2013ACM-ICPC杭州赛区全国邀请赛

 

 题意：给你一个n，m，l，r。n个格子围成一个圈，一个机器人从第一个格子开始，经过m个步，每个步知道走几格，但是不知道方向，问你最后在l到r之间的概率。

#include<iostream>

#include<cstdio>

#include<cstring>



using namespace std;



const int N=220;



int n,m,l,r;

double dp[2][N];



int main(){



    //freopen("input.txt","r",stdin);



    while(~scanf("%d%d%d%d",&n,&m,&l,&r)){

        if(n==0 && m==0 && l==0 && r==0)

            break;

        for(int i=0;i<=n;i++)

            dp[0][i]=0;

        dp[0][0]=1;

        int x,cur=0;

        while(m--){

            scanf("%d",&x);

            for(int i=0;i<n;i++)

                dp[cur^1][i]=0;

            for(int i=0;i<n;i++){

                if(dp[cur][i]==0)

                    continue;

                dp[cur^1][((i-x)%n+n)%n]+=0.5*dp[cur][i];

                dp[cur^1][(i+x)%n]+=0.5*dp[cur][i];

            }

            cur^=1;

        }

        double ans=0;

        for(int i=l-1;i<r;i++)  //注意这里的范围

            ans+=dp[cur][i];

        printf("%.4lf\n",ans);

    }

    return 0;

}