卡特兰数 题解

bsoj7107 ，来源不明。

卡特兰数

题目描述

今天，接触信息学不久的小A刚刚学习了卡特兰数。

卡特兰数的一个经典定义是，将 (n​) 个数依次入栈，合法的出栈序列个数。

小A觉得这样的情况太平凡了。于是，他给出了 (m) 组限制，每个限制形如 ((f_i,g_i)) ，表示 (f_i) 不能在 (g_i) 之后出栈。

他想求出：在满足了这 (m) 组限制的前提下，共有多少个合法的出栈序列。他不喜欢大数，你只需要求出答案在模 (998244353) 意义下的值即可。

输入格式

输入第一行为两个非负整数 (n) 、 (m) ，含义题面已给出。

接下来 (m​) 行，每行两个正整数， ((f_i,g_i)​) 表示一组限制。

输出格式

输出一行，为一个非负整数，表示你求得的答案 (mod 998244353​) 。

样例输入

3 1
2 3

样例输出

3

样例解释

可以验证 ({1,2,3 }​) ， ({2,1,3 }​) ， ({2,3,1 }​) 都是合乎条件的。

数据规模

(n le 300) , (m le frac{n(n-1)}{2}) , (f_i,g_i le n) 。

部分数据的 (m) 较小。

考前一天写题解会不会涨RP

套路但是还是感觉很巧妙的转换。关键在于把限制形式化的写出来方便观察。

考虑Catalan数的递推公式推导方式——枚举最后出栈的元素 (k​) 。其递推式为

[egin{aligned} c_0 &= 1 \ c_n &= sum_{k=1}^n c_{k-1} c_{n-k} end{aligned} ]

加入限制后无法使用排列的相对顺序dp，故需改造原来的递推。定义 (h_{l,r}) 为排列 ([l,r]) 不同的出栈序列方案数，若暂不考虑限制，有

[h_{l,r} = sum_{k=l}^r h_{l,k-1} h_{k+1,r} ]

显然 (c_n = h_{1,n}​) 。

递推中，出栈序列的相对顺序为 (([l,k-1],[k+1,r],{ k })) 。限制 ((f_i, g_i)) 要求 (f_i) 在 (g_i) 之前出栈。故满足 (g_i in [l,k-1], f_i in [k,r]) 或者 (g_i in [k+1,r], f_i = k) 的 (k) 都不合法，不应对 (h_{l,r}) 产生贡献。故在此dp中枚举所有限制判断合法性即可获得 (O(n^3 m)​) 的算法。

这种方法看着就很naive，考虑优化合法性判断。容易发现把 ((f_i,g_i)) 投到二维平面上后，判断合法性只需询问 ((k,l) - (r,k-1)) 、 ((k,k+1) - (k,r)) 两个矩形中是否有点。二维前缀和处理即可 (O(1)) 判断。于是这样就可以 (O(n^3)) 通过本题了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
using namespace std;
typedef long long ll;
ll Rd(){
	ll ans=0;bool fh=0;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-') fh=1; c=getchar();}
	while(c>='0'&&c<='9') ans=ans*10+c-'0',c=getchar();
	if(fh) ans=-ans;
	return ans;
}
const ll MOD=998244353;
#define _ %MOD
ll PMod(ll x){
	if(x>=MOD) return x-MOD;
	else if(x<0) return x+MOD;
	else return x;
}

const ll MXN=305;
ll N,M;
ll A[MXN][MXN];
ll S[MXN][MXN];
void SpawnSum(){
	for(ll i=1;i<=N;i++)
		for(ll j=1;j<=N;j++)
			S[i][j]=S[i-1][j]+S[i][j-1]-S[i-1][j-1]+A[i][j];
}
bool Check(ll x1,ll y1,ll x2,ll y2){
	return S[x2][y2]-S[x1-1][y2]-S[x2][y1-1]+S[x1-1][y1-1];
}
ll f[MXN][MXN];
void Solve(){
	for(ll i=0;i<=N;i++) f[i+1][i]=1;
	for(ll len=1;len<=N;len++){
		for(ll l=1;l+len-1<=N;l++){
			ll r=l+len-1;
			for(ll k=l;k<=r;k++){
				if(Check(k,l,r,k-1)) continue;
				if(Check(k,k+1,k,r)) continue;
				f[l][r]=(f[l][r]+f[l][k-1]*f[k+1][r])_;
			}
		}
	}
	cout<<f[1][N];
}
int main(){
	N=Rd();M=Rd();
	for(ll i=0;i<=N;i++) for(ll j=0;j<=N;j++) A[i][j]=S[i][j]=0;
	for(ll i=1;i<=M;i++){
		ll x=Rd(),y=Rd();
		A[x][y]=1;
	}
	SpawnSum();
	Solve();
	return 0;
}

2020/12/04