Race to 1 Again LightOJ - 1038（概率 dp）

最新推荐文章于 2023-03-29 10:28:27 发布

做一只大熊猫

最新推荐文章于 2023-03-29 10:28:27 发布

阅读量80

点赞数 1

分类专栏：概率 DP

本文链接：https://blog.csdn.net/qq_34261446/article/details/117395226

版权

概率 DP 专栏收录该内容

7 篇文章 0 订阅

订阅专栏

该博客探讨了一道涉及概率和动态规划的数学问题。题目描述了如何通过等概率选择因子的方式将一个数转化为1，求解期望的操作次数。博主通过初始化因子数组、预处理因子信息以及建立期望DP状态转移方程来解决这个问题。代码中展示了如何计算从每个数到1的期望操作次数，并给出了完整的C++实现。

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

题目链接

题意

给一个数 n ，每次操作这个 n 会等概率的变成它的一个因子，问经过多少次操作，这个数 n 变成 1，求期望操作次数

思路

一道很明显，有经典的期望 dp，
期望 dp 是从后往前递推而来，对于这题来说就是：从最终状态 1 -> 递推出 n 这个值期望的次数，
总感觉这样说不太对，其实我的理解就是，我们对当前的数 x 进行一次操作后 x 会变成的它的因子，那么当前的 dp 状态就要从它的子状态（即：因子的 dp 状态）转移而来。
我们设从 x 变到 1 期望次数为 dp [x],
举个转移的例子：对于 x = 6 的时候，6 的因子有 6、3、2、1 共 4 个，那状态转方程为：
因此我们可以预处出 1~1e5 中每个数的因子，然后从后往前进行 dp 状态转移

代码

#include <bits/stdc++.h>
using namespace std;
#define db  double
#define ll  long long
#define sc  scanf
#define pr  printf
#define fi  first
#define se  second
#define pb  push_back
#define m_p make_pair
#define Pir pair<int, int>
#define inf 0x3f3f3f3f
#define INF 0x3f3f3f3f3f3f3f3f
/*==========ACMer===========*/
const int N = 1e5 + 10;
vector<int> fac[N];

void init()
{
    for (int i = 1; i < N; i ++)
    {
        for (int j = 1; j * j <= i; j ++)
        {
            if (i % j == 0)
            {
                fac[i].pb(j);
                if (i / j != j) fac[i].pb(i / j);
            }
        }
    }
}

db dp[N];
void pre_do()
{
    init();
    for (int i = 2; i < N; i ++)
    {
        int m = fac[i].size();
        db p = 1.0 / (m - 1);
        dp[i] = m * p;
        for (auto x : fac[i])
        {
            if (x == i) continue;
            dp[i] += dp[x] * p;
        }
    }
}

int main()
{
    pre_do();
    int T, cas = 1; sc("%d", &T);
    while (T --)
    {
        int n; sc("%d", &n);
        pr("Case %d: %f\n", cas ++, dp[n]);
    }

    return 0;
}