期望与概率DP专题/高斯消元

高斯消元

---

总结：概率和期望问题的求解都要依靠递推式，概率问题按照递推式进行dp，期望问题根据递推式化简的难易程度通过化简再dp或高斯消元求解，该类问题大部分仍属于各个类型的dp求解问题

2018.9.3-概率与期望dp学习

Discovering Gold

一排1到n的格子，每个格子上有黄金 ai ,你最开始在 1 号,每一次投骰子决定到哪一个格子，超出1~n范围则重新投掷，你到了哪个格子就得到哪个格子的金币，问最终在n 能得到金币的期望

#include<stdio.h>
#include<string.h>
const int MAX=1e2+5;
int a[MAX];
double dp[MAX];
int t,n;
int main()
{
    scanf("%d",&t);
    for(int tc=1;tc<=t;++tc)
    {
        memset(dp,0,sizeof(dp));
        scanf("%d",&n);
        for(int i=1;i<=n;++i) scanf("%d",&a[i]);
        dp[1]=1.0;
        for(int i=1;i<n;++i)
        {
            for(int j=1;j<=6&&i+j<=n;++j)
            {
                int tmp=6;
                if(i+6>n) tmp=n-i;
                dp[i+j]+=dp[i]/tmp;
            }
        }
        double ans=0;
        for(int i=1;i<=n;++i) ans+=dp[i]*a[i];
        printf("Case %d: %f\n",tc,ans);
    }
    return 0;
}

HDU3366 Passage

有n扇门和m枚金币，每扇门后有3种可能：1.出路概率p，2.另一扇门概率1-p-q，3.失去一枚金币并且出现另一扇门概率q，门可选择，求在最优选择情况下找到出路得概率

根据p/q从大到小对每扇门排序即为最优策略，注意q可能的值为0，需eps处理

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int MAX=1e3+5;
const double eps=1e-16;
double dp[MAX][15],ans;
int t,n,m;
struct P
{
    double p,q;
}a[MAX];
inline bool cmp(P& x,P& y)
{
    return (x.p+eps)/(x.q+eps)>(y.p+eps)/(y.q+eps);
}
int main()
{
    scanf("%d",&t);
    for(int tc=1;tc<=t;++tc)
    {
        memset(dp,0,sizeof(dp));
        scanf("%d%d",&n,&m);
        for(int i=0;i<n;++i) scanf("%lf%lf",&a[i].p,&a[i].q);
        sort(a,a+n,cmp);
        dp[0][m]=1.0;ans=0;
        for(int i=0;i<n;++i)
        {
            for(int j=m;j>=0;--j)
            {
                ans+=dp[i][j]*a[i].p;
                if(i<n-1&&j>=1)dp[i+1][j-1]+=dp[i][j]*a[i].q;
                if(i<n-1)dp[i+1][j]+=dp[i][j]*(1.0-a[i].p-a[i].q);
            }
        }
        printf("Case %d: %.5f\n",tc,ans);
    }
    return 0;
}

Island of Survival

小岛上有man\tiger\deer三种生物，每天都有两种相遇：1.tiger-man  2.tiger-tiger  3.tiger-deer  4.man-deer/deer  5.deer-deer

求man活下来的概率，当且仅当tiger==0。：本题求的是man能活下来的最终概率而不是某种情况下仍存活的概率，题中3\4\5种情况只能延长man的存活时间，而不能影响man的存活概率，因此只需要对1\2情况进行讨论。

#include<stdio.h>
#include<string.h>
const int MAX=1e3+5;
double dp[MAX];
int t,n,m;
int main()
{
    scanf("%d",&t);
    for(int tc=1;tc<=t;++tc)
    {
        memset(dp,0,sizeof(dp));
        scanf("%d%d",&n,&m);
        if(n==0) {printf("Case %d: 1\n",tc);continue;}
        if(n%2==1) {printf("Case %d: 0\n",tc);continue;}
        dp[n]=1;
        for(int i=n;i>=2;i-=2) dp[i-2]=dp[i]*(i-1)/(i+1);
        printf("Case %d: %.6f\n",tc,dp[0]);
    }
    return 0;
}

Bag of mice

有n只白鼠m只黑鼠，小D和小P二人随机各抓1只并且小D先手，小P抓完后会有1只小鼠出逃，求小D先抓到白鼠的概率

与上题相似，抓到白鼠则游戏结束，只需讨论2人都抓不到的概率和2种小鼠随机出逃的概率，答案为各种情况下抓到白鼠的概率之和：  本题我的解法中用eps忽略不会到达的情况，需要注意eps精度的设置，题中精度偏差不超过1e-9，eps=1e-11可AC

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int MAX=1e3+1;
const double eps=1e-11;
double dp[MAX][MAX];
int n,m;
int main()
{
    scanf("%d%d",&n,&m);
    memset(dp,0,sizeof(dp));
    dp[n][m]=1;
    double ans=0,tmp;
    for(int i=n;i>0;--i)
    {
        for(int