TYVJ 2002 扑克牌

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

---

【题目分析】
lyd留的一道作业题。莫名其妙被秀了一脸（划去），被难了半天。其实是一道数学期望DP。由于大王和小王的存在，使问题复杂了。找不到好的递推的方法。只能使用记忆化搜索。然后分别考虑普通扑克和大小王。a,b,c,d,e,f分别表示四种牌，以及大王小王的状态。0，1，2，3代表四种花色。4表示还没有出现。然后每种状态初值为1（这种状态本身就需要翻一次）.然后加上各种情况的期望值。如果期望值大于54，就是-1.
感觉tyvj好坑啊，用printf输出实数，无论怎么都是-0.000，不得已用了cout。各种复杂的流输出。看起来很难，代码不长。状态转移详见代码。

---

【代码】

#include <cstdio>
#include <cstring>
#include <cmath>
#include<iomanip>
#include <algorithm>
#include <iostream>
using namespace std;
double dp[15][15][15][15][5][5];
bool vis[15][15][15][15][5][5];
int t[4];
double dfs(int a,int b,int c,int d,int e,int f)
{
    if (vis[a][b][c][d][e][f]) return dp[a][b][c][d][e][f];
    if (((a+(e==0)+(f==0))>=t[0])&&((b+(e==1)+(f==1))>=t[1])&&((c+(e==2)+(f==2))>=t[2])&&((d+(e==3)+(f==3))>=t[3])) return dp[a][b][c][d][e][f]=0.0;
    int sum=a+b+c+d+(e!=4)+(f!=4);
    double F=1;
    if (a<13) F+=dfs(a+1,b,c,d,e,f)*(13-a)/(54-sum);
    if (b<13) F+=dfs(a,b+1,c,d,e,f)*(13-b)/(54-sum);
    if (c<13) F+=dfs(a,b,c+1,d,e,f)*(13-c)/(54-sum);
    if (d<13) F+=dfs(a,b,c,d+1,e,f)*(13-d)/(54-sum);
    double md=100; if (e==4) {for (int i=0;i<4;++i) md=min(dfs(a,b,c,d,i,f)/(54-sum),md); F+=md;}
           md=100; if (f==4) {for (int i=0;i<4;++i) md=min(dfs(a,b,c,d,e,i)/(54-sum),md); F+=md;}
    vis[a][b][c][d][e][f]=true;
    return dp[a][b][c][d][e][f]=F;
}
int main()
{
    memset(vis,false,sizeof vis);
    int a,b,c,d;
    scanf("%d%d%d%d",&t[0],&t[1],&t[2],&t[3]);
    a=t[0];b=t[1];c=t[2];d=t[3];
    double ans=dfs(0,0,0,0,4,4);
    if (ans>54.0) printf("-1.000\n");
    else cout<<setiosflags(ios::fixed)<<setprecision(3)<<ans;
}