组合数

n!中质因子p的个数

版本一：

int cal(int n,int p){
 int ans = 0;
 for(int i=2;i <= n;i++){//遍历2~n
    int temp = i;
    while(temp%p==0){//只要temp还是p的倍数
       ans++;  //p的个数加1
       temp/=p; //temp除以p
    }
  }
  return ans;
}

版本二：

int cal(int n,int p){
  int ans = 0;
  while(n){
    ans+=n/p; //累加n/p^k
    n/=p;     //相当于分母多乘一个p
   }
  return ans;
}

版本三：

int cal(int n,int p){
  if(n < p) return 0;//n<p时1~n中不可能有质因子p
  return n / p + cal(n / p, p);//返回n/p加上（n/p）!中质因子p的个数
 }

组合数的计算

方法一：

long long C(long long n,long long m){
 long long ans = 1;
 for(long long i = 1;i <= n;i++){
    ans*=i;
  }
 for(long long i = 1;i <= m;i++){
   ans/=i;
  }
 for(long long i = 1;i <= n-m;i++){
   ans/=i; 
 }
 return ans;
}

方法二：

long long res[67][67]={0};
long long C(long long n,long long m){
 if(m == 0||m == n) return 1;
 if(res[n][m]!=0) return res[n][m];
 return res[n][m]=C(n-1,m) + C(n-1,m-1);
}

const int n = 60;
void calC(){
 for(int i=0;i<=n;i++){
    res[i][0] = res[i][i] = 1;
  }
 for(int i=2;i<= n;i++){
   for(int j=0;j<=i/2;j++){
     res[i][j] = res[i-1][j] + res[i-1][j-1];
     res[i][i-j] = res[i][j];
   }
  }
}

方法三：

long long C(long long n, long long m){
 long long ans = 1;
 for(long long i=1;i<= m;i++){
     ans = ans*(n-m+i) / i; //注意一定先乘再除
 }
 return ans;
}

方法一：

递归：

int res[1010][1010]={0};
int C(int n,int m,int p){
  if(m == 0 || m == n) return 1;
  if(res[n][m] != 0) return res[n][m];
  return res[n][m] = (C(n-1,m)+C(n-1,m-1))%p;
}

递推：

void calC(){
 for(int i=0;i<= n;i++){
   res[i][0]=res[i][i]=1;
  }
 for(int i=2;i<=n;i++){
  for(int j=0;j<=i/2;j++){
      res[i][j]=(res[i-1][j]+res[i-1][j-1])%p;
      res[i][i-j]=res[i][j];
   }
 }
}

方法二：

int prime[maxn];
int C(int n,int m,int p){
  int ans =1;
  for(int i=0;prime[i]<=n;i++){
    int c=cal(n,prime[i]) - cal(m,prime[i]) - cal(n-m,prime[i]); 
    //快速幂计算prime[i]^c%p
    ans = ans * binaryPow(prime[i],c,p)%p;
   }
  return ans;
}

方法三：

1，m < p,且p时素数

int C(int n,int m,int p){
 int ans =1;
 for(int i=1;i<=m;i++){
   ans=ans*(n-m+i)%p;
   ans=ans*inverse(i,p)%p;//求i模p的逆元
 }
  return res;
}

2，m任意，且p是素数

int C(int n, int m,int p){
 int ans=1,numP=0; //ans存放计算结果，numP统计分子中的p比分母中的p多几个
 for(int i=1;i<=m;i++){
   int temp = n -m +i;//分子
   while(temp%p==0){  //去掉分子中的所有p，同时累计numP
    numP++;
    temp/=p;
   }
   ans=ans*temp%p;
   temp =i;
  while(temp%p==0){
    numP--;
    temp/=p;
   }
  ans=ans*inverse(temp,p)%p;
  }
 if(numP>0) return 0;
 else return ans;
}

方法四：Lucas定理

int Lucas(int n,int m){
 if(m == 0) return 0;
 return C(n%p,m%p)*Lucas(n/p,m/p)%p; 
}