点击打开链接
题意:自己读
公式
∑i=0n(p(n+1)∗(1−p)(n+1)∗C(n)(2n−1)∗i+(1−p)(n+1)∗p(n+1)∗C(n)(2n−1)∗i)
直接计算的话,会溢出 。
取一个ln之后 就可以计算了。
ln(Cmn)=∑ni=m+1(lni)−∑(n−m)i=2(lni)
Cmn=exp(ln(Cmn))
#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;
}