概率dp

#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
);
}
return 0;
}