【题解】[牛客网NOIP赛前集训营-提高组（第五场）]C.串串 组合数

题目链接

---

#include<cstdio>
const int N=4e3+10,mod=1e9+7;
int a,b,c,d,fac[N],inv[N],ans;
int C(int n,int m)
{
	if(n<0||m<0)return 0;
	return 1ll*fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int main()
{
	//freopen("in.txt","r",stdin);
	fac[0]=inv[0]=inv[1]=1;
	for(int i=1;i<N;i++)fac[i]=(1ll*fac[i-1]*i)%mod;
	for(int i=2;i<N;i++)inv[i]=(mod-1ll*mod/i*inv[mod%i]%mod)%mod;
	for(int i=2;i<N;i++)inv[i]=(1ll*inv[i-1]*inv[i])%mod;
	scanf("%d%d%d%d",&a,&b,&c,&d);
	if(!c||!d)
	{
		printf("%d\n",C(a+b,a));return 0;
	}
	a-=c;b-=d;
	for(int i=0;i<=a;i++)
	    for(int j=0;j<=b;j++)
	        ans=(ans+1ll*C(i+j,i)*C(a-i+d-1,d-1)%mod*C(b-j+c-1,c-1)%mod)%mod;
	printf("%d\n",1ll*ans*C(c+d,c)%mod);
	return 0;
}

总结

避免重复很关键，以前也有一道类似的