Codeforces Round 900 (Div. 3) F 证明+质因数分解-CSDN博客

本文链接：https://blog.csdn.net/qq_63548094/article/details/133352943

文章讲述了如何利用质因子分解和唯一分解定理证明函数d(n)对于大整数n，其值必定能整除n。通过维护质因子及其指数的集合，对n和d(n)进行质因数分解，检查n中因子的指数是否满足d(n)的指数条件来判断整除关系。

摘要由CSDN通过智能技术生成

结论：若 $a$ 存在，则 $d(n)$ 必定能整除 $n$

由此可证明当 $gcd(a,n)=1$ 时， $d(n)$ 为积性函数。

进一步推导： $d(n\cdot a)=n\rightarrow d(n)\cdot d(a)=n$

即 $d(n)$ 必定能整除 $n$ ，证明完毕。

由于 $n\cdot x$ 的值非常大，不能直接表示，考虑对因子下手。注意到 $d(n)\leq 10^{9^{^{}}}$ ，由唯一分解定理可知，一个大于1的正整数 $N$ ，如果它的标准分解式为 $N=A^{a}\cdot B^{b}...Z^{z}$ ，那么它的正因子个数，即 $d(n)=(1+a)(1+b)...(1+z)$ 。

因此，可以维护一个质因子及其指数的集合，每当乘以一个 $x$ 时，再将 $x$ 的质因子分解加入 $n$ 中，即可得到当前 $n\cdot x$ 的质因子分解式。计算 $d(n)$ 时同理，对每一个 $n$ 的分解式中每一个指数+1进行质因子分解。若 $d(n)$ 能整除 $n$ ，则对于 $n$ 分解式中的任意一个底数，其指数必须大于 $d(n)$ 分解式的指数。

//对n或x质因数分解
void divide(int n)
{
	for(int i=2;i*i<=n;i++)
		while(n%i==0)  mp[i]++,n/=i;
	if(n>1)  mp[n]++;
}
//检验d(n)能否整除n
bool check()
{
	d.clear();
	for(auto [l,r]:mp)//对指数质因数分解
	{
		r++;
		for(int i=2;i*i<=r;i++)
			while(r%i==0)  d[i]++,r/=i;
		if(r>1)  d[r]++;
	}
	for(auto [l,r]:d)
		if(mp[l]<r)  return 0;//若对于同一底数，n中指数小于d中指数，则不能整除
	return 1;
}