LightOJ1038 Race to 1 Again（概率期望）

题意：

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

数据范围：n<=1e5

解法：

期望推导裸题

令d[i]为i变成1的期望步数
显然d[i]由d[i的因子]转移而来

推式子:
d[i]=(d[f1]+d[f2]+...d[fm])/m+1
其中f1表示i的第一个因子,
假设共m个因子,fm表示i的最后一个因子，
显然fm=i,因为i本身就是自己最大的因子,
式子移项得:
d[i]-d[fm]/m=(d[f1]+d[f2]+...d[f(m-1)])/m+1
式子左边等于d[i]*(m-1)/m
将(m-1)/m乘到右边得:
d[i]=((d[f1]+d[f2]+...d[f(m-1)])/m+1)*(m/(m-1))

预处理每个数的期望,对于询问O(1)输出

code：

#include<bits/stdc++.h>
using namespace std;
const int maxm=1e5+5;
double d[maxm];
signed main(){
    d[1]=0;
    for(int i=2;i<maxm;i++){
        int cnt=0;
        double sum=0;
        for(int j=1;j*j<=i;j++){
            if(i%j==0){
                sum+=d[j],cnt++;
                if(i/j!=j)sum+=d[i/j],cnt++;
            }
        }
        d[i]=(sum/cnt+1)*cnt/(cnt-1);
    }
    int T;scanf("%d",&T);
    int cas=1;
    while(T--){
        int n;scanf("%d",&n);
        printf("Case %d: ",cas++);
        printf("%.6f\n",d[n]);
    }
    return 0;
}

---