容器

题目描述

DP

用括号表示一段区间。
设f[i,j,k]表示做到第i个位置，目前未被匹配的左括号有j个，还有k只跳蚤没有左括号。
转移可以直接枚举这个位置放多少右括号以及下一个位置放多少左括号，乘上对应组合数。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
using namespace std;
typedef long long ll;
const int maxn=40+10,mo=1011110011;
int f[maxn][maxn][maxn],fac[maxn],inv[maxn];
int i,j,k,x,y,l,t,n,m,p,ans;
int qsm(int x,int y){
    if (!y) return 1;
    int t=qsm(x,y/2);
    t=(ll)t*t%mo;
    if (y%2) t=(ll)t*x%mo;
    return t;
}
int C(int n,int m){
    if (n<m||m<0) return 0;
    return (ll)fac[n]*inv[m]%mo*inv[n-m]%mo;
}
int main(){
    freopen("container.in","r",stdin);freopen("container.out","w",stdout);
    scanf("%d%d%d",&n,&m,&p);
    fac[0]=1;
    fo(i,1,40) fac[i]=(ll)fac[i-1]*i%mo;
    inv[40]=qsm(fac[40],mo-2);
    fd(i,39,0) inv[i]=(ll)inv[i+1]*(i+1)%mo;
    f[0][0][m]=1;
    fo(i,0,n)
        fo(j,0,p)
            fo(k,0,m)
                if (f[i][j][k])
                    fo(x,0,j)
                        fo(y,0,k)
                            (f[i+1][j-x+y][k-y]+=(ll)f[i][j][k]*C(j,x)%mo*C(k,y)%mo)%=mo;
    ans=f[n+1][0][0];
    (ans+=mo)%=mo;
    printf("%d\n",ans);
}