P4780-Phi的反函数【dfs】

1|0正题

题目链接:https://www.luogu.com.cn/problem/P4780

---

1|1题目大意

给出$n$，求一个最小的$x$满足$\varphi(x)=n$。
若不存在或者大于$2^{31}$则输出$-1$。

$1\leq n\leq 2^{31}$

---

1|2解题思路

考虑用$\varphi$比较常用的公式，把$n$拆成若干个$\prod (p_i-1)*p_i^{c_i}$的形式。因为这个不会超过$log$个所以可以暴力搜索比较小的质数，然后直到$n$剩下一个$p_i+1$时或$1$时再暴力判断。

---

1|3code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
using namespace std;
const ll N=46360;
ll n,ans,cnt,pri[N/10];
bool v[N];
void Prime(){
	for(ll i=2;i<N;i++){
		if(!v[i])v[i]=1,pri[++cnt]=i;
		for(ll j=1;j<=cnt&&i*pri[j]<N;j++){
			v[i*pri[j]]=1;
			if(i%pri[j]==0)break; 
		}
	}
	return;
}
bool IsPri(ll x){
	if(x%2==0)return 0;
	for(ll i=3;i*i<=x;i+=2)
		if(x%i==0)return 0;
	return 1;
}
void dfs(ll phi,ll x,ll k){
	if(phi>(1ll<<31))return;
	if(x==1){
		ans=min(ans,phi);
		return;
	}
	if(x>sqrt(n)&&IsPri(x+1))
		ans=min(ans,phi*(x+1));
	if(pri[k]>x)return;
	for(ll i=k;i<=cnt;i++){
		if(x%(pri[i]-1)==0){
			ll z=x/(pri[i]-1),p=phi*pri[i];
			dfs(p,z,i+1); 
			while(z%pri[i]==0){
				p*=pri[i];z/=pri[i];
				dfs(p,z,i+1);
			}
		}
	}
	return;
}
signed main()
{
	scanf("%lld",&n);
	Prime();ans=(1ll<<32);
	dfs(1,n,1);
	if(ans==(1ll<<32))puts("-1");
	else printf("%lld\n",ans);
	return 0;
}

__EOF__

本文作者：QuantAsk
本文链接：https://www.cnblogs.com/QuantAsk/p/15114327.html
关于博主：退役OIer，GD划水选手
版权声明：本博客所有文章除特别声明外，均采用 BY-NC-SA 许可协议。转载请注明出处！
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！