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

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

#include<bits/stdc++.h>
using namespace std; const int maxn = ;
const double esp = 1e-; 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=;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+;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>=;i--){
if(fabs(E[i][i])<esp)continue;
for(int j=i+;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=;i<=m;i++){
int u,v;
scanf("%d%d",&u,&v);
mp[u][v]=mp[v][u]=;
d[u]++;d[v]++;
e[i].u=u;e[i].v=v;
} //建立矩阵
E[n][n]=b[]=;
for(int i=;i<n;i++){
E[i][i]=;
for(int j=;j<=n;j++)
if(mp[i][j])E[i][j]-=1.0/d[j];
}
guass(); for(int i=;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+,e++m,cmp);
double ans=;
for(int i=;i<=m;i++)
ans+=e[i].E*i;
printf("%.3lf\n",ans);
}

最新文章

  1. BZOJ 4576: [Usaco2016 Open]262144
  2. php 教程列表
  3. 交换芯片收发包的 DMA 实现原理
  4. RequireJS 循环依赖报 模块undefined 处理方案
  5. Linux 下安装 jdk步骤:
  6. dirname和basename命令
  7. jqmobi 子菜单 高亮效果
  8. cc.RepeatForever和cc.Spawn冲突
  9. max Sum(简单动态规划)
  10. JDNI
  11. 关于Webapp的注意事项
  12. springMVC和spring的集成
  13. Hadoop Mapreduce中wordcount 过程解析
  14. 在n个任意不相同的数中,输出r个数的组合,并且n和r由键盘输入。
  15. Django ORM中常用字段和参数
  16. python+selenium+unnitest写一个完整的登陆的验证
  17. Git(介绍和安装)
  18. Win10系列:C#应用控件进阶1
  19. Log4j 把不同包的日志打印到不同位置
  20. elasticsearch 请求体查询方式整理

热门文章

  1. 单源最短路径问题1 (Bellman-Ford算法)
  2. vue-为什么子组件中的data选项必须是函数?
  3. 配置Tomcat-8.5 JVM内存参数
  4. ncurses库的介绍与安装
  5. String类的substring()方法
  6. 解决ios10及以上Safari双击和双指缩放无法禁止的问题
  7. MyBatis mappers元素标签及其属性、配置
  8. 解决:Map的area属性标签鼠标Hover可以给area加背景
  9. Parallels Desktop Centos 设置IP
  10. websokect的原理