2018.09.28 牛客网contest/197/A因子（唯一分解定理）

SC.ldxcaicai

于 2018-09-28 23:34:03 发布

阅读量155

点赞数 1

CC 4.0 BY-SA版权

分类专栏： # 唯一分解定理 # 数论文章标签：数论

本文链接：https://blog.csdn.net/dreaming__ldx/article/details/82891563

数论同时被 2 个专栏收录

28 篇文章

订阅专栏

唯一分解定理

1 篇文章

订阅专栏

本文介绍了一种在算法竞赛中应用唯一分解定理解决特定数学问题的方法。通过计算质因子在n!中出现的次数，找到最小比率以求解问题。使用C++实现并详细解释了代码逻辑。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门
比赛的时候由于变量名打错了调了很久啊。
这道题显然是唯一分解定理的应用。
我们令 $P=a_1^{p_1}*a_2^{p_2}*...*a_k^{p_k}$ 。
那么对于每一个质因子 $a_i$ 我们都可以计算出它在 $n!$ 中出现 $b_i$ 次，每次计算是 $O(log_{a_i} n)$ 的。
由于对于每个 $P$ 每个 $a_i$ 会出现 $p_i$ 次因此我们只需要取 $m i n$ { $b_i/k_i$ }就行了。
代码：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll n,p,ans=1e18;
inline ll calc(ll tmp){
	ll mul=tmp,cnt=0;
	while(mul<=n&&mul>0){ 
		cnt+=n/mul;
		mul*=tmp;
	}
	return cnt;
}
ll pri[100005],vis[10005],tot=0;
inline void init(){
	for(int i=2;i<=p;++i){
		if(!vis[i])pri[++tot]=i;
		for(int j=1;j<=tot;++j){
			int k=pri[j]*i;
			if(k>p)break;
			vis[k]=1;
			if(i%pri[j]==0)break;
		}
	}
}
int main(){
	cin>>n>>p,init();
	ll len=sqrt(p);
	for(ll i=1;pri[i]<=len;++i){
		if(p%pri[i])continue;
		ll siz=0;
		while(p%pri[i]==0)++siz,p/=pri[i];
		ans=min(ans,calc(pri[i])/siz);
	}
	if(p!=1)ans=min(ans,calc(p));
	cout<<ans;
	return 0;
}