Description

在米国有一所大学,名叫万国歌剧与信息大学(UniversalOperaandInformaticasUniversity)。简称UOI大学。UO
I大学的建筑与道路分布很有趣,每对建筑之间有且仅有一条直接或者间接的路径相连,更加明确的说,就是成树
形分布。其实在设计时,对于大学的N个建筑,总共有M条道路可以修建,每条道路都有一个距离值Disti和一个美
学值Valuei。一个设计方案的距离值和美学值定义为该设计方案内包含的道路的距离值与美学值之和。投资方的要
求只有设计方案的距离值最小。大神出于对树的喜爱所以决定设计方案必须是一棵树。因为要参加UOI,所以当时
大神就急急忙忙地随机选择了一个合法的方案。但其实存在很多合法的方案,假设每种设计方案取的概率是均等的
,那么设计方案的美学值期望是多少?

Input

第一行两个整数,N和M,意义如上所述。
第二行到第M+1行,每行4个整数,Xi,Yi,Disti,Valuei,分别表示这条道路连接的
两个建筑的编号,距离值以及美学值。
输入保证至少有一种合法方案。
100%的数据保证N<=10000M<=200000
100%的数据保证距离值相同的道路数小于30,同时不保证没有重边。

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: 期望的更多相关文章

  1. 【BZOJ4637】期望 Kruskal+矩阵树定理

    [BZOJ4637]期望 Description 在米国有一所大学,名叫万国歌剧与信息大学(UniversalOperaandInformaticasUniversity).简称UOI大学.UOI大学 ...

  2. bzoj4637:期望

    思路:最小生成树计数只不过加了一个期望,由于期望具有线性性质,就可以转化为每条边的期望之和,那么一条边的期望如何求呢,在最小生成树记数中,是把相同边权的一起处理,之后把属于连通块内的点缩点,也就是说, ...

  3. 【BZOJ-3143】游走 高斯消元 + 概率期望

    3143: [Hnoi2013]游走 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2264  Solved: 987[Submit][Status] ...

  4. bzoj1415[NOI2005]聪聪和可可-期望的线性性

    这道题之前我写过一个巨逗比的写法(传送门:http://www.cnblogs.com/liu-runda/p/6220381.html) 当时的原因是这道题可以抽象出和"绿豆蛙的归宿&qu ...

  5. hdu 4481 Time travel(高斯求期望)(转)

    (转)http://blog.csdn.net/u013081425/article/details/39240021 http://acm.hdu.edu.cn/showproblem.php?pi ...

  6. 【BZOJ3036】绿豆蛙的归宿 概率与期望

    最水的概率期望,推荐算法合集之<浅析竞赛中一类数学期望问题的解决方法> #include <iostream> #include <cstdio> using na ...

  7. UVA&&POJ离散概率与数学期望入门练习[4]

    POJ3869 Headshot 题意:给出左轮手枪的子弹序列,打了一枪没子弹,要使下一枪也没子弹概率最大应该rotate还是shoot 条件概率,|00|/(|00|+|01|)和|0|/n谁大的问 ...

  8. 【BZOJ-1426】收集邮票 概率与期望DP

    1426: 收集邮票 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 261  Solved: 209[Submit][Status][Discuss] ...

  9. 【BZOJ-1419】Red is good 概率期望DP

    1419: Red is good Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 660  Solved: 257[Submit][Status][Di ...

随机推荐

  1. C. Polycarpus' Dice

    在每个位置讨论一下最大值最小值的取值范围就行 #include<cstdio> #include<iostream> #define maxn 200003 using nam ...

  2. BIOS设置

    主界面:  标准BLOS设置: 高级BLOS设置: 1.设置从光盘驱动 电脑启动时默认是从硬盘启动,但是在安装操作系统时或者是使用某些特殊软件时,可能需要从光盘.软盘或者U盘启动,这时就需要设置第一启 ...

  3. asp.net 中json字符串转换

    List<ATTVal> Replys = JsonParser.FromJson<List<ATTVal>>(attrValueStr);

  4. String的常规使用集合

    今天先附上代码君: package com.jacob.javase; import java.io.UnsupportedEncodingException; /* *探讨String:  *  * ...

  5. SQL注入测试平台 SQLol -6.CHALLENGES挑战

    SQLol上面的挑战共有14关,接下来我们一关一关来突破. Challenge 0 目的是让查询返回所有的用户名,而不是只有一个. SELECT username FROM users WHERE u ...

  6. 安装arbotix simulator仿真环境--9

    原创博客:转载请表明出处:http://www.cnblogs.com/zxouxuewei/  周学伟 安装之前:首先确保已经正常制作了ros工作空间并且安装了rbx1功能包: cd ~/catki ...

  7. Hibernate--Enum类型的set集合映射到数据库(xml配置文件实现方式)

    使用enum 存储Permission的值 package demo; public enum Permission {    CREATE,DELETE,UPDATE;} Role与Permissi ...

  8. nginx和apache下的url rewrite

    将服务器上面的数据同步到本地之后,发现打开首页显示不正常,本地服务器是apache,经过打开url rewrite之后本地首页正常显示. 原因是phpwind本身支持了url rewrite的功能,但 ...

  9. 作业:用HTML制作邮箱登陆界面

    <body leftmargin="200" rightmargin="200"> <font size="45" > ...

  10. (转) 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”, ...