【jzoj3760】【BJOI2014】【Euler】

题目大意

欧拉函数 φ(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;
}