蓝书——217. 绿豆蛙的归宿 拓扑排序/dfs 两种方法

概率dp经典套路：设f[i]为：第i个点到终点的期望步数。

对一个边u -> v.

u走到v点的概率是  1/deg[u].

走到v点再到终点的步数是：dis(u,v)+f[v].

到v点这个事件带来的期望的贡献是：1/deg[u]*(dis(u,v)+f[v]).

f[u]={任取v属于u的直接可达点}

初状态：f[n]=0;

求f[1].

任意的x->y   f[y]必须在f[x]之前就算出结果。

要么是反向图DAG的拓扑

要么直接dfs

 

dfs

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
//typedef __int128 LL;
//typedef unsigned long long ull;
//#define F first
//#define S second
typedef long double ld;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef pair<ld,ld> pdd;
const ld PI=acos(-1);
const ld eps=1e-9;
//unordered_map<int,int>mp;
#define ls (o<<1)
#define rs (o<<1|1)
#define pb push_back
//#define a(i,j) a[(i)*(m+2)+(j)]  //m是矩阵的列数
//pop_back()
const int seed=131;
const int M = 1e5+7;

int head[M],cnt;
void init(){cnt=0,memset(head,0,sizeof(head));}
struct EDGE{int to,nxt, val;}ee[M*2];
void add(int x,int y,int z){ee[++cnt].nxt=head[x],ee[cnt].to=y,ee[cnt].val=z,head[x]=cnt;}

double f[M];//每个点到终点的期望长度 输出f[1] 
int deg[M];//度数 
bool vs[M];
int n;
void dfs(int u)
{
	if(u==n)return ;//终点的f[u]必须 
	if(vs[u])return ;
	vs[u]=true;
	for(int i=head[u];i;i=ee[i].nxt)
	{
		int v=ee[i].to,w=ee[i].val;
		//cout<<u<<" "<<v<<" "<<w<<endl;
		dfs(v);
		f[u]+=1.0/deg[u]*(w+f[v]);
	}
}
int main()
{
//	ios::sync_with_stdio(false);
//  	cin.tie(0);
  	int m,u,v;
	int w;
  	cin>>n>>m;
  	for(int i=1;i<=m;i++)
  	scanf("%d%d%d",&u,&v,&w),add(u,v,w),deg[u]++;
  	dfs(1);
  //	for(int i=1;i<=n;i++)
  //	cout<<i<<"      "<<f[i]<<endl;
  	printf("%.2f\n",f[1]);
	return 0;
}

拓扑排序

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
//typedef __int128 LL;
//typedef unsigned long long ull;
//#define F first
//#define S second
typedef long double ld;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef pair<ld,ld> pdd;
const ld PI=acos(-1);
const ld eps=1e-9;
//unordered_map<int,int>mp;
#define ls (o<<1)
#define rs (o<<1|1)
#define pb push_back
//#define a(i,j) a[(i)*(m+2)+(j)]  //m是矩阵的列数
//pop_back()
const int seed=131;
const int M = 1e5+7;

int head[M],cnt;
void init(){cnt=0,memset(head,0,sizeof(head));}
struct EDGE{int to,nxt, val;}ee[M*2];
void add(int x,int y,int z){ee[++cnt].nxt=head[x],ee[cnt].to=y,ee[cnt].val=z,head[x]=cnt;}

double f[M];//每个点到终点的期望长度 输出f[1] 
int deg[M],du[M];//度数 
bool vs[M];

int main()
{

  	int n,m,u,v;
	int w;
  	cin>>n>>m;
  	for(int i=1;i<=m;i++)
  	scanf("%d%d%d",&u,&v,&w),add(v,u,w),deg[u]++,du[u]++;
  	queue<int>q;
  	q.push(n);
  	while(!q.empty()){
  		int u=q.front();q.pop();
  		for(int i=head[u];i;i=ee[i].nxt)
  		{
  			int v=ee[i].to,w=ee[i].val;
  			f[v]+=(f[u]+w)/deg[v];
  			du[v]--;
  			if(du[v]==0)q.push(v);
		}
	}
  	printf("%.2f\n",f[1]);
	return 0;
}