TYVJ P2002 扑克牌

背景

Admin生日那天，Rainbow来找Admin玩扑克牌……
玩着玩着Rainbow觉得太没意思了,于是决定给Admin一个考验~~~

描述

Rainbow把一副扑克牌(54张)随机洗开，倒扣着放成一摞。然后Admin从上往下依次翻开每张牌，每翻开一张黑桃、红桃、梅花或者方块，就把它放到对应花色的堆里去。
Rainbow想问问Admin，得到A张黑桃、B张红桃、C张梅花、D张方块需要翻开的牌的张数的期望值E是多少？
特殊地，如果翻开的牌是大王或者小王，Admin将会把它作为某种花色的牌放入对应堆中,使得放入之后E的值尽可能小。
由于Admin和Rainbow还在玩扑克，所以这个程序就交给你来写了~

输入格式

输入仅由一行，包含四个用空格隔开的整数，A,B,C,D。

输出格式

输出需要翻开的牌数的期望值E，四舍五入保留3位小数。
如果不可能达到输入的状态，输出-1.000。

测试样例1

输入

样例输入1
1 2 3 4

样例输入2
15 15 15 15

输出

样例输出1
16.393

样例输出2
-1.000

备注

对于100%的数据，0<=A,B,C,D<=15

lydrainbowcat - "Admin生日"杯NOIP模拟赛 第三题
Blog: www.lydrainbowcat.tk  Email: lyd@tyvj.cn

决策挺复杂的。

现在不想写解析，日后补。

/*by SilverN*/
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
const int mxn=15;
double f[mxn][mxn][mxn][mxn][5][5];
bool vis[mxn][mxn][mxn][mxn][5][5];
int x1,x2,x3,x4;
double dp(int a,int b,int c,int d,int e,int g){
    if(vis[a][b][c][d][e][g])return f[a][b][c][d][e][g];
    int i;

    if(((a+(e==0)+(g==0))>=x1) && ((b+(e==1)+(g==1))>=x2) && ((c+(e==2)+(g==2))>=x3) && ((d+(e==3)+(g==3))>=x4))
        return f[a][b][c][d][e][g]=0;
    int sum=a+b+c+d+(e!=4)+(g!=4);
    double tmp=1;
    //普通 
    if(a<13)tmp+=dp(a+1,b,c,d,e,g)*(13-a)/(54-sum);
    if(b<13)tmp+=dp(a,b+1,c,d,e,g)*(13-b)/(54-sum);
    if(c<13)tmp+=dp(a,b,c+1,d,e,g)*(13-c)/(54-sum);
    if(d<13)tmp+=dp(a,b,c,d+1,e,g)*(13-d)/(54-sum);
    //王
    double es;
    if(e==4){
        es=100;
        for(i=0;i<=3;i++) es=min(dp(a,b,c,d,i,g)/(54-sum),es);
        tmp+=es;
    }
    if(g==4){
        es=100;
        for(i=0;i<=3;i++) es=min(dp(a,b,c,d,e,i)/(54-sum),es);
        tmp+=es;
    }
    
    vis[a][b][c][d][e][g]=1;
    return f[a][b][c][d][e][g]=tmp;
}
int main(){
    memset(vis,0,sizeof(vis));
    int a,b,c,d;
    scanf("%d%d%d%d",&x1,&x2,&x3,&x4);
    double ans=dp(0,0,0,0,4,4);
    if(ans>54)ans=-1;
    printf("%.3lf
",ans);
    return 0;
}