LightOJ 1038 Race to 1 Again 概率DP

Race to 1 Again LightOJ - 1038

对于一个整数 $D$，每次从他的因数中选一个 $d$，然后令 $D=D/d$，直到 $D=1$ 为止，问从 $D$ 到 $1$ 执行这样操作的期望次数。

用 $E[i]$ 表示整数 $i$ 变为 $1$ 的期望次数，$n$ 表示因子的总个数，$a[k]$ 表示第 $k$ 个因子，第 $1$ 个因子 $a[1]=1$，第 $n$ 个因子 $a[n]=i$，则：

$$\begin{array}{rl}E[i]& =\frac{1}{n}\sum _{k=1}^n(E[a[k]]+1)\\ E[i]& =\frac{E[i]}{n}+\frac{1}{n}\sum _{k=2}^{n-1}E[a[k]]+1\\ \frac{n-1}{n}E[i]& =\frac{1}{n}\sum _{k=2}^{n-1}E[a[k]]+1\\ E[i]& =\frac{1}{n-1}(\sum _{k=2}^{n-1}E[a[k]]+n)\end{array}$$

代码如下：

#include<iostream>
#include<cstdio>
//#define WINE
#define MAXN 100010
using namespace std;
double d[MAXN];
int n,j,T,iCase;
int main(){
#ifdef WINE
    freopen("data.in","r",stdin);
#endif
    d[0]=d[1]=0;
    for(int i=2;i<MAXN;i++){
        n=0;
        for(j=2;j*j<i;j++){
            if(i%j==0){
                n+=2;
                d[i]+=d[j];
                d[i]+=d[i/j];
            }
        }
        if(j*j==i){
            n+=1;
            d[i]+=d[j];
        }
        n+=2;// 1 and i
        d[i]+=n;
        d[i]/=(n-1);
    }
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        printf("Case %d: %.6lf\n",++iCase,d[n]);
    }
    return 0;
}