马尔科夫过程（概率dp）（北理16校赛）

时间限制1秒 内存限制64M

题目描述：

        给一个正整数N，每次可以在不超过N的素数中随机选择一个P，如果P是N的约数，则把N变成N/P，否则N不变。问把N最终变成1的期望次数是多少？

 

输入格式：

        输入第一行为数据组数T ( T <= 100 )，以下T行每行一个正整数N ( 1 <= N <= 1000000 )。

 

输出格式：

        对于每组数据，输出“Case #x: y”,其中x是指第x个用例，y是把N最终变成1的期望次数，输出保留三位小数。

 

样例输入：

        2

        3

        13

 

样例输出：

        Case #1: 2.000

        Case #2: 6.000

 

样例解释：

        对于第二组样例，小于等于13的素数有2、3、5、7、11、13，一共6个，其中是13的约数的只有13，随机选中13的概率是1/6，有5/6概率保持不变，不论选中的是否是N的约数都将计次。

题目大意：

              一个数N，每次有相同概率选择比他的素数，如果选的数字是他的因数，那么就把N/这个数，问最后得到1的期望选择次数

解题思路：

             假设E(n)为把n变为1选择次数的期望，他就等于从那些素数选择一个的概率*（选完这个素数之后n变为的那个数它的期望值+1）

             如果把n写为p1^r1+p2^r2+...+pk^rk，小于等于n的数除去p1,p2,..pk得到的数共有cnt个，那么E(n)=1/k*（(1+E(n/p1)）+（1+E(n/p2))+...+(1+E(n/pk))+cnt*(1+E(n)))。

             通过化简和移项可以得到左边是E(n)右边是E(n/pi)的式子，那么就可以用素数筛法筛出素数，之后每次求解取小于等于n的素数代入公式即可。

            可以用dp或者记忆化搜索来做，dp时间复杂度O(sqrt(n)*n)时间太长了，记忆化搜索可以避免求解没用的值，所以这道题用记忆化搜索做的。

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=1e6+200;
int p[maxn],ip[maxn],no=0;
double dp[maxn];
void init()
{
	for(int i=0;i<=1000006;i++)
		ip[i]=1;
	ip[0]=ip[1]=0;
	for(int i=2;i<=1000006;i++)
	{
		if(ip[i])
		{
			p[no++]=i;
			for(int j=2*i;j<=1e6;j+=i)
				ip[j]=0;
		}
	}
}
double f(int n)
{
	if(n==1)return dp[1]=0;
	if(dp[n]!=0)return dp[n];
	int cnt=0;
	double tem=0;
	int k;
	for( k=0;p[k]<=n;k++)
	{
		if(n%p[k]==0)tem+=1+f(n/p[k]);
		else 
		{
			tem+=1;
			cnt++;
		}
	}
	return dp[n]=tem/(k-cnt);
}
int main()
{
	int t,ca=1,n;
	scanf("%d",&t);
	init();
	dp[2]=1;dp[3]=2;
	dp[1]=0;
	while(t--)
	{
		scanf("%d",&n);
		printf("Case #%d: %.3lf\n",ca++,f(n));
	}
	return 0;
}