分析:
这个题很坑,注意乘除要同时进行,否则精度不够。
三分查找写挫了,贴一个非三分代码。
代码:
#include <bits/stdc++.h>
using namespace std;
long long n,p;
int main()
{
cin>>n>>p;
int x=0,i,t1,t2;
while((x+1)*(n+x+2)<(x+2)*(n+x+2-p)) x++;
long double tmp=(long double)p*(long double)(x+1);
t1=n-p+2;
t2=n+x-p+2;
for(i=0;;i++)
{
if(t1+i<=n) tmp*=(long double)(t1+i);
if(t2+i<=n+x+1) tmp/=(long double)(t2+i);
if((t1+i>n)&&(t2+i>n+x+1)) break;
}
printf("%.12Lf",tmp);
return 0;
}