hdu 新生晚会

新生晚会
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 12054 Accepted Submission(s): 4169


Problem Description
开学了,杭电又迎来了好多新生。ACMer想为新生准备一个节目。来报名要表演节目的人很多,多达N个,但是只需要从这N个人中选M个就够了,一共有多少种选择方法?


Input
数据的第一行包括一个正整数T,接下来有T组数据,每组数据占一行。
每组数据包含两个整数N(来报名的人数,1<=N<=30),M(节目需要的人数0<=M<=30)



Output
每组数据输出一个整数,每个输出占一行


Sample Input
5
3 2
5 3
4 4
3 6
8 0


Sample Output
3
10
1
0

1


这一倒要用gcd即最大公约数对每一位进行约分,这样就不会产生溢出了


# include <iostream>
using namespace std;
long long int ans_a;
long long int ans_b;
int a[100];//数组a和数组b分别用来存分子的阶乘和分母的阶乘
int b[100];
int gcd(int a, int b){
	if(a%b!=0){
		return gcd(b, a%b);
	}
	return b;
}
int main(){
	int n, t, i, j, k, temp, m;
	int cnt1, cnt2;
	cin>>t;
	for(i=1; i<=t; i++){
		cin>>n>>m;
		if(m>n){
			printf("0\n");
			continue;
		}
		if(m==n||m==0){
			printf("1\n");
			continue;
		}
		cnt1=0;cnt2=0;
		ans_a=1;ans_b=1;
		for(j=1; j<=n; j++){
			a[cnt1++]=j;
		}
		for(j=1; j<=m; j++){
			b[cnt2++]=j;
		}
		for(j=1; j<=n-m; j++){
			b[cnt2++]=j;
		}
		for(j=0; j<=cnt1-1; j++){
			for(k=0; k<=cnt2-1; k++){
				temp=gcd(a[j], b[k]);
				a[j]=a[j]/temp;
				b[k]=b[k]/temp;
			} 
		}
		for(j=0; j<=cnt1-1; j++){
			ans_a=ans_a*a[j];
		}
		for(j=0; j<=cnt2-1; j++){
			ans_b=ans_b*b[j];
		}
		cout<<(ans_a/ans_b)<<endl;
	}
	return 0;
}


你可能感兴趣的:(算法,阶乘,模拟,ACM,gcd)