山东省第四届ACM大学生程序设计竞赛 The number of steps 概率dp

The number of steps

Time Limit: 1000ms   Memory limit: 65536K  有疑问？点这里^_^

题目描述

    Mary stands in a strange maze, the maze looks like a triangle（the first layer have one room，the second layer have two rooms，the third layer have three rooms …）. Now she stands at the top point(the first layer), and the KEY of this maze is in the lowest layer’s leftmost room. Known that each room can only access to its left room and lower left and lower right rooms .If a room doesn’t have its left room, the probability of going to the lower left room and lower right room are a and b (a + b = 1 ). If a room only has it’s left room, the probability of going to the room is 1. If a room has its lower left, lower right rooms and its left room, the probability of going to each room are c, d, e (c + d + e = 1). Now , Mary wants to know how many steps she needs to reach the KEY. Dear friend, can you tell Mary the expected number of steps required to reach the KEY?

输入

There are no more than 70 test cases. 

 

In each case , first Input a positive integer n(0

The input is terminated with 0. This test case is not to be processed.

输出

Please calculate the expected number of steps required to reach the KEY room, there are 2 digits after the decimal point.

示例输入

30.3 0.70.1 0.3 0.60 

示例输出

3.41

提示

 

来源

2013年山东省第四届ACM大学生程序设计竞赛

示例程序

主人公站在一个金字塔顶端往下面走

规定只能往左 下 左下 走然后对应相对的概率

求从1,1到n，n的期望

有点像数塔从底往上推求得

ACcode：

#include <map>#include <queue>#include <cmath>#include <cstdio>#include <cstring>#include <stdlib.h>#include <iostream>#include <algorithm>#define maxn 100using namespace std;double dp[maxn][maxn];double a,b,c,d,e;int n;int main(){    while(scanf("%d",&n),n){        scanf("%lf%lf%lf%lf%lf",&a,&b,&c,&d,&e);        memset(dp,0,sizeof(dp));        for(int i=n-1;i>0;--i)            dp[n][i]+=dp[n][i+1]+1;        for(int i=n-1;i>0;--i){            dp[i][i]+=a*(dp[i+1][i+1]+1)+b*(dp[i+1][i]+1);            for(int j=i-1;j>0;--j)                dp[i][j]+=c*(dp[i+1][j+1]+1)+d*(dp[i+1][j]+1)+e*(dp[i][j+1]+1);        }        printf("%.2lf\n",dp[1][1]);    }    return 0;}