组合数判断溢出问题

组合数溢出

对于以下一段有问题的组合数代码：

#include<stdio.h>
#include<stdlib.h>
long long factorial(int n) {
	long long m = 1;
	for (int i = 1; i <= n; i++)
	{
		m *= i;
	}
	return m;
}
long long c(int n,int m) {
	return factorial(n) / (factorial(m) * factorial(n - m));
}
int main() {
	int n = 21;
	int m = 1;
    long long d = 0;
	d=c(n,m);
	printf("%11d",d);
	return 0;
}

问题：
在测试的时候如果规定n=21,m=1的话会发现：

程序的返回值竟然是-1，本来应该是21,但是结果不符合，显然在这个过程中出现了溢出。
其实，如果当n=20,m=1的时候结果还是正确的，是20，从21开始就不行了

解决方法

其实最好的解决办法就是在这个过程中进行“约分”
像 n!/m!=(m+1)(m+2)…(n-1)n 虽然在过程中不能说完全就不存在中间结果溢出，以下给出一种办法，虽然说这个方法不是说一定不存在缺陷，但是至少比上面的能“坚持”的范围更广了。
将c函数改为：

long long c(int n, int m)
{
	if (m < n - m) m = n - m;//数学问题，自行思考
	long long ans = 1;
	for (int i = m + 1; i <= n; i++) ans *= i;
	for (int i = 1; i <= n - m; i++) ans /= i;
	return ans;
}

去测试会发现，结果已经正确：