hdu 5001 从任意点出发任意走d步不经过某点概率

来源：互联网发布：乔丹生涯数据总汇编辑：程序博客网时间：2024/05/01 01:41

http://acm.hdu.edu.cn/showproblem.php?pid=5001

给定n个点m条边的无向图问从任意点出发任意走d步，从不经过某个点的概率

本想先算路过每个点的概率然后用1减去这个概率但是由于可以重复路过所以无法判断是不是第一次经过这个点

所以我们应该直接做不路过的概率即类似bfs的一步步走如果走到了要计算的点就停下来意思就是除了要计算概率的那个点不能走以外其他都能走最后统计下概率的和就是不路过这个点的概率

#include <cstdio>#include <cstdlib>#include <cmath>#include <cstring>#include <string>#include <queue>#include <vector>#include<map>#include <iostream>#include <algorithm>using namespace std;#define RD(x) scanf("%d",&x)#define RD2(x,y) scanf("%d%d",&x,&y)#define clr0(x) memset(x,0,sizeof(x))typedef long long LL;typedef pair<int,int> pii;const double INF=2100000000;int n,m,D;vector<int> g[55];double dp[55][10050];int main(){    int _,u,v;    RD(_);    while(_--){        RD(n);RD2(m,D);        for(int i = 1;i <= n;++i)   g[i].clear();        for(int i = 1;i <= m;++i){            RD2(u,v);            g[u].push_back(v);            g[v].push_back(u);        }        double ans;        for(int i = 1;i <= n;++i){            clr0(dp);            ans = 0;            for(int d = 0;d <= D;++d){                if(!d){                    for(int j = 1;j <= n;++j)if(i != j){                        dp[j][0] = 1.0/(double)n;                    }                }                else{                    for(int j = 1;j <= n;++j)if(i != j){                        for(int k = 0;k < g[j].size();++k){                            int m = g[j][k];                            if(m != i)                                dp[j][d] += dp[m][d-1]/(double)g[m].size();                        }                    }                }            }            for(int j = 1;j <= n;++j)if(i != j){                ans += dp[j][D];            }            printf("%.9lf\n",ans);        }    }    return 0 ;}

0 0