P5279

这个故事告诉我们，要全心全意赞美太阳。
题解

#include<bits/stdc++.h>
using namespace std;
const unsigned long long mod=998244353;
struct data {
	int dp[3][3];
	inline int* operator [](const int& x) {
		return dp[x];
	}
	data() {for(int i=0; i<3; i++)for(int j=0; j<3; j++)dp[i][j]=-1;}
	friend bool operator <(data a,data b) {
		for(int i=0; i<3; i++)for(int j=0; j<3; j++)if(a[i][j]!=b[i][j])return a[i][j]<b[i][j];
		return false;
	}
	friend bool operator !=(data a,data b) {
		for(int i=0; i<3; i++)for(int j=0; j<3; j++)if(a[i][j]!=b[i][j])return true;
		return false;
	}
	inline void trs(data& c,int del) {
		for(int i=0; i<3; i++)for(int j=0; j<3; j++)if(dp[i][j]!=-1)for(int k=0,rem=del-i-j; k<3&&i+j+k<=del; k++,rem--)c[j][k]=max(c[j][k],min(i+dp[i][j]+rem/3,4));
	}
	inline bool ckhu() {
		for(int i=0; i<3; i++)for(int j=0; j<3; j++)if(dp[i][j]>=4)return true;
		return false;
	}
};
struct nod {
	data is_pair[2];
	int cnt_pair;
	friend bool operator <(nod a,nod b) {
		if(a.cnt_pair!=b.cnt_pair)return a.cnt_pair<b.cnt_pair;
		if(a.is_pair[0]!=b.is_pair[0])return a.is_pair[0]<b.is_pair[0];
		if(a.is_pair[1]!=b.is_pair[1])return a.is_pair[1]<b.is_pair[1];
		return false;
	}
	inline void clear() {
		is_pair[0]=data();
		is_pair[1]=data();
		cnt_pair=-2333;
	}
	inline bool ckhu() {
		if(cnt_pair>=7) {
			clear();
			return true;
		}
		if(is_pair[1].ckhu()) {
			clear();
			return true;
		}
		return false;
	}
	friend nod operator +(nod a,int b) {
		if(a.cnt_pair==-2333)return a;
		nod c;
		if(b>=2)a.is_pair[0].trs(c.is_pair[1],b-2);
		a.is_pair[0].trs(c.is_pair[0],b);
		a.is_pair[1].trs(c.is_pair[1],b);
		c.cnt_pair=a.cnt_pair+(b>=2);
		c.ckhu();
		return c;
	}

};
inline nod stat() {
	nod res;
	res.clear();
	res.cnt_pair=0;
	res.is_pair[0][0][0]=0;
	return res;
}
map <nod,int> mrk;
int tot;
int mp[4100][5];
queue <nod> q;
int ed;
inline void bfs() {
	nod st=stat();
	q.push(st);
	mrk[st]=++tot;
	while(!q.empty()) {
		st=q.front();
		q.pop();
		int nu=mrk[st];
		for(int i=0; i<=4; i++) {
			nod tw=st+i;
			if(mrk.count(tw))mp[nu][i]=mrk[tw];
			else mp[nu][i]=mrk[tw]=++tot,q.push(tw);
		}
	}
	st.clear();
	ed=mrk[st];
}
unsigned long long fac[4100],inv[4100],ifac[4100],dp[2][440][4100],xs[4100];
unsigned long long* p1[440];
unsigned long long* p2[440];
int sum[4100];
int n;
inline void pre() {
	fac[0]=1;
	for(int i=1; i<4100; i++)fac[i]=fac[i-1]*i%mod;
	inv[0]=inv[1]=1;
	for(int i=2; i<4100; i++)inv[i]=(mod-mod/i)*inv[mod%i]%mod;
	ifac[0]=1;
	for(int i=1; i<4100; i++)ifac[i]=ifac[i-1]*inv[i]%mod;
}
inline void calc() {
	int cap=n*4-13;
	for(int i=0; i<=cap; i++)p1[i]=dp[0][i];
	for(int i=0; i<=cap; i++)p2[i]=dp[1][i];
	p1[0][mrk[stat()]]=1;
	for(int z=1; z<=n; z++) {
		for(int j=0; j<=4-sum[z]; j++)xs[j]=fac[4-sum[z]]*ifac[4-sum[z]-j]%mod*ifac[j]%mod;
		for(int i=cap; i>=0; i--) {
			for(int j=1; j<=tot; j++)p2[i][j]=0;
			for(int j=1; j<=tot; j++)if(p1[i][j]!=0)for(int tmp=0; tmp<=4-sum[z]&&i+tmp<=cap; tmp++)(p2[i+tmp][mp[j][tmp+sum[z]]]+=p1[i][j]*xs[tmp])%=mod;
		}
		for(int i=0; i<=cap; i++)swap(p1[i],p2[i]);
	}
}
int main() {
	bfs();
	pre();
	scanf("%d",&n);
	for(int i=1,w,t; i<=13; i++)scanf("%d%d",&w,&t),sum[w]++;
	calc();
	unsigned long long ans=0;
	for(int p=1,q=4*n-14; p<=4*n-13; p++,q--)for(int j=1; j<=tot; j++)if(j!=ed)(ans+=p1[p][j]*fac[p]%mod*fac[q])%=mod;
	printf("%lld",ans*ifac[4*n-13]%mod+1);
	return 0;
}