概率

条件概率: P(B|A) = P(AB) / P(A) 为事件A发生条件下，事件B发生的条 件概率

全概率公式：

 

概率期望 https://blog.csdn.net/consciousman/article/details/73730834
概率DP https://www.cnblogs.com/kuangbin/archive/2012/10/02/2710606.html
概率DP https://blog.csdn.net/myjs999/article/details/81022546

https://ac.nowcoder.com/acm/contest/358/B

条件概率简答题

设A,B是两个事件，且P(B)>0,则在事件B发生的条件下，事件A发生的条件概率（conditional probability)为：

P(A|B)=P(AB)/P(B)

p(第二次穿女装)=p(两次穿女装)/p(第一次穿女装)

ac:

#include<bits/stdc++.h>
using namespace std;
#define ll long long
ll n,m,x,y,t;
 
ll gcd(ll a,ll b)
{
    return b==0?a:gcd(b,a%b);
}
 
int main()
{
    scanf("%lld%lld%lld%lld%lld",&n,&m,&x,&y,&t);
    ll a=x*m*(x-1)*(m-1)+y*n*(y-1)*(n-1);
    ll b=(n-1)*(m-1)*(x*m+y*n);
    if(t==0)
        printf("%.3lf\n",1.0*a/b);
    else
        printf("%lld/%lld\n",a/gcd(a,b),b/gcd(a,b));
    return 0;
}

1.https://vjudge.net/problem/POJ-2151

题意:

有t个队伍,m道题

给出每个队做出每题的概率p[i][j]

问每个队都做出至少1题并且冠军队做出n题的概率

解析:

冠军队是做题最多的,转化为存在一个队做出n题及以上的概率

求解转化为:每队至少做1题-每队仅仅做(1~n-1)题的概率

dp[i][j][k]为第i队在前j题做出k道题的概率

初始化:dp[i][0][0]=1.0;

显然:dp[i][j][k]=dp[i][j-1][k-1]*p[i][j]+dp[i][j-1][k]*(1-p[i][j]);//做出第j+未做出第j

注意dp[i][j][0]的情况

ac:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define MAXN 1005
using namespace std;
double p[MAXN][MAXN];
double dp[MAXN][31][31];//第i队在前j题解出k题的概率

int m,t,n;

void init()
{
    memset(dp,0,sizeof(dp));
    for(int i=1;i<=t;i++)
    {
        dp[i][0][0]=1.0;
        for(int j=1;j<=m;j++)
            dp[i][j][0]=dp[i][j-1][0]*(1.0-p[i][j]);

        for(int j=1;j<=m;j++)
            for(int k=1;k<=j;k++)
                dp[i][j][k]=dp[i][j-1][k-1]*p[i][j]+dp[i][j-1][k]*(1-p[i][j]);
    }
}

int main()
{
    while(cin>>m>>t>>n)
    {
        if(m==0&&n==0&&t==0)
            break;
        for(int i=1;i<=t;i++)
            for(int j=1;j<=m;j++)
                cin>>p[i][j];
        init();

        double p1=1.0,p2=1.0;
        for(int i=1;i<=t;i++)//每队至少做出1题
            p1=p1*(1-dp[i][m][0]);

        for(int i=1;i<=t;i++)
        {
            double f=0;//第i队做出1~n-1题的概率
            for(int j=1;j<=n-1;j++)
                f+=dp[i][m][j];
            p2*=f;//每个队做出1~n-1题的概率
        }
        printf("%.3f\n",p1-p2);
    }
    return 0;
}

2.https://vjudge.net/problem/LightOJ-1079

题意:

给定n个银行,每个银行有a[i]元,如果去抢的话有p[i]的概率被抓

求被抓概率不大于v的情况下,最多抢多少钱?

解析:

概率01背包

一般01背包是以体积为背包容量,但是这里概率是浮点数,不是很方便

呢么我们可以将背包容量设置为钱数,值为抢到该值钱下逃脱的概率

得抢一定钱的情况下,最大的逃生概率

ac:

#include<bits/stdc++.h>
#define MAXN 10005
using namespace std;
double dp[MAXN],p[101];
int a[101];

int main()
{
    int t,n,cas=1;
    double v;
    scanf("%d",&t);
    while(t--)
    {
        memset(dp,0,sizeof(dp));
        scanf("%lf%d",&v,&n);
        int sum=0;
        for(int i=1;i<=n;i++)
            scanf("%d%lf",&a[i],&p[i]),sum+=a[i];

        dp[0]=1.0;//不抢逃脱的概率
        for(int i=1;i<=n;i++)
            for(int j=sum;j>=a[i];j--)
                    dp[j]=max(dp[j],dp[j-a[i]]*(1.0-p[i]));

        int maxs=0;
        for(int i=sum;i>=0;i--)
        {
            if(dp[i]>=1.0-v)
            {
                printf("Case %d: %d\n",cas++,i);
                break;
            }
        }
    }
    return 0;
}

3.Robot

题意:

一个环形的数字圈,起始点为1,每次走x步

随机向前向后,概率各位0.5,询问y点的概率

解析:

简单的一道概率dp题

概率dp一般从前往后推

这题数据非常大,要用滚动数组优化

注意顺序,最优子结构,处理好边界

ac:

#include<bits/stdc++.h>
#define ll long long
#define MAXN 505
using namespace std;
double dp[2][MAXN]={0};
//MAXN非常大,但是计算的时候i只会对i+1产生影响
//所以我们用滚动数组优化空间

int main()
{
    int n,m,l,r,x;
    while(scanf("%d%d%d%d",&n,&m,&l,&r)!=EOF)
    {
        if(n==0&&m==0&&l==0&&r==0)
            break;
        for(int i=0;i<n;i++)
            dp[0][i]=0;
        dp[0][0]=1.00;//从0出发
        int now=0;
        while(m--)
        {
            scanf("%d",&x);
            now=now^1;//滚动数组优化
            for(int i=0;i<n;i++)//滚动数组,重新计算,要清零
                dp[now][i]=0;
            for(int i=0;i<n;i++)//每个点的机器人都移动一次
            {
                if(dp[now^1][i]==0)
                    continue;
                dp[now][(i+x)%n]+=0.5*dp[now^1][i];//正向,取模后构成循环
                dp[now][(i-x+n)%n]+=0.5*dp[now^1][i];//反向
            }
        }
        double ans=0;
        for(int i=l-1;i<r;i++)
            ans+=dp[now][i];
        printf("%.4f\n",ans);
    }
    return 0;
}

4.https://vjudge.net/problem/UVALive-6672#author=0

题意:

你可以选择买icpc的票或者acm的票,每一轮选出一个成功者

icpc的票被选中的概率是acm票的两倍

但是这个概率是根据持票人数变化的

现在你知道买两种票的人数a,b和会选出n个成功的人

求选icpc票最终成功的概率和选择acm票最终成功的概率

解析:

我在选之前,一定是已经选定了买了一种票.

我们要算两种清空下的概率,这道题卡时间卡的非常紧,我们遍历一次完成.

dp[i][j]为我选择第icpc票,第i轮有j个人成功

dp2[i][j]为我选择第acm票,第i轮有j个人成功

ans1+=算出每种情况下的出现的概率*(我选择icpc后成功的概率)

ans2+=算出每种在情况下出现的概率*(我选择acm后成功的概率)

ac:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<set>
#include<map>
#include<vector>
#include<algorithm>
#include<cmath>
#include<queue>
#include<sstream>
#include<cstdlib>
using namespace std;

const int maxn=3000+5;
int n,a,b;
double dp[maxn][maxn];///dp[i][j]表示到第i轮时有j个人竞选成功
double dp2[maxn][maxn];///dp2[i][j]表示到第i轮有j个人竞选成功
double ans1,ans2;
int ua,ub,la,lb;
int main()
{
    while(~scanf("%d%d%d",&n,&a,&b))
    {
        if(n>a+b)//必定成功
        {
            printf("1.0000000000000000\n");
            printf("1.0000000000000000\n");
            continue;
        }
        ans1=0.0;ans2=0.0;
        memset(dp,0,sizeof(dp));
        memset(dp2,0,sizeof(dp2));
        dp[0][0]=1.0;
        dp2[0][0]=1.0;
        for(int i=0;i<n;i++)//第i轮
        {
            for(int j=0;j<=i&&j<=a;j++)//枚举选icpc成功人的个数,i-a就是选acm成功的人数了
            {
                if(fabs(dp[i][j])<1e-15) continue;
                ua=j;//成功的a
                ub=i-j;//成功的b
                la=a-ua;//原a人里现剩下的人数
                lb=b-ub;//原b人里现剩下的人数

                double fenmu=1.0/(la*2.0+lb+2.0);//自己该轮选icpc的成功概率
                ans1+=dp[i][j]*2.0*fenmu;//累加每种情况的下,我选择icpc后的成功概率
                if(la>0)
                    dp[i+1][j+1]+=dp[i][j]*la*2.0*fenmu;//la里有个人以icpc竞选成功
                if(lb>0)
                    dp[i+1][j]+=dp[i][j]*lb*1.0*fenmu;//lb里有个人以acm竞选成功

                fenmu=1.0/(la*2.0+lb+1.0);//自己该轮选择acm的成功概率
                ans2+=dp2[i][j]*fenmu;//累加每种情况的下,我选择acm后的成功概率
                if(la>0)
                    dp2[i+1][j+1]+=dp2[i][j]*la*2.0*fenmu;//同上
                if(lb>0)
                    dp2[i+1][j]+=dp2[i][j]*lb*fenmu;//同上
            }
        }
        printf("%.12f\n",ans1);
        printf("%.12f\n",ans2);
    }
    return 0;
}

http://acm.hdu.edu.cn/showproblem.php?pid=6693

题意:

有n个物品,每个物品有一定概率成功,让你选择若干件使得成功一次的概率最大

解析:

取两件概率:A*(1-B)+B*(1-A)=A+B-2AB,如果A大,呢么A超过0.5,呢么就只能取一件,取多件的概率和不会超过0.5

贪心取,每次取大的,用前缀和处理,一次算出结果

ac:

#include<bits/stdc++.h>
using namespace std;
double p[10001]={0};

int main()
{
    int t,n;
    scanf("%d",&t);
    while(t--)
    {
        double kk=0;
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            scanf("%lf",&p[i]),kk=max(kk,p[i]);
        if(kk>=0.500000000)
            printf("%.12f\n",kk);
        else{
            sort(p+1,p+n+1);
            double ans=p[n],sum=1-p[n];
            for(int i=n-1;i>=1;i--)
            {
                double tmp=ans*(1-p[i])+sum*p[i];
                if(tmp>ans)
                {
                    ans=tmp;
                    sum=sum*(1-p[i]);
                }
            }
            printf("%.12f\n",ans);
        }
    }
    return 0;
}

http://acm.hdu.edu.cn/showproblem.php?pid=6656

题意:

#include<bits/stdc++.h>
#define ll long long
#define mod 1000000007
#define MAXN 500005
using namespace std;
ll r[MAXN],s[MAXN],c[MAXN];
ll dp[MAXN];
ll sum[MAXN];
ll x[MAXN];

ll finv(ll a,ll b)
{
    ll res=1;
    b=b-2;
    while(b)
    {
        if(b&1)
            res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}

void init(ll n)
{
    memset(dp,0,sizeof(dp));
    memset(sum,0,sizeof(sum));
    for(int i=1;i<=n;i++)
    {
        ll kk=finv(r[i],mod);
        ll a=(((s[i]*kk)%mod)*c[i])%mod;
        ll b=(((((s[i]-r[i])%mod)*kk)%mod)*((sum[i-1]-sum[x[i]-1]+mod)%mod))%mod;
        dp[i]=(a+b+mod)%mod;
        sum[i]=(sum[i-1]+dp[i]+mod)%mod;
    }
}

int main()
{
    ll t,n,m,l,rr;
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld%lld",&n,&m);
        for(int i=1;i<=n;i++)
            scanf("%lld%lld%lld%lld",&r[i],&s[i],&x[i],&c[i]);
        init(n);
        while(m--)
        {
            scanf("%lld%lld",&l,&rr);
            printf("%lld\n",(sum[rr-1]-sum[l-1]+mod)%mod);
        }
    }
    return 0;
}