【LightOj 1038】 Race to 1 Again （概率+dp）

题意

有一个数 $D$ 每次除以他的一个因数得到一个新的 $D'$，问 $D$ 变成1的次数的期望

题解

$dp[i]=\frac{dp[k]}{约数个数}+1 \quad k$为 $i$ 的约数
$(1-\frac{1}{约数个数})dp[i]=\frac{dp[k]}{约数个数}+1 \quad k$为 $i$ 的约数($k ≠ i$)
$\frac{约数个数-1}{约数个数}dp[i]=\frac{dp[k]+约数个数}{约数个数} \quad k$为 $i$ 的约数($k ≠ i$)
$dp[i]=\frac{dp[k]+约数个数}{约数个数-1} \quad k$为 $i$ 的约数($k ≠ i$)

代码

#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std; 
#define N 100010
int n,T;
double dp[N];
int main(){
    scanf("%d",&T);
    dp[1]=0;
    for(int i=2;i<=100000;i++){
        int num=0;
        for(int j=1;j<=sqrt(i);j++){
            if(i%j==0){
                dp[i]+=dp[j];num++;
                if(i/j!=j) dp[i]+=dp[i/j],num++;
            }
        }
        dp[i]=(dp[i]+num)/(num-1);
    }
    for(int tt=1;tt<=T;tt++){
        scanf("%d",&n);
        printf("Case %d: %lf\n",tt,dp[n]);
    }
    return 0;
}