HDU--4405（概率DP）

2014-12-30 11:35:57

思路：入门级概率DP。（虽然也写了挺久orz...）

　　其实就是个很简单的记忆化搜索 或者 递推，用dp[i]表示从 i 出发直至完成飞行需要的期望步数。

　　如果该点有flight line，那么dp[i] = dp[flight[i]]，否则 dp[i] = (dp[i + k] + 1) / 6.0 （1<=k<=6)

 1 #pragma comment(linker,"/STACK:1024000000,1024000000")
 2 #include <cstdio>
 3 #include <cstring>
 4 #include <cstdlib>
 5 #include <cmath>
 6 #include <vector>
 7 #include <map>
 8 #include <set>
 9 #include <stack>
10 #include <queue>
11 #include <iostream>
12 #include <algorithm>
13 using namespace std;
14 #define lp (p << 1)
15 #define rp (p << 1|1)
16 #define getmid(l,r) (l + (r - l) / 2)
17 #define MP(a,b) make_pair(a,b)
18 typedef long long ll;
19 typedef unsigned long long ull;
20 typedef pair<int,int> pii;
21 const int INF = 1 << 30;
22 
23 int N,M;
24 double dp[100020];
25 int vis[100020];
26 
27 double Solve(int p){
28     if(dp[p]) return dp[p];
29     if(p >= N) return 0.0;
30     if(vis[p]){
31         dp[p] += Solve(vis[p]);
32         return dp[p];
33     }
34     for(int i = 1; i <= 6; ++i)
35         dp[p] += (Solve(p + i) + 1.0) / 6.0;
36     return dp[p];
37 }
38 
39 int main(){
40     int a,b;
41     while(scanf("%d%d",&N,&M) != EOF){
42         if(N == 0 && M == 0) break;
43         memset(dp,0,sizeof(dp));
44         memset(vis,0,sizeof(vis));
45         for(int i = 1; i <= M; ++i){
46             scanf("%d%d",&a,&b);
47             vis[a] = b;
48         }
49         Solve(0);
50         printf("%.4f\n",dp[0]);
51     }
52     return 0;
53 }

 

转载于:https://www.cnblogs.com/naturepengchen/articles/4193252.html