BZOJ3143 游走 [期望][高斯消元]

本人蒟蒻，反正本只会用高斯消元做这道题，求出每个情况的系数，上期望算法乱搞就行了。
注意高斯消元的for循环，起止点很乱，最好手动推一推。

#include<bits/stdc++.h>using namespace std;const int N=505,M=250500;const int eps=1e-9;int mp[N][N],d[N],a[M],b[M],m,n;double C[N][N],val[M],num[N],ans;template<class T>inline void read(T &res){    static char ch;T flag=1;    while((ch=getchar())<'0'||ch>'9')if(ch=='-')flag=-1;res=ch-48;    while((ch=getchar())>='0'&&ch<='9')res=res*10+ch-48;res*=flag;}double ABS(double x){return x<0?-x:x;}void swap(double &x,double &y){double t=x;x=y;y=t;}void Gauss(){    for(register int i=1;i<=n-1;i++){        int tmp=i;        for(register int j=i+1;j<=n;++j)if(ABS(C[j][i])>ABS(C[tmp][i])+eps)tmp=j;        if(tmp-i)for(register int j=i;j<=n+1;j++)swap(C[i][j],C[tmp][j]);        for(register int j=i+1;j<=n;j++)            for(register int l=n+1;l>=i;l--)                C[j][l]=C[i][l]*C[j][i]/C[i][i]-C[j][l];    }    for(register int i=n;i>=1;i--){        num[i]=C[i][n+1];        for(register int j=i+1;j<=n;j++)num[i]-=num[j]*C[i][j];        (C[i][i])?num[i]/=C[i][i]:num[i]=0;    }}int main(){    read(n),read(m);    for(register int i=1;i<=m;i++){        read(a[i]),read(b[i]);        mp[a[i]][b[i]]++,mp[b[i]][a[i]]++,d[a[i]]++,d[b[i]]++;    }    for(register int i=1;i<=n-2;++i)        for(register int j=1;j<=n-1;++j)            C[i][j]=(i-j)?(1.0*mp[i][j])/d[j]:-1;    for(register int i=1;i<=n-2;++i)C[i][n]=(i-1)?0:-1;    C[n-1][n]=1;    for(register int i=1;i<=n-1;i++)C[n-1][i]=1.0*mp[n][i]/d[i];    num[n--]=0;    Gauss();    for(register int i=1;i<=m;i++)      val[i]=num[a[i]]/(1.0*d[a[i]])+num[b[i]]/(1.0*d[b[i]]);    sort(val+1,val+1+m);    for(register int i=1;i<=m;i++)ans+=val[i]*(m-i+1);    printf("%.3lf\n",ans+eps);    return 0;}