传送门

Solution

对于条件一：记录一个\(cnt\)表示牌个数\(≥2\)的个数

设\(dp_{i,0/1,j,k}\)表示考虑了\(1...i\)，当前是否有对子，以\(i-1\),\(i\)开始的顺子数为\(j\),\(k\)个，的最大面子数

其中\(dp_i\)是一个大小为\(18\)的状态，我们通过搜索可以发现这样的状态很少

不胡的状态有\(2091\)个，可以直接开一个\(map\)来记录，这里要记得重载小于号

将期望看成，如果摸了\(i(i≥13)\)张牌，如果仍然不胡则贡献\(1\)

答案需要求出有\(i\)张牌时仍然不胡的方案数，并乘一个排列数\((4n-i)!\)

仍然考虑dp，设\(f_{i,j,k}\)表示考虑\(1...i\)，当前状态为\(j\)，一共选了\(k\)张牌的方案数

转移时枚举\(i+1\)面值选多少张即可

总复杂度\(O(n^2S)\)，\(S\)表示状态数

Code

#include<bits/stdc++.h>
using namespace std;
#define reg register
inline int read()
{int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}return x*f;
}
const int P=998244353,C[5][5]={{1},{1,1},{1,2,1},{1,3,3,1},{1,4,6,4,1}};
int Mul(int x,int y){return (1ll*x*y)%P;}
int Add(int x,int y){return (x+y)%P;}
int n,w[15],t[105];
void rw(int &x,int y){if(min(y,4)>x)x=y;}
struct State
{int f[18],cnt;State(){memset(f,-1,sizeof f);f[0]=cnt=0;}bool operator <(const State&o)const{if(cnt!=o.cnt) return (cnt<o.cnt);for(int i=0;i<18;++i)if(f[i]!=o.f[i]) return (f[i]<o.f[i]);return false;}bool HU(){return cnt>=7||f[9]>=4;}friend State Trans(State a,int b){register int i,j,k;State c;c.cnt=a.cnt+(b>=2);for(i=0;i<3;++i)for(j=0;j<3;++j)for(k=0;k<3&&i+j+k<=b;++k){if(~a.f[i*3+j]) rw(c.f[j*3+k],a.f[i*3+j]+i+(b-i-j-k)/3);if(~a.f[i*3+j+9])rw(c.f[9+j*3+k],a.f[9+i*3+j]+i+(b-i-j-k)/3);if(i+j+k+2<=b&&~a.f[i*3+j]) rw(c.f[9+j*3+k],a.f[i*3+j]+i);}return c;}
}st[2100];int tot;int tr[2100][5];
std::map<State,int> mp;
int dfs(State cur)
{if(cur.HU()) return 0;if(mp.find(cur)!=mp.end()) return mp[cur];st[mp[cur]=++tot]=cur;int now=tot;for(int i=0;i<=4;++i) tr[now][i]=dfs(Trans(cur,i));return now;
}
int f[2][2100][405],fac[450],inv[450],ans;
int main()
{dfs(State());register int sum,i,j,k,l;for(inv[0]=inv[1]=fac[0]=fac[1]=1,i=2;i<=420;++i)fac[i]=Mul(fac[i-1],i),inv[i]=Mul((P-P/i),inv[P%i]);for(i=2;i<=420;++i) inv[i]=Mul(inv[i],inv[i-1]);n=read();for(i=1;i<=13;++i) w[i]=read(),read(),++t[w[i]];f[0][1][0]=1;for(sum=0,i=1;i<=n;sum+=t[i],++i){memset(f[i&1],0,sizeof f[i&1]);for(j=1;j<=tot;++j)for(k=sum;k<=(i-1)<<2;++k)if(f[(i&1)^1][j][k])for(l=t[i];l<=4;++l)if(tr[j][l]){int nxt=tr[j][l];f[i&1][nxt][k+l]=Add(f[i&1][nxt][k+l],Mul(f[(i&1)^1][j][k],Mul(C[4-t[i]][l-t[i]],Mul(inv[k-sum],fac[k+l-sum-t[i]]))));}}ans=0;for(i=0;i<=(n<<2);ans=Add(ans,Mul(l,fac[(n<<2)-i])),++i)for(l=0,j=1;j<=tot;l=Add(l,f[n&1][j][i]),++j);printf("%d\n",Mul(ans,inv[(n<<2)-13]));return 0;
}

Blog来自PaperCloud，未经允许，请勿转载，TKS！