HDU 4465 - Candy（概率与数学优化）

2012成都Regional的B题，花了一个小时推出了式子，但是搞了好久发现都控制不了精度，后来突然想到组合数可以用log优化，改了之后就AC了

比较水的概率题

 

#include <stdio.h>
#include <math.h>
#define maxn 200005

double f[2*maxn];
double c(int a,int b)
{
    return f[a]-f[a-b]-f[b];
}
int main()
{
    f[0]=0;
    for(int i=1;i<2*maxn;i++)
    {
        f[i]=f[i-1]+log(1.0*i);
    }
    double ans;
    double p;
    int n,cas=0;
    while(scanf("%d%lf",&n,&p)==2)
    {
        ans=0;
        double u=log(p),v=log(1-p);
        for(int i=1;i<=n;i++)
        {
            double k1=c(2*n-i,n-i)+(n-i)*u+(n+1)*v;
            double k2=c(2*n-i,n-i)+(n-i)*v+(n+1)*u;
            ans+=i*(exp(k1)+exp(k2));
        }
        printf("Case %d: %lf\n",++cas,ans);
    }
    return 0;
}