hdu5001（概率dp）

题目连接：http://acm.hdu.edu.cn/showproblem.php?pid=5001

题意：一个人随即从一个点出发，到达邻接点的概率相同，求出走d步都不会到达1~n点的每一点i的概率。

分析：每次删掉i点去计算d步后到达其他点的概率，加起来的都是没有到达i点的概率。

一开始算出x步到达每点的概率，然后d步后（1-dp[x][i]）^d，错了，因为算每点到达的概率时已经算了到达i点时的概率了。

网络赛时这道水题做不出，好殇。。。。。

#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>
#include <cstdlib>
#include <stack>
#include <vector>
#include <set>
#include <map>
#define LL long long
#define mod 100000000
#define inf 0x3f3f3f3f
#define eps 1e-9
#define N 100010
#define FILL(a,b) (memset(a,b,sizeof(a)))
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
using namespace std;
double dp[10010][60];
vector<int>g[100];
int main()
{
    int t,n,m,d;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d%d",&n,&m,&d);
        for(int i=0;i<=n;i++)g[i].clear();
        for(int i=1;i<=m;i++)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            g[u].push_back(v);
            g[v].push_back(u);
        }
        for(int i=1;i<=n;i++)
        {
            FILL(dp,0);
            for(int j=1;j<=n;j++)
                if(j!=i)dp[0][j]=1.0/n;
            for(int x=1;x<=d;x++)
            {
                for(int j=1;j<=n;j++)
                {
                    if(i==j)continue;
                    for(int k=0,sz=g[j].size();k<sz;k++)
                    {
                        int v=g[j][k];
                        if(v!=i)dp[x][v]+=dp[x-1][j]/sz;
                    }
                }
            }
            double res=0;
            for(int j=1;j<=n;j++)
                if(i!=j)res+=dp[d][j];
            printf("%.10lf
",res);
        }
    }
}