题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4465
这个题目公式很好推,但是得出正确结果还有很多问题要解决
首先说公式一共是两盒,我们只考虑一盒,另外一盒同理
从剩下0个到剩下n个的概率都能算出来,所以利用组合数公式可以直接得到
num=(n-i)*(C(n+i,i)*p^i*(1-p)^n*(1-p));
其实前面的都好理解,但是最后一个1-p怎么来的,其实很简单,就是最后等
这盒没了之后还要选到这盒的意思
懂了上面的公式之后就是求了
直接求组合受不了肯定超范围或者精度损失的一塌糊涂
对于组合数,那就取对数呗
C(n,m)=log(n!/(m!*(n-m)!)=logn!-logm!-log((n-m))
现在就是怎么求logn!的问题了
logn!=log1+log2+log3+......logn;
滴推一下就好了!
#include <iostream>
#include <math.h>
#include <string.h>
#include <algorithm>
#include <stdio.h>
using namespace std;
#define maxn 400005
int n;
double p,q;
double f[maxn];
int init(){
f[0]=0;
for(int i=1;i<maxn;i++)
f[i]=f[i-1]+log(i*1.0);
return 0;
}
double C(int n,int m){
return f[n]-f[m]-f[n-m];
}
int main(){
int i,j,k=0;
double ans;
init();
while(scanf("%d%lf",&n,&p)!=EOF){
ans=0;
q=1-p;
for(i=0;i<=n;i++)
ans+=(n-i)*(exp(C(n+i,i)+(n+1)*log(q*1.0)+i*log(p*1.0))+exp(C(n+i,i)+(n+1)*log(p*1.0)+i*log(q*1.0)));
printf("Case %d: %lf\n",++k,ans);
}
return 0;
}