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;}