概率与期望 ———— T1 [LOJ 1083]

素影·流年

已于 2023-01-16 16:46:41 修改

阅读量106

点赞数

分类专栏：信息学竞赛文章标签：动态规划算法

于 2023-01-05 21:58:25 首次发布

本文链接：https://blog.csdn.net/Brain_Gym/article/details/128555444

版权

信息学竞赛专栏收录该内容

11 篇文章 2 订阅

订阅专栏

该问题讨论的是在整数N随机变成其因子的过程下，期望需要多少次才能变成1。通过动态规划的方法，计算每个整数变成1的期望步数，最终得出N的期望步数。代码示例给出了C++的解决方案，使用了动态规划数组dp来存储每个数的期望值。

摘要由CSDN通过智能技术生成

题目

给你一个整数 $N(1 <= N <= 10 ^ 5)$ ，每次 $N$ 会随机变成 $N$ 的某个因子，问期望几次会变成 $1$

分析

1. 画图理解期望

在这里插入图片描述

平均值

期望就是所有可能取值的平均值，即
$\sum_{i = 1}^{\infty} X_i$

概率

期望也可以理解为概率，即
$P_i 表示取到X_i的概率，则 \\ E(x) = \sum_{i = 1}^{\infty}P_i \times X_i$
$i\in{\pm1}$

2.数学式的表达（dp）

$dp_x表示整数x变成1的期望，n 表示 x 的因数个数则\\ dp_x = \frac{\sum_{i\in{i \equiv 0\pmod{x}}}dp_i}{n} + 1$
但是 $\in{i \equiv 0\pmod{x} }$ 所以不是dp,我们进行移向
得
$dp_x = \frac{(\sum_{i\in{i \equiv 0\pmod{x},i \neq x}}{dp_i}) + n}{n - 1}$
初始化 $dp_1 = 0$ 即可转移

#include <bits/stdc++.h>
using namespace std;
const int N = 1e5;
int n;
double dp[N];
int main() {
	scanf("%d", &n);
	dp[1] = 0;
	for (int i = 2; i <= n; i++) {
		double sum = 0;
		int cnt = 0;
		for (int j = 1; j * j <= i; j++) {
			if ( i % j == 0 ) {
				cnt++;
				sum += dp[j];
				if (j * j != i) {
					cnt ++;
					sum += dp[i / j];
				}
			}
		}
		sum += cnt;
		dp[i] = sum / (cnt - 1);
	}
	printf("%.2lf",dp[n]);
	return 0;
}