╰( ̄▽ ̄)╭

Morenan 被困在了一个迷宫里。

迷宫可以视为 N 个点 M 条边的有向图,其中 Morena n处于起点 S , 迷宫的终点设为 T 。

可惜的是 , Morenan 非常的脑小 , 他只会从一个点出发随机沿着一条从该点出发的有向边 , 到达另一个点 。

这样 , Morenan 走的步数可能很长 , 也可能是无限,更可能到不了终点。

若到不了终点,则步数视为无穷大。

但你必须想方设法求出 Morenan 所走步数的期望值。

(⊙ ▽ ⊙)

一开始看着道题,就觉得是tarjan缩点后,转化成DAG上的问题。

当原图是DAG时

设fi表示第i个点走到终点的距离。

容易有fi=∑j∈next(i)1|next(i)|∗(fj+1),

其中next(i)是i的后继集合。

很容易使用拓扑排序来完成动态规划。

当原图是一般的有向图时

利用tarjan算法可以把图中的强连通分量找出来。

对于任意一个强连通分量,我们利用高斯消元来求解出强连通分量中的每个点的f值。

套上拓扑排序,就能够解决。


时间复杂度为O(n∗L3),其中L为最大强连通分量的大小。

实际时间复杂度则远远不到。

( ̄~ ̄)

高斯(gauss)消元:

1.目标

对n条n元一次方程求解。

2.步骤

首先:

把某条方程的第一个元的系数化为一,利用它可以消去其他方程的第一元。

以此类推,可以消去第二个元、第三个元……

最后,会只剩下一条只关于最后一元的方程,实际上就是这一元的根。

然后:

得出最后一元的根后,就可以代入先前的方程。

于是可以算出所有元的根。

(⊙v⊙)

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#include<math.h>
#define ll long long
using namespace std;
const char* fin="jzoj2758.in";
const char* fout="jzoj2758.out";
const int Inf=0x7fffffff;
const int maxn=20007,maxm=2000007,maxk=107;
int n,m,St,En,i,j,k,N,ti;
int fi[maxn],ne[maxm],la[maxm],tot;
int Fi[maxn],Ne[maxm],La[maxm],Tot;
int dfn[maxn],low[maxn],num,st[maxn],be[maxn];
bool bz[maxn],End;
int inf[maxn];
int az[maxn];
int ha[maxn][maxk],ru[maxn],id[maxn],tow[maxn];
double f[maxn],gs[maxk][maxk];
int b[maxn],head,tail;
void add_line(int a,int b){
tot++;
ne[tot]=fi[a];
la[tot]=b;
fi[a]=tot;
}
void Add_line(int a,int b){
Tot++;
Ne[Tot]=Fi[a];
La[Tot]=b;
Fi[a]=Tot;
}
void tarjan(int v){
int i,j,k;
dfn[v]=low[v]=++num;
bz[st[j=++st[0]]=v]=true;
for (k=fi[v];k;k=ne[k])
if (!dfn[la[k]]){
tarjan(la[k]);
low[v]=min(low[v],low[la[k]]);
}else if (bz[la[k]]) low[v]=min(low[la[k]],low[v]);
if (low[v]==dfn[v]){
N++;
be[N]=0;
while (st[0]>=j){
ha[N][++ha[N][0]]=st[st[0]];
be[st[st[0]]]=N;
bz[st[st[0]--]]=false;
}
}
}
double ABS(double x){
return x>0?x:-x;
}
void add(int v){
b[++tail]=v;
}
void count(int v,int n){
int i,j,k,l;
double tmp,tmd;
ti++;
for (i=1;i<=n-1;i++){
for (j=1;j<n;j++){
if (az[j]<ti && (gs[j][i]>10e-13 || gs[j][i]<-10e-13)){
az[j]=ti;
tow[i]=j;
tmp=1.0/gs[j][i];
for (k=1;k<=n;k++) gs[j][k]*=tmp;
for (k=1;k<n;k++){
if (k==j) continue;
tmd=gs[k][i];
for (l=1;l<=n;l++)
gs[k][l]-=tmd*gs[j][l];
}
break;
}
}
}
for (i=n-1;i>=1;i--){
f[ha[v][tow[i]]]=-gs[tow[i]][n];
for (j=1;j<=n-1;j++){
gs[j][n]-=gs[j][i]*gs[tow[i]][n];
}
}
}
void work(int v){
int i,j,k;
for (i=1;i<=ha[v][0];i++) id[ha[v][i]]=i;
for (i=1;i<=ha[v][0];i++){
for (j=1;j<=ha[v][0]+1;j++) gs[i][j]=0;
j=0;
for (k=Fi[ha[v][i]];k;k=Ne[k]){
gs[i][ha[v][0]+1]+=1;
gs[i][i]-=1;
if (f[La[k]]>10e-7 || La[k]==En) gs[i][ha[v][0]+1]+=f[La[k]];
else gs[i][id[La[k]]]+=1;
}
}
count(v,ha[v][0]+1);
}
void update(int v){
int i,j,k;
for (i=1;i<=ha[v][0];i++)
for (k=fi[ha[v][i]];k;k=ne[k]){
inf[be[la[k]]]=max(inf[be[la[k]]],inf[v]);
if (!--ru[be[la[k]]]) add(be[la[k]]);
}
}
void topsort(){
int i,j,k;
head=tail=0;
memset(f,0,sizeof(f));
f[En]=0;
for (i=n+1;i<=N;i++)
if (!ru[i]){
add(i);
}
while (head++<tail){
if (be[En]==b[head]) inf[b[head]]=1;
if (inf[b[head]]==0) inf[b[head]]=2;
if (inf[b[head]]==1 && b[head]!=be[En]) work(b[head]);
update(b[head]);
}
}
int main(){
scanf("%d%d%d%d",&n,&m,&St,&En);
N=n;
for (i=1;i<=m;i++){
scanf("%d%d",&j,&k);
if (j!=En){
Add_line(j,k);
add_line(k,j);
}
}
for (i=1;i<=n;i++) if (!dfn[i]) tarjan(i);
for (i=1;i<=n;i++)
for (k=fi[i];k;k=ne[k])
if (be[la[k]]!=be[i]) ru[be[la[k]]]++;
topsort();
//printf("%3.lf",f[St]);
if (inf[be[St]]!=1) printf("INF");
else printf("%.3lf",f[St]);
return 0;
}

【JZOJ2758】【SDOI2012】走迷宫(labyrinth)的更多相关文章

  1. BZOJ 2707: [SDOI2012]走迷宫( tarjan + 高斯消元 )

    数据范围太大不能直接高斯消元, tarjan缩点然后按拓扑逆序对每个强连通分量高斯消元就可以了. E(u) = 1 + Σ E(v) / degree(u) 对拍时发现网上2个程序的INF判断和我不一 ...

  2. BZOJ 2707: [SDOI2012]走迷宫 [高斯消元 scc缩点]

    2707: [SDOI2012]走迷宫 题意:求s走到t期望步数,\(n \le 10^4\),保证\(|SCC| \le 100\) 求scc缩点,每个scc高斯消元,scc之间直接DP 注意每次清 ...

  3. 【BZOJ2707】[SDOI2012]走迷宫 Tarjan+拓扑排序+高斯消元+期望

    [BZOJ2707][SDOI2012]走迷宫 Description Morenan被困在了一个迷宫里.迷宫可以视为N个点M条边的有向图,其中Morenan处于起点S,迷宫的终点设为T.可惜的是,M ...

  4. SDOI2012 走迷宫

    走迷宫 Morenan被困在了一个迷宫里.迷宫可以视为N个点M条边的有向图,其中Morenan处于起点S,迷宫的终点设为T.可惜的是,Morenan非常的脑小,他只会从一个点出发随机沿着一条从该点出发 ...

  5. bzoj 2707 [SDOI2012]走迷宫(SCC+高斯消元)

    Description Morenan被困在了一个迷宫里.迷宫可以视为N个点M条边的有向图,其中Morenan处于起点S,迷宫的终点设为T.可惜的是,Morenan非常的脑小,他只会从一个点出发随机沿 ...

  6. BZOJ2707 [SDOI2012]走迷宫 【概率dp + tarjan + 高斯消元】

    题目 Morenan被困在了一个迷宫里.迷宫可以视为N个点M条边的有向图,其中Morenan处于起点S,迷宫的终点设为T.可惜的是,Morenan非常的脑小,他只会从一个点出发随机沿着一条从该点出发的 ...

  7. [SDOI2012]走迷宫 (强连通分量缩点,动态规划,高斯消元)

    题面 Morenan被困在了一个迷宫里.迷宫可以视为N个点M条边的有向图,其中Morenan处于起点S,迷宫的终点设为T.可惜的是,Morenan非常的脑小,他只会从一个点出发随机沿着一条从该点出发的 ...

  8. bzoj2702[SDOI2012]走迷宫

    题意:给你一个有向图,点数10000,边数1000000,SCC大小不超过100(按数据范围的写法只有第三部分数据满足这个条件,不过第二部分数据并没有出现大小大于100个点的SCC,我是用数组大小为1 ...

  9. BZOJ.2707.[SDOI2012]走迷宫(期望 Tarjan 高斯消元)

    题目链接 一个点到达终点的期望步数 \(E_i=\sum_{(i,j)\in G}\frac{E_j+1}{out[i]}\),\(out[i]\)为点\(i\)的出度. 那么对于一个DAG可以直接在 ...

随机推荐

  1. Python读写文件学习笔记

    一. 基础 1.创建文件夹 import os os.makedirs('I:\\pythonWorkPace') # 创建文件夹 2. 获取文件夹里面文件列表 import os # os.make ...

  2. MaxCompute问答整理之8月

    本文是基于对MaxCompute产品的学习进度,再结合开发者社区里面的一些问题,进而整理成文.希望对大家有所帮助. 问题一.通过数据源数据增量同步后,如何查看某一条数据具体被同步到MaxCompute ...

  3. 线条之美,玩转SVG线条动画

    线条之美,玩转SVG线条动画 作者:AlloyTeam www.alloyteam.com/2017/02/the-beauty-of-the-lines-break-lines-svg-animat ...

  4. Python玩转人工智能最火框架 TensorFlow应用实践

    Python玩转人工智能最火框架 TensorFlow应用实践 整个课程都看完了,这个课程的分享可以往下看,下面有链接,之前做java开发也做了一些年头,也分享下自己看这个视频的感受,单论单个知识点课 ...

  5. ajax无刷新上传文件

    网页上传文件最简单的方式就是通过表单上传了,但是在提交表单的时候会导致网页刷新,但有的时候我们不想网页刷新,有什么办法呢,我们可以使用ajax上传文件来做到这一点.只有ajax还不行,还需要JavaS ...

  6. 一个网页登陆界面写30多个测试Case——测试之道

    转自博文:http://www.cnblogs.com/I-am-Betty/p/3566411.html 具体需求: 有一个登陆页面, (假如上面有2个textbox, 一个提交按钮. 请针对这个页 ...

  7. TZ_10_spring-sucrity 服务器和页面的权限控制

    1.在服务器端我们可以通过Spring security提供的注解对方法来进行权限控制. Spring Security在方法的权限控制上支持三种类型的注解,JSR-250注解.@Secured注解和 ...

  8. SSM3-SVN的安装和搭建环境

    1.安装svn 2.创建仓库 3.设置用户 . 4.eclipse和svn的集成 eclipse里安装SVN插件,一般来说,有两种方式: 直接下载SVN插件,将其解压到eclipse的对应目录里 使用 ...

  9. iview 小问题笔记

    总结一下用 iview 组件期间遇到的小问题.一边做项目一边做笔记. 官方文档:iview API 按需引入: 1,iview 分页组件 <Page :total="aboutPage ...

  10. Http中GET和POST请求的区别

    https://mp.weixin.qq.com/s/_oLnkDQn3TE7XdFWCT5Img GET请求 GET /books/?sex=man&name=Professional HT ...