BZOJ4832: [Lydsy2017年4月月赛]抵制克苏恩

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

int t,a,b,c,k;

double f[55][8][8][8]; 

int main(){
    scanf("%d",&t);
    f[i][j][p][q]//攻击了i次，一滴血的j个，二滴血的p个，三滴血的q个 
    while(t--){
        scanf("%d%d%d%d",&k,&a,&b,&c);
        memset(f,0,sizeof(f));
        f[0][a][b][c]=1.;
        double ans=0.;
        for(int i=0;i<k;i++){
            for(int j=0;j<=7;j++){
                for(int p=0;p<=7;p++){
                    for(int q=0;q<=7;q++){
                        if(j+p+q>7)break;
                        f[i+1][j][p][q]+=f[i][j][p][q]/(j+p+q+1);
                        ans+=f[i][j][p][q]/(j+p+q+1);
                        if(j+p+q==7){
                            if(j)f[i+1][j-1][p][q]+=f[i][j][p][q]*j/(j+p+q+1);  
                            if(p)f[i+1][j+1][p-1][q]+=f[i][j][p][q]*p/(j+p+q+1);//攻击二滴血的，二滴血的-1，一滴血的+1 
                            if(q)f[i+1][j][p+1][q-1]+=f[i][j][p][q]*q/(j+p+q+1);
                        }else{
                            if(j)f[i+1][j-1][p][q]+=f[i][j][p][q]*j/(j+p+q+1);
                            if(p)f[i+1][j+1][p-1][q+1]+=f[i][j][p][q]*p/(j+p+q+1);
                            if(q)f[i+1][j][p+1][q]+=f[i][j][p][q]*q/(j+p+q+1);//攻击三滴血的-1，二滴血+1，三滴血+1 
                        }
                    }
                }
            }
        }
        printf("%.2f\n",ans);
    }
    return 0;
}