BZOJ 1778 [Usaco2010 Hol] Dotp 驱逐猪猡
Description
Input
Output
Sample Input
1 2
Sample Output
0.666666667
0.333333333
HINT
Source
Gold
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
高斯消元+拆点~
拆点,把每个点拆成两个点,i+n是在i点爆炸了,i是没有爆炸,其中i+n是不能用来转移其他点的。
p/q要强制转化;p/q不能记成一个double再调用,否则会被卡精度……
#include<cstdio>int n,m,x,y,fi[301],w[89701],ne[897021],cnt,du[301];
double p,q,a[601][602];int read()
{int totnum=0,f=1;char ch=getchar();while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}while(ch>='0' && ch<='9') {totnum=(totnum<<1)+(totnum<<3)+ch-'0';ch=getchar();}return totnum*f;
}void add(int u,int v)
{w[++cnt]=v;ne[cnt]=fi[u];fi[u]=cnt;w[++cnt]=u;ne[cnt]=fi[v];fi[v]=cnt;
}void guass(int n)
{for(int i=1;i<n;i++){for(int j=i+1;j<=n+1;j++) a[i][j]/=a[i][i];a[i][i]=1;for(int j=i+1;j<=n;j++){for(int k=i+1;k<=n+1;k++) a[j][k]-=a[j][i]*a[i][k];a[j][i]=0;}}for(int i=n;i>1;i--){for(int j=i+1;j<=n+1;j++) a[i][j]/=a[i][i];a[i][i]=1;for(int j=i-1;j;j--){for(int k=i+1;k<=n+1;k++) a[j][k]-=a[j][i]*a[i][k];a[j][i]=0;}}
}int main()
{n=read();m=read();p=read();q=read();for(int i=1;i<=m;i++) x=read(),y=read(),add(x,y),du[x]++,du[y]++;for(int i=1;i<=n;i++) a[i][i]=-1.0,a[i+n][i+n]=-1.0,a[n+i][i]=(double)p/q;for(int i=1;i<=n;i++)for(int j=fi[i];j;j=ne[j]) a[i][w[j]]=(1.0-(double)p/q)/(double)du[w[j]];a[1][n*2+1]=-1;guass(n*2);for(int i=1;i<=n;i++) printf("%.9lf\n",a[i+n][2*n+1]/a[i+n][i+n]);return 0;
}
BZOJ 1778 [Usaco2010 Hol] Dotp 驱逐猪猡相关推荐
- [BZOJ 1778][Usaco2010 Hol]Dotp 驱逐猪猡
1778: [Usaco2010 Hol]Dotp 驱逐猪猡 Time Limit: 10 Sec Memory Limit: 64 MB Submit: 690 Solved: 269 [Sub ...
- BZOJ 1778: [Usaco2010 Hol]Dotp 驱逐猪猡 [高斯消元 概率DP]
1778: [Usaco2010 Hol]Dotp 驱逐猪猡 题意:一个炸弹从1出发p/q的概率爆炸,否则等概率走向相邻的点.求在每个点爆炸的概率 高斯消元求不爆炸到达每个点的概率,然后在一个点爆炸就 ...
- bzoj 1778: [Usaco2010 Hol]Dotp 驱逐猪猡
Description 奶牛们建立了一个随机化的臭气炸弹来驱逐猪猡.猪猡的文明包含1到N (2 <= N <= 300)一共N个猪城.这些城市由M (1 <= M <= 44, ...
- bzoj 1778 [Usaco2010 Hol]Dotp 驱逐猪猡
http://www.elijahqi.win/archives/3299 Description 奶牛们建立了一个随机化的臭气炸弹来驱逐猪猡.猪猡的文明包含1到N (2 <= N <= ...
- BZOJ 1778 [Usaco2010 Hol]Dotp 驱逐猪猡 ——期望DP
思路和BZOJ 博物馆很像. 同样是高斯消元 #include <map> #include <ctime> #include <cmath> #include & ...
- BZOJ 1778 Usaco2010 Hol Dotp 驱逐猪猡 期望DP+高斯消元
题目大意:给定一个无向图,炸弹从1号节点出发,每个时刻有P/Q的概率爆炸,如果某个时刻没有爆炸,就会等概率沿着随机一条出边走到下一个城市,求最终每个城市的爆炸概率 #include <cmath ...
- BZOJ 1778: [Usaco2010 Hol]Dotp 驱逐猪猡 概率与期望+高斯消元
这个还挺友好的,自己相对轻松能想出来~ 令 $f[i]$ 表示起点到点 $i$ 的期望次数,则 $ans[i]=f[i]\times \frac{p}{q}$ #include <cmath&g ...
- BZOJ 1778: [Usaco2010 Hol]Dotp 驱逐猪猡 (高斯消元)
题面 题目传送门 分析 令爆炸概率为PPP.设 f(i)=∑k=0∞pk(i)\large f(i)=\sum_{k=0}^{\infty}p_k(i)f(i)=∑k=0∞pk(i),pk(i)p ...
- bzoj 1778: [Usaco2010 Hol]Dotp 驱逐猪猡【dp+高斯消元】
算是比较经典的高斯消元应用了 设f[i]为i点答案,那么dp转移为f[u]=Σf[v]*(1-p/q)/d[v],意思是在u点爆炸可以从与u相连的v点转移过来 然后因为所有f都是未知数,高斯消元即可( ...
最新文章
- 高性能千万级定时任务管理服务forsun使用详解
- 5G 标准 — R16
- MSSQLServer基础01(数据类型)
- idea 编译显示source1.3不支持泛型(请使用source5或更高版本)
- UVALive 3942 Remember the Word(字典树+DP)
- dart系列之:dart类的扩展
- 【Nutch2.3基础教程】集成Nutch/Hadoop/Hbase/Solr构建搜索引擎:安装及运行【集群环境】
- C/C++高级算法之绘制曼德布洛特集
- 1711: [Usaco2007 Open]Dingin吃饭
- 程序员的自我修养——读《软技能-代码之外的生存指南》笔记
- 数学与计算机科学学院宣传词,【榜样人物】筑梦前行 不负韶华——数学与计算机科学学院姚彬彬同学的考研心路...
- 第五章 调优案列分析与实战(待续)
- 集合涉及到的排序方式
- 树莓派HDMI转VGA线有无源
- Kobe -接小球游戏
- P1740 Ink on paper
- java 转成数字,java汉字转成数字
- Executors 源码解析(JDK8)
- 微软识花 android版 (非官方)
- HCNP学习笔记-路由交换网络(2)-IS-IS协议原理及配置