题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5001
【解题报告】
概率DP。令DP[i][ map[j][k] ] 表示第i步从第j个城市走到j的第k个子节点的概率。
转移方程: dp[i][ map[j][k] ]+=dp[i-1][j]*1.0/map[j].size();
显然如果dp[i-1][j]不为0,说明第i-1步可以到达j,那么在j点就可以由第i-1步递推至第i步,每一层i之间不会相互干扰,所以不会有重复或遗漏。
那么此时对第i个节点,忽略它,然后进行dp即可保证所有可行路径中不包含i,那么累加dp[d][1..n|i]即为第i组答案。
时间复杂度为O(d*n^2),时限15000ms,可能我姿势比较丑或者T过大,实测2100ms过了这道题。
【参考代码】
#include<cstdio>
#include<iostream>
#include<stack>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
vector<int>mp[60];
double dp[10000+10][60];
int n,m,d;
void solve( int x )
{
memset(dp,0,sizeof(dp));
for(int i=1;i<=n;i++)dp[0][i]=1.0/n;
for(int i=0;i<d;i++)
{
for(int j=1;j<=n;j++)if(j!=x)
{
for(int k=0;k<mp[j].size();k++)
{
dp[i+1][ mp[j][k] ]+= dp[i][j] *1.0/(double)mp[j].size();
}
}
}
double ans=0.0;
for(int i=1;i<=n;i++)if(i!=x)
{
ans+=dp[d][i];
}
printf("%.5f\n",ans);
}
int main()
{
int t;
cin>>t;
while( t-- )
{
scanf("%d%d%d",&n,&m,&d);
for(int i=0;i<60;i++)mp[i].clear();
for(int i=1;i<=m;i++)
{
int u,v;
scanf("%d%d",&u,&v);
mp[u].push_back(v);
mp[v].push_back(u);
}
for(int i=1;i<=n;i++)solve(i);
}
return 0;
}