题目大意
欧拉函数 φ(n) 定义为不超过正整数 n 并且与 n 互素的整数的数目。
可以证明 φ(n) = n ∗ ∏ (1 − 1 / pi). 其中 pi(1 <= i <= k)是 n 的全部素因子。
已知 y,求最小的自然数 x 使得 φ(x) = y.
多组询问。
解题思路
观察可知y的约数-1包含了所有x的质因子,所以我们可以求出x的质因子。通过变形可以用y表示x,通过选取质因子可以求出x,同时选取质因子也会构成y,我们通过枚举质因子选取情况达到同时满足x和y的限制。观察可知可能的质因子的个数不会太多,构成x的质因子的个数更少,所以可以在时限范围内通过,另外有一个结论x<=10y。
code
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LF double
#define LL long long
#define min(a,b) ((a<b)?a:b)
#define max(a,b) ((a>b)?a:b)
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
using namespace std;
int const MxCnt=1e4;
LL Ans,Cnt,Pri[MxCnt+10];
bool CPri(LL N){
int Mx=sqrt(N);
fo(i,2,Mx)if(N%i==0)return 0;
return 1;
}
void Prep(LL N){
Cnt=0;int Mx=sqrt(N);
fo(i,1,Mx)if(N%i==0){
if(CPri(i+1))Pri[++Cnt]=i+1;
if(CPri(N/i+1))Pri[++Cnt]=N/i+1;
}
}
bool Cmp(LL X,LL Y){return X>Y;}
void Dfs(int Pos,LL Y,LL X){
if(Y==1)Ans=min(Ans,X);
if((Pos>Cnt)||(Y==1)||(X>Ans))return;
for(LL i=Pri[Pos]-1;Y%i==0;i*=Pri[Pos])Dfs(Pos+1,Y/i,X/(Pri[Pos]-1)*Pri[Pos]);
Dfs(Pos+1,Y,X);
}
int main(){
freopen("d.in","r",stdin);
freopen("d.out","w",stdout);
int T;scanf("%d",&T);
fo(cas,1,T){
LL Y;scanf("%lld",&Y);
Prep(Y);
sort(Pri+1,Pri+Cnt+1,Cmp);
Ans=(Y==1)?Y:Y*10;Dfs(1,Y,Y);
printf("%lld\n",Ans);
}
return 0;
}