bzoj2337

为什么一道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);
}