组合数

一、关于 n! 的一个问题
求解 n! 中有几个质因子
n! 中有 (n/p + n/p² + n/p³+……）个质因子p，其中除法均为向下整除

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

利用算法何以很快计算出 n! 的末尾有几个零：由于末尾0的个数等于 n！中因子10的个数，而这又等于 n! 中质因子5的个数，因此只需带入 cal(n,5) 就可获得结果。
递归版：
n!中质因子p的个数，实际上等于1~n中 p 的倍数的个数 n/p加上(n/p)!中质因子p的个数。

int cal(int n,int p){
	if(n < p) return 0;
	return n / p + cal(n / p,p);
}

二、组合数的计算
组合数C(n,m)是指从 n 个不同元素中选取m个元素的方案（m<=n),
C(n,m) = n! / (m!(n-m)!)

后文直接用 LL 代替 long long

#define LL long long

方法一：通过定义式直接计算（n<=20）

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

方法二：通过递推公式计算
一般在 n=67、m=33时溢出
C(n,m) = C(n-1,m) + C(n-1,m-1)
简洁的递归：

LL C(LL n,LL m){
	if(m==0||m==n) return 1;
	else return C(n-1,m)+C(n-1,m-1);
}

去除重复运算的递归：

LL res[67][67] = {0};
LL C(LL n,LL 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);	//赋值res[n][m]并返回
}

递推打表：

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];	//递归计算C（i，j)
			res[i][i-j] = res[i][j];		//C(i,i-j) = C(i,j)
		}
	}
}

方法三：通过定义式的变形来计算
一般在n=62,m=31时开始溢出

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