【BZOJ 1415】 [Noi2005]聪聪和可可

最新推荐文章于 2021-09-18 17:33:27 发布

suxb201

最新推荐文章于 2021-09-18 17:33:27 发布

阅读量370

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/sxb_201/article/details/50849274

本文介绍了一种使用SPFA（Shortest Path Faster Algorithm）进行最短路径搜索，并结合DFS（Depth-First Search）算法求解特定问题的方法。通过实例展示了如何建立图结构、运用SPFA算法更新节点间的关系，以及如何通过DFS算法进一步计算节点间的期望距离。

1 #include <cstdio>
2 #include <iostream>
3 #include <cmath>
4 #include <algorithm>
5 #include <cstring>
6 using namespace std;
7
8 int p[ 1001][ 1001];
9 double ans[ 1001][ 1001];
10 int n,m;
11 int tot,g[ 1001],nnext[ 1001* 2],num[ 1001* 2];
12 int team[ 1000000];
13 int d[ 1001];
14 int head,tail;
15 bool b[ 1001];
16 void Add( int x, int y)
17 {
18     tot++;
19     nnext[tot]=g[x];
20     g[x]=tot;
21     num[tot]=y;
22 }
23 void SPFA( int t)
24 {
25      while(head<tail)
26     {
27          int x=team[++head]; b[x]= false;
28          for( int i=g[x];i;i=nnext[i])
29         {
30              int tmp=num[i];
31              if(d[tmp]>d[x]+ 1||(d[tmp]==d[x]+ 1&&x<p[tmp][t]))
32             {
33                 d[tmp]=d[x]+ 1;
34                 p[tmp][t]=x;
35                  if(!b[tmp])
36                 {
37                     team[++tail]=tmp;
38                     b[tmp]= true;
39                 }
40             }
41         }
42     }
43 }
44
45 double DFS( int x, int y)
46 {
47 //     cout<< x<<' '<<y<<endl;
48
49      if(x==y) return 0.0;
50      if(ans[x][y]) return ans[x][y];
51      if(p[x][y]==y) return 1.0; if(p[p[x][y]][y]==y) return 1.0;
52
53      double aa= 0;
54      int cnt= 0;
55      for( int i=g[y];i;i=nnext[i])
56         aa+=DFS(p[p[x][y]][y],num[i])+ 1.0,cnt++;
57     aa+=DFS(p[p[x][y]][y],y)+ 1.0;
58     aa/=( double)(cnt+ 1);
59     ans[x][y]=aa;
60      return ans[x][y];
61 }
62
63 int main()
64 {
65     cin>>n>>m; int s,t;cin>>s>>t;
66      for( int i= 1;i<=m;i++)
67     {
68          int x,y;
69         cin>>x>>y;
70         Add(x,y);
71         Add(y,x);
72     }
73
74 //     cout<<"fasfa";
75
76      for( int i= 1;i<=n;i++)
77     {
78         memset(d, 63, sizeof(d));
79         d[i]= 0;
80         head=tail= 0;
81         team[++tail]=i;b[i]= true;
82         SPFA(i);
83     }
84
85      //     cout<<endl;
86     printf( " %.3lf\n ",DFS(s,t));
87      return 0;
88 }