uva 1635 唯一分解定理

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1E5 + 10;
bool prim[maxn];
vector<pair<int, int> > Fta;
vector<int>primes, ans;
int n, m, kase, Irrelevant[maxn];
void init()
{
	for (int i = 2; i <= maxn; i++)
		if (!prim[i])
		{
			primes.push_back(i);
			for (int j = i + i; j <= maxn; j++)
				prim[j] = 1;
		}
}
void prime_factors(int n)
{
	int m = floor(sqrt(n) + 0.5);
	for (int i = 2; i  <= m; i++)
	{
		int cnt = 0, p = i;
		if (n % p == 0)
		{
			while (n % p == 0) n /= p, cnt++;
			Fta.push_back(make_pair(p, cnt));
		}
	}
	if (n > 1) Fta.push_back(make_pair(n, 1));
}
int main(int argc, char const *argv[])
{
	init();
	while (cin >> n >> m)
	{
		n--;//计算的是C(0~n-1,n-1)!;
		ans.clear();
		Fta.clear();
		prime_factors(m);
		memset(Irrelevant, 0, sizeof(Irrelevant));
		for (int i = 0; i < Fta.size(); i++)
		{
			int p = Fta[i].first, cnt = 0, x;
			for (int k = 1; k < n; k++)
			{
				x = n - k + 1;
				while (x % p == 0) {x /= p; cnt++;}
				x = k;
				while (x % p == 0) {x /= p; cnt--;}
				if (cnt < Fta[i].second) Irrelevant[k] = 1;
			}
		}
		for (int i = 1; i < n; i++) if (!Irrelevant[i]) ans.push_back(i + 1);
		cout  << ans.size() << endl;
		for (int i = 0; i < ans.size(); i++)
			cout << ans[i] << (i == ans.size() - 1 ? "\n" : " ");
		if (!ans.size()) cout << endl;//千万别忘了没有答案，也要换行！
	}
	return 0;
}

唯一分解定理的模板题，顺带附带线性打表求素数

虽然是模板题，但是解法不好想，可以利用公式C(k,n) = (n-k+1)/k*C(n,k-1) 来计算

但还是绕不过去大数（显然上不了大树） 于是就有了很有意思的观点，对于每一个m的质数考察,考察对于每一个c（），是否能被整除，对于下一个c（）由公式可得，这样上一个计算出来的指数可以被下一个继承，而本来需要高精度表示的数字，就只需要指数来表示就可以了。十分精妙