TZOJ 5463: 学习之数论欧拉筛+分解质因数_c++欧拉shai分解质因数-CSDN博客

本文链接：https://blog.csdn.net/weixin_51423794/article/details/120156171

题目：

请添加图片描述

题解：

第一眼看到后就想到了分解质因数。
先把正整数K给分解质因数了。然后利用约数和他的次数，求出最大的n
但这样还不够（可能是crq在折腾负载平衡的问题），算法t了。于是我在求质因数的时候用了欧拉筛把复杂度降了下来，最后ac代码也只跑了287ms。难道欧拉分解质因数这么快的吗？

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair<LL,int>pii;
const int N=1e6+5;
int prime[N],factor[N],cnt;
void init()
{
	for(int i=2;i<N;i++)
	{
		if(!factor[i])
		{
			prime[++cnt]=i;
		}
		for(int j=1;j<=cnt&&prime[j]<N/i;j++)
		{
			factor[prime[j]*i]=prime[j];
			if(i%prime[j]==0) break;
		}
	}
}
int main()
{
	init();
	LL k;
	while(~scanf("%lld",&k))
	{
		vector<pii>vec;
		for(int i=1;prime[i]<=k/prime[i];i++)
		{
			LL p=(LL)prime[i];
			if(k%p==0)
			{
				int ans=0;
				while(k%p==0)
				{
					ans++;
					k/=p;
				}
				vec.push_back({p,ans});
			}
			
		}
		if(k>1)vec.push_back({k,1});
		LL ans=0;
		for(int j=0;j<vec.size();j++)
		{
			LL p=vec[j].first;int n=vec[j].second;
			//cout<<p<<" "<<n<<endl;
			for(LL i=p;1;i+=p)
			{
				int cnt=0;
				LL res=i,ii=i;
				while(ii%p==0)
				{
					cnt++;
					ii/=p;
				}
				//cout<<res<<" "<<cnt<<" "<<n<<endl;
				n-=cnt;
				if(n<=0)
				{
					ans=max(ans,res);
					break;
				}
			}
		}
		cout<<ans<<endl;
	}
}