xjoj171 水晶球

Description

很久很久以前，有一个叫小猫猫的女巫，她用自己的法力制作了 n (1 ≤ n ≤ 1,000,000,000) 个水晶球。一旦这些水晶球被连起来，它们将会产生一股令人难以想象的力量。为防止别有用心的人使用这些水晶球，小猫猫又制作了 f (1 ≤ f ≤ 1,000,000,000)个外表和之前的一样，但是没有任何法力的水晶球，并把这些水晶球放在了一起。

　　大魔头飞猪打算找出哪些水晶球是有法力的，并用它们来统治世界。为达到这一邪恶企图，飞猪召集了一批地精工人来尝试水晶球的各种组合，一共 Cn+fn 　种可能（即便这样，找到真正的水晶球也是遥遥无期……）。飞猪打算把这些任务分配给每个地精工人，而且每种组合只做一次试验。令她头疼的是，这些地精们还成立了一个工会，它们需要保证所有参加工作的地精的工作量是完全一样的。这样，就需要选择合适数量的地精来为飞猪大魔头工作，以保证这些工作能被平均地分配下去。

　　为了节省时间，飞猪需要召集尽可能多的地精。目前一共有 g (1 ≤ g ≤ 100000) 个地精愿意为她工作。大魔头飞猪重金聘请你来为她写个程序，要你计算出实际为她工作的地精数（小于等于 g），这个数必须能整除 Cn+fn　，并且尽可能大。

Input

第一行，一个整数 C (1≤ C ≤100)，表示一共有 C 组数据。

后面 C 行，每行三个整数 n, f, g 用空格分开。

Output

C 行，每行对应输出最多能为大魔头工作的地精数。

Sample Input

4
3 1 3
3 3 50
4 3 4
15634 456 5000

Sample Output

2
20
1
4990

分析，其实就是求C(n,m)是否能够被k整除，这个问题可以用素因子分解的方法解决：即对于k的每一个素因子p，看C(n,m)中是否包含p的次数大于等于k包含p的次数。由于C(n,m)=n!/(m!*(n-m)!)它们所包含的p的个数就是阶乘中每一个因子所包含的p的个数的和。而这个和可以很容易计算出来。代码如下：

/*
 * main.cpp
 *
 *  Created on: 2011-11-19
 *      Author: wwf
 *      Accepted
 *
 *      思路，求解C(n+f,n)能否被k整除，可以按照整数分解为素因子的幂的乘积的形式，来比较
 *      而C(n+f,n)中含有的素因子p的次数，则可以比较容易算出来，即(n+f)!中含p的次数，减去n和f中含有p的次数的和。
 *      从而k整除C(n+f,n)就等价于k中含有的每一个素因子p的次数都不大于C(n+f,n)中含有的p的次数。
 *
 */

#include<iostream>
#include<vector>
//#include<cmath>
#include<ctime>
using namespace std;

long factor_num(long end, long factor);

void init_numbers(void);

int solveOne(int n,int f,int g);

int brutalSolveOne(int n, int f, int g);

vector<vector<int> > numbers(100001);

int main(void) {


	init_numbers();

	long C;
	cin >> C;
	for (; C > 0; C--) {
		int n, f, g;
		cin >> n >> f >> g;
		
		int a=solveOne(n,f,g);
		cout<<a<<endl;

	}

}

long factor_num(long end, long factor) {
	long result = 0;
	while (end > 0) {
		end /= factor;
		result += end;
	}
	return result;
}

void init_numbers(void){
	for(int number=2;number<numbers.size();number++){
		if(numbers[number].size()==0){
			for(int j=number;j<numbers.size();j+=number){
				numbers[j].push_back(number);
			}
		}
	}
}

int solveOne(int n,int f,int g){
	int result=1;
	vector<int> powers(g+1);
	for(int number=2;number<=g;number++){

		//judge if number can divide C(n+f,n)
		bool candivide=true;
		{//number is not prime
			for(vector<int>::iterator it=numbers[number].begin();it!=numbers[number].end();it++){
				if(powers[(*it)]==0)powers[(*it)]=factor_num(n+f,(*it))-factor_num(n,(*it))-factor_num(f,(*it));

				int tmp=number;
				int power=0;
				while(tmp%(*it)==0){
					tmp/=(*it);
					power++;
				}
				if(power>powers[(*it)]){
					candivide=false;
					break;
				}
			}
		}
		if(candivide){
			result=number;
		}
	}
	return result;
}

/*
原始方法，只是为了测试用 

*/
int brutalSolveOne(int n, int f, int g){
	int q=1;
	int n1=1,f1=1,f2=1;
	while(true){
		if(f1<=f&&q%f1==0){
			q/=f1;
			f1++;
		}else if(f2<=n&&q%f2==0){
			q/=f2;
			f2++;
		}else if(n1<=n+f){
			q*=n1;
			n1++;
		}else{
			break;
		}
	}
	int result=1;
	for(int i=2;i<=g;i++){
		if(q%i==0)result=i;
	}
	return result;
}