20200205B组T1组合数计算

组合数计算

题目描述
计算组合数，其中n、m和C(n, m)均满足64位无符号整数取值范围。
输入
第一行为一个整数T (1 <= T <= 10)，表示测试数据组数。
接下来T行，每行有两个数，分别为n和m。
输出
对于每一对n和m，输出C(n, m)的值。每组数据各占一行。
样例输入
2
4 2
5 3
样例输出
6
10
思路
对于小于2000的数据，直接杨辉三角
而对于剩下的数据，用暴力解决
先判断一下m和n-m哪个大，就直接一次性约分
C（5,2）就可以直接变成5* 4 / 2 *1了
然后在乘分子的时候乘一个就去跟分母的因数匹配，找到能整除的先除，最后输出

#include<cstdio>
#include<cstring>
#define rg register
#define ull unsigned long long
using namespace std;
int tot;
ull n,m,t,h,end,s,yh[2005][2005];
bool b[1005];
int main()
{
	for (rg int i=1;i<=2000;i++)
		yh[i][1]=1,yh[i][i]=1;
	for (rg int i=2;i<=2000;i++)
		for (rg int j=2;j<i;j++)
			yh[i][j]=yh[i-1][j]+yh[i-1][j-1];
	scanf("%d",&tot);
	for (;tot;tot--)
	{
		scanf("%llu%llu",&n,&m);
		if (n<2000&&m<2000)
		{
			printf("%llu\n",yh[n+1][m+1]);
			continue;
		}
		if (n==m||n-m==n)
		{
			printf("1\n");
			continue;
		}
		memset(b,true,sizeof(b));
		t=n;
		if (m>n-m) h=m;
		else h=n-m;
		end=n-h;
		h++;
		s=1;
		if (end==1) 
		{
			printf("%llu\n",n);
			continue;
		}
		for (rg ull i=h;i<=t;i++)
		{
			s*=i;
			for (rg ull j=2;j<=end;j++)
			{
				if (s%j==0&&b[j]) s/=j,b[j]=false;
			}
		}
		printf("%llu\n",s);
	}	
	return 0;
}