HDU Candy(数学 组合)

题目链接：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;
}