[ICPC2019 南京网络赛]-D.Robots-拓扑期望递推

说在前面

感觉自己好菜呀
推了一个多小时都没出来，最后还是看的题解emmm

---

题目

计蒜客T41301传送门

题目大意

给出一张N个点M条边的有向无环图，1号点为全图唯一的入度为0的点，N号点为全图唯一的出度为0的点
现在你从1号点出发，每单位时间，你有相同的概率 走向相邻节点或原地不动。第 $i$ 单位时间内你的消耗为 $i$
询问你走到 N 点的期望消耗

数据范围

有T组数据，$1\le T\le 10$
对于每组数据，$2\le N\le 10^5,1\le M\le 2\times 10^5$
保证$\sum n\le 4\ast 10^5,\sum m\le 5\ast 10^5$

---

解法

概率顺推，期望逆推
显然，这道题不能一次性得出答案，因为走的每一步消耗不一样，需要对这样的消耗预处理

定义$t[u]$ 表示从 $u$ 走到 $N$ 的期望时间，转移很简单：
$t[u]=\sum \frac{t[v]}{outdeg[u]+1}+\frac{t[u]}{outdeg[u]+1}+1$，移项化简一下，左边就只有$t[u]$了

定义$f[u]$ 表示从 $u$ 走到 $N$ 的期望花费，用 $t[u]$ 来辅助转移
$f[u]=\sum \frac{f[v]}{outdeg[u]+1}+\frac{f[u]}{outdeg[u]+1}+t[u]$，化简一下，左边只剩$f[u]$，就可以转移了

注意这个 $f$ 的转移，其实我们把每单位时间的消耗换了个计算时间。假如已知一条路径用时 k，我们可以把代价计算反过来：在时间1计入k的代价，时间2计入k-1的代价…在时间k时计入1的代价
因为期望倒推，我们每次倒推的时候，实际上是在路径上加了一个起点u，这样操作，增加的代价就是$t[u]$

---

下面是代码

​#include <queue>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;

int Cas , head[100005] , tp , N , M , cnt_out[100005] , out[100005] ;
double t[100005] , f[100005] ;
vector<int> fr[100005] ;
queue<int> que ;

struct Path{
	int pre , to ;
} p[200005] ;

void In( int t1 , int t2 ){
	p[++tp] = ( Path ){ head[t1] , t2 } , head[t1] = tp ;
	cnt_out[t1] ++ , out[t1] ++ ;
}

void clear(){
	tp = 0 ;
	while( !que.empty() ) que.pop() ;
	memset( head , 0 , sizeof( head ) ) ;
	memset( t , 0 , sizeof( t ) ) ;
	memset( f , 0 , sizeof( f ) ) ;
	memset( cnt_out , 0 , sizeof( cnt_out ) ) ;
	memset( out , 0 , sizeof( out ) ) ;
}

void preWork(){
	scanf( "%d%d" , &N , &M ) ;
	for( int i = 1 ; i <= N ; i ++ ) fr[i].clear() ;
	for( int i = 1 , u , v ; i <= M ; i ++ ){
		scanf( "%d%d" , &u , &v ) ;
		In( u , v ) , fr[v].push_back( u ) ;
	}
}

void get_t( int u ){
	if( u == N ) return ;
	double tmp = 0 ;
	for( int i = head[u] ; i ; i = p[i].pre ) tmp += t[ p[i].to ] ;
	t[u] = ( tmp + out[u] + 1 ) / out[u] ;
}

void get_f( int u ){
	if( u == N ) return ;
	double tmp = 0 ;
	for( int i = head[u] ; i ; i = p[i].pre ) tmp += f[ p[i].to ] + t[u] ;
	f[u] = ( tmp + t[u] ) / out[u] ;
}

void solve(){
	que.push( N ) ;
	while( !que.empty() ){
		int u = que.front() ; que.pop() ;
		get_t( u ) ;
		get_f( u ) ;
	//	printf( "%d : %lf %lf\n" , u , t[u] , f[u] ) ;
		for( int i = fr[u].size() - 1 ; i >= 0 ; i -- ){
			if( !( -- cnt_out[ fr[u][i] ] ) ) que.push( fr[u][i] ) ;
		}
	}
	printf( "%.2lf\n" , f[1] ) ;
}

int main(){
	scanf( "%d" , &Cas ) ;
	for( int i = 1 ; i <= Cas ; i ++ ){
		clear() ;
		preWork() ;
		solve() ;
	}
}