题意: 任何一个大于1的整数, 经过若干次除以自己的因子之后可以变为1, 求该变换字数的数学期望值。
思路: 若n的约数有1, a[0], a[1], ... , a[i], n; 即共有i + 2个约数;
则 E(n) = (E(1) / (i + 2) + E(a[0]) / (i + 2) + E(a[1]) / i + 2) + ... + E(a[i]) / (i + 2) + E(n) / (i + 2) + 1);
移项可得 E(n) = (i + 2) / (i + 1) * (E(1) / (i + 2) + E(a[0]) / (i + 2) + E(a[1]) / i + 2) + ... + E(a[i]) / (i + 2) + 1));
其中E(1) = 0; 递推即可求解
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<queue>
#include<vector>
#include<cmath>
#include<functional>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 10;
const double eps = 1e-9;
const int inf = 1 << 30;
vector<int> v[maxn];
int vis[maxn];
double dp[maxn];
void init() {
for(int i = 2; i < maxn; i++) {
for(int j = i + i; j < maxn; j += i) {
v[j].push_back(i);
//vis[j] = 1;
}
}
memset(dp, 0, sizeof dp);
for(int i = 2; i < maxn; i++) {
for(int j = 0; j < v[i].size(); j++) {
dp[i] += 1.0 * dp[v[i][j]] / (v[i].size() + 2);
}
dp[i] += 1;
dp[i] = dp[i] * (v[i].size() + 2.0) / (v[i].size() + 1.0);
}
// for(int j = 0; j < v[100000].size(); j++) {
// printf("%d ", v[100000][j]);
// }
// for(int i = 2; i < 20; i++) {
// printf("i = %d ", i);
// for(int j = 0; j < v[i].size(); j++) {
// printf("%d ", v[i][j]);
// }
// printf("\n");
// }
}
int main() {
init();
int t, v = 1, n;
scanf("%d", &t);
while(t--) {
scanf("%d", &n);
printf("Case %d: %.8f\n", v++, dp[n]);
}
return 0;
}