UVa 1639 (期望) Candy

题意：

两个盒子里各有n颗糖，每天有p的概率从第一个盒子里取一颗糖，1-p的概率从第二个盒子里去一颗糖。直到某一天打开某个盒子忽然发现没糖了，求另一个盒子里剩余糖果数的期望。

分析：

紫书上面已经分析的很清楚了，而且也给出了解决精度损失问题的方法，就是先取对数然后再乘幂。

 #include <cstdio>
 #include <cmath>

 const int maxn =  + ;
 long double logF[maxn *  + ];

 long double logC(int n, int m)
 {
     return logF[n] - logF[m] - logF[n-m];
 }

 int main()
 {
     //freopen("in.txt", "r", stdin);
     for(int i = ; i <= maxn*; ++i)
         logF[i] = logF[i-] + log(i);

     int n, kase = ;
     double p;
     while(scanf("%d%lf", &n, &p) == )
     {
         double ans = 0.0;
         for(int i = ; i <= n; ++i)
         {
             long double c = logC(*n-i, n);
             long double v1 = c + (n+)*log(p) + (n-i)*log(-p);
             long double v2 = c + (n+)*log(-p) + (n-i)*log(p);
             ans += i * (exp(v1) + exp(v2));
         }
         printf("Case %d: %.6f\n", ++kase, ans);
     }

     return ;
 }

代码君