#include typedef uint64_t u64; typedef int64_t i64; using namespace std; template struct modint{ u64 val; modint(i64 val_=0):val((val_%i64(mod)+i64(mod))%i64(mod)){} modint operator-(){ return (val==0)?0:mod-val; } modint operator+(modint rhs){ return modint(*this)+=rhs; } modint operator-(modint rhs){ return modint(*this)-=rhs; } modint operator*(modint rhs){ return modint(*this)*=rhs; } modint operator/(modint rhs){ return modint(*this)/=rhs; } modint operator^(i64 rhs){ return modint(*this)^=rhs; } modint &operator+=(modint rhs){ val+=rhs.val,val-=((val>=mod)?mod:0); return (*this); } modint &operator-=(modint rhs){ val+=((val>=1; } return (*this)=res; } bool operator==(modint rhs){ return val==rhs.val; } bool operator!=(modint rhs){ return val!=rhs.val; } friend std::ostream &operator<<(std::ostream& os,modint x){ return os<<(x.val); } friend std::istream &operator>>(std::istream& is,modint& x){ u64 t; is>>t,x=t; return is; } }; template struct matrix{ vector> val; int height,width; matrix(int H,int W,T init=0):val(H,vector (W,init)),height(H),width(W){} vector &operator[](int i){ return val[i]; } matrix operator+(matrix rhs){ return matrix(*this)+=rhs; } matrix operator*(matrix rhs){ return matrix(*this)*=rhs; } matrix operator+=(matrix rhs){ for(int i=0;i res(height,rhs.width,0); for(int i=0;iheight;++i){ for(int j=0;jwidth;++k){ res[i][j]+=(*this)[i][k]*rhs[k][j]; } } } return (*this)=res; } matrix identity(int n){ matrix res(n,n); for(int i=0;i0){ if(k%2==1){ res*=now; } now*=now; k/=2; } return (*this)=res; } }; template struct comb{ vector dat,idat; comb(int mx=3000000):dat(mx+1,1),idat(mx+1,1){ for(int i=1;i<=mx;++i){ dat[i]=dat[i-1]*i; } idat[mx]/=dat[mx]; for(int i=mx;i>0;--i){ idat[i-1]=idat[i]*i; } } T operator()(int n,int k){ if(n<0||k<0||n mint; string S; int K; cin>>S>>K; int N=S.size()/2; vector>> dp(2*N+1,vector> (2*N+5,vector (2*N+5))); dp[0][0][0]=1; for(int i=0;i<2*N;++i){ for(int j=0;j<=2*N;++j){ for(int k=0;k<=2*N;++k){ if(S[i]=='('){ dp[i+1][j+1][k]+=dp[i][j][k]; if(0 A(N+1,N+1),B(N+1,1); B[0][0]=1; for(int i=0;i<=N;++i){ A[i][i]=N*(2*N-1)-i*i-(N-i)*(N-i); if(0 C; for(int i=0;i<=N;++i){ ans+=dp[2*N][0][2*i]*B[i][0]/(C(N,i)^2); } cout<