BZOJ3253 : 改编
设$f[x][y]$表示从x和y出发相遇的期望长度,则$f[x][x]=0$,且$f[x][y]$对称,共$C(n,2)$个未知量。
列出方程组$G$,得到$G\times F=B$。
高斯消元求出$G$的逆矩阵$G'$,则$F=G'\times B$,对于每个询问代入计算即可。
时间复杂度$O(n^6+tn^4)$。
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=235;
int n,m,T,i,j,k,x,y,o,id[25][25],cnt,d[25],e[25][25];
double p[25],f[N][N],g[N][N],B[N][N],len[N],t,ans;
int main(){
scanf("%d%d%d",&n,&m,&T);
for(i=1;i<=n;i++)for(j=1;j<=n;j++)e[i][j]=-1;
for(i=1;i<=n;i++)e[i][i]=0;
for(i=1;i<=m;i++){
scanf("%d%d",&x,&y);
d[x]++,d[y]++;
e[x][y]=e[y][x]=i;
}
for(i=1;i<=n;i++)scanf("%lf",&p[i]);
for(i=1;i<=n;i++)for(j=i+1;j<=n;j++)id[i][j]=++cnt;
for(i=1;i<=n;i++)for(j=i+1;j<=n;j++){
o=id[i][j];
f[o][o]=1;
for(x=1;x<=n;x++)if(~e[i][x])for(y=1;y<=n;y++)if(~e[j][y]){
t=1;
if(x==i)t*=p[i];else t*=(1.0-p[i])/d[i];
if(y==j)t*=p[j];else t*=(1.0-p[j])/d[j];
B[o][e[i][x]]+=t;
B[o][e[j][y]]+=t;
if(x!=y)f[o][x<y?id[x][y]:id[y][x]]-=t;
}
}
for(i=1;i<=cnt;i++)g[i][i]=1;
for(i=1;i<=cnt;i++){
for(k=i,j=i+1;j<=cnt;j++)if(fabs(f[j][i])>fabs(f[k][i]))k=j;
if(k!=i)for(j=1;j<=cnt;j++)swap(f[i][j],f[k][j]),swap(g[i][j],g[k][j]);
for(j=i+1;j<=cnt;j++)for(t=f[j][i]/f[i][i],k=1;k<=cnt;k++)f[j][k]-=f[i][k]*t,g[j][k]-=g[i][k]*t;
}
for(i=cnt;i;i--){
for(j=cnt;j>i;j--)for(t=f[i][j],k=1;k<=cnt;k++)f[i][k]-=f[j][k]*t,g[i][k]-=g[j][k]*t;
for(t=f[i][i],j=1;j<=cnt;j++)f[i][j]/=t,g[i][j]/=t;
}
while(T--){
for(i=1;i<=m;i++)scanf("%lf",&len[i]);
scanf("%d%d",&x,&y);
if(x==y){puts("0.00");continue;}
o=x<y?id[x][y]:id[y][x];
ans=0;
for(i=1;i<=cnt;i++)for(j=1;j<=m;j++)ans+=g[o][i]*B[i][j]*len[j];
printf("%.2f\n",ans);
}
return 0;
}
BZOJ3253 : 改编的更多相关文章
- 【转】监听按钮除OnClick外其他事件的方法,附简易改编的UIButton类
http://lib.csdn.net/article/unity3d/38463 作者:IceFantasyLcj 大家好,我是雨中祈雨.一直以来,CSDN都是我最好的编程助手.这是我在CSDN的第 ...
- 【转】 NGUI 监听按钮除OnClick外其他事件的方法,附简易改编的UIButton类
http://blog.csdn.net/icefantasylcj/article/details/49450555 大家好,我是雨中祈雨.一直以来,CSDN都是我最好的编程助手.这是我在CSDN的 ...
- NOIP2014无线网络发射器选址改编1
问题描述 随着智能手机的日益普及,人们对无线网的需求日益增大.某城市决定对城市内的公共场所覆盖无线网. 假设该城市的布局为由严格平行的129条东西向街道和129条南北向街道所形成的网格状,并且相邻的平 ...
- APK改之理 手游修改改编安卓程序工具安装使用教程
APK改之理 手游修改改编安卓程序工具安装使用教程 --APK破解付费程序 apk改之理是pc平台上一款非常好用的apk反编译工具,他将反编译以及签名等功能集中在一起,并且拥有非常人性化的操作界面,如 ...
- C语言 · 8皇后问题改编
8皇后问题(改编) 问题描述 规则同8皇后问题,但是棋盘上每格都有一个数字,要求八皇后所在格子数字之和最大. 输入格式 一个8*8的棋盘. 输出格式 所能得到的最大数字和 样例输入 1 2 3 4 5 ...
- 2017Facebook面试题改编“一面砖墙 ”
题目:一面砖墙 这道题改编自网上Facebook去年的一道面试题,是hihoCoder的1494题(https://hihocoder.com/problemset/problem/1494) 这道题 ...
- Atcoder Regular Contest 092 A 的改编
原题地址 题目大意 给定平面上的 $n$ 个点 $p_1, \dots, p_n$ .第 $i$ 点的坐标为 $(x_i, y_i)$ .$x_i$ 各不相同,$y_i$ 也各不相同.若两点 $p_i ...
- Cocos2d-x 3.1 一步一步地做改编
本文并不想谈论的屏幕改编或真理的概念.假设不知道cocos2d-x的,请先看这篇文章:http://www.cocoachina.com/gamedev/cocos/2014/0516/8451.ht ...
- C++ ABI之名字改编(以Qt为例)
在C++中,由于重载等技术的存在,编译器要将函数.结构体.类等等的信息传递给链接器,就不能像C语言那样简单地通过函数名来完成,它需要提供额外的参数信息,而还要和C语言共用链接器,这就需要用到名字改编( ...
随机推荐
- 论文阅读笔记四:CTPN: Detecting Text in Natural Image with Connectionist Text Proposal Network(ECCV2016)
前面曾提到过CTPN,这里就学习一下,首先还是老套路,从论文学起吧.这里给出英文原文论文网址供大家阅读:https://arxiv.org/abs/1609.03605. CTPN,以前一直认为缩写一 ...
- [转]PO BO VO DTO POJO DAO概念及其作用(附转换图)
来源:http://www.blogjava.net/vip01/archive/2013/05/25/92430.html J2EE开发中大量的专业缩略语很是让人迷惑,尤其是跟一些高手讨论问题的时候 ...
- net core体系-web应用程序-4net core2.0大白话带你入门-1目录
asp.net core2.0大白话带你入门 本系列包括: 1.新建asp.net core项目2.web项目目录解读3.配置访问地址4.环境变量详解5.配置文件6.日志7.DI容器8.服务的生命周期 ...
- selenium课程笔记
selenium课程笔记第一天(2017-7-1) 一.配置火狐浏览器 运行:firefox.exe -p -no -remote selenium课程笔记第二天 用Eclipse+java+sele ...
- Java03动手动脑
1.当JAVA里定义的函数中去掉static后,怎么办? static代表静态,由于main函数是静态的,如果自己定义的函数方法加了static则在类加载时就一起加载了.但如果不写static,就必须 ...
- Kafka-API
ctrl+Hnew 它的实现类ctrl+r替换格式化ctrl+alt+l ctrl+fctrl+alt+v 替换 < " < < > > Ka ...
- 01++ Bookshelf 2
http://poj.org/problem?id=3628 就是比原题多了一个要求,输出>=m的最小值 kisang~独立做出来的都开心<( ̄︶ ̄)> #include<cs ...
- shiro学习总结(一)----初识shiro
本系列内容大多总结自官网和张开涛的<跟我学Shiro> 一.shiro简介 1.1.shiro有什么用? shiro是一个功能强大使用简单的java安全框架,主要提供了五大功能: 1.认证 ...
- Spring里的applicationContext.getBeansWithAnnotation
getBeansWithAnnotation Map<String,Object> getBeansWithAnnotation(Class<? extends Annotation ...
- ORA-01440 要减小精度或者标度.则要修改的列必须为空.
--1 创建备份表 drop table contract_kangjia_back; create table contract_kangjia_back as select * from cont ...