hdoj 5001 概率dp

hdoj 5001

题意：一个人可以在任意点出发，以相同概率经过相邻点，问不经过各点的概率。

思路：dp[i][j]是第i步到达j点，枚举所有的点Vj累加所有的dp[i][Vj]即可。

#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
double dp[10020][60], ans[60];
vector<int>edge[60];

main() {
    int n, m, d, t;
    scanf("%d", &t);
    while(t--) {
        scanf("%d %d %d", &n, &m, &d);
        for(int i = 0; i < 60; i++) edge[i].clear();
        memset(ans, 0, sizeof ans);
        for(int i = 0; i < m; i++) {
            int a, b;
            scanf("%d %d", &a, &b);
            edge[a].push_back(b);
            edge[b].push_back(a);
        }
        for(int v = 1; v <= n; v++) {
            memset(dp, 0, sizeof dp);
            for(int i = 1; i <= n; i++) dp[0][i] = 1.0 / n;
            ans[v] += dp[0][v];
            for(int i = 0; i < d; i++){
                for(int j = 1; j <= n; j++){
                    if(j == v) continue;
                    for(int k = 0; k < edge[j].size(); k++){
                        if(j == edge[j][k]) continue;
                        dp[i + 1][edge[j][k]] += (1.0 / edge[j].size()) * dp[i][j];
                    }
                }
                ans[v] += dp[i + 1][v];
            }
        }
        for(int i = 1; i <= n; i++) {
            printf("%.10f\n", 1.0 - ans[i]);
        }
    }
}