BZOJ 3143 游走(贪心+期望+高斯消元)

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

总分的期望值=每条边的期望经过次数*边的编号 之和。

不论我们如何编号，每条边的期望经过次数是不会变的，要使得边权和的期望最小，只需要贪心地使期望次数和边权倒序对应即可。

考虑如何求每条边的经过次数，记每个点度数为di，期望通过次数为xi,每条边期望通过次数为yi则

变形一下高斯消元即可。

# include <cstdio>
# include <cstring>
# include <cstdlib>
# include <iostream>
# include <vector>
# include <queue>
# include <stack>
# include <map>
# include <bitset>
# include <set>
# include <cmath>
# include <algorithm>
using namespace std;
# define lowbit(x) ((x)&(-x))
# define pi acos(-1.0)
# define eps 1e-
# define MOD
# define INF
# define mem(a,b) memset(a,b,sizeof(a))
# define FOR(i,a,n) for(int i=a; i<=n; ++i)
# define FDR(i,a,n) for(int i=a; i>=n; --i)
# define bug puts("H");
# define lch p<<,l,mid
# define rch p<<|,mid+,r
# define mp make_pair
# define pb push_back
typedef pair<int,int> PII;
typedef vector<int> VI;
# pragma comment(linker, "/STACK:1024000000,1024000000")
typedef long long LL;
inline int Scan() {
    int x=,f=; char ch=getchar();
    while(ch<''||ch>''){if(ch=='-') f=-; ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-''; ch=getchar();}
    return x*f;
}
inline void Out(int a) {
    if(a<) {putchar('-'); a=-a;}
    if(a>=) Out(a/);
    putchar(a%+'');
}
const int N=;
//Code begin...

struct Edge{int u, v; double w;}edge[N*N];
int dee[N], equ, var;
bool G[N][N];
double a[N][N], x[N];

int Guass(){
    int i, j, k, col, max_r;
    for (k=,col=; k<equ&&col<var; ++k,++col) {
        max_r=k;
        for (i=k+; i<equ; ++i) if (fabs(a[i][col])>fabs(a[max_r][col])) max_r=i;
        if (fabs(a[max_r][col])<eps) return ;
        if (k!=max_r) {
            for (j=col; j<var; ++j) swap(a[k][j],a[max_r][j]);
            swap(x[k],x[max_r]);
        }
        x[k]/=a[k][col];
        for (j=col+; j<var; ++j) a[k][j]/=a[k][col];
        a[k][col]=;
        for (i=; i<equ; ++i) if (i!=k) {
            x[i]-=x[k]*a[i][col];
            for (j=col+; j<var; ++j) a[i][j]-=a[k][j]*a[i][col];
            a[i][col]=;
        }
    }
    return ;
}
bool comp(Edge a, Edge b){return a.w>b.w;}
int main ()
{
    int n, m, u, v;
    scanf("%d%d",&n,&m);
    FOR(i,,m) scanf("%d%d",&u,&v), G[u][v]=G[v][u]=true, ++dee[u], ++dee[v], edge[i].u=u, edge[i].v=v;
    equ=var=n-;
    FOR(i,,n-) {
        if (i==) x[i-]=;
        a[i-][i-]=;
        FOR(j,,n-) {
            if (!G[i][j]) continue;
            a[i-][j-]=(-1.0)/dee[j];
        }
    }
    Guass();
    FOR(i,,m) {
        u=edge[i].u; v=edge[i].v;
        edge[i].w=x[u-]/dee[u]+x[v-]/dee[v];
    }
    sort(edge+,edge+m+,comp);
    double ans=;
    FOR(i,,m) ans+=edge[i].w*i;
    printf("%.3f\n",ans);
    return ;
}