HDU 4405 Aeroplane chess 概率DP

题目大意：

就是现在有一个飞行棋的游戏，起点是0， 每次掷骰子的点数(0~6)为走的步数,当走到n点及其之后时游戏结束，其中有一些点的位置在走到上面的时候会飞行岛其他点，飞到的位置如果又可以飞行则继续飞行至下一处一直到不能飞行为止再掷骰子，问需要掷骰子多少次完成游戏，求这个期望

大致思路：

首先由于飞行的路线不会出现往回走的情况，所以是个很简单的概率DP, 方程等见代码

代码如下：

Result  :  Accepted     Memory  :  1436 KB     Time  :  31 ms

/*
 * Author: Gatevin
 * Created Time:  2014/11/30 19:08:57
 * File Name: Asuna.cpp
 */
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#include<iomanip>
using namespace std;
const double eps(1e-8);
typedef long long lint;

int n,m;
int next[100010];
/*
 * 用next[i] 表示在点i处会飞行到哪里停下
 */

double E[100010];
/*
 * 用E[i]表示从点i处出发，完成游戏的期望掷骰子的次数，那么
 * E[i] = sigma(E[next[i + j]]) + 1, (1 <= j <= 6);
 */

int main()
{
    int x,y;
    while(scanf("%d %d", &n, &m), n || m)
    {
        for(int i = 0; i <= n + 6; i++) next[i] = i;
        while(m--)
        {
            scanf("%d %d", &x, &y);
            next[x] = y;
        }
        for(int i = 0; i <= n + 6; i++)
        {
            int tmp = i;
            while(tmp != next[tmp])//由于之多只会有M <= 1000条飞行路径，不会超时
                tmp = next[tmp];
            next[i] = tmp;
        }
        memset(E, 0, sizeof(E));
        for(int i = n - 1; i >= 0; i--)
        {
            for(int j = 1; j <= 6; j++)
                E[i] += E[next[i + j]]/6;
            E[i] += 1;
        }
        printf("%.4f\n", E[0]);
    }
    return 0;
}