うまくいけば、この問題は、あなたが消える麻雀を愛することはありません
私たちは、各ブランド柯が表示された回数をカウントし、何の統計は、いくつかのカードではありません
カードのデッキを宣告し、次のようなことができない、この質問
この質問の場合:
セット\([I] [J F ] [K] [0/1] \) 前者表す\(Iは\)カードの種類、ジュンコ\は、((I-1、I 、I + 1)\) が登場\ (J \)回、ジュンコ\は、((I、I + 1 、I + 2)\) が登場\(K \)と/鳥の頭の数まで顔せず、回。転送は比較的簡単です
私たちは、柯見つけるために(私は\)\この寸法はそれほど重要ではありませんDP値が超えていない強制\(4 \) (以上\(4 \)無駄にする)、鳥の頭のカウントが超えていない(7 \)\(類似)を、検索Keが、状態の異なる性質の発作への総バースト\ (2091 \) A
KE KEは、各状態の後ろに追加する([0,4] \でX \ \) KE自動機械を形成するために張+1ポイントカード、私たちは彼女と無免許の自動機と呼ばれます
私たちは、それぞれの状態、柯は上に行くと、ライセンスのない自動機してもらう、とするかどうかを決定することができます
私たちは、セット\(DP [I] [J ] [k]を\) の前に見て示している(私は\)\カードの種類は、無免許の自動機械にとに\(J \)状態、触れた\(k個\)張を胡、柯の期待に最終カウントの種類のブランド、数ではありません
我々柯スクロールする配列を使用するには、\(私は\)最適化されたスペースをロールアウト
#include <bits/stdc++.h>
#define mod 998244353
#define N 405
#define getchar nc
using namespace std;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int read()
{
register int x=0,f=1;register 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;
}
inline void write(register int x)
{
if(!x)putchar('0');if(x<0)x=-x,putchar('-');
static int sta[20];register int tot=0;
while(x)sta[tot++]=x%10,x/=10;
while(tot)putchar(sta[--tot]+48);
}
inline int Min(register int a,register int b)
{
return a<b?a:b;
}
inline int Max(register int a,register int b)
{
return a>b?a:b;
}
struct node{
int f[3][3][2],cnt;
inline void init()
{
memset(f,-1,sizeof(f));
f[0][0][0]=cnt=0;
}
inline int hu()
{
if(cnt>=7)
return 1;
for(register int i=0;i<3;++i)
for(register int j=0;j<3;++j)
if(f[i][j][1]>=4)
return 1;
return 0;
}
}rt,S[2100];
bool operator < (node a,node b){
if(a.cnt!=b.cnt)
return a.cnt<b.cnt;
for(register int i=0;i<3;++i)
for(register int j=0;j<3;++j)
for(register int k=0;k<2;++k)
if(a.f[i][j][k]!=b.f[i][j][k])
return a.f[i][j][k]<b.f[i][j][k];
return 0;
}
int tot=0;
map<node,int> ma;
inline node trans(register node u,register int w)
{
node v;
v.init();
v.cnt=Min(u.cnt+(w>=2),7);
for(register int i=0;i<3;++i)
for(register int j=0;j<3;++j)
{
if(~u.f[i][j][0])
{
for(register int k=0;k<3&&i+j+k<=w;++k)
v.f[j][k][0]=Max(v.f[j][k][0],Min(u.f[i][j][0]+i+(w-i-j-k>=3),4));
if(w>=2)
for(register int k=0;k<3&&i+j+k<=w-2;++k)
v.f[j][k][1]=Max(v.f[j][k][1],Min(u.f[i][j][0]+i,4));
}
if(~u.f[i][j][1])
for(register int k=0;k<3&&i+j+k<=w;++k)
v.f[j][k][1]=Max(v.f[j][k][1],Min(u.f[i][j][1]+i+(w-i-j-k>=3),4));
}
return v;
}
inline void build(register node u)
{
if(u.hu())
return;
if(ma.find(u)!=ma.end())
return;
ma[u]=++tot;
S[tot]=u;
for(register int i=0;i<=4;++i)
build(trans(u,i));
}
int n,s[N],ans;
int fac[N],inv[N],invf[N],tr[2100][5],f[2][2100][N];
inline int C(register int n,register int m)
{
return 1ll*fac[n]*invf[m]%mod*invf[n-m]%mod;
}
int main()
{
rt.init();
build(rt);
invf[0]=inv[0]=inv[1]=fac[0]=1;
for(register int i=1;i<N;++i)
fac[i]=1ll*fac[i-1]*i%mod;
for(register int i=2;i<N;++i)
inv[i]=1ll*inv[mod%i]*(mod-mod/i)%mod;
for(register int i=1;i<N;++i)
invf[i]=1ll*invf[i-1]*inv[i]%mod;
n=read();
for(register int i=1;i<=13;++i)
++s[read()],read();
for(register int i=1;i<=tot;++i)
for(register int j=0;j<=4;++j)
tr[i][j]=ma[trans(S[i],j)];
f[0][1][0]=1;
for(register int i=1,sum=0;i<=n;++i)
{
int now=i&1,pre=now^1;
memset(f[now],0,sizeof(f[now]));
for(register int j=1;j<=tot;++j)
for(register int k=s[i];k<=4;++k)
{
if(!tr[j][k])
continue;
int w=1ll*C(4-s[i],k-s[i])*fac[k-s[i]]%mod;
for(register int l=0;l<=4*n-k;++l)
{
if(!f[pre][j][l])
continue;
f[now][tr[j][k]][l+k]=(0ll+f[now][tr[j][k]][l+k]+1ll*f[pre][j][l]*w%mod*C(k+l-sum-s[i],k-s[i])%mod)%mod;
}
}
sum+=s[i];
}
for(register int i=13,w=1;i<=4*n;++i)
{
int now=0;
for(register int j=1;j<=tot;++j)
now=(0ll+now+f[n&1][j][i])%mod;
ans=(0ll+ans+1ll*now*w%mod)%mod;
w=1ll*w*inv[4*n-i]%mod;
}
write(ans);
return 0;
}