UVA11427+UVA11762——概率递推

觉得这两个蓝书上的题很有代表性，特意记录一下，蓝书P142-143页。

UVA11427——玩纸牌

题目链接：https://vjudge.net/problem/UVA-11427

题目大意：每天 晚上都会玩纸牌，假如每盘游戏获胜概率为p，各盘输赢是独立的。你只有在当晚赢的次数严格大于p时才会去睡觉，并且每晚最多玩n盘。如果获胜比率在最后都没有超过p，你只能垂头丧气睡觉并且之后都不会玩纸牌。问你玩纸牌的期望天数为多少。

很经典的一个题，你只有在一个晚上没有达标才会停止游戏，而每晚上游戏都是独立的，所以这是一个马尔科夫过程（不知道标书对不对。。。），所以我们只需要计算出每天获胜的概率，即可完成此题。

那么每天不获胜的概率为多少呢？可以发现这是一个组合计数的过程，满足所有比率小于p的过程都可以。所以可以根据这个递推，设dp[i][j]表示前i局中枚举结束后获胜比率均不超过p，且前i局共获胜j局的概率。则转移方程为

当j/i<=p时，dp[i][j]=dp[i-1][j]*(1-p)+dp[i-1][j-1]*p；其他dp[i][j]=0。

是不是感觉很像组合数的递推过程？其实就是。

最后，dp[n][0]+dp[n][1]+.....就是所求答案Q；

最后，假设第k天结束，那么概率就是Q*(1-Q)^(k-1)，那么只需要推一个无穷级数的和即可。

最后答案就是1/Q。

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

const int maxn=100+50;
double dp[maxn][maxn];
int n,a,b;
double p;
const double epx=1e-6;

int main(){
    int T;
    scanf("%d",&T);
    int index=1;
    while(T--){
        int n,a,b;
        double p;
        scanf("%d/%d%d",&a,&b,&n);
        p=(double)a/b;
        memset(dp,0,sizeof(dp));
        dp[0][0]=1.0;
        dp[0][1]=0.0;
        for(int i=1;i<=n;i++){
            for(int j=0;j*b<=a*i;j++){
                dp[i][j]=dp[i-1][j]*(1-p);
                if(j)dp[i][j]+=dp[i-1][j-1]*p;
            }
        }
        double Q=0.0;
        for(int j=0;j*b<=a*n;j++)Q+=dp[n][j];
        printf("Case #%d: %d\n",index++,(int)(1/Q));
    }
    return 0;
}

UVA11762——得到一

题目链接：https://vjudge.net/problem/UVA-11762

题目大意：给一个整数N，每次可以在不超过N的素数中随机选择一个P，如果P是N的约束，N/=P；否则N不变，问你将N变为1的期望操作次数为多少？

思路：蓝书上的思路：这题可以看做一个随机转移的状态机，每个状态的出发的所有概率之和为1，那么就可以算出所有转移的概率，这是一个马尔科夫的过程：转移出去之后的状态与之前无关。

设f[i]表示当前的数为i时接下来需要选择的次数，则可以列出以下转移方程：

             f[6]=1+f[6]*1/3+f[3]*1/3+f[2]*1/3;

边界f[1]=0。用全期望公式展开后可以得到一般公式：

设不超过x的素数个数为p[x]个，其中g[x]个是x的因子，则

因为x/y<x，可以用记忆化搜索推出f(x)，具体看代码。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <cmath>
#include <queue>
#include <vector>
using namespace std;

const int maxn=1000000+10;
vector<int>prime;
int vis[maxn];

void get_prime(){
    memset(vis,0,sizeof(vis));
    for(int i=2;i<maxn;i++){
        if(!vis[i]){
            prime.push_back(i);
        }
        for(int j=0;j<prime.size();j++){
            if(i*prime[j]>=maxn)break;
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
}

int vis2[maxn];
double f[maxn];
double dp(int x){
    if(x==1)return 0.0;
    if(vis2[x])return f[x];
    vis2[x]=1;
    double &ans=f[x];
    ans=0.0;
    int g=0,p=0;//累加g(x)与p(x)
    for(int i=0;i<prime.size()&&prime[i]<=x;i++){
        p++;
        if(x%prime[i]==0){
            g++;
            ans+=dp(x/prime[i]);
        }
    }
    ans=(ans+p)/g;
    return ans;
}
int main(){
    int T;
    scanf("%d",&T);
    int index=1;
    get_prime();
    memset(vis2,0,sizeof(vis2));
    while(T--){
        int n;
        scanf("%d",&n);
        printf("Case %d: %.10lf\n",index++,dp(n));
    }
    return 0;
}