CodeForces 785 D.Anton and School - 2（组合数学）

Description

定义一个长度为$n$的括号序列$s$是一个$RSBS$当且仅当其满足以下四个条件：

1.长度非空($n\neq 0$)

2.长度为偶数

3.前$\frac{n}{2}$个括号是左括号

4.后$\frac{n}{2}$个括号是右括号

现在给出一个括号序列，问其有多少个是$RSBS$的子序列

Input

输入一个长度不超过$200000$的括号序列

Output

输出该括号序列的是$RSBS$的子序列个数，结果模$10^9+7$

Sample Input

$)(()()$

Sample Output

6

Solution

设以第$i$个左括号为子序列中最后一个括号，那么如果要形成一个$RSBS$，则需要从该左括号前面选取$i$个左括号，从该左括号后面选取$i+1$个右括号，假设该左括号前共有$l$个左括号，后面共有$r$个右括号，则方案数为$\sum\limits_{i=0}^{min(l,r-1)}C_l^i\cdot C_r^{i+1}=\sum\limits_{i=1}^{min(l,r-1)}C_l^i\cdot C_r^{r-1-i}=C_{l+r}^{r-1}$

故只要统计前$i$个括号中左括号数量$l[i]$及第$i$个括号到最后一个括号中右括号数量$r[i]$，答案即为$\sum\limits_{i=1}^{|s|}C_{l[i]+r[i]-1}^{r[i]-1}[s[i]='(']$

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=200005;
#define mod 1000000007
int fact[maxn],inv[maxn];
void init(int n=2e5)
{
    fact[0]=1;
    for(int i=1;i<=n;i++)fact[i]=(ll)i*fact[i-1]%mod;
    inv[1]=1;
    for(int i=2;i<=n;i++)inv[i]=mod-(ll)(mod/i)*inv[mod%i]%mod;
    inv[0]=1;
    for(int i=1;i<=n;i++)inv[i]=(ll)inv[i-1]*inv[i]%mod; 
}
int C(int n,int m)
{
    if(m<0||m>n)return 0;
    return (ll)fact[n]*inv[m]%mod*inv[n-m]%mod;
}
void add(int &x,int y)
{
    x=x+y>=mod?x+y-mod:x+y;
}
int l[maxn],r[maxn];
char s[maxn]; 
int main()
{
    init();
    scanf("%s",s+1);
    int n=strlen(s+1);
    l[0]=r[n+1]=0;
    for(int i=1;i<=n;i++)l[i]=l[i-1]+(s[i]=='(');
    for(int i=n;i>=1;i--)r[i]=r[i+1]+(s[i]==')');
    int ans=0;
    for(int i=1;i<=n;i++)
        if(s[i]=='(')add(ans,C(l[i]-1+r[i],r[i]-1));
    printf("%d\n",ans);
    return 0;
}