【JoyOI 2002】扑克牌

传送门

---

Problem

Rainbow 把一副扑克牌$(54$张$)$随机洗开，倒扣着放成一摞。然后 Admin 从上往下依次翻开每张牌，每翻开一张黑桃、红桃、梅花或者方块，就把它放到对应花色的堆里去。

Rainbow 想问问 Admin，得到 $\mathrm{A}$ 张黑桃、$\mathrm{B}$ 张红桃、$\mathrm{C}$ 张梅花、$\mathrm{D}$ 张方块需要翻开的牌的张数的期望值 $\mathrm{E}$ 是多少？

特殊地，如果翻开的牌是大王或者小王，Admin 将会把它作为某种花色的牌放入对应堆中，使得放入之后 $\mathrm{E}$ 的值尽可能小。

由于 Admin 和 Rainbow 还在玩扑克，所以这个程序就交给你来写了。

数据范围：$0\le \mathrm{A},\mathrm{B},\mathrm{C},\mathrm{D}\le 15$。

---

Solution

期望概率DP简单题。

我们设 $f[a,b,c,d,x,y]$ 表示已经翻了 $a$ 张黑桃、$b$ 张红桃、$c$ 张梅花、$d$ 张方块，并且小王状态为 $x$，大王状态为 $y$ 的最小期望。

小王的状态指，当 $x=4$ 时，没翻到小王；当 $x=0,1,2,3$ 时翻到了小王，且小王作为对应的花色（$0$ 表示黑桃、$1$ 表示红桃、$2$ 表示梅花、$3$ 表示张方块），大王和小王同理。

此时翻开的牌的总数为 $sum=a+b+c+d+[x$ ≠ $4]+[y$ ≠ $4]$。目前还剩 $54-sum$ 张牌，$13-a$ 张黑桃，那再翻一张黑桃的概率为 $\frac{13-a}{54-num}$。翻开这张牌后，还需翻开的期望张数就是 $f[a+1,b,c,d,x,y]$，其他牌同理。

特别的，对于大小王，它们有 $\frac{1}{54-num}$ 的概率被翻到，此时枚举它最为那种花色的牌，使期望最小，即 $mi{n}_{0\le x^{\prime }\le 3}f[a,b,c,d,x^{\prime },y]$。

那么综上，有：$$\begin{gathered} f[a,b,c,d,x,y] \\ =\frac{13-a}{54-num}\times f[a+1,b,c,d,x,y]+\frac{13-b}{54-num}\times f[a,b+1,c,d,x,y] \\ +\frac{13-c}{54-num}\times f[a,b,c+1,d,x,y]+\frac{13-d}{54-num}\times f[a,b,c,d+1,x,y] \\ +\frac{1}{54-num}\times mi{n}_{0\le x^{\prime }\le 3}f[a,b,c,d,x^{\prime },y](ifx=4) \\ +\frac{1}{54-num}\times mi{n}_{0\le y^{\prime }\le 3}f[a,b,c,d,x,y^{\prime }](ify=4) \end{gathered}$$

初始化：把翻开的牌数达到要求的期望赋为 $0$。

目标：$f[0,0,0,0,4,4]$。

用记忆化搜索DP就可以了。

---

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int A,B,C,D;
double f[15][15][15][15][5][5];
bool vis[15][15][15][15][5][5];
double dp(int a,int b,int c,int d,int x,int y){
	if(vis[a][b][c][d][x][y])  return f[a][b][c][d][x][y];
	vis[a][b][c][d][x][y]=true;
	int i=a+(x==0)+(y==0),j=b+(x==1)+(y==1),k=c+(x==2)+(y==2),l=d+(x==3)+(y==3);
	if(i>=A&&j>=B&&k>=C&&l>=D)  return f[a][b][c][d][x][y]=0;
	int sum=i+j+k+l;double ans=1,Min;
	if(a<13)  ans+=dp(a+1,b,c,d,x,y)*(13-a)/(54-sum);
	if(b<13)  ans+=dp(a,b+1,c,d,x,y)*(13-b)/(54-sum);
	if(c<13)  ans+=dp(a,b,c+1,d,x,y)*(13-c)/(54-sum);
	if(d<13)  ans+=dp(a,b,c,d+1,x,y)*(13-d)/(54-sum);
	if(x==4){
		Min=10000;
		for(i=0;i<=3;++i)  Min=min(Min,dp(a,b,c,d,i,y));
		ans+=Min/(54-sum);
	}
	if(y==4){
		Min=10000;
		for(i=0;i<=3;++i)  Min=min(Min,dp(a,b,c,d,x,i));
		ans+=Min/(54-sum);
	}
	return f[a][b][c][d][x][y]=ans;
}
int main(){
	scanf("%d%d%d%d",&A,&B,&C,&D);
	double ans=dp(0,0,0,0,4,4);
	printf("%.3lf",(ans>54)?-1:ans);
	return 0;
}