【題解】[牛客網NOIP賽前集訓營-提高組(第五場)]C.串串 組合數
阿新 • • 發佈:2018-11-02
#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; }
總結
避免重複很關鍵,以前也有一道類似的