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

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 }