LightOJ - 1038 Race to 1 Again(概率DP)

Forwhy_

已于 2024-01-10 14:34:40 修改

阅读量188

点赞数 1

文章标签：数据结构动态规划算法 c++ 概率dp

于 2023-03-29 10:28:27 首次发布

本文链接：https://blog.csdn.net/weixin_62802134/article/details/129831249

版权

Race to 1 Again - LightOJ 1038 - Virtual Judge (vjudge.net)

题意

给定一个整数n,每次操作可以对当前数进行除以它的某个因子。除以哪个因子是随机的，求把n变成1的期望步数。

思路

每个数分解成它的因子，它的因子又能继续往下分解，最后分解为一，那么我们可设：

$p_{j}$ 代表j的因子个数
$dp_{j}$ 代表j分解为1的期望次数

则有状态转移方程

设j的因子为 $a_{i}$ ，其中最大因子为 $a_{k}$

$f_{j} = \frac{f_{a_{1}} + 1 + f_{a_{2}}+1+f_{a_{3}}+1+...+f_{a_{k}}+1}{p_{j}}$

代码

#include <bits/stdc++.h>
#define ll long long
namespace IO{
  inline int read(){
    int X=0, W=0;char ch=getchar();
    while(!isdigit(ch)) W|=ch=='-', ch=getchar();
    while(isdigit(ch)) X=(X<<1)+(X<<3)+(ch^48), ch=getchar();
    return W?-X:X;
  }
  inline void write(ll x){
    if(x<0) x=-x, putchar('-');
    if(x>9) write(x/10);
    putchar(x%10+'0');
  }
  inline void sprint(ll x){write(x);putchar(32);}
  inline void eprint(ll x){write(x);putchar(10);}
}
using namespace IO;
using namespace std;
double f[100100];
int p[100100];

signed main(){
  f[1] = p[1] = 0;
  for(int i = 1; i < 100100; i++){
    if(i > 1) f[i] = (f[i] + 1) / p[i] + 1;
    for(int j = i * 2; j < 100100; j += i) f[j] += f[i], p[j]++;
  }
  int t = read();
  for(int c = 1; c <= t; c++){
    int n = read();
    printf("Case %d: %.15lf\n", c, f[n]);
  }
}