概率与数学期望

期望的线性性质：

        X Y为两个事件，a b为各自发生的概率，E为期望，则

基本做法：

        用图表示事件之间的转换关系，表示事件u的期望，表示事件u发生的概率，事件j表示事件u可以转移到的事件，那么

        采用记忆化搜索或者递推的方式进行求解，注意如果使用递推，需要从终点开始按拓扑序倒推，使用记忆化搜索从起点开始递归即可（推荐使用记忆化搜索，比较方便）

---

绿豆蛙的归宿 

217. 绿豆蛙的归宿 - AcWing题库

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int N=1e5+5;
int n,m;
int h[N],e[N<<1],ne[N<<1],w[N<<1],idx;
void add(int a,int b,int c){
	e[idx]=b;
	w[idx]=c;
	ne[idx]=h[a];
	h[a]=idx++;
}
int son[N];
double f[N];
double dp(int u){
	if(f[u]>=0) return f[u];
	f[u]=0;
	for(int i=h[u];~i;i=ne[i]){
		int j=e[i];
		f[u]+=(w[i]+dp(j))/son[u];
	}
	return f[u];
}
int main(){
	cin>>n>>m;
	memset(h,-1,sizeof h);
	while(m--){
		int a,b,c;cin>>a>>b>>c;
		add(a,b,c);
		son[a]++;
	}
	memset(f,-1,sizeof f);
	cout<<fixed<<setprecision(2)<<dp(1);

} 

---

扑克牌

218. 扑克牌 - AcWing题库

#include <bits/stdc++.h>
using namespace std;

const int N = 20;
const double INF = 1e20;

int A, B, C, D;
double f[N][N][N][N][5][5];

double dfs(int a, int b, int c, int d, int x, int y)
{
    double &v = f[a][b][c][d][x][y];
    if (v >= 0) return v;
    // 每堆已经使用的张数 
    int na = a + (x == 0) + (y == 0);
    int nb = b + (x == 1) + (y == 1);
    int nc = c + (x == 2) + (y == 2);
    int nd = d + (x == 3) + (y == 3);
    // 已经满足就不需要再搜了 
    if (na >= A && nb >= B && nc >= C && nd >= D) return v = 0;
    // 当前剩余牌数 
	int sum = a + b + c + d + (x != 4) + (y != 4);
    sum = 54 - sum;
    if (sum <= 0) return v = INF;
    v = 1;
    /*
    为什么这里初始化1？
    可以把 DP 数组放到图里，因为摸一张牌会到下一个状态，所以每条边长度为 1
    那么 f[i] = p[1] * (f[1] + 1) + p[2] * (f[2] + 1) + ... + p[k] * (f[k] + 1)
              = p[1] + p[2] + ... + p[k] + p[1] * f[1] + p[2] * f[2] + ... + p[n] * f[n]
    显然，p[1] + p[2] + ... + p[k] = 1，所以 v 要初始化 1
    */
    if (a < 13) v += (13.0 - a) / sum * dfs(a + 1, b, c, d, x, y);
    if (b < 13) v += (13.0 - b) / sum * dfs(a, b + 1, c, d, x, y);
    if (c < 13) v += (13.0 - c) / sum * dfs(a, b, c + 1, d, x, y);
    if (d < 13) v += (13.0 - d) / sum * dfs(a, b, c, d + 1, x, y);
    if (x == 4)
    {
        double mn = INF;
        for (int i = 0; i < 4; i ++ ) mn = min(mn, 1.0 / sum * dfs(a, b, c, d, i, y));
        v += mn;
    }
    if (y == 4)
    {
        double mn = INF;
        for (int i = 0; i < 4; i ++ ) mn = min(mn, 1.0 / sum * dfs(a, b, c, d, x, i));
        v += mn;
    }
    return v;
}

int main()
{
    scanf("%d%d%d%d", &A, &B, &C, &D);
    memset(f, -1, sizeof f);
    double res = dfs(0, 0, 0, 0, 4, 4);
    if (res > 60) res = -1.000;
    printf("%.3lf\n", res);
    return 0;
}