比较经典的题,题解看网上的。。https://www.cnblogs.com/GXZlegend/p/7054536.html

自己sort弄错了。。还以为是高斯消元写歪了。。

#include<bits/stdc++.h>
using namespace std;const int maxn = 505;
const double esp = 1e-10;struct Edge{int u,v;double E;}e[maxn*maxn];
int mp[maxn][maxn],d[maxn],n,m;
double E[maxn][maxn],b[maxn];
int cmp(Edge a,Edge b){return a.E>b.E;}void guass(){for(int i=1;i<=n;i++){int maxx=i;for(int j=i;j<=n;j++){if(fabs(E[j][i])>esp&&fabs(E[j][i])>fabs(E[maxx][i]))maxx=j;}if(maxx!=i){swap(E[maxx],E[i]);swap(b[maxx],b[i]);}if(fabs(E[i][i])<esp)continue;for(int j=i+1;j<=n;j++){if(fabs(E[j][i])<esp)continue;double rate=E[j][i]/E[i][i];for(int k=i;k<=n;k++)E[j][k]-=rate*E[i][k];b[j]-=rate*b[i];}}for(int i=n;i>=1;i--){if(fabs(E[i][i])<esp)continue;for(int j=i+1;j<=n;j++)b[i]-=E[i][j]*b[j];b[i]/=E[i][i];}
}int main(){scanf("%d%d",&n,&m);for(int i=1;i<=m;i++){int u,v;scanf("%d%d",&u,&v);mp[u][v]=mp[v][u]=1;d[u]++;d[v]++;e[i].u=u;e[i].v=v;}//建立矩阵E[n][n]=b[1]=1;for(int i=1;i<n;i++){E[i][i]=1;for(int j=1;j<=n;j++)if(mp[i][j])E[i][j]-=1.0/d[j];}guass();for(int i=1;i<=m;i++){int u=e[i].u,v=e[i].v;e[i].E=b[u]/d[u] + b[v]/d[v];}sort(e+1,e+1+m,cmp);double ans=0;for(int i=1;i<=m;i++)ans+=e[i].E*i;printf("%.3lf\n",ans);
}

转载于:https://www.cnblogs.com/zsben991126/p/11058012.html

期望dp+高斯消元——bzoj3143相关推荐

  1. LightOJ 1151 Snakes and Ladders (期望DP + 高斯消元)

    Description 'Snakes and Ladders' or 'Shap-Ludu' is a game commonly played in Bangladesh. The game is ...

  2. ICPC 2005 hangzhou Generator (UVA1358)KMP + 期望DP / 高斯消元

    整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 Generator Weblink https://www.luogu.com.cn/problem/ ...

  3. P4457-[BJOI2018]治疗之雨【期望dp,高斯消元】

    正题 题目链接:https://www.luogu.com.cn/problem/P4457 题目大意 开始一个人最大生命值为nnn,剩余hphphp点生命,然后每个时刻如果生命值没有满那么有1m+1 ...

  4. 【BZOJ1778】[Usaco2010 Hol]Dotp 驱逐猪猡 期望DP+高斯消元

    [BZOJ1778][Usaco2010 Hol]Dotp 驱逐猪猡 Description 奶牛们建立了一个随机化的臭气炸弹来驱逐猪猡.猪猡的文明包含1到N (2 <= N <= 300 ...

  5. 【BZOJ2337】XOR和路径,概率期望DP+高斯消元

    Time:2016.08.27 Author:xiaoyimi 转载注明出处谢谢 传送门 思路: 与游走思路有一定相似的地方 对答案的每一位进行判断 通过高斯消元解出每个点到n xor路径为1的概率 ...

  6. BZOJ 1778 Usaco2010 Hol Dotp 驱逐猪猡 期望DP+高斯消元

    题目大意:给定一个无向图,炸弹从1号节点出发,每个时刻有P/Q的概率爆炸,如果某个时刻没有爆炸,就会等概率沿着随机一条出边走到下一个城市,求最终每个城市的爆炸概率 #include <cmath ...

  7. 【BZOJ3640】JC的小苹果 概率DP+高斯消元

    [BZOJ3640]JC的小苹果 Description 让我们继续JC和DZY的故事. "你是我的小丫小苹果,怎么爱你都不嫌多!" "点亮我生命的火,火火火火火!&qu ...

  8. BZOJ 3270: 博物馆 1778: 驱逐猪猡 【概率DP+高斯消元】

    题目描述: 中文题面,不多解释.1778传送门 3270 传送门 (博物馆)题目分析: 也许很多人做概率题的时候都有种虚幻感..感觉莫名其妙就得出一个期望.概率,一知半解... 所以我在这里仔细地剖析 ...

  9. BZOJ 3270: 博物馆 [概率DP 高斯消元]

    http://www.lydsy.com/JudgeOnline/problem.php?id=3270 题意:一张无向图,一开始两人分别在$x$和$y$,每一分钟在点$i$不走的概率为$p[i]$, ...

  10. BZOJ_1778_[Usaco2010 Hol]Dotp 驱逐猪猡_概率DP+高斯消元

    BZOJ_1778_[Usaco2010 Hol]Dotp 驱逐猪猡_概率DP+高斯消元 题意: 奶牛们建立了一个随机化的臭气炸弹来驱逐猪猡.猪猡的文明包含1到N (2 <= N <= 3 ...

最新文章

  1. 网易云音乐的消息队列改造之路
  2. XGBOOST带试验源码
  3. c# 指定打开某个路径下的CMD_(win10下sublime通过配置JSON调用MATLAB直接运行程序)(转载)...
  4. 【博客】csdn搬家到wordpress
  5. C#面试题——附答案
  6. 十大常用个人linux,Linux—常用十大命令
  7. PHP笔记-平润年判断例子
  8. 赛玛共享按摩椅前端蓝色UI+分销返利+去除短信注册验证
  9. [LeetCode] Majority Element II
  10. composer 安装dev包_「PHP编程」如何搭建私有Composer包仓库?
  11. Tensorflow 实战 Google 深度学习框架(第2版)---- 10.2.2节 P272 代码
  12. 系统服务器ping偶尔超时,云服务器ping超时原因
  13. LINUX下,ffmpeg增加NVIDIA硬件编解码的步骤及解决办法
  14. 维智WSD-A2系列伺服驱动器用户手册(MECHATROLINK-Ⅱ总线通信型)
  15. python如何求平方_python中求平方
  16. 学妹跑过来问我为啥Xshell 打不开了,让我帮她处理下【手把手讲解】
  17. Appium+夜神模拟器模拟人为操作
  18. 开源财务会计软件(搬运)
  19. 《Web安全渗透全套教程(40集)》学习笔记 | SQL注入攻击及防御
  20. 智联招聘如何设置公司屏蔽

热门文章

  1. 拓端tecdat|电力消耗模型构建、分析和预测
  2. 拓端tecdat|R语言中的多类别问题的绩效衡量:F1-score 和广义AUC
  3. Java EE实战教程 servlet (一)
  4. android 找不到类文件,Android库:当“实现项目”用于库的模块依赖时,找不到类文件...
  5. MatCaffe问题: make matcaffe 'matlab/+caffe/private/caffe_.mexa64' failed
  6. Caffe学习:使用pycaffe读取mean.binaryproto文件参数
  7. Ubuntu16.04 Jupyter安装(Ipython Notebook安装)
  8. Tf-idf推导及理解
  9. 用java判断x奇或偶_改善java程序——用偶判断,不用奇判断
  10. python中不同类型的数据不能相互运算_python语法入门(数据类型、运算)