【CF1924D】Balanced Subsequences 【卡特兰数,组合数,数学】

C F 1924 D . CF1924D. CF1924D.Balanced Subsequences

题意为:给定 n , m , k n,m,k n,m,k,求有多少个由 n n n( m m m) 组成的序列满足最长的合法括号子序列的长度恰为 2 k 2k 2k(对 1 0 9 + 7 10^9+7 109+7 取模)。

思路

我们很容易想到这跟 C a t e l a n Catelan Catelan 数有关。

首先我们知道:长度为 2 k 2k 2k 的合法括号序列包含 k k k 个左括号和 k k k 右括号。

首先有 k > min ⁡ { n , m } k>\min\{n,m\} k>min{n,m},则不合法。

我们现在由 n − k n-k nk 个废物左括号,和 m − k m-k mk 个废物右括号。

不难想出,肯定把废物右括号放在废物左括号前,以免它们再从产生新的贡献。

将括号序列抽象成折线,初始位置在 ( 0 , 0 ) (0,0) (0,0) x x x 轴表示当前括号的总数,而 y y y 轴表示左括号 − - 右括号,也就是说,折线只能往右上或者右下走。

而我们有 m − k m-k mk 个废物右括号会失配,所以最终在某一横坐标下,纵坐标会到达 k − m k-m km,因为当废物右括号失配时,前面的括号肯定全部匹配完毕了。

而我们最终要到达的坐标为 ( n + m , n − m ) . (n+m,n-m). (n+m,nm).

题意转换为:从 ( 0 , 0 ) (0,0) (0,0) ( n + m , n − m ) (n+m,n-m) (n+m,nm) 的合法路径(合法判定为:路径不越过 y = k − m y = k - m y=km 的路径(为什么是 y = k − m y=k-m y=km,因为后面 m − k m-k mk 都是没用的,所以不能合用,此为达到 k k k 的基础))。

我们不难想到用容斥原理:而减去部分就为 ( 0 , 0 ) → ( n + m , 2 × k − n − m ) (0,0)\rightarrow(n+m,2 \times k - n - m) (0,0)(n+m,2×knm) ,其方案数不难得出为: C n + m k − 1 C_{n+m}^{k-1} Cn+mk1

所以 a n s = C n + m k − C n + m k − 1 . ans=C_{n+m}^{k}-C_{n+m}^{k-1}. ans=Cn+mkCn+mk1.

我们就有了一下代码:

#include <iostream>
#include <cstdio>
#include <algorithm>
#define int long long
#define N 5006
using namespace std;
const int mod = 1e9+7;
int n,m,jc[N],k;
int qpow(int a,int b) {
	int res = 1;
	while (b) {
		if (b & 1) res = res * a % mod;
		a = a * a % mod;
		b >>= 1;
	}
	return res;
}
int C(int a,int b) {
	if (a < 0 || b < 0 || a < b) return 0;
	return jc[a] * qpow(jc[b] * jc[a - b] % mod,mod - 2) % mod;
}
signed main(){
	int T;
	jc[0] = 1;
	for (int i = 1;i <= 5000;i ++) jc[i] = jc[i - 1] * i % mod;
	cin >> T;
	for (;T--;) {
		scanf("%lld%lld%lld",&n,&m,&k);
		if (min(n,m) < k) printf("0\n");
		else printf("%lld\n",(C(n + m,k) - C(n + m,k - 1) + 2 * mod) % mod);
	} 
//	cout << (C(n + m,n) - C(n + m,n + 1) + 2 * mod) % mod;
	return 0;
}
  • 24
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值