X-因子链(factor)

$$\mathrm{X}-因子链\left(\mathrm{f}\mathrm{a}\mathrm{c}\mathrm{t}\mathrm{o}\mathrm{r}\right)$$

题目链接：$\mathrm{S}\mathrm{S}\mathrm{L}比赛1502$

题目

给一个正整数 $X$ ，一个长度为 $m$ 的X-因子链是指这样一个序列： $X_0=1$ ， $X_1$ ， $X_2$ ， $。。。$ ， $X_m=X$ 满足： $X_i<X_i+1$ 同时 $X_i|X_i+1$ ( $X_i+1$ 能被 $X_i$ 整除)
要求X-因子链的最大长度 $Len$ 和长度为 $Len$ 的X-因子链的数量。

输入

一个正整数 $X$

输出

一行，两个整数，分别表示最大长度和该长度链的种数。

样例输入

100

样例输出

4 6

数据范围

对于 $20\%$ 的数据： $X<=20,000;$
对于 $100\%$ 的数据： $X<=2^{31};$ 且保证因子链最大长度小于等于 $20;$

思路

这道题是一道数学题。

我们很容易可以发现，其实因子链的长度就是质因子的数量。
那接着就看方案数了，那其实就是一个全排列问题，可以得出公式：
$$ans=\frac{num!}{a_1!\ast a_2!\ast ...\ast a_x!}$$
（ $num$ 为质因子的数量， $a_i$ 为第 $i$ 中质因子有多少个）

那就暴力算出来，就完事了。

代码

#include<cstdio>
#define ll long long

using namespace std;

ll x, su[50001], num, zhi[50001], ans = 1;
bool susu[50001];

int main() {
	susu[0] = susu[1] = 1;//求出素数
	for (ll i = 2; i <= 50000; i++)
		if (!susu[i]) {
			su[++su[0]] = i;
			for (ll j = i + i; j <= 50000; j += i)
				susu[j] = 1;
		}
	
	scanf("%lld", &x);//读入
	for (ll i = 1; i <= su[0]; i++)
		if (x % su[i] == 0) {
			++zhi[0];
			while (x % su[i] == 0) {//分解质因数
				x /= su[i];
				num++;
				zhi[zhi[0]]++;
			}
		}
	
	printf("%lld ", num);//输出
	
	for (ll i = 2; i <= num; i++)//求出组合方案
		ans *= i;
	for (ll i = 1; i <= zhi[0]; i++) {
		for (int j = 2; j <= zhi[i]; j++)
			ans /= j; 
	}
	
	printf("%lld", ans);//输出
	
	return 0;
}