【BZOJ 2878】 [Noi2012]迷失游乐园

Description

放假了，小Z觉得呆在家里特别无聊，于是决定一个人去游乐园玩。进入游乐园后，小Z看了看游乐园的地图，发现可以将游乐园抽象成有n个景点、m条道路的无向连通图，且该图中至多有一个环（即m只可能等于n或者n-1）。小Z现在所在的大门也正好是一个景点。小Z不知道什么好玩，于是他决定，从当前位置出发，每次随机去一个和当前景点有道路相连的景点，并且同一个景点不去两次（包括起始景点）。贪玩的小Z会一直游玩，直到当前景点的相邻景点都已经访问过为止。小Z所有经过的景点按顺序构成一条非重复路径，他想知道这条路径的期望长度是多少？小Z把游乐园的抽象地图画下来带回了家，可是忘了标哪个点是大门，他只好假设每个景点都可能是大门（即每个景点作为起始点的概率是一样的）。同时，他每次在选择下一个景点时会等概率地随机选择一个还没去过的相邻景点。

Input

第一行是两个整数n和m，分别表示景点数和道路数。 接下来行，每行三个整数Xi, Yi, Wi，分别表示第i条路径的两个景点为Xi, Yi，路径长Wi。所有景点的编号从1至n，两个景点之间至多只有一条道路。

Output

共一行，包含一个实数，即路径的期望长度，保留五位小数

Sample Input

4 3
1 2 3
2 3 1
3 4 4

Sample Output

6.00000

【样例解释】样例数据中共有6条不同的路径: 路径 长度 概率
1-->4 8 1/4
2-->1 3 1/8
2-->4 5 1/8
3-->1 4 1/8
3-->4 4 1/8
4-->1 8 1/4
因此期望长度 = 8/4 + 3/8 + 5/8 + 4/8 + 4/8 + 8/4 = 6.00
【评分方法】本题没有部分分，你程序的输出只有和标准答案的差距不超过0.01时，才能获得该测试点的满分，否则不得分。
【数据规模和约定】对于100%的数据，1 <= Wi <= 100。 测试点编号 n m 备注
1 n=10 m = n-1 保证图是链状
2 n=100 只有节点1的度数大于2
3 n=1000 /
4 n=100000 /
5 n=100000 /
6 n=10 m = n /
7 n=100 环中节点个数<=5
8 n=1000 环中节点个数<=10
9 n=100000 环中节点个数<=15
10 n=100000 环中节点个数<=20

 

这是我迄今为止做过的最麻烦的一道树形DP，还是那句话，我的DP就是弱，你打我啊

对于树来说

利用F[i]表示每一个点到达它的子树的期望长度

利用D[i]=Σ（F[j]+e[i,j]）(j为i的儿子)

利用du[i]表示一个点的度数

则可求得F[i]=D[i]/du[i]

如此进行一遍dfs即可求出以一个点为根的数值

此时我们再进行一遍dfs，求出以任一点为根的D数组

主要利用D[j]+=(D[i]-f[j]-e[i,j])/(du[i]-1)+e[i,j](j为i的儿子)

这样就可以解出所有的D值

至于M==N的情况，我们可以暴力处理环，毕竟环只有最大20

首先是判断去掉环的情况，那么图就变成了无根森林，和树的处理方式相同，进行一遍dfs

再是找出环，利用gg[i]数组表示从i点出发的期望长度，g数组表示从i点出发单向走的长度

gg[i]=g[j]+g[k]+e[i,j]+e[j,k](i,k为i左右相连的边)

暴力求解出环上的所有gg

最后环上点d[i]+=gg[i]，du[i]+=2;

 #include<cstdio>
 #define ld long double
 #include<cstring>
 #include<iostream>
 using namespace std;
 const int N=;
 struct ee{int to,next,w;}e[N*];
 int n,m,u,v,w,cnt,root,timer;
 int head[N],du[N],vis[N],fa[N];
 ld f[N],gg[N],g[N],d[N];
 bool c[N];
 void ins(int u,int v,int w){
     e[++cnt].to=v;e[cnt].next=head[u];e[cnt].w=w;head[u]=cnt;
     e[++cnt].to=u;e[cnt].next=head[v];e[cnt].w=w;head[v]=cnt;
 }
 void dfs1(int x){
     vis[x]=;
     for (int i=head[x];i;i=e[i].next){
         int v=e[i].to;
         if (!vis[v]&&!c[v]){
             du[x]++;
             dfs1(v);
             d[x]+=f[v]+e[i].w;
         }
     }
     if (du[x]) f[x]=(d[x]+0.0)/du[x];
     if(x!=root) du[x]++;
 }

 void dfs2(int x){
     vis[x]=;
     for(int i=head[x];i;i=e[i].next){
         int v=e[i].to;
         if(!vis[v]&&!c[v]){
             d[v]+=(d[x]-f[v]-e[i].w)/max(,(du[x]-))+e[i].w;
             dfs2(v);
         }
     }
 }

 void find(int x){
     vis[x]=++timer;
     for (int i=head[x];i;i=e[i].next){
         int v=e[i].to;
         if(!vis[v]) fa[v]=x,find(v);
         else if(fa[x]!=v&&vis[v]<vis[x]){
             c[v]=;
             for(;x!=v;)
             c[x]=,x=fa[x];
         }
     }
 }

 void solve_cycle(int x,int fa){
     bool flag=;
     g[x]=;//1 wa
     for (int i=head[x];i;i=e[i].next){
         int v=e[i].to;
         if(v!=root&&v!=fa&&c[v]){
             flag=;
             solve_cycle(v,x);
             g[x]+=g[v]+e[i].w;
         }
     }
     if(x==root) return;
     int k=du[x];if(!k)k++;
     if(flag)g[x]=(d[x]+0.0)/k;
     else k=du[x]+,g[x]=(g[x]+d[x])/k;
 }

 int main(){
     scanf("%d%d",&n,&m);
     for (int i=;i<=m;i++){
         scanf("%d%d%d",&u,&v,&w);
         ins(u,v,w);
     }
     if(m==n-){
         root=;
         dfs1();
         memset(vis,,sizeof(vis));
         dfs2();
     }else{
         find();
         memset(vis,,sizeof(vis));
         for (int i=;i<=n;i++) if(c[i]) {
             root=i;
             dfs1(i);
         }
         for (int i=;i<=n;i++) if(c[i]) {
             root=i;
             solve_cycle(i,);
             gg[i]=g[i];
         }
         memset(vis,,sizeof(vis));
         for(int i=;i<=n;i++) if(c[i]){
             du[i]+=;
             d[i]+=gg[i];
         }
         for(int i=;i<=n;i++) if(c[i]) dfs2(i);
     }
     double ans=;
     for (int i=;i<=n;i++) ans+=(d[i]+0.0)/du[i];
     printf("%.5lf\n",ans/n);
 }