bzoj4637: 期望
Description
Input
Output
一行一个整数,即满足总道路长度最小的情况下,设计方案的美学值期望。要求保留5位小数
按边权升序加边同时缩点,忽略缩点产生的自环,对同一权值且 加入此权值的边后在同个联通块内 的一组边,用矩阵树定理计算出生成树个数以及删去每条边后的生成树个数,于是可得一条边在最小生成树中的概率,统计答案。似乎要用long double才能过,但计算行列式时无论用long double还是模意义下的整数计算都能过。
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long double ld;
#define double ld
int _(){
int x=,c=getchar(),f=;
while(c<)c=='-'&&(f=-),c=getchar();
while(c>)x=x*+c-,c=getchar();
return x*f;
}
int es[],enx[],e0[],ep=,ev[];
struct edge{
int a,b,c,d;
void init(){
a=_();b=_();c=_();d=_();
}
void adde(){
es[ep]=b;enx[ep]=e0[a];ev[ep]=d;e0[a]=ep++;
es[ep]=a;enx[ep]=e0[b];ev[ep]=d;e0[b]=ep++;
}
}e[];
bool operator<(edge a,edge b){
return a.c<b.c;
}
double ans=;
int n,m,f[],id[],idp=,idt[],tk=,ed[],ID[],IDP;
int v[][],v1,ee[][],eep,os[][],op=;
int get(int x){
int a=x,c;
while(x!=f[x])x=f[x];
while(x!=f[a])c=f[a],f[a]=x,a=c;
return x;
}
void gid(int x){
if(idt[x]!=tk)idt[x]=tk,id[x]=++idp;
}
void dfs(int w){
if(ed[w]!=tk)ed[w]=tk,ID[w]=++IDP;
for(int i=e0[w],u;i;i=enx[i]){
u=es[i];
if(ed[u]!=tk)dfs(u);
if(ID[w]<ID[u]){
++v[ID[w]][ID[w]];
++v[ID[u]][ID[u]];
--v[ID[w]][ID[u]];
--v[ID[u]][ID[w]];
ee[eep][]=ID[w];
ee[eep][]=ID[u];
ee[eep++][]=ev[i];
}
}
}
const double _0=1e-;
double solve(int n){
static double a[][];
double s=;
for(int i=;i<=n;++i)
for(int j=;j<=n;++j)
a[i][j]=v[i][j];
for(int i=;i<=n;++i){
if(fabs(a[i][i])<_0){
int t=i;
for(int j=i+;j<=n;++j)if(fabs(a[j][i])>fabs(a[t][i]))t=j;
if(fabs(a[t][i])<_0)return ;
for(int j=i;j<=n;++j)swap(a[i][j],a[t][j]);
}
for(int j=i+;j<=n;++j)if(a[j][i]){
double x=a[j][i]/a[i][i];
for(int k=i;k<=n;++k)a[j][k]-=x*a[i][k];
}
}
for(int i=;i<=n;++i)s*=a[i][i];
return s;
}
void chk(int x){
if(ed[x]==tk)return;
IDP=;++tk;eep=;
dfs(x);
double v0=solve(IDP-);
for(int i=;i<eep;++i){
int x=ee[i][],y=ee[i][];
--v[x][x],--v[y][y],++v[x][y],++v[y][x];
double v1=v0-solve(IDP-);
ans+=v1/v0*ee[i][];
++v[x][x],++v[y][y],--v[x][y],--v[y][x];
}
for(int i=;i<=IDP;++i){
for(int j=;j<=IDP;++j)v[i][j]=;
}
}
int main(){
n=_();m=_();
for(int i=;i<=n;++i)f[i]=i;
for(int i=;i<m;++i)e[i].init();
std::sort(e,e+m);
for(int i=,j=;i<m;){
for(++tk,idp=;j<m&&e[i].c==e[j].c;++j);
for(int k=i;k<j;++k){
int x=get(e[k].a),y=get(e[k].b);
if(x==y){
e[k].a=-;
continue;
}
os[op][]=e[k].a;os[op++][]=e[k].b;
gid(x);gid(y);
e[k].a=id[x];
e[k].b=id[y];
e[k].adde();
}
while(op)--op,f[get(os[op][])]=get(os[op][]);
for(;i<j;++i)if(~e[i].a){
chk(e[i].a);
chk(e[i].b);
}
for(int t=;t<=idp;++t)e0[t]=;
ep=;
}
printf("%.5Lf",ans);
return ;
}
bzoj4637: 期望的更多相关文章
- 【BZOJ4637】期望 Kruskal+矩阵树定理
[BZOJ4637]期望 Description 在米国有一所大学,名叫万国歌剧与信息大学(UniversalOperaandInformaticasUniversity).简称UOI大学.UOI大学 ...
- bzoj4637:期望
思路:最小生成树计数只不过加了一个期望,由于期望具有线性性质,就可以转化为每条边的期望之和,那么一条边的期望如何求呢,在最小生成树记数中,是把相同边权的一起处理,之后把属于连通块内的点缩点,也就是说, ...
- 【BZOJ-3143】游走 高斯消元 + 概率期望
3143: [Hnoi2013]游走 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2264 Solved: 987[Submit][Status] ...
- bzoj1415[NOI2005]聪聪和可可-期望的线性性
这道题之前我写过一个巨逗比的写法(传送门:http://www.cnblogs.com/liu-runda/p/6220381.html) 当时的原因是这道题可以抽象出和"绿豆蛙的归宿&qu ...
- hdu 4481 Time travel(高斯求期望)(转)
(转)http://blog.csdn.net/u013081425/article/details/39240021 http://acm.hdu.edu.cn/showproblem.php?pi ...
- 【BZOJ3036】绿豆蛙的归宿 概率与期望
最水的概率期望,推荐算法合集之<浅析竞赛中一类数学期望问题的解决方法> #include <iostream> #include <cstdio> using na ...
- UVA&&POJ离散概率与数学期望入门练习[4]
POJ3869 Headshot 题意:给出左轮手枪的子弹序列,打了一枪没子弹,要使下一枪也没子弹概率最大应该rotate还是shoot 条件概率,|00|/(|00|+|01|)和|0|/n谁大的问 ...
- 【BZOJ-1426】收集邮票 概率与期望DP
1426: 收集邮票 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 261 Solved: 209[Submit][Status][Discuss] ...
- 【BZOJ-1419】Red is good 概率期望DP
1419: Red is good Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 660 Solved: 257[Submit][Status][Di ...
随机推荐
- C. Polycarpus' Dice
在每个位置讨论一下最大值最小值的取值范围就行 #include<cstdio> #include<iostream> #define maxn 200003 using nam ...
- BIOS设置
主界面: 标准BLOS设置: 高级BLOS设置: 1.设置从光盘驱动 电脑启动时默认是从硬盘启动,但是在安装操作系统时或者是使用某些特殊软件时,可能需要从光盘.软盘或者U盘启动,这时就需要设置第一启 ...
- asp.net 中json字符串转换
List<ATTVal> Replys = JsonParser.FromJson<List<ATTVal>>(attrValueStr);
- String的常规使用集合
今天先附上代码君: package com.jacob.javase; import java.io.UnsupportedEncodingException; /* *探讨String: * * ...
- SQL注入测试平台 SQLol -6.CHALLENGES挑战
SQLol上面的挑战共有14关,接下来我们一关一关来突破. Challenge 0 目的是让查询返回所有的用户名,而不是只有一个. SELECT username FROM users WHERE u ...
- 安装arbotix simulator仿真环境--9
原创博客:转载请表明出处:http://www.cnblogs.com/zxouxuewei/ 周学伟 安装之前:首先确保已经正常制作了ros工作空间并且安装了rbx1功能包: cd ~/catki ...
- Hibernate--Enum类型的set集合映射到数据库(xml配置文件实现方式)
使用enum 存储Permission的值 package demo; public enum Permission { CREATE,DELETE,UPDATE;} Role与Permissi ...
- nginx和apache下的url rewrite
将服务器上面的数据同步到本地之后,发现打开首页显示不正常,本地服务器是apache,经过打开url rewrite之后本地首页正常显示. 原因是phpwind本身支持了url rewrite的功能,但 ...
- 作业:用HTML制作邮箱登陆界面
<body leftmargin="200" rightmargin="200"> <font size="45" > ...
- (转) A Survival Guide to a PhD
A Survival Guide to a PhD Sep 7, 2016 This guide is patterned after my “Doing well in your courses”, ...