HDU 5001 Walk 概率DP

题目大意：

一个人随即从一个点出发，到达邻接点的概率相同，求出走d步都不会到达1~n点的每一项的不可能概率（这里第一次随即取的点是要求的点也算到达过了）

 

这道题开始一直在计算到达那一点的可能性，最后用1-ans[i]，但到最后还是没有找到自己哪里错了，有机会再看看

 

后来直接计算不可能概率，通过dp找到一直走d步后还能到达非自己所要的num的那些点，每次访问到num点，都直接跳过不理它，使其概率不能下传就行了

然后在solve函数中，把所有除num外的dp[j][d]相加就成功了

 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <vector>
using namespace std;
#define N 52
vector<int> G[N];
int n,m,d;
double dp[N][10005],ans[N];
void bfs(int num)
{
    memset(dp,0,sizeof(dp));
    for(int i=1;i<=n;i++)
        if(i!=num)
            dp[i][0] = 1.0/n;
        else dp[i][0]=0; //对于已经到达的点不往下继续传值

    for(int i=0;i<d;i++){
        for(int j=1;j<=n;j++){
            for(int m=0;m<(int)G[j].size();m++){
                //因为求不到的可能，所以遇到自己需要到达的点就不计算，使其保持空值，无法向下传值
                if(G[j][m]!=num){
                    dp[G[j][m]][i+1]+=dp[j][i]*1.0/(int)G[j].size();
                    //cout<<dp[j][i]<<endl;
                }
            }
        }
    }
}

void solve()
{
    memset(ans,0,sizeof(ans));
    for(int i=1;i<=n;i++){
        bfs(i);
        for(int j=1;j<=n;j++){
            if(j!=i)
                ans[i]+=dp[j][d];
                //将所有能走到最后一步且不是自己所要终点的可能性相加即是不到i的概率
        }
    }
}

int main()
{
    int T,a,b;
    scanf("%d",&T);
    while(T--){

        scanf("%d%d%d",&n,&m,&d);

        for(int i=1;i<=n;i++) G[i].clear();

        for(int i=0;i<m;i++){
            scanf("%d%d",&a,&b);
            G[a].push_back(b);
            G[b].push_back(a);
        }

        solve();

        for(int i=1;i<=n;i++)
            printf("%.10f\n",ans[i]);
    }
    return 0;
}