Race to 1 Again LightOJ - 1038

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[];
 int main()
 {
     int i,j,t;
     for(i=;i<=;i++)
     {
         k=;
         anss=;
         for(j=;j<sqrt(i);j++)
         {
             if(i%j!=)    continue;
             anss=anss+ans[j]+ans[i/j];
             k+=;
         }
         t=sqrt(i);
         if(t*t==i)
         {
             anss+=ans[t];
             k++;
         }
         ans[i]=anss/(k-)+(double)k/(k-);
     }
     scanf("%d",&T);
     for(i=;i<=T;i++)
     {
         scanf("%d",&n);
         printf("Case %d: %.9f\n",i,ans[n]);
     }
     return ;
 }