【组合数学】 HDOJ 5184 Brackets

附上BC的官方题解：

当n为奇数的时候答案是0。
先判断字符串的前面是否符合括号匹配，即对于任何前缀左括号个数>=右括号个数。
设左括号个数为a右括号个数为b, m=n/2，问题可以转化为在平面中从座标(a,b)沿网格走到(m,m) 且不跨过x=y这一条直线的方法数。数据太大，普通DP和搜索都不行的。
问题可以进一步转化为从(a-n,b-n)到(0,0)且不跨过x=y的方法数。再对称一下，转化到(0,0)到(n-b,n-a)不跨过x=y的方法数。对于从(0,0)点走到(p,q)点不跨过x=y的方法数是

p−q+1p+1Cqp+q

证明如下：
我们可以通过总的数目来减掉非法的数目即可。
把(0,0)和(p,q)都往下移一格，非法数目即为(0,-1)到(p,q-1)且路径中至少有一点和x=y相交的方法数。记(d,d)为从(0,-1)到(p,q-1)路径中最先和x=y相交的点。则由于对称性(-1,0)到(d,d)的方法数和(0,-1)到(d,d)的方法数是相同的。所以(0,-1)到(p,q-1)且与x=y相交的方法数和(-1,0)到(p,q-1)的方法数是相同的。
所以答案是

Cqp+q−Cq−1p+q=p−q+1p+1Cqp+q

然后对100W以内的数字进行一个阶乘处理，就可以O(1)得出答案了。

#include <iostream>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <climits>
#include <cstdlib>
#include <cmath>
#include <time.h>
#define maxn 1000005
#define maxm 100005
#define eps 1e-7
#define mod 1000000007
#define INF 0x3f3f3f3f
#define PI (acos(-1.0))
#define lowbit(x) (x&(-x))
#define mp make_pair
#define ls o<<1
#define rs o<<1 | 1
#define lson o<<1, L, mid
#define rson o<<1 | 1, mid+1, R
#define pii pair<int, int>
#pragma comment(linker, "/STACK:16777216")
typedef long long LL;
typedef unsigned long long ULL;
//typedef int LL;
using namespace std;
LL qpow(LL a, LL b){LL res=1,base=a;while(b){if(b%2)res=res*base;base=base*base;b/=2;}return res;}
LL powmod(LL a, LL b){LL res=1,base=a;while(b){if(b%2)res=res*base%mod;base=base*base%mod;b/=2;}return res;}
//headLL f[maxn];
LL g[maxn];
char s[maxn];
int n;void init()
{f[0] = 1;for(int i = 1; i <= 1000000; i++) f[i] = f[i-1] * i % mod;g[1000000] = powmod(f[1000000], mod - 2);for(int i = 999999; i >= 0; i--) g[i] = g[i+1] * (i+1) % mod;
}inline LL calc(int a, int b)
{return f[a+b] * (b + 1 - a) % mod * g[a] % mod * g[b+1] % mod;
}void work()
{if(n % 2) {printf("0\n");return;}int a = 0, b = 0, ok = 1;for(int i = 0; s[i]; i++) {if(s[i] == '(') a++;else b++;if(b > a) ok = 0;}if(!ok) {printf("0\n");return;}n /= 2;a = n - a;b = n - b;if(a < 0 || b < 0) {printf("0\n");return;}printf("%I64d\n", calc(a, b));
}int main()
{init();while(scanf("%d%s", &n, s)!=EOF) work();return 0;
}

【组合数学】 HDOJ 5184 Brackets相关推荐

HDOJ 5184 Brackets 卡特兰数扩展
既求从点(0,0)仅仅能向上或者向右而且不穿越y=x到达点(a,b)有多少总走法... 有公式: C(a+b,min(a,b))-C(a+b,min(a,b)-1) /// 折纸法证明卡特兰数: h ...
卡特兰数相关问题 hdu 5184 Brackets
题解: 当n为奇数的时候答案是0. 先判断字符串的前面是否符合括号匹配,即对于任何前缀左括号个数>=右括号个数. 设左括号个数为a右括号个数为b, m=n/2,问题可以转化为在平面中从座标(a, ...
BestCoder Round #32
问题描述目前,我们用PM2.5的含量来描述空气质量的好坏.一个城市的PM2.5含量越低,它的空气质量就越好.所以我们经常按照PM2.5的含量从小到大对城市排序.一些时候某个城市的排名可能上升,但是他 ...
【HDU 5184】 Brackets （卡特兰数）
Brackets Problem Description We give the following inductive definition of a "regular brackets& ...
【HDOJ】1261 字串数【组合数学--排列+代数】
传送门:字串数题意: sum(A1,A2,...,An)的阶乘/(A1!)/(A2!)/.../(An!) 因为数字可能很大,所以要用到大数类,Java的大数用起来比较爽,就用Java写了. 好不容 ...
HDOJ ACM 题目
转载 HDOJ 题目分类(转) 1001 整数求和水题 1002 C语言实验题--两个数比较水题 1003 1.2.3.4.5... 简单题 1004 渊子赛马排序+贪心的方法归并 1005 H ...
(区间dp 或记忆化搜素 )Brackets -- POJ -- 2955
http://poj.org/problem?id=2955 Description We give the following inductive definition of a "reg ...
不畏浮云遮望眼--离散数学和组合数学
不畏浮云遮望眼,基础很重要! 离散数学是算法和数据结构的基础,而算法和数据结构又是什么的基础?不解释了. 1.<离散数学及其应用> 作者: (美)Kenneth H. R 出版社: 机械工 ...
并查集 HDOJ 1232 畅通工程
题目传送门 1 /* 2 并查集(Union-Find)裸题 3 并查集三个函数:初始化Init,寻找根节点Find,连通Union 4 考察:连通边数问题 5 */ 6 #include <c ...

【组合数学】 HDOJ 5184 Brackets

【组合数学】 HDOJ 5184 Brackets相关推荐

最新文章

热门文章