设为首页 加入收藏

TOP

hdu 1204 糖果大战 (Markov Chains求期望)
2015-11-21 01:03:53 来源: 作者: 【 】 浏览:1
Tags:hdu 1204 糖果 大战 Markov Chains 期望

糖果大战

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 2336 Accepted Submission(s): 797


Problem Description 生日Party结束的那天晚上,剩下了一些糖果,Gandon想把所有的都统统拿走,Speakless于是说:“可以是可以,不过我们来玩24点,你不是已经拿到了一些糖果了吗?这样,如果谁赢一局,就拿走对方一颗糖,直到拿完对方所有的糖为止。”如果谁能算出来而对方算不出来,谁就赢,但是如果双方都能算出或者都不能,就算平局,不会有任何糖果的得失。
Speakless是个喜欢提前想问题的人,既然他发起了这场糖果大战,就自然很想赢啦(不然可就要精光了-_-)。现在他需要你的帮忙,给你他每局赢的概率和Gardon每局赢的概率,请你给出他可能获得这场大战胜利的概率。

Input 每行有四个数,Speakless手上的糖果数N、Gardon手上的糖果数M(0<=N,M<=50)、一局Speakless能解答出来的概率p、一个问题Gardon能解答出来的概率q(0<=p,q<=1)。

Output 每行一个数,表示Speakless能赢的概率(用百分比计算,保留到小数点后2位)。

Sample Input
50 50 0.5 0.5
10 10 0.51 0.5
50 50 0.51 0.5

Sample Output
0.50
0.60
0.88


用f[i]表示Speakless有i颗糖果的时赢的概率,则f[i]=p(1-q)*f[i+1]+q(1-p)*f[i-1]+(1-p*(1-q)-q*(1-p))*f[i]; 这个方程表示当前状态赢的期望等于当前状态能到达的所有状态赢的期望的和。化简得, (f[i+1]-f[i-1])*p*(1-q)=f[i]-f[i-1]*(q*(1-p));
令k=(q*(1-p))/(p*(1-q));
设f[i+1]-f[i]= (f[i]-f[i-1])*k;

已知 f[0]=0;f[n+m]=1; (必输和必赢)

f[1]-f[0]=f[1]; f[2]-f[1]=(f[1])*k; f[3]-f[2]=(f[2]-f[1])*k=f[1]*k^2; ... f[n]-f[n-1]=(f[n-1]-f[n-1])*k=f[1]*k^(n-1); 对数列求和,得 f[n]=(1-k^n)/(1-k); f[n+m]=(1-k^(n+m))/(1-k)=1;
则f[n]=f[n]/f[n+m]=(1-k^n)/(1-k^(n+m));

#include
  
   
#include
   
     #include
    
      #include
     
       #include
      
        #include
       
         using namespace std; #define ll long long #define N 25 const int M=256; const double eps=1e-4; int main() { int n,m,n1,n2,i; double k,s,p,q; while(~scanf("%d%d%lf%lf",&n,&m,&p,&q)) { if(m==0) s=1; else if(n==0||p==0||q==1) s=0; else { k=q*(1-p)/p/(1-q); if(fabs(k-1.0)
        
         









】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
分享到: 
上一篇ZOJ 3872 Beauty of Array(数学) 下一篇NYOJ 737 石子合并(一) (区间D..

评论

帐  号: 密码: (新用户注册)
验 证 码:
表  情:
内  容: