Race to 1 Again LightOJ

Race to 1 Again LightOJ - 1038

题意：有一个数字D，每次把D变为它的一个因数（变到所有因数的概率相等，可能是本身），变到1后停止。求对于某个初始的D变到1的期望步数。

x的因子有p[1],...,p[k]

那么ans[x]=1/k*(ans[p[1]]+1)+...+1/k*(ans[p[k]]+1)
=1/k*(ans[p[1]]+...+ans[p[k-1]])+1/k*ans[p[k]]+1
(k-1)/k*ans[x]=1/k*(ans[p[1]]+...+ans[p[k-1]])+1
ans[x]=1/(k-1)*(ans[p[1]]+...+ans[p[k-1]])+k/(k-1)

#include<cstdio>
#include<cmath>
using namespace std;
int k,T,n;
double anss,ans[100100];
int main()
{
    int i,j,t;
    for(i=2;i<=100000;i++)
    {
        k=2;
        anss=0;
        for(j=2;j<sqrt(i);j++)
        {
            if(i%j!=0)    continue;
            anss=anss+ans[j]+ans[i/j];
            k+=2;
        }
        t=sqrt(i);
        if(t*t==i)
        {
            anss+=ans[t];
            k++;
        }
        ans[i]=anss/(k-1)+(double)k/(k-1);
    }
    scanf("%d",&T);
    for(i=1;i<=T;i++)
    {
        scanf("%d",&n);
        printf("Case %d: %.9f
",i,ans[n]);
    }
    return 0;
}