bzoj 1272: [BeiJingWc2008]Gate Of Babylon

题意：

传送门

题解：

先考虑没有限制的情况，用隔板法显然就是$C_{n+m}^n$
那么因为t很小，所以容斥下就好了。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
LL n,m,t,p,B[20];
LL f[100010],inv[100010];
LL C(LL n,LL m)
{
    if(m>n) return 0;
    return f[n]*inv[n-m]%p*inv[m]%p;
}
LL lucas(LL n,LL m)
{
    if (m>n) return 0;
    LL ans=1LL;
    for (;m;n/=p,m/=p)
        ans=ans*C(n%p,m%p)%p;
    return ans;
}
void pre()
{
    f[1]=inv[1]=inv[0]=1;
    for(LL i=2;i<=p;i++) f[i]=f[i-1]*i%p,inv[i]=inv[p%i]*(p-p/i)%p;
    for(LL i=2;i<=p;i++) inv[i]=inv[i]*inv[i-1]%p;
}
LL ans=0;
void dfs(LL x,LL k,LL num)
{
    if(k<0) return;
    if(x>t)
    {
        LL tmp=lucas(k+n,n);
        if(num&1) ans=(ans-tmp)%p;
        else ans=(ans+tmp)%p;
        return;
    }
    dfs(x+1,k-B[x]-1,num+1);
    dfs(x+1,k,num);
}
int main()
{
    scanf("%lld %lld %lld %lld",&n,&t,&m,&p);
    for(LL i=1;i<=t;i++) scanf("%lld",&B[i]);
    pre();
    dfs(1,m,0);
    printf("%lld",(ans+p)%p);
}