糖果大战 - HDU 1204 非公式递推版本

糖果大战

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

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
  
  
  
  


  
  

题意：好像不用说大家也明白……【毕竟是中文的题目。

思路：在网上查题解看到的基本上都是推公式的，但这道题其实递推也能做。首先设A赢的概率为k，那么f[n]=f[n-1]*(1-k)+f[n+1]*k（不懂的考虑下n m 都是1的情况），但是这个公式里面基本上有两个是不知道的，所以设f[1]=1,因为f[0]=0,f[n+m]=1，当推到f[n+m]时，这个结果应该大于1，那么这个结果就是因为f[1]多了f[n+m]倍所造成的，那么f[1]就应该等于1/f[n+m]，同理，剩下的所有的都应该除以f[n+m]。另外考虑几个特殊情况。

AC代码如下：

#include<cstdio>
#include<cstring>
using namespace std;
double f[1010];
int main()
{ int n,m,i,j,len;
  double p,q,k,num;
  while(~scanf("%d%d%lf%lf",&n,&m,&p,&q))
  { if(n==0)
    { printf("0.00\n");continue;}
    if(m==0)
    { printf("1.00\n");continue;}
    if(p==1 || q==0)
    { printf("1.00\n");continue;}
    if(q==1 || p==0)
    { printf("0.00\n");continue;}
    k=p*(1-q)/(q*(1-p)+p*(1-q));
    len=n+m;
    f[0]=0;
    f[1]=1;
    f[len]=1;
    for(i=2;i<=len;i++)
    f[i]=(f[i-1]-f[i-2]*(1-k))/k;
    num=f[len];
    printf("%.2f\n",f[n]/num);
  }
}