UVA-11761-马尔可夫/记忆化搜索

https://vjudge.net/problem/UVA-11762

　　　　给出一个整数n，每次随机挑选一个小于等于n的素数，如果是n的因子，n变为n/x ,否则不变，问n变为1的期望挑选次数。

　　f[i]=1/(m1+m2)*(SUM{ f[i] } + SUM{ f[i/x] }) +1,化简后记忆化搜索就好了，筛素数的时候要保存所有的素数所以不能根号优化。

　　

#include<iostream>
#include<cstring>
#include<queue>
#include<cstdio>
#include<stack>
#include<set>
#include<map>
#include<cmath>
#include<ctime>
#include<time.h> 
#include<algorithm>
using namespace std;
#define mp make_pair
#define pb push_back
#define debug puts("debug")
#define LL long long 
#define pii pair<int,int>
#define eps 1e-10
bool is[1001000];
int prime[100000],tot;
double f[1001000];
void init(){
    is[0]=is[1]=1;
    for(LL i=2;i<=1000000;++i){
        if(!is[i]){
            prime[tot++]=i;
            for(LL j=i*i;j<=(LL)1000000;j+=i)
               is[j]=1;
        }
    }
}
double dfs(int u){
    if(f[u]) return f[u];
    if(u==1) return f[u]=0;
    int g=0,p=0;
    f[u]=0;
    for(int i=0;i<tot&&prime[i]<=u;++i){
        p++;
        if(u%prime[i]==0){
            f[u]+=dfs(u/prime[i]);    
            g++;
        }
    }
f[u]=(f[u]+p)/g;
    return f[u];
}
int main()
{
    int n,m,i,j,k,t;
    int cas=0;
    init();
    cin>>t;
    while(t--){
        memset(f,0,sizeof(f));
        scanf("%d",&n);
        printf("Case %d: %.11f
",++cas,dfs(n));
    }
    return 0; 
}