概率专题·期望

最新推荐文章于 2022-05-14 09:12:08 发布
最新推荐文章于 2022-05-14 09:12:08 发布
阅读量655 收藏
点赞数
分类专栏： # 系列 期望
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/nike0good/article/details/50446471
版权

uva 11021 Tribbles

k只麻球，每活一天就会死亡，但第二天可能会生一些麻球，具体说来，生i个麻球的概率为pi ,求m天后所有麻球都死亡的概率
fi=p0+p1fi−1+p2fi−22+...+pn−1fi−1n−1 

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
#include<iomanip> 
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=next[p])  
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define vi vector<int> 
typedef long long ll;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define MAXN (1000+10)
int T,n,k,m;
double p[MAXN],f[MAXN];
int main()
{
//  freopen("uva11021.in","r",stdin);
//  freopen(".out","w",stdout);

    T=read();
    For(kcase,T) {
        n=read(),k=read(),m=read();
        Rep(i,n) cin>>p[i];
        f[0]=0; f[1]=p[0];
        Fork(i,2,m) {
            f[i]=0;
            Rep(j,n) f[i]+=p[j]*pow(f[i-1],j);                      
        }
        printf("Case #%d: %.7lf\n",kcase,pow(f[m],k));      
    }
    return 0;
}

UVA 11722 Joining with Friend

题意：假设你会在时间区间[t1,t2]中的任意时刻以相同的概率到达火车站，你朋友则在[s1,s2]中的任意时刻以相同的概率到达火车站,你们在到达后停留w分钟，问你们会面的概率。若存在一个时刻你们同时在火车站视为一定会面。
题解：概率空间是坐标系中的一个矩形，分情况讨论算面积。

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
ld t1,t2,s1,s2,w;
int main()
{
//  freopen("uva11722.in","r",stdin);
//  freopen(".out","w",stdout);

    int T=read();
    For(kcase,T) {
        cin>>t1>>t2>>s1>>s2>>w;
        if (s2-t2>s1-t1) swap(t1,s1),swap(t2,s2);
        ld a[5]={0,s1-t2,s2-t2,s1-t1,s2-t1};
        ld S=(t2-t1)*(s2-s1),p=0;

        if (-w<a[1]) {
        } else if (-w<a[2]) {
            p+=pow(t2-w-s1,2)/2;    
        } else if (-w<a[3]) {
            p+=(t2-s1-w+t2-s2-w)*(s2-s1)/2;
        } else if (-w<a[4]) {
            p+=S-pow(s2+w-t1,2)/2;
        } else p+=S;
        if (w<a[1]) {
            p+=S;
        } else if (w<a[2]) {
            p+=S-(t2+w-s1)*(t2+w-s1) /2;    
        } else if (w<a[3]) {
            p+=S- (t2-s1+w+t2-s2+w)*(s2-s1)/2;
        } else if (w<a[4]) {
            p+=pow(s2-w-t1,2)/2;
        } 
        printf("Case #%d: ",kcase);
        cout<<(S-p)/S<<endl;
    }

    return 0;
}

UVA 11427 Expect the Expected

题意：每天晚上你会玩牌。第一把赢了就去睡觉，输了继续。你每一把赢的概率是p，当你赢的比例严格大于p时，才会去睡觉。但是你一个晚上最多玩n把游戏，如果哪天获胜比例一直没有大于p，你以后再也不玩了。求玩游戏天数的期望。

假设单独一天获胜概率为Q
则 Ex=Q+2(1−Q)Q+3(1−Q)2Q+...
高中数学作差法可得 Ex=1/Q

计算Q可以dp， dpi,j表示前i天赢j局时，最后‘垂头丧气′的概率
dp[i,j]=dp[i−1,j](1−p)+dp[i−1,j−1]p(j/i≤p) 

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
typedef long long ll;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
const int maxn = 100 + 10;
int n;
double dp[maxn][maxn];
int main()
{
//  freopen("uva11427.in","r",stdin);
//  freopen(".out","w",stdout);

    int T=read();
    For(kcase,T) {
        int a,b;
        scanf("%d/%d %d",&a,&b,&n);
        double p=(double)a/b;

        MEM(dp) dp[0][0]=1;
        For(i,n)
            Rep(j,i) 
            {
                if ((double)j/i>p) break;
                dp[i][j]=dp[i-1][j]*(1-p)+( j ? dp[i-1][j-1]*p : 0 ); 
            } 
        double ans=0;
        Rep(j,n) ans+=dp[n][j];
        printf("Case #%d: %d\n",kcase,(int)(1/ans)); 


    }


    return 0;
}

UVA 11762 Race to 1

其实直接用马尔可夫过程就好

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
typedef long long ll;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
const int MAXN = (1000000+10);
int primes[MAXN],size=0;
bool b[MAXN]={0};
void make_prime() {
    int n=MAXN-10;
    Fork(i,2,n) {
        if (!b[i]) primes[++size]=i;
        For(j,size) {
            if (i*primes[j]>n) break;
            b[i*primes[j]]=1;
            if (i%primes[j]==0) break;
        }
    }
}

double dp[MAXN]={0};
bool vis[MAXN]={0,1};
double calc(int x) {
    if (vis[x]) return dp[x];
    dp[x]=1; 
    int siz=0,t=0;
    double c=0;
    for(int i=1;primes[i]<=x&&i<=size;++i) {
        ++siz;
        if (x%primes[i]==0) ++t,c+=calc(x/primes[i]);
    }
    return dp[x] = (c+siz)/t;   
}


int main()
{
//  freopen("uva11762.in","r",stdin);
//  freopen(".out","w",stdout);
    make_prime();
    int T=read();
    For(kcase,T) {
        double ans=calc(read());
        printf("Case %d: %.8lf\n",kcase,ans);
    }


    return 0;
}

UVA 10491 So you want to be a 2n-aire?

很容易搞错的是，每次拿到题目后’我’马上知道答对概率，要根据这个概率决策的。 

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define eps (1e-9)
#define SI(a) ((a).size())
typedef long long ll;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
int n;
double t,p,dp[50];
void work() {
    dp[n]=(ll)1<<(ll)n;
    if (fabs(1-t)<eps) {
        cout<<dp[n]<<endl; 
        return;
    }
    RepD(i,n-1) {
        double p2=((ll)(1<<(ll)i));
        double f=p2 / dp[i+1];
        if (f<=t) dp[i]=p*dp[i+1];
        else dp[i] = ( (f-t)*p2 + (1-f)*(f+1) /2*dp[i+1] ) / (1-t); 
    }
    cout<<dp[0]<<endl;
}
int main()
{
//  freopen("uva10900.in","r",stdin);
//  freopen(".out","w",stdout);
    cout<<setiosflags(ios::fixed)<<setprecision(3);
    while(cin>>n>>t && n) { 
        p=(1+t)/2;
        work();
    }
    return 0;
}

UVA 11346 Probability

积分求面积，SB题

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
typedef long long ll;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
int main()
{
//  freopen("uva11346.in","r",stdin);
    int T=read();
    double a,b,S;
    while(cin>>a>>b>>S) {
        if (a*b<=S) {
            puts("0.000000%");
        } else if (S<1e-7) {
            puts("100.000000%");
        }else {
            printf("%.6lf%\n",100-(1+log(a)-log(S/b))*S/a/b*100);
        }
    }   


    return 0;
}

UVA 11637 Garbage Remembering Exam

期望具有线性性质，原题可以拆成每个单词无效无效概率的和，取反计算有效概率，线状排列中与第i个有效的概率为与在线排列中其所有距离不超过k的在都在环中超过k个的概率，暴力统计

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
typedef long long ll;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
const int MAXN = 100000 + 10;
int k,n;
double p[MAXN];
double solve(int n,int k) {
//  {2*k+1..n-1-2k}
    MEM(p) 
    int x=2*k+1;
    p[n]=1;
    For(i,2*k)
    {
        p[n]=p[n]*(n-2*k+i)/(n-2*k-1+i);
    }
    ForkD(X,x,n) {
        p[X-1]=p[X]*(X-2*k)/X; 
    }   
    double ans=0;
    For(i,n) {
        int len=n-(min(i+k,n)-max(i-k,1)+1);
        ans+=1-p[len];
    } 
    return ans;
}
int main()
{
//  freopen("uva11637.in","r",stdin);
//  freopen(".out","w",stdout);

    int kcase=1;
    while (cin>>n>>k && n && k) {
        double ans=0;
        if (n==1) ans=0;
        else if (n-1-2*k>0) {
            ans=solve(n,k);
        } else {
            ans=n;
        }   
        printf("Case %d: %.4lf\n",kcase++,ans);
    }
    return 0;
}

UVA 11605 Lights inside a 3d Grid

用线性期望转成所有灯开的概率，由于这题是用xor和
根据
{(1−p+p)n=C0npn+C1n−1(1−p)pn−1+...+Cnn(1−p)n(1−p−p)n=(−1)nC0npn+(−1)n−1C1n−1(1−p)pn−1+...+Cnn(1−p)n
可以求出 C1n(1−p)n−1p+C3n(1−p)n−3p3+...+Ckn(1−p)n−kpk∣k是奇数 

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
typedef long long ll;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
int M,N,P,K;
void calc(int i,int j,int k) {

}
int main()
{
//  freopen("uva11605.in","r",stdin);
//  freopen(".out","w",stdout);

    int T=read();
    For(kcase,T) {
        double ans=0;
        cin>>N>>M>>P>>K;
        For(i,N)
            For(j,M)
                For(k,P) {
                    double p=(2.0*i*(N-i+1)-1)*(2*j*(M-j+1)-1)*(2*k*(P-k+1)-1)/(double)N/N/M/M/P/P;
                    ans+=1-pow(1.0-2*p,K);
                }
        printf("Case %d: %.10lf\n",kcase,ans/2);
    }

    return 0;
}

UVA 10491 Cows and Cars

注意所有门是一次打开的

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
typedef long long ll;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
int main()
{
//  freopen("uva11181.in","r",stdin);
//  freopen(".out","w",stdout);

    double a,b,c;
    while(cin>>a>>b>>c) {
        cout<<setiosflags(ios::fixed)<<setprecision(5);
        cout<< (a*b+b*(b-1))/(a+b)/(a+b-c-1)<<endl;
    }

    return 0;
}

UVA 1413 Game

卡精度问题，首先可以用区间dp算出一个某个区间的某侧人拿球时传左，右的概率，然后考虑第k个人拿球时，出现某个区间某侧人拿球的期望，只要n-1个人拿球了，后面无论如何都是n赢

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
const int MAXN = 50 +10;
int n,k;
double p[MAXN]; 
const int L = 0, R = 1;
ld f[MAXN][MAXN][2][2]; // [left,Right,Ball,Excepted Direction]

ld dp[MAXN][MAXN][2]; // P[Left,Right,Excepted Direction]

int main()
{
//  freopen("uva1413.in","r",stdin);
//  freopen(".out","w",stdout);

    while(scanf("%d%d",&n,&k)==2) {
        MEM(p) MEM(f) MEM(dp)

        For(i,n-1) scanf("%lf",&p[i]);

        For(i,n-1) {
            f[i][i][L][L]=f[i][i][R][L]= 1 - p[i];
            f[i][i][L][R]=f[i][i][R][R]= p[i];
        }
        // 用等比数列算出的公式 ：在区间[I,J]，球在k端时，从l端出去的概率  
        For(d,n) {
            For(i,n-d-1) {
                int j=i+d;
                f[i][j][L][L]= (1-p[i]) / ( 1 - p[i] + p[i]*f[i+1][j][L][R] );
                f[i][j][L][R]=  p[i]*f[i+1][j][L][R] / ( 1 - p[i] + p[i]*f[i+1][j][L][R] ); //如果直接用上面的概率减会卡精度
                f[i][j][R][R]= (p[j]) / ( p[j] + (1 - p[j]) * f[i][j-1][R][L] );
                f[i][j][R][L]=  ( (1 - p[j]) * f[i][j-1][R][L]) / ( p[j] + (1 - p[j]) * f[i][j-1][R][L] );
            }
        }
        // 在区间[I,J]，球在k端的概率 
        dp[k][k][R]=1;

        For(d,n-1) {
            For(i,n-d-1) {

                int j=i+d;
                dp[i][j][L] = dp[i+1][j][L] * f[i+1][j][L][L] + dp[i+1][j][R] * f[i+1][j][R][L];
                dp[i][j][R] = dp[i][j-1][L] * f[i][j-1][L][R] + dp[i][j-1][R] * f[i][j-1][R][R];
            }
        }

        double ans = dp[1][n-1][L] + dp[1][n-1][R];
        cout<<setiosflags(ios::fixed)<<setprecision(10);
        cout<<ans<<endl;


    }

    return 0;
}

UVA 1498 Activation

排队，注意各种特判， fi,j 表示队列有i人排第j时“倒霉”的机率，先解方程求出 fi,i 递推。

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define vi vector<int> 
#define pi pair<int,int>
#define SI(a) ((a).size())
typedef long long ll;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
const int MAXN = (2000+10); 
int N,M,K;
double p1,p2,p3,p4;
double f[MAXN][MAXN];  
void prework() {
    For(i,N) {
        double k=p2/(1-p1);
        f[i][i]=0;
        double k2=pow(k,i-1),b2=0;
        Fork(j,2,i) b2=b2*k+(p3*f[i-1][j-1]+p4 * (bool)(j<=K) )/ (1-p1);
        f[i][i]= (k2 * p4 / (1-p1) + b2 ) / ( 1 - k2 * p2 / (1 - p1) ) ;
        f[i][1]=( p2 * f[i][i] + p4 ) / ( 1 - p1 );
        Fork(j,2,i-1) {
            f[i][j] = (p2 * f[i][j-1] + p3 * f[i-1][j-1] + p4 * (j<=K) ) / (1-p1) ;
        }

    }
}

int main()
{
//  freopen("uva1498.in","r",stdin);
//  freopen(".out","w",stdout);

    while(cin>>N>>M>>K) {
        cin>>p1>>p2>>p3>>p4;
        if (p3 == 0 && p4 == 0 || p1 == 1 || p3 == 1 || p2 == 1 ) {
            puts("0.00000"); continue;
        } 
        MEM(f)
        prework();
        cout<<setiosflags(ios::fixed)<<setprecision(5);
        cout<<f[N][M]<<endl; 
    }


    return 0;
}
nike0good
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
高三概率重要复习专题2.doc
09-10
10. **随机变量的期望与概率的乘积**：如`P(X=k) * k`表示随机变量X取值k时的期望贡献。 11. **概率的加法原理**：多个事件中至少有一个发生的概率等于各个事件发生的概率之和减去它们共同发生的概率。 12. **古典...
必知C++算法之概率基本问题
Dxs4396的博客
01-23 4808
概率 概率、期望计算（期望是概率和随机变量乘积的总和） 往往利用古典概率进行计算（组合数学） 概率应用 1.利用随机来改进著名算法（快速排序） 2.随机数的发生器（用给定的随机数发生器构造另外一个） 8只球队，有3个强队，其余都是弱队，随机把他们分成4组比赛，每组两个队，问两强不相遇的概率是多大 1.首先求出8只球队分成4组比赛的方法数 7x5x3x1 = 105种 2.没有两强相遇的方法数 在5...
概率与期望
zhn的博客
05-14 162
1.多重集排列 多重集概念 多重集排列 例题-擅长解密的小红同学
期望&概率dp总结
十分残念的博客
01-25 4857
总算刷完kuangbin期望&amp;概率专题了，下面总结一下心得和题解！1.期望dp期望dp通常逆推，即从结果推向初始状态，也可以用记忆化搜索进行dp；E=Σp1*(E1+X1)+Σp2*(E+X2)其中E为当前状态的期望，E1为下一个状态的期望，p1和X1分别为将当前状态转移到下一个状态的概率和花费，p2和X2分别为保持当前状态的概率和花费。最后化简为E=(Σp1*(E1+X1)+Σp2*X2...
ACM 概率&期望
我的ACM，我的梦！！！
06-26 5744
概率&期望好久没写博客了，最近刷完了概率期望的专题，特地总结一下。概率(1) 基本概率知识 学过概率论课程的话，下面的都是基础了。 ①条件概率：p(A|B) = p(AB) / p(B)。 P(A|B)指B发生的条件下A发生的概率。 ②全概率：p(A)=∑ni=1p(A|Bi)p(A) = \sum_{i=1}^n p(A|B_i)。 全概率公式的关键在于划分样本空间，需要把所有可能情况不
高三数学专题复习概率 教案.doc
10-29
离散型随机变量的期望值表示所有可能结果的平均值，计算公式为E(X) = ΣXi*P(Xi)，其中Xi是每个可能的结果，P(Xi)是对应的结果发生的概率。方差则衡量随机变量分布的离散程度，公式为Var(X) = E[(X - E(X))^2]。 ...
专题资料（2021-2022年）保险学课件例题.doc
10-08
1. 风险衡量：例题1介绍了如何计算概率分布中的标准差和离散系数，这是衡量风险的重要指标。标准差反映数据的波动程度，离散系数则是相对风险度量，通过将标准差除以预期损失，表示风险相对于平均损失的大小。 2. ...
专题资料（2021-2022年）北邮网络学院工程数学阶段作业四.docx
10-08
北邮网络学院的工程数学阶段作业四主要涵盖了概率论与数理统计的相关知识，特别是随机变量的性质、独立性、相关系数、协方差、概率密度函数以及随机变量的期望与方差等内容。以下是这些知识点的详细解释： 1. **...
专题资料（2021-2022年）二叉树定价模型.doc
10-08
在这个案例中，我们还需要确定股价上升的概率p，使得无风险利率r和期望股票价格的折现值匹配。计算出p后，期权的价值即可得出。 两步二叉树模型则进一步增加了计算的精确度。通过将总时间T分为两个较小的时间步长，...
概率和期望的小总结
yiqzq的博客
08-13 364
发现了一个在求期望上面很有用的公式 基本的表示形式是这样子的 Ex=p+(1−p)∗(Ex+1)E_x=p+(1-p)*(E_x+1)Ex​=p+(1−p)∗(Ex​+1), 其中ExE_xEx​表示数学期望,ppp表示成功的概率,所以这个式子是用来求某一个事件成功的期望次数.把它理解成几何分布的形式的话,那么ExE_xEx​的值就等于1p\frac{1}{p}p1​,与实际情况相符。 下面是几个...
概率&期望
wangqianqianya的博客
07-13 440
概率DP f[i]=∑p[i→j]*f[j]+w[i→j] 其中p[]表示转移的概率，w[]表示转移对答案的贡献。一般来说，初始状态确定时可用顺推，终止状态确定时可用逆推。 如： 涂格子1 n个格子，每次随机涂一个，求涂满m个格子的期望次数。 如概述所说，因为最终状态确定，使用逆推。 设计状态f[i]表示涂了i 个格子，到涂满m 个格子还要涂的期望次数。初始状态是f[m]=0 。转移时...
UVA - 11637 Garbage Remembering Exam （组合+概率）
那些年
09-03 1143
Little Tim is now a graduate,and is thinking about higher studies. However, he first needs to appear in anexam whose preparation alone includes memorizing the meanings of over 3500words! After going
概率期望
wangtao的博客
08-05 247
数学期望 E(x)=Σ每一种状态*对应的概率 状态很多,有时候不能完全枚举. 大多数题都是找公式,规律,动态规划计算,数学期望的dp一般要逆推,处理好边界 一般的dp是到dp[x]状态有多少,答案为dp[n] 而数学期望的dp[x]已经在状态x,距离n还差多少,答案为dp[0]. 学习:https://www.cnblogs.com/kuangbin/archive/2012/10/02/...
专题·数学概率与期望【including 条件概率，贝叶斯定理， 全概率公式，数学期望， 绿豆蛙的归宿
樱狸✨愿携清风归来处
04-02 3300
初见安~~~又开启数论的探索啦~~：） 一。概率 1.基本定义 在概率论中，我们把一个随机事件的一个可能结果称为其样本点，其所有样本点构成的集合称之为样本空间。（注意，随机事件并不一定只有一种可能结果）在样本空间中，我们称事件所包含的子集为随机事件。 概率的定义就很简单了，我们也都知道样本空间中的任意随机事件的概率不会超过1不会小于0. 就比如我们抛硬币连续扔三次（不考虑侧面稳落地），有...
概率与期望的学习
FSYo 的博客
04-05 631
目录 一些公式 P1291 [SHOI2002]百事世界杯之旅 [WOJ3010] 骰子 [WOJ3007] Dumb Bones [WOJ2276]挑战 [WOJ3806]收集邮票 [WOJ3083] 收集宝石 [WOJ3080]迷宫 [WOJ2658]单机游戏 [WOJ1111] 钉子和小球 简单游戏 P3211 [HNOI2011]XOR和路径 P2473 [SCO...
浅谈概率与数学期望
beautiful_CXW的博客
10-12 7280
由于鄙人水平有限，故只能给出一些基本的概率入门知识，以及做一些解释，并对一些基础题目诸如三门问题及醉鬼问题的推理。 浅谈概率 先由一道十分经典之题目承题： 三门问题 即有三扇门，有一扇门后由一辆车，另外两扇门分别有一只羊，此时你打开哪扇门，门里的东西就是你的（当然车的价值比羊大且你并不知道每扇门的后面有什么）。此时你已选择了一扇门，但是有知道门后面分别是什么的人给你打开了一扇门后面是羊的门，并给你...
概率期望小结
Destiny
03-19 1374
从一开始就觉得期望和概率很难，现在也是。应该是因为这个玩意儿比较抽象吧，来几道简单例题。 T1 在陆地花费的时间是固定的，那么它对期望的贡献当然也是本身啦。 考虑在通过河的时候，最好的情况是船正好向你走来，那么你花费的时间是Lv\frac{L}{v}vL​，最坏的呢？此时船刚向对面走去，那么花费的时间为3Lv\frac{3L}{v}v3L​。我们求期望，一般有两种方法：1、直接定义期望，递推求...
UVA 11637 - Garbage Remembering Exam(组合概率)
Remilia's
08-06 1625
UVA 11637 - Garbage Remembering Exam 题目链接 题意：大概意思是，有n个单词，分别打乱放在一个环形的，一个非环形里面，环形的两个单词距离为顺时针逆时针的最小值，非环形的就是位置的差的绝对值，如果有一对单词，在两个里面的距离都是不大于k，那么这单词为无效单词，问平均会出现多少个无效单词 思路：组合概率，假设在非环形形成了一个随机序列，那么我们给它
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值