HDU 4465 HDOJ Candy 2012ACM亚洲赛成都赛区B题

http://acm.hdu.edu.cn/showproblem.php?pid=4465

本题题解非博主原创。

是队里的神牛YCL在比赛的时候灵光一闪想到的，目前依旧只能无尽的膜拜OOOOrz

首先公式是很好推的

　　for(int i=0;i<=n;i++) ans+=C(n+i,n)*(p^(n+1) * (1-p)^i + (1-p)^(n+1) * p^i)

目前依然不知道官方是怎样做到不溢出。。基本网上看到的题解都是在想办法解决中间过程溢出的问题，这里提供下队里YCL神牛比赛时灵光一闪的做法。。

看代码应该就能看懂了。

 1 #include<cstdio>

 2 #include<cmath>

 3 #include<cstdlib>

 4 #include<algorithm>

 5 using namespace std;

 6 int main()

 7 {

 8     int n,C=0;

 9     double p,c,p1,p2,s1,s2,ans;

10     while(~scanf("%d%lf",&n,&p)){

11         ans=c=0.0;

12         p1=log(p);

13         p2=log(1-p);

14         s1=(n+1)*p1;

15         s2=(n+1)*p2;

16         for(int i=0;i<n;i++){

17             if(c+s1 > -30 || c+s2 > -30)ans+=(exp(c+s1)+exp(c+s2))*(n-i);

18             c+=log(n+i+1)-log(i+1);

19             s1+=p2;

20             s2+=p1;

21         }

22         printf("Case %d: %lf\n",++C,ans);

23     }

24     return 0;

25 }