为什么一道SB题纠结了这么久。。。
设f[i]为i到n的1的概率(或期望)
由于正着推会出锅(f[n]=0)所以要倒着推
则由i到它的son都是等概率的,即1/degree[i]
PS:为什么bzoj1877系数又是1/degree[son]?是因为一个是由i到son,一个是由son到i啊。。。
列方程即可
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define maxn 110
#define maxm 21011
using namespace std;
struct edge{
int r,nxt,w;
}e[maxm];
double ans,f[maxn][maxn],res,r[maxn];
int n,m,esz,head[maxn];
double d[maxn];
void addedge(int u,int v,int w){
e[++esz].r=v;e[esz].nxt=head[u];
d[u]++;head[u]=esz;e[esz].w=w;
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1,u,v,w;i<=m;++i){
scanf("%d%d%d",&u,&v,&w);
addedge(u,v,w);
if(u!=v)addedge(v,u,w);
}
for(int p=0;p<=30;++p){
for(int i=1;i<=n;++i)
for(int j=1;j<=n+1;++j)f[i][j]=0;
for(int i=1;i<=n;++i)f[i][i]=1;
for(int u=1;u<n;++u){
for(int t=head[u];t;t=e[t].nxt){
int v=e[t].r;
if((e[t].w>>p)&1)f[u][v]+=1/d[u],f[u][n+1]+=1/d[u];
else f[u][v]-=1/d[u];
}
}
/* for(int i=1;i<=n;++i,putchar(10))
for(int j=1;j<=n+1;j++)
printf("%.2lf ",f[i][j]);
printf("---------------\n");
*/ for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j)if(i!=j){
double t=f[j][i]/f[i][i];
for(int k=1;k<=n+1;++k)
f[j][k]-=f[i][k]*t;
}
}
// printf("[%.3lf]",f[1][n+1]/f[1][1]);
ans+=(f[1][n+1]/f[1][1])*(1<<p);
}
printf("%.3lf",ans);
}