【BZOJ-3270】博物馆 高斯消元 + 概率期望
3270: 博物馆
Time Limit: 30 Sec Memory Limit: 128 MB
Submit: 292 Solved: 158
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
1 2
0.5
0.5
Sample Output
HINT
对于100%的数据有 n <= 20,n-1 <= m <= n(n-1)/2
Source
Solution
这个题想了一会,转化后会好想很多。
首先直接列方程组显然什么都解不出来,所以转化一下,令二元组$(x,y)$表示第一个人在$x$位置,第二个人在$y$位置的概率;
上述状态当做dp来做转移很显然$$dp[x][y]=p[x][y]*dp[x][y]+\frac {1-p[fx]}{d[fx]}*p[y]*dp[fx][y]+\frac {1-p[fy]}{d[fy]}*p[x]*dp[x][fy]+\frac {1-p[fx]}{d[fx]}*\frac {1-p[fy]}{d[fy]}*dp[fx][fy] \quad \quad \quad (f[A][B]=1)$$
把所有这样的状态设成未知数$X_{i}$,然后列方程,高斯消元得解。
列方程的复杂度$O(N^{4})$,解方程的复杂度是$O((N^{2})^{3})$的,总的复杂度是$O(N^{6})$
Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
inline int read()
{
int x=,f=; char ch=getchar();
while (ch<'' || ch>'') {if (ch=='-') f=-; ch=getchar();}
while (ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}
return x*f;
}
#define MAXN 500
#define eps 1e-5
int N,M,A,B,id[MAXN][MAXN],ID,mp[MAXN][MAXN],d[MAXN],flag;
double a[MAXN][MAXN],X[MAXN],p[MAXN];
inline void Debug()
{
puts("=========");
for (int i=; i<=ID; i++,puts(""))
for (int j=; j<=ID+; j++) printf("%.2lf ",a[i][j]);
puts("");
}
inline void Gauss()
{
flag=;
for (int i=; i<=ID; i++)
{
int mx=i;
for (int j=i+; j<=ID; j++)
if (abs(a[j][i])>abs(a[mx][i])) mx=j;
swap(a[i],a[mx]);
if (abs(a[i][i])<eps) {flag=-; continue;}
for (int j=i+; j<=ID+; j++) if (abs(a[i][j])>) a[i][j]/=a[i][i];
a[i][i]=; for (int j=i+; j<=ID; j++)
{
for (int k=i+; k<=ID+; k++)
a[j][k]-=a[j][i]*a[i][k];
a[j][i]=;
}
//Debug();
}
for (int i=,f=; i<=ID; i++,f=)
{
for (int j=; j<=ID && f; j++)
if (abs(a[i][j])>eps) f=;
if (abs(a[i][M+])>eps && f) flag=;
}
if (flag==) return;
for (int i=ID; i>=; i--)
{
X[i]=a[i][ID+];
for (int j=i+; j<=ID; j++) X[i]-=X[j]*a[i][j];
}
}
int main()
{
N=read(),M=read(),A=read(),B=read();
for (int i=,x,y; i<=M; i++) x=read(),y=read(),mp[x][y]=mp[y][x]=,d[x]++,d[y]++;
for (int i=; i<=N; i++) mp[i][i]=;
for (int i=; i<=N; i++)
for (int j=; j<=N; j++) id[i][j]=++ID;
for (int i=; i<=N; i++) scanf("%lf",&p[i]);
for (int i=; i<=N; i++)
for (int j=; j<=N; j++)
{
a[id[i][j]][id[i][j]]=-1.0;
for (int fi=; fi<=N; fi++)
for (int fj=; fj<=N; fj++)
if (fi!=fj)
{
int x=id[i][j],y=id[fi][fj];
if (mp[fi][i] && mp[fj][j])
{
if (fi==i && fj==j) a[x][y]+=p[fi]*p[fj];
if (fi!=i && fj==j) a[x][y]+=((1.0-p[fi])/d[fi]) * (p[fj]);
if (fi==i && fj!=j) a[x][y]+=((1.0-p[fj])/d[fj]) * (p[fi]);
if (fi!=i && fj!=j) a[x][y]+=((1.0-p[fi])/d[fi]) * ((1.0-p[fj])/d[fj]);
}
}
}
a[id[A][B]][ID+]=-1.0;
// Debug();
Gauss();
// Debug();
for (int i=; i<=N; i++) printf("%.6lf ",X[id[i][i]]);
return ;
}
【BZOJ-3270】博物馆 高斯消元 + 概率期望的更多相关文章
- 【BZOJ3143】【HNOI2013】游走 && 【BZOJ3270】博物馆 【高斯消元+概率期望】
刚学完 高斯消元,我们来做几道题吧! T1:[BZOJ3143][HNOI2013]游走 Description 一个无向连通图,顶点从1编号到N,边从1编号到M. 小Z在该图上进行随机游走,初始时小 ...
- BZOJ3270 博物馆(高斯消元+概率期望)
将两个人各自所在点视为状态,新建一个图.到达某个终点的概率等于其期望次数.那么高斯消元即可. #include<iostream> #include<cstdio> #incl ...
- 【BZOJ-3143】游走 高斯消元 + 概率期望
3143: [Hnoi2013]游走 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2264 Solved: 987[Submit][Status] ...
- [JLOI2012]时间流逝 树上高斯消元 概率期望
题面 题意:(感觉题面写的题意是错的?)有\(n\)种能量不同的圈,设当前拥有的圈的集合为\(S\),则: 1,每天有\(p\)概率失去一个能量最小的圈.特别的,如果\(S = \varnothing ...
- 洛谷P4457/loj#2513 [BJOI2018]治疗之雨(高斯消元+概率期望)
题面 传送门(loj) 传送门(洛谷) 题解 模拟赛的时候只想出了高斯消元然后死活不知道怎么继续--结果正解居然就是高斯消元卡常? 首先有个比较难受的地方是它一个回合可能不止扣一滴血--我们得算出\( ...
- 高斯消元与期望DP
高斯消元可以解决一系列DP序混乱的无向图上(期望)DP DP序 DP序是一道DP的所有状态的一个排列,使状态x所需的所有前置状态都位于状态x前: (通俗的说,在一个状态转移方程中‘=’左侧的状态应该在 ...
- HDU4870_Rating_双号从零单排_高斯消元求期望
原题链接:http://acm.hdu.edu.cn/showproblem.php?pid=4870 原题: Rating Time Limit: 10000/5000 MS (Java/Other ...
- hdu 4870 rating(高斯消元求期望)
Rating Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Sub ...
- BZOJ 2337: [HNOI2011]XOR和路径 [高斯消元 概率DP]
2337: [HNOI2011]XOR和路径 题意:一个边权无向连通图,每次等概率走向相连的点,求1到n的边权期望异或和 这道题和之前做过的高斯消元解方程组DP的题目不一样的是要求期望异或和,期望之间 ...
随机推荐
- Oracle安装
1.根据自己的操作系统,到Oracle官网下载相应的安装包 下载地址:http://download.oracle.com/otn/nt/oracle11g/112010/win32_11gR2_cl ...
- SQL中varchar和nvarchar的区别
varchar(n)长度为 n 个字节的可变长度且非 Unicode 的字符数据.n 必须是一个介于 1 和 8,000 之间的数值.存储大小为输入数据的字节的实际长度,而不是 n 个字节. nvar ...
- UITabBarController 升级定制
UITabBarController 定制 特点 用法 1.准备工作: 加入你的相关图片,放入了Assets.xcassets; 导入Categroy文件夹(这个里面的文件,在这里不详细说明了,有疑问 ...
- [Android]使用Dagger 2依赖注入 - API(翻译)
以下内容为原创,欢迎转载,转载请注明 来自天天博客:http://www.cnblogs.com/tiantianbyconan/p/5092525.html 使用Dagger 2依赖注入 - API ...
- Mybatis环境
第一步:下载jar包并导入 1.mysql驱动包 2.mybatis环境包 第二步:创建MYSQL数据库 由于这是用于测试,只创建了test-usreinfo数据表 第三步:在src文件夹中创建myb ...
- SlidingMenu 侧滑菜单的用法
很多APP都有侧滑菜单的功能,部分APP左右都是侧滑菜单~SlidingMenu 这个开源项目可以很好帮助我们实现侧滑功能,将侧滑菜单引入项目的方式很多中,先通过例子介绍各种引入方式,然后给大家展示个 ...
- 架构从最简单的数据访问框架(ORM)到资源调度和治理中心(SOA)说起
随着互联网的发展,网站应用的规模不断扩大,常规的垂直应用架构已无法应对,分布式服务架构以及流动计算架构势在必行,亟需一个治理系统确保架构有条不紊的演进. 单一应用架构当网站流量很小时,只需一个应用,将 ...
- Oracle forall bulk collect批量数据更新
对于数据量较大的插入操作可采用此种方法操作,注意: limit减少内存占用,如果数据量较大一次性全部加载到内存中,对PGA来说压力太大,可采用limit的方法一次加载一定数量的数据,建议值通常为100 ...
- 关于Oracle表连接
表连接注意left join on与where的区别: select * from dept; select * from emp; select * from emp a right outer j ...
- Django模板的继承
一.extend 1.extend继承模板 2.一个文件中只能继承一个模板 3.extend继承模板中的所有内容,模板的内容包括:html的head和body ,eg: