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[];
 int prime[],tot;
 double f[];
 void init(){
     is[]=is[]=;
     for(LL i=;i<=;++i){
         if(!is[i]){
             prime[tot++]=i;
             for(LL j=i*i;j<=(LL);j+=i)
                is[j]=;
         }
     }
 }
 double dfs(int u){
     if(f[u]) return f[u];
     if(u==) return f[u]=;
     int g=,p=;
     f[u]=;
     for(int i=;i<tot&&prime[i]<=u;++i){
         p++;
         if(u%prime[i]==){
             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=;
     init();
     cin>>t;
     while(t--){
         memset(f,,sizeof(f));
         scanf("%d",&n);
         printf("Case %d: %.11f\n",++cas,dfs(n));
     }
     return ;
 }