这道题就是要求你一直抽卡,直到抽到粉卡(然后歪了)。再抽卡的过程中,每一种卡被抽到的概率还会有所变化,具体变化见原题。
我们发现, v v v 的最小值为 0.1 0.1 0.1,说明整个抽卡的流程中,抽卡次数最多也只有 15 15 15 次左右,完全可以爆搜。
我们可以在将以下五个参数传入搜索中:三种卡被抽中的概率 a , b , c a,b,c a,b,c,抽到这一步时一共抽了几次卡 s t e p step step 以及抽到这一步时的概率 s u m sum sum,每一次将 c × s t e p × s u m c\times step\times sum c×step×sum 累加给答案,然后按照题意更新每一个参数即可。注意精度。
代码:
#include<cmath>
#include<cstdio>
#define int long long
using namespace std;
double a,b,c,v,ans;
const double lep=1e-10;
int T;
void dfs(double x,double y,double z,int s,double sum){
ans+=sum*z*s;
double xx=x,yy=y,zz=z,tot=sum;
sum*=x;//假设抽到了第一种卡
if(x>=v){//按照题意更新概率即可
if(y>=lep){
y+=v/2;
}
if(y>=lep){
z+=v/2;
}else{
z+=v;
}
x-=v;
}else if(x>=lep){
if(y>=lep){
y+=x/2;
}
if(y>=lep){
z+=x/2;
}else{
z+=x;
}
x=0;
}
/*
判断:如果“抽到第一张卡”这个操作不合法,也就是抽到第一张卡的概率降为了0,那么,抽到第三张卡的概率是绝对不会变的。
而如果抽到第三张卡的概率变了,说明抽到第一张卡的概率不为0,也就是这个说:这个操作是合法的,可以进行下一步搜索。
*/
if(zz!=z){
dfs(x,y,z,s+1,sum);
}
x=xx,y=yy,z=zz,sum=tot;
//同理,假设抽到了第二种卡
sum*=y;
if(y>=v){
if(x>=lep){
x+=v/2;
}
if(x>=lep){
z+=v/2;
}else{
z+=v;
}
y-=v;
}else if(y>=lep){
if(x>=lep){
x+=y/2;
}
if(x>=lep){
z+=y/2;
}else{
z+=y;
}
y=0;
}
if(zz!=z){
dfs(x,y,z,s+1,sum);
}
}
signed main(){
scanf("%d",&T);
while(T--){
scanf("%lf%lf%lf%lf",&a,&b,&c,&v);
dfs(a,b,c,1,1);
printf("%.12lf\n",ans);
ans=0;
}
return 0;
}