lightoj 1038 概率dp

最新推荐文章于 2020-07-15 15:50:24 发布

cadongllas

最新推荐文章于 2020-07-15 15:50:24 发布

阅读量208

点赞数

分类专栏：数学动态规划

本文链接：https://blog.csdn.net/cadongllas/article/details/52422749

版权

数学同时被 2 个专栏收录

9 篇文章 0 订阅

订阅专栏

动态规划

9 篇文章 0 订阅

订阅专栏

题意：
给出一个数字Ｄ
我们可以选择１－Ｄ中可以被Ｄ整除的数字，然后用Ｄ出得到一个新的数字Ｄ‘；
然后在找所有Ｄ’的因子，用Ｄ‘除，直到得到１；

问除的次数的期望值

设一个数的约数有m个，E[n] = E[a[1]]/m+E[a[2]]/m+...+E[a[m]]/m+1 (因为又除一次，所以+1)
整理得：E[n]=(E[a[1]]+E[a[2]]+...+E[a[m-1]]+m)/(m-1)

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <vector>
#define rep(i, j, k) for(int i = j; i <= k; i++)
#define maxn 100009

using namespace std;

int vis[maxn], n;
double dp[maxn];

void dfs (int x)
{
	if (x == 1)
	{
		dp[1] = 0.0;
		vis[1] = 1;
		return;
	}
	if (vis[x])
		return;
	vis[x] = 1;
	double tot = 1.0;
	int num = 0;
	for (int i = 1; i * i <= x; i++)
	if (x % i == 0)
	{
		num += i * i == x ? 1 : 2;
		if (vis[i] == 0)
			dfs (i);
		if (i * i != x)
		if (vis[x / i] == 0)
			dfs (x / i);
		if (i * i != x)
		tot += dp[x / i];
		tot += dp[i];
	}
	//printf ("dfs %d %d\n", x, num);
	dp[x] = tot / (1.0 * (num - 1) ) + 1.0;
	//printf ("%.3lf\n", dp[x]);
}

int main ()
{
	int ti;
	cin >> ti;
	rep (Time, 1, ti)
	{
		memset (vis, 0, sizeof (vis));
		memset (dp, 0, sizeof (dp));
		scanf ("%d", &n);
		dfs (n);
		printf ("Case %d: %.10f\n", Time, dp[n]);
	}
	return 0;
}