Description###

一个无向连通图,顶点从1编号到N,边从1编号到M。 小Z在该图上进行随机游走,初始时小Z在1号顶点,每一步小Z以相等的概率随机选 择当前顶点的某条边,沿着这条边走到下一个顶点,获得等于这条边的编号的分数。当小Z 到达N号顶点时游走结束,总分为所有获得的分数之和。 现在,请你对这M条边进行编号,使得小Z获得的总分的期望值最小。

Input###

第一行是正整数N和M,分别表示该图的顶点数 和边数,接下来M行每行是整数u,v(1<=u,v<=N),表示顶点u与顶点v之间存在一条边。 输入保证30%的数据满足N<=10,100%的数据满足2<=N<=500且是一个无向简单连通图。

Output###

仅包含一个实数,表示最小的期望值,保留3位小数。

Sample Input###

3 3

2 3

1 2

1 3

Sample Output###

3.333

HINT###

边(1,2)编号为1,边(1,3)编号2,边(2,3)编号为3。


想法##

果断高斯消元解期望方程啊。

如果我们知道了每条边期望被经过的次数,那么贪心按期望值从小到大、编号从大到小编就好了

但是每条边期望被经过的次数并不好搞,每个点的期望经过次数是可以搞出来的

列方程组的时候注意,1号点最初的期望值就有1,n号点只能进不能出

假设某条边两端的点为u,v,每个点期望被经过次数为e[i],每个点的度数为d[i]

则这条边期望被经过的次数为 \(e[u]/d[u] + e[v]/d[v]\)

注意要特判一下,如果u、v中某一个为n的话,这条边只能从那个不是n的点过来,所以被经过的期望为\(e[u]/d[u]\)或\(e[v]/d[v]\)

然后排个序贪心就行啦


代码##

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath> #define eps 0.001 using namespace std; const int N = 505;
typedef double db[N][N]; int n,m;
void gauss(db A){
for(int i=1;i<=n;i++){
int r=i;
for(int j=i+1;j<=n;j++)
if(fabs(A[j][i])-fabs(A[r][i])>eps) r=j;
for(int j=i;j<=n;j++) swap(A[i][j],A[r][j]);
for(int j=i+1;j<=n;j++){
double f=A[j][i]/A[i][i];
for(int k=i;k<=n+1;k++) A[j][k]-=f*A[i][k];
}
}
for(int i=n;i>0;i--){
for(int j=i+1;j<=n;j++)
A[i][n+1]-=A[i][j]*A[j][n+1];
A[i][n+1]/=A[i][i];
}
} db a;
struct edge{
int u,v;
double f;
edge(int u=0,int v=0,double f=0.0):u(u),v(v),f(f) {}
bool operator < (const edge &b) const { return f>b.f; }
}d[N*N/2];
int s[N]; int main()
{
int x,y;
scanf("%d%d",&n,&m);
for(int i=0;i<m;i++){
scanf("%d%d",&x,&y);
if(x>y) swap(x,y);
d[i]=edge(x,y);
if(y==n) a[y][x]=1.0;
else a[x][y]=a[y][x]=1.0;
s[x]++; s[y]++;
}
for(int i=1;i<=n;i++) {
a[i][i]=1.0;
if(i==1) a[i][n+1]=1.0;
for(int j=1;j<=n;j++){
if(j==i) continue;
if(a[i][j]) a[i][j]*=-1.0/s[j];
}
} gauss(a);
for(int i=0;i<m;i++){
if(d[i].v==n) d[i].f=a[d[i].u][n+1]*1.0/s[d[i].u];
else d[i].f=a[d[i].u][n+1]*1.0/s[d[i].u]+a[d[i].v][n+1]*1.0/s[d[i].v];
}
sort(d,d+m); double ans=0.0;
for(int i=0;i<m;i++) ans+=d[i].f*(i+1);
printf("%.3lf\n",ans); return 0;
}

[bzoj3143] [洛谷P3232] [HNOI2013] 游走的更多相关文章

  1. 洛谷P3232[HNOI2013]游走

    有一个无向简单连通图,顶点从 \(1\) 编号到 \(n\),边从 \(1\) 编号到 \(m\) 小Z在该图上进行随机游走,初始时小Z在\(1\)号顶点,每一步小Z以相等的概率随机选 择当前顶点的某 ...

  2. 洛谷 P3232 [HNOI2013]游走

    链接: P3232 题意: 和上次考试 T4 的简化且无修改一样,经典图上高斯消元求期望. 分析: 要求出每个点的期望出发次数 \(f_i\),每个点度数为 \(d_i\),有 \[f1=\sum\d ...

  3. 洛谷P3232 [HNOI2013]游走(高斯消元+期望)

    传送门 所以说我讨厌数学……期望不会高斯消元也不会……好不容易抄好了高斯消元板子被精度卡成琪露诺了…… 首先,我们先算出走每一条边的期望次数,那么为了最小化期望,就让大的期望次数乘上小编号 边的期望次 ...

  4. 题解 P3232 [HNOI2013]游走

    洛谷P3232[NOI2013]游走 题目描述 给定一个 n 个点 m 条边的无向连通图,顶点从 1 编号到 n,边从 1 编号到 m. 小 Z 在该图上进行随机游走,初始时小 Z 在 1 号顶点,每 ...

  5. P3232 [HNOI2013]游走 解题报告

    P3232 [HNOI2013]游走 题目描述 一个无向连通图,顶点从\(1\)编号到\(N\),边从\(1\)编号到\(M\). 小Z在该图上进行随机游走,初始时小Z在1号顶点,每一步小Z以相等的概 ...

  6. P3232 [HNOI2013]游走——无向连通图&&高斯消元

    题意 一个无向连通图,顶点从1编号到N,边从1编号到M. 小Z在该图上进行随机游走,初始时小Z在1号顶点,每一步小Z以相等的概率随机选 择当前顶点的某条边,沿着这条边走到下一个顶点,获得等于这条边的编 ...

  7. P3232 [HNOI2013]游走

    吐槽 傻了傻了,对着题解改了好长时间最后发现是自己忘了调用高斯消元了... 思路 期望题,分配编号,显然编号大的分给贡献次数小的,所以需要知道每个边被经过次数的期望 然后边被经过的次数的期望就是连接的 ...

  8. BZOJ 3143 Luogu P3232 [HNOI2013]游走 (DP、高斯消元)

    题目链接: (bzoj) https://www.lydsy.com/JudgeOnline/problem.php?id=3143 (luogu) https://www.luogu.org/pro ...

  9. 【BZOJ3143】[Hnoi2013]游走 期望DP+高斯消元

    [BZOJ3143][Hnoi2013]游走 Description 一个无向连通图,顶点从1编号到N,边从1编号到M. 小Z在该图上进行随机游走,初始时小Z在1号顶点,每一步小Z以相等的概率随机选 ...

随机推荐

  1. java的四种代码块

    用{}括起来的称为代码块: 普通代码块:类中方法的方法体 构造代码块:类中{}直接括起来的语句,每次创建对象都会被调用,先于构造函数执行 静态代码块:类中static{}括起来的语句,只执行一次,先于 ...

  2. [板子]Kruskal

    众所周知求最小生成树的两种方法: 1.Kruskal 2.Prim 这里只挂第一种,因为noip掌握第一种就够了. 两种做法的区别可以参考这个博客:http://blog.csdn.net/molln ...

  3. Java并发机制的底层实现原理之volatile应用,初学者误看!

    volatile的介绍: Java代码在编译后会变成Java字节码,字节码被类加载器加载到JVM里,JVM执行字节码,最终需要转化为汇编指令在CPU上执行,Java中所使用的并发机制依赖于JVM的实现 ...

  4. nginx负载均衡的几种模式

    nginx 的 upstream目前支持 4 种方式的分配 ).轮询(默认) 每个请求按时间顺序逐一分配到不同的后端服务器,如果后端服务器down掉,能自动剔除. ).weight 指定轮询几率,we ...

  5. 【Kubernetes】容器集群管理常用命令笔记

    一.集群部署-查询集群状态 ①查询k8s master各组件健康状态: kubectl get componentstatus ②查询k8s node健康状态: kubectl get node 二. ...

  6. Channel 9视频整理【3】

    Will 保哥 微软mvp https://channel9.msdn.com/Niners/Will_Huang 繁体中文视频 Visual Studio 2017 新功能探索 https://ch ...

  7. 用mingw静态编译Qt4.8.2和Qt5.1.1(需要修改不少源码)

    因为一些乱七八糟的原因,我需要用mingw静态编译Qt4.8.2和Qt5.1.1.经历了一天的折腾之后,自觉编译一下Qt还是件颇为麻烦的事情,故将过程略作总结,以备不时之需. 首先,在编译之前,我需要 ...

  8. Jenkins配置QQ邮箱发送邮件

    1.登陆QQ邮箱 2. 在“帐户”里开启“POP3/SMTP”并获取授权码 3. 发送短信验证验证后得到下面验证码 aeoygabszxfecbdj #验证吗 点击确定之后,服务已经开启 4. Jen ...

  9. C++简单实现Log日志类轻量级支持格式化输出变量

    CLog 头 代码很简单 如果需要的直接Ctrl+C  ----Ctrl+V 即可 #ifndef __CLOG__ #define __CLOG__ #include <windows.h&g ...

  10. vps批量管理,服务器批量管理,3389批量管理工具

    注册的第一天,把我的工具弄上来.如果有一个软件,如果你有vps虚拟机,服务器,3389中的任意一种两种或全部,而且还是批量的,需要管理.那么,你可能会想找个软件来管理吧,毕竟,一个人忙不过来管理这么多 ...