TYVJ2002 扑克牌

卢克生日那天，汉来找卢克玩扑克牌，玩着玩着汉觉得太没意思了，于是决定给卢克一个考验
汉把一副扑克牌(54张)随机洗匀，倒扣着放成一摞。然后卢克从上往下一次翻开每张牌，每翻开一张黑桃，红桃，梅花或方块，就把它放到对应的花色的堆里去。
汉想问问卢克，得到A张黑桃，B张红桃，C张梅花，D张方块需要翻开的牌的张数的期望值E是多少
特殊的，如果翻开的牌是大小王，那么卢克可以把它作为某种花色的牌放入对应堆中。卢克会采取最优策略，使得放入之后E的值尽可能的小
显然卢克靠原力的感知可以无视这个问题，汉为给卢克丢了一个没意义的问题而懊恼，于是决定把问题丢给还是绝地学徒的你
虽然你并不能善用原力，但是好在你会c++，于是你可以靠程序解决这个问题。
卢克和汉还在玩牌，所以你可以在计算出结果后直接将答案汇报给他们
0<=A,B,C,D<=15 四舍五入保留三位小数输出

f[a,b,c,d,x,y]表示黑桃取了a张，红桃取了b张，梅花取了c张，方块取了d张，小王状态为x，大王状态为y时的期望值
具体来说，x = 4表示没有用过小王，x = 0~3表示用过小王，且把小王作为对应的花色(0代表黑桃，1代表红桃，2代表梅花，3代表方块)
可以知道的是，随着翻开的牌数的增加，得到某种花色的概率是在变化的。
当前已经翻开的牌总数sum = (a + b + c + d + (x == 4) + (y == 4)）。到目前为止，还剩下54-sum张牌，其中13-a张黑桃，13-b张红桃，13-c张梅花，13-d张方块
以黑桃为例：
翻开一张黑桃的概率为(13 - a)/(54 - sum)。还需要翻开的牌的期望张数为f[a+1,b,c,d,x,y]。对于红桃，梅花，方块，情况类似
特别地，当x = 4时，有1/(54 - sum)的概率翻开小王。根据题意，应选择把小王看做某种花色，是期望值尽量小，即min{f[a,b,c,d,x',y]}(0<=x'<=3).对于大王，情况类似
对于大王，情况类似

方程显然且省略，建议直接看代码

初值：若已经翻开的牌数达到了题目要求的数量，则期望值位0.例如已经翻开的黑桃张数为a+(x==0)+(y==0)，其余同理
目标f[0,0,0,0,4,4]

在数学期望递推，数学期望Dp中，通常把终止状态作为初值，把起始状作为目标状态，倒着进行计算。
原因：以本题为例：
根据数学期望的定义，若我们正着计算，则还需求出从起始状态到达每个终止状态的概率，与F值相乘求和才能得到答案，增加了难度，且易出错
而倒着计算，因为起始状态F[0,0,0,0,4,4]唯一，所以它的概率一定为1，最后直接输出f[0,0,0,0,4,4]即为所求

 #include<bits/stdc++.h>
 using namespace std;
 int A, B, C, D;
 double f[][][][][][];
 bool vis[][][][][][];

 inline int read() {
     int x = , y = ;
     char ch = getchar();
     while(!isdigit(ch)) {
         if(ch == '-') y = -;
         ch = getchar();
     }
     while(isdigit(ch)) {
         x = (x << ) + (x << ) + ch - '';
         ch = getchar();
     }
     return x * y;
 }

 double Dp_to_make_ans(int a, int b, int c, int d, int p, int q) {
     if(vis[a][b][c][d][p][q]) return f[a][b][c][d][p][q];
     vis[a][b][c][d][p][q] = ;
     int h = a, j = b, k = c, l = d;
     if(p == ) h++; if(p == ) j++; if(p == ) k++; if(p == ) l++;
     if(q == ) h++; if(q == ) j++; if(q == ) k++; if(q == ) l++;
     if(h >= A && j >= B && k >= C && l >= D) return f[a][b][c][d][p][q] = ;
     int sum = h + j + k + l;
     double f_ans = ;
     if(a < ) f_ans += Dp_to_make_ans(a + , b, c, d, p, q) * ( - a) / ( - sum);
     if(b < ) f_ans += Dp_to_make_ans(a, b + , c, d, p, q) * ( - b) / ( - sum);
     if(c < ) f_ans += Dp_to_make_ans(a, b, c + , d, p, q) * ( - c) / ( - sum);
     if(d < ) f_ans += Dp_to_make_ans(a, b, c, d + , p, q) * ( - d) / ( - sum);
     double shiki;
     if(p == ) {
         shiki = ;
         for(int i = ; i <= ; ++i) shiki = min(shiki, Dp_to_make_ans(a, b, c, d, i, q));
         f_ans = f_ans + shiki / ( - sum);
     }
     if(q == ) {
         shiki = ;
         for(int i = ; i <= ; ++i) shiki = min(shiki, Dp_to_make_ans(a, b, c, d, p, i));
         f_ans = f_ans + shiki / ( - sum);
     }
     return f[a][b][c][d][p][q] = f_ans;
 }

 int main() {
 //    freopen("test.in", "r", stdin);
 //    freopen("test.out", "w", stdout);
     A = read(), B = read(), C = read(), D = read();
     double ans = Dp_to_make_ans(, , , , , );
     if(ans > ) ans = -;
     printf("%0.3lf", ans);
     return ;
 }