Robots 2019南京网络赛 （概率dp）

Robots

[Time Limit: 1000 ms quad Memory Limit: 262144 kB ]

题意

有一个机器人要从 (1) 点走到 (n) 点，每走一步都需要一天，并且这一步的代价是已经经过的天数。

思路

将问题拆成两个问题，(dp1[i]) 表示从 (i) 到 (n) 的天数期望，(dp2[i]) 表示从 (i) 到 (n) 的代价期望。
那么就很容易得到

[dp1[u] = frac{sum_{u->v}(dp1[v]) + dp1[u]}{u.size()+1} + 1\ dp1[u] = frac{sum_{u->v}(dp1[v]) + u.size()+1}{u.size()} \ dp2[u] = frac{sum_{u->v}(dp2[v]) + dp2[u]}{u.size()+1} + dp1[u]\ dp2[u] = frac{sum_{u->v}(dp2[v]) + (u.size()+1)*dp1[u]}{u.size()} ]

/*************************************************************** 
    > File Name    : a.cpp
    > Author       : Jiaaaaaaaqi
    > Created Time : Mon 09 Sep 2019 10:07:54 PM CST
 ***************************************************************/

#include <map>
#include <set>
#include <list>
#include <ctime>
#include <cmath>
#include <stack>
#include <queue>
#include <cfloat>
#include <string>
#include <vector>
#include <cstdio>
#include <bitset>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define  lowbit(x)  x & (-x)
#define  mes(a, b)  memset(a, b, sizeof a)
#define  fi         first
#define  se         second
#define  pb         push_back
#define  pii        pair<int, int>

typedef unsigned long long int ull;
typedef long long int ll;
const int    maxn = 1e5 + 10;
const int    maxm = 1e5 + 10;
const ll     mod  = 1e9 + 7;
const ll     INF  = 1e18 + 100;
const int    inf  = 0x3f3f3f3f;
const double pi   = acos(-1.0);
const double eps  = 1e-8;
using namespace std;

int n, m;
int cas, tol, T;

vector<int> vv[maxn], rv[maxn];
int ind[maxn];
double dp1[maxn], dp2[maxn];

void solve() {
	queue<int> q;
	dp1[n] = dp2[n] = 0;
	q.push(n);
	while(!q.empty()) {
		int u = q.front();
		q.pop();
		if(u != n) {
			int sz = rv[u].size();
			{
				double ans = 0;
				for(auto v : rv[u]) {
					ans += dp1[v];
				}
				dp1[u] = 1.0*(ans+sz+1)/sz;
			}
			{
				double ans = 0;
				for(auto v : rv[u]) {
					ans += dp2[v];
				}
				dp2[u] = 1.0*(ans+1.0*(sz+1)*dp1[u])/sz;
			}
		}
		for(auto v : vv[u]) {
			ind[v]--;
			if(ind[v] == 0)	q.push(v);
		}
	}
}

int main() {
	// freopen("in", "r", stdin);
	scanf("%d", &T);
	while(T--) {
		scanf("%d%d", &n, &m);
		for(int i=1; i<=n; i++) {
			vv[i].clear();
			rv[i].clear();
			dp1[i] = dp2[i] = 0;
		}
		for(int i=1, x, y; i<=m; i++) {
			scanf("%d%d", &x, &y);
			rv[x].pb(y);
			vv[y].pb(x);
			ind[x]++;
		}
		solve();
		printf("%.2f
", dp2[1]);
	}
	return 0;
}