tyvj2002 扑克牌

背景 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。

动态规划。思维难度不大，但是细节比较多，容易写错。
dp[i][j][k][l][p][q]表示四种牌分别翻出i,j,k,l张，大小王的状态分别为p,q的期望次数。
记忆化搜索即可，对于翻出四种花色的情况累加，翻出大小王的情况考虑当成的颜色取最小值累加。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<iomanip>
using namespace std;
double dp[17][17][17][17][7][7];
bool vis[17][17][17][17][7][7];
int a,b,c,d;
double min(double x,double y)
{
    return x<y?x:y;
}
double dfs(int aa,int bb,int cc,int dd,int p,int q)
{
    int i,j,k,x,y,z,sum;
    double ans,temp;
    if (vis[aa][bb][cc][dd][p][q]) return dp[aa][bb][cc][dd][p][q];
    vis[aa][bb][cc][dd][p][q]=1;
    if (aa+(p==1)+(q==1)>=a&&bb+(p==2)+(q==2)>=b
      &&cc+(p==3)+(q==3)>=c&&dd+(p==4)+(q==4)>=d) return 0;
    sum=54-(aa+bb+cc+dd+(p>0)+(q>0));
    ans=1;
    if (aa<13) ans+=dfs(aa+1,bb,cc,dd,p,q)*(13-aa)/sum;
    if (bb<13) ans+=dfs(aa,bb+1,cc,dd,p,q)*(13-bb)/sum;
    if (cc<13) ans+=dfs(aa,bb,cc+1,dd,p,q)*(13-cc)/sum;
    if (dd<13) ans+=dfs(aa,bb,cc,dd+1,p,q)*(13-dd)/sum;
    if (!p)
    {
        temp=10000;
        /*if (aa<13)*/ temp=min(temp,dfs(aa,bb,cc,dd,1,q)/sum);
        /*if (bb<13)*/ temp=min(temp,dfs(aa,bb,cc,dd,2,q)/sum);
        /*if (cc<13)*/ temp=min(temp,dfs(aa,bb,cc,dd,3,q)/sum);
        /*if (dd<13)*/ temp=min(temp,dfs(aa,bb,cc,dd,4,q)/sum);
        ans+=temp;
    }
    if (!q)
    {
        temp=10000;
        /*if (aa<13)*/ temp=min(temp,dfs(aa,bb,cc,dd,p,1)/sum);
        /*if (bb<13)*/ temp=min(temp,dfs(aa,bb,cc,dd,p,2)/sum);
        /*if (cc<13)*/ temp=min(temp,dfs(aa,bb,cc,dd,p,3)/sum);
        /*if (dd<13)*/ temp=min(temp,dfs(aa,bb,cc,dd,p,4)/sum);
        ans+=temp;
    }
    dp[aa][bb][cc][dd][p][q]=ans;
    return ans;
}
int main()
{
    int i,j,k,m,n,p,q,x,y,z,tot;
    scanf("%d%d%d%d",&a,&b,&c,&d);
    tot=0;
    if (a>13) tot+=a-13;
    if (b>13) tot+=b-13;
    if (c>13) tot+=c-13;
    if (d>13) tot+=d-13;
    if (tot>2)
      printf("-1.000\n");
    else
      cout<<setiosflags(ios::fixed)<<setprecision(3)<<dfs(0,0,0,0,0,0)<<endl;
}