首先可以各位分开求和

定义$f(i)$表示从i到n的期望值,然后经过一些常识,发现$f(n)=1$的时候的转移,然后直接转移,也可以找到$f(n)=0$的转移。

然后高斯消元31次就可以了。

#include <map>
#include <ctime>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)
#define maxn 20005 int h[maxn],to[maxn],ne[maxn],du[maxn],w[maxn],en=0,n,m;
double a[105][105],ans[105]; void add(int a,int b,int c)
{du[a]++;to[en]=b;ne[en]=h[a];w[en]=c;h[a]=en++;} void Gauss()
{
F(i,1,n-1)
{
int tmp=i;
F(j,i,n-1)if(fabs(a[j][i])>fabs(a[tmp][i]))tmp=j;
F(j,i,n+1) swap(a[i][j],a[tmp][j]);
F(j,i+1,n-1)
{
double tmp=a[j][i]/a[i][i];
F(k,i,n+1) a[j][k]-=a[i][k]*tmp;
}
}
D(i,n-1,1)
{
F(j,i+1,n)
a[i][n+1]-=a[i][j]*ans[j];
ans[i]=a[i][n+1]/a[i][i];
}
} double ret=0.0; void solve(int x)
{
memset(a,0,sizeof a);
memset(ans,0,sizeof ans);
F(i,1,n-1)
{
for (int j=h[i];j>=0;j=ne[j])
if ((w[j]>>x)&1)
a[i][to[j]]+=1,a[i][n+1]+=1;
else a[i][to[j]]-=1;
a[i][i]+=du[i];
}
// F(i,1,n) F(j,1,n+1) printf("%.6f%c",a[i][j],j==n+1?'\n':' ');printf("\n\n");
Gauss();
ret+=ans[1]*(1<<x);
} int main()
{
memset(h,-1,sizeof h);
scanf("%d%d",&n,&m);
F(i,1,m)
{
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
add(a,b,c);if (a!=b) add(b,a,c);
}
D(i,30,0) solve(i);
printf("%.3f\n",ret);
}

  

最新文章

  1. js_多个引号的用法
  2. Js设置及获取Cookie的方法
  3. java中的包以及内部类的介绍
  4. Docker私有仓库 Registry中的镜像管理
  5. bzoj4080
  6. appium desktop 版本发布
  7. struts2+hibernate+spring注解版框架搭建以及简单测试(方便脑补)
  8. 鼠标跟随效果 vue或者js通用
  9. java 调用 api接口
  10. git 入门教程之协同开发
  11. Java中构造函数传参数在基本数据类型和引用类型之间的区别
  12. 关于Unity中的小案例之运动的小船以及摄像机跟随技术(专题五)
  13. android 系统广播
  14. POJ2506:Tiling(递推+大数斐波那契)
  15. java利用poi读取excel异常问题
  16. Codeforces 954H Path Counting 【DP计数】*
  17. 浅谈对MVC的理解
  18. 2018牛客多校2 - J farm 随机乱搞/二进制分组
  19. 【JS设计模式】装饰者模式
  20. Django 之 JsonResponse 对象

热门文章

  1. Azure资源管理工具Azure PowerShell介绍
  2. Resize a UIImage the right way
  3. html备忘录
  4. RMQ求区间最大最小值
  5. (5)JSTL的xml标签库
  6. OpenCV Haar AdaBoost源码改进据说是比EMCV快6倍
  7. Educational Codeforces Round 40千名记
  8. (54)zabbix链接及解除模板链接
  9. shell进阶
  10. Perl学习之四:语句