LightOJ 1038 Race to 1 Again [概率dp]

点击打开题目

题意: 任何一个大于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;
}