牛客网 KY105 整除问题

参考下述解法：
1.素数筛
2.对n!质因子分解算法
3.对a质因子分解算法
int nFac[MAXN]={0};//记录n!的质因数分解后的结果，n!存在质因数i，则aFac[i]=指数
int aFac[MAXN]={0};//记录a的质因数分解的结果，a存在质因数i，则aFac[i]=指数
比较n！和a的相同质因数的指数相除的最小值为所求。
例如:
6!=2⁴⁺³2+5^1
10=2^1*51
有公共质因数：
2 指数相除为4/1=4
5 指数相除为1/1=1
所以答案为1

这里是将a质因子分解后在n！的质因子中逐个将其消去，每成功消去a的全部质因子后就代表可以除一个a，就能得到k的值。

#include <iostream>
#include <cstring>
#include <vector>
#include <cmath>
using namespace std;

int num[1001];											//存放各因子的个数

int main()
{
	int n, a;
	while (cin >> n >> a) {
		vector<int> yin_zi;
		int k = 0;

		memset(num, 0, 1001);

		for (int i = 2; i <= n; ++i) {
			int t = i;
			for (int j = 2; j <= sqrt(t); ++j) {			//所求阶乘的全部因子
				while (t % j == 0) {
					t /= j;
					num[j]++;
				}

				if (t == 1) break;
			}
			if (t > 1) num[t]++;
		}

		for (int j = 2; j <= sqrt(a); ++j) {			//所给数的全部因子
			while (a % j == 0) {
				a /= j;
				yin_zi.push_back(j);
			}

			if (a == 1) break;
		}
		if (a > 1) yin_zi.push_back(a);

		bool test = true;
		while (test) {
			for (int i = 0; i < yin_zi.size(); ++i) {
				num[yin_zi[i]]--;							//对应因子消去
				if (num[yin_zi[i]] < 0) test = false;
			}
			if (!test) break;
			k++;
		}

		cout << k << endl;

	}

}