【基础数学】 HDU 4465 Candy Log搞

点击打开链接
题意：自己读
公式

$$\sum_{i=0}^n \left(p^\left(n+1\right)*(1-p)^\left(n+1\right)*C{^\left(n\right)_\left(2n-1\right)}*i+（1-p）^\left(n+1\right)*p^\left(n+1\right)*C{^\left(n\right)_\left(2n-1\right)}*i\right)$$
直接计算的话，会溢出 。
取一个ln之后 就可以计算了。
$ln(C{^m_n）}=\sum_{i=m+1}^n(lni)-\sum_{i=2}^\left(n-m\right)(lni)$
$C{^m_n}=exp(ln(C{^m_n）})$

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;
#include <queue>
#include <stack>
#include <vector>
#include <bitset>
#include <deque>
#include <set>
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <map>
typedef long long LL;
typedef unsigned long long ULL;
const int INF = 1<<29;
const LL mod = 1e9+7;
const int MAXN = 200100;//点数的最大值
const int MAXM = 520;//边数的最大值
const double pi = acos(-1);
double In[MAXN];
int main()
{
    In[0]=0;
    for(int i=1;i<MAXN;i++)
    {
        In[i]=In[i-1]+log(i);
    }
    int cas=0,n;
    double p;
    while(scanf("%d%lf",&n,&p)!=EOF)
    {
        double ans=0;
        for(int i=1;i<=n;i++)
        {
            double go1=(n+1)*log(p)+(n-i)*log(1-p)+log(i)+In[2*n-i]-In[n];
            if(n-i>=2)
                go1-=(In[n-i]-In[1]);
            double go2=(n+1)*log(1-p)+(n-i)*log(p)+log(i)+In[2*n-i]-In[n];
            if(n-i>=2)
                go2-=(In[n-i]-In[1]);
            ans+=exp(go1)+exp(go2);
        }
        printf("Case %d: %.6lf\n",++cas,ans);
    }
    return 0;
}