UVA - 1639 Candy （概率，精度）

X表示剩下的糖数量，如果最后打开的是p对应的盒子。划分：Xi表示剩下i个糖，最后一次选的概率为p，

前面的服从二项分布。根据全概率公式和期望的线性性，求和就好了。

精度处理要小心，n很大，组合数会很大，p的部分很小，要取对数，而且中间计算精度也要用long double才够。

组合数的对数预处理一下或者递推一下就好了。

/*********************************************************
*      --------------Tyrannosaurus---------              *
*   author AbyssalFish                                   *
**********************************************************/
#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef long double ld;

const int maxn2 = 4e5+;
ld sum_fac_ln[maxn2];

//inline double ln_(int x) { return sum_fac_ln[x]-sum_fac_ln[x-1]; }
inline ld ln_C(int m,int n)
{
    return sum_fac_ln[m]-sum_fac_ln[n]-sum_fac_ln[m-n];
}

const ld eps = 1e-;

double solve(int n, ld p)
{
    ld q = -p;
    if(p < eps || q < eps) return n;
    ld E = ;
    int _2n = n<<;
    ld ln_p = log(p), ln_q = log(q);
    ld nln_p = (n+)*ln_p;
    ld nln_q = (n+)*ln_q;
    for(int x = ; x <= n; x++){
        ld lnC = ln_C(_2n-x,n);
        E += x*( exp( lnC + nln_p + (n-x)*ln_q ) + exp(lnC + nln_q + (n-x)*ln_p ) );
    }
    return E;
}

//#define LOCAL
int main()
{
#ifdef LOCAL
    freopen("in.txt","r",stdin);
#endif
    for(int i = ; i < maxn2; i++) sum_fac_ln[i] = log((ld)i);
    for(int i = ; i < maxn2; i++) sum_fac_ln[i] += sum_fac_ln[i-];
    int kas = , n;
    double p;

    while(~scanf("%d%lf",&n,&p)){
        printf("Case %d: %.6lf\n", ++kas, solve(n,p));
    }
    return ;
}