概率dp

题目：将一个数每次除以它的因子，求除到1的期望

 思路：dp[1]=0;dp[i] = （dp[i] + sum(dp[j])(所有因子dp和）+ sum（因子个数和）+1）/sum（因子个数）；

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <map>
#include <set>
#include <iomanip>

using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
#define maxn 100005
#define MOD 1000000007
#define mem(a , b) memset(a , b , sizeof(a))
#define LL long long
#define INF 100000000
int n;
double dp[maxn];

void init()
{
    for(int i = 3 ; i < maxn ; i ++)
    {
        double num = 0;   //记录i因子个数
        for(int j = 1 ; j*j <= i ; j ++ )            //求i因子个数
        {
            if(i % j == 0)
            {
                dp[i] += dp[j];      //期望次数必须累加
                num++;
                if(i / j != j && i / j != i)   //保证平方数不重复
                {
                    dp[i] += dp[i/j];
                    num++;
                }
            }
        }
        dp[i] = (dp[i] + num + 1) / num;   //累加了所有因子期望，每次的概率相等，取平均值
    }
}

int main()
{
    mem(dp , 0);
    dp[1] = 0.0;dp[2] = 2.0;
    init();
    int t;
    cin >> t;
    while(t-- )
    {
        scanf("%d" , &n);
        printf("%.6lf\n" , dp[n]);
    }
    return 0;
}