博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
CodeForces 605 E. Intergalaxy Trips
阅读量:5327 次
发布时间:2019-06-14

本文共 3764 字,大约阅读时间需要 12 分钟。

E. Intergalaxy Trips
time limit per test:2 seconds
memory limit per test:256 megabytes
input:standard input
output:standard output

The scientists have recently discovered wormholes — objects in space that allow to travel very long distances between galaxies and star systems.

The scientists know that there are n galaxies within reach. You are in the galaxy number 1 and you need to get to the galaxy number n. To get from galaxy i to galaxy j, you need to fly onto a wormhole (i, j) and in exactly one galaxy day you will find yourself in galaxy j.

Unfortunately, the required wormhole is not always available. Every galaxy day they disappear and appear at random. However, the state of wormholes does not change within one galaxy day. A wormhole from galaxy i to galaxy j exists during each galaxy day taken separately with probability pij. You can always find out what wormholes exist at the given moment. At each moment you can either travel to another galaxy through one of wormholes that exist at this moment or you can simply wait for one galaxy day to see which wormholes will lead from your current position at the next day.

Your task is to find the expected value of time needed to travel from galaxy 1 to galaxy n, if you act in the optimal way. It is guaranteed that this expected value exists.

Input

The first line of the input contains a single integer n (1 ≤ n ≤ 1000) — the number of galaxies within reach.

Then follows a matrix of n rows and n columns. Each element pij represents the probability that there is a wormhole from galaxy i to galaxy j. All the probabilities are given in percents and are integers. It is guaranteed that all the elements on the main diagonal are equal to 100.

Output

Print a single real value — the expected value of the time needed to travel from galaxy 1 to galaxy n if one acts in an optimal way. Your answer will be considered correct if its absolute or relative error does not exceed 10 - 6.

Namely: let's assume that your answer is a, and the answer of the jury is b. The checker program will consider your answer correct, if .

Examples
input
3 100 50 50 0 100 80 0 0 100
output
1.750000000000000
input
2 100 30 40 100
output
3.333333333333333
Note

In the second sample the wormhole from galaxy 1 to galaxy 2 appears every day with probability equal to 0.3. The expected value of days one needs to wait before this event occurs is .

题意:

给出一张$n$个点$n^{2}$条边的有向图,每条边每天的出现概率为$p[i][j]$,求从$1$到$n$的期望天数...

分析:

首先我们考虑两个点的情况,也就是第二个样例,从$1$到$2$的边的出现概率为$0.3$,所以我们此时求的期望天数就是期望第几天会出现这条边:$ans=\sum _{i=0}^{+∞}0.7^{i}$,收敛一下就是$\frac {1}{0.3}$...

然后我们再考虑多个点的情况,如果我们要从$i$走到$j$,必须满足的是走到$j$之后的结果要比$i$优,否则就不走...所以我们是每次选取一个最优的点去更新其他的点,这就是一个$dijkstra$的过程,更新的方式就是$f[i]=\frac {(p[i][j_{1}]*f[j_{1}]+(1-p[i][j_{1}])*p[i][j_{2}]*f[j_{2}]+……+1)}{1-(1-p[i][j_{1}])(1-p[i][j_{2}])……}$...

代码:

#include
#include
#include
#include
#include
//by NeighThornusing namespace std;const int maxn=1000+5;int n,vis[maxn];double f[maxn],a[maxn],b[maxn],p[maxn][maxn];struct M{ int x; double y; friend bool operator < (M a,M b){ return a.y>b.y; } M(int a=0,double b=0.0){ x=a;y=b; } };priority_queue
q;signed main(void){ scanf("%d",&n); for(int i=1,x;i<=n;i++) for(int j=1;j<=n;j++) scanf("%d",&x),p[i][j]=x/100.0; for(int i=1;i<=n;i++) a[i]=b[i]=1.0,f[i]=1e30; f[n]=0;q.push(M(n,0)); while(!q.empty()){ int top=q.top().x;q.pop(); if(vis[top]) continue; vis[top]=1; for(int i=1;i<=n;i++) if(p[i][top]>0&&!vis[i]){ a[i]+=b[i]*p[i][top]*f[top]; b[i]*=1.0-p[i][top]; f[i]=a[i]/(1.0-b[i]); q.push(M(i,f[i])); } } printf("%.15f\n",f[1]); return 0;}

  


By NeighThorn

转载于:https://www.cnblogs.com/neighthorn/p/6439160.html

你可能感兴趣的文章
Linux设置环境变量的方法
查看>>
构建自己的项目管理方案
查看>>
利用pca分析fmri的生理噪声
查看>>
div水平居中且垂直居中
查看>>
epoll使用具体解释(精髓)
查看>>
AndroidArchitecture
查看>>
安装Endnote X6,但Word插件显示的总是Endnote Web"解决办法
查看>>
python全栈 计算机硬件管理 —— 硬件
查看>>
大数据学习
查看>>
简单工厂模式
查看>>
Delphi7编译的程序自动中Win32.Induc.a病毒的解决办法
查看>>
Objective-C 【关于导入类(@class 和 #import的区别)】
查看>>
倍福TwinCAT(贝福Beckhoff)常见问题(FAQ)-点击运行按钮进入到运行状态报错Error starting TwinCAT System怎么办 AdsWarning1823怎么办...
查看>>
【转】javascript 中的很多有用的东西
查看>>
Centos7.2正常启动关闭CDH5.16.1
查看>>
Android 监听返回键、HOME键
查看>>
Android ContentProvider的实现
查看>>
sqlserver 各种判断是否存在(表名、函数、存储过程等)
查看>>
给C#学习者的建议 - CLR Via C# 读后感
查看>>
Recover Binary Search Tree
查看>>