code:
//预处理fac
int C(int n,int m){
if(m<0||m>n)return 0;
return fac[n]*ppow(fac[m]*fac[n-m]%mod,mod-2,mod)%mod;
}
//预处理fac和 fac的inv
int C(int n,int m){
if(m<0||m>n)return 0;
return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
..
int fac[maxm],inv[maxm];
int C(int n,int m){
if(m<0||m>n)return 0;
return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int A(int n,int m){
if(n<0||n-m<0)return 0;
return fac[n]*inv[n-m]%mod;
}
void init(){
fac[0]=1;
for(int i=1;i<maxm;i++)fac[i]=fac[i-1]*i%mod;
inv[maxm-1]=ppow(fac[maxm-1],mod-2,mod);
for(int i=maxm-2;i>=0;i--)inv[i]=(i+1)*inv[i+1]%mod;
}
//打表
int c[maxm][maxm];
void init(){
c[0][0]=1;
for(int i=1;i<maxm;i++){
c[i][0]=1;
for(int j=1;j<=i;j++){
c[i][j]=c[i-1][j-1]+c[i-1][j];
//c[i][j]%=mod;
}
}
}
球&盒问题结论:
以下均为n球m盒
1.球同盒不同非空:C(n-1,m-1)
2.球同盒不同可空:C(n+m-1,m-1)或C(n+m-1,n)