得到2颗星的期望次数可以解方程算但得到1颗星的算不了。
因为每关都要通过,我们按Yi从小到大处理,f[i][j],pi[i][j]表示玩了i关,获得了j颗星的期望次数,概率
如果(n-i)*2=m-j,那接下来每关都要玩2星,也只会在这些关弄2星因为他们的Yi最大
否则我们只要求通过就行了
通过的期望次数是能算的,通过了后获得1/2颗星的期望是Xi:Yi,就可以dp了
code:
#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn = 2010;
const int maxm = 4010;
int n,m;
struct node{int x,y;}a[maxn];
inline bool cmp(const node x,const node y){return x.y<y.y;}
double f[maxn][maxm],pi[maxn][maxm];
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) scanf("%d%d",&a[i].x,&a[i].y);
sort(a+1,a+n+1,cmp);
f[0][0]=0.0; pi[0][0]=1.0;
for(int i=1;i<=n;i++)
{
double x=(double)a[i].x/1000.0,y=(double)a[i].y/1000.0;
double p=1.0/(x+y),p2=1.0/y;
double px=x/(x+y),py=y/(x+y);
for(int j=i-1;j<=2*n;j++)
{
double lp=pi[i-1][j],tmp=f[i-1][j];
if(2*(n-i+1)==m-j) f[i][j+2]+=tmp+p2*lp,pi[i][j+2]+=lp;
else
{
f[i][j+1]+=tmp*px+p*px*lp,pi[i][j+1]+=lp*px;
f[i][j+2]+=tmp*py+p*py*lp,pi[i][j+2]+=lp*py;
}
}
}
double ans=0.0;
for(int i=m;i<=2*n;i++) ans+=f[n][i];
printf("%.8lf\n",ans);
return 0;
}