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

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

比较水的概率题

#include <stdio.h>#include <math.h>#define maxn 200005double 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;}