bzoj 3143(期望与DP+高斯消元)

3143: [Hnoi2013]游走

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 2156   Solved: 951
[ Submit][ Status][ Discuss]

Description

一个无向连通图，顶点从1编号到N，边从1编号到M。 
小Z在该图上进行随机游走，初始时小Z在1号顶点，每一步小Z以相等的概率随机选 择当前顶点的某条边，沿着这条边走到下一个顶点，获得等于这条边的编号的分数。当小Z 到达N号顶点时游走结束，总分为所有获得的分数之和。 
现在，请你对这M条边进行编号，使得小Z获得的总分的期望值最小。

Input

第一行是正整数N和M，分别表示该图的顶点数 和边数，接下来M行每行是整数u，v(1≤u,v≤N)，表示顶点u与顶点v之间存在一条边。 输入保证30%的数据满足N≤10，100%的数据满足2≤N≤500且是一个无向简单连通图。

Output

仅包含一个实数，表示最小的期望值，保留3位小数。

Sample Input

3 3
2 3
1 2
1 3

Sample Output

3.333

解题思路：首先肯定是期望次数最少的编号最大，然后就是求每条

边的期望次数，可以转化到求点的期望次数。高斯消元。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
int n,m,len;
long double f[510][510];
int to[500000],next[500000],h[500000];
int xg[500000],yg[500000];
long double zhi[500000]; 
long double du[500000]; 
 
inline int read()
{
char y; int x=0,f=1; y=getchar();
while (y<'0' || y>'9') {if (y=='-') f=-1; y=getchar();}
while (y>='0' && y<='9') {x=x*10+int(y)-48; y=getchar();}
return x*f;
}


void insert(int x,int y)
 {
  ++len; to[len]=y; next[len]=h[x]; h[x]=len;
 }


void work()
 {
  for (int i=1;i<=n;++i)
  {
  double ogg=f[i][i];
  for (int j=i;j<=n+1;++j)
  f[i][j]/=ogg;
  for (int j=i+1;j<=n;++j)
  {
    ogg=f[j][i];
    for (int k=i;k<=n+1;++k)
     {
       f[j][k]-=ogg*f[i][k];
}
}
 }
for (int i=n-1;i>=1;--i)
{
for (int j=i+1;j<=n;++j)
{
double ogg=f[i][j];
f[i][n+1]-=ogg*f[j][n+1];
f[i][j]=0;
 }
}
 }


int main()
{
n=read(); m=read();
for (int i=1;i<=m;++i)
{
int x,y; x=read(); y=read(); ++du[x]; ++du[y];
insert(x,y); insert(y,x);
xg[i]=x; yg[i]=y;
}
for (int i=1;i<=n;++i)
{
int u=h[i];
f[i][i]=1;
if (i==1) f[i][n+1]=1;
        while (u!=0)
         {
          if (to[u]!=n)
           {
            f[i][to[u]]=-1/du[to[u]];
  }
          u=next[u];
}
     }
    work();
    f[n][n+1]=0;
for (int i=1;i<=m;++i)
{
int x=xg[i]; int y=yg[i];
zhi[i]=f[x][n+1]/du[x]+f[y][n+1]/du[y];
}
long double ans=0;
sort(zhi+1,zhi+m+1);
for (int i=m;i>=1;--i)
{
ans+=zhi[i]*(m-i+1);
}
printf("%.3Lf",ans);
}