【BZOJ-3643】Phi的反函数 数论 + 搜索

3643: Phi的反函数

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 141  Solved: 96
[Submit][Status][Discuss]

Description

Input

Output

Sample Input

4

Sample Output

5

HINT

Source

By Zky

Solution

首先答案和N一定是同阶的，所以，可以很暴力的线筛扫一遍求解。

然后根据欧拉函数的式子，我们实际上是可以爆搜的。

爆搜他的质因子然后去凑答案，加最优性剪枝就可以跑过。

最关键的是依据欧拉函数的定义式找到规律！

Code

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
#define LL long long
inline int read()
{
    int x=,f=; char ch=getchar();
    while (ch<'' || ch>'') {if (ch=='-') f=-; ch=getchar();}
    while (ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}
    return x*f;
}
#define MAXN 5100000
int X,N;
int prime[MAXN],flag[MAXN],cnt;
void Getprime()
{
    flag[]=; cnt=;
    for (int i=; i<=N; i++)
        {
            if (!flag[i]) prime[++cnt]=i;
            for (int j=; j<=cnt && prime[j]*i<=N; j++)
                {
                    flag[i*prime[j]]=;
                    if (prime[j]%i==) break;
                }
        }
}
LL ans=(1LL<<);
LL sqr(LL x) {return (LL)x*x;}
bool check(int x)
{
//  if (flag[x]) return 0;
    for (int i=; sqr(prime[i])<=x && i<=cnt; i++) if (x%prime[i]==) return ;
    return ;
}
void DFS(int dep,LL sum,LL x,int last)
{
    if (sum>=ans) return;
    if (x==) {ans=sum; return;}
    if (check(x+) && sqr(x)>X) ans=min(ans,sum*(x+));
    for (int i=last+; (prime[i]-)<=x && sqr(prime[i]-)<=X; i++)
        if (!(x%(prime[i]-)))
            {
                LL xx=(LL)x/(prime[i]-),summ=(LL)sum*prime[i];
                DFS(dep+,summ,xx,i);
                while (!(xx%prime[i])) xx=(LL)xx/prime[i],summ=(LL)summ*prime[i],DFS(dep+,summ,xx,i);
            }
}
int main()
{
    X=read();  N=int(sqrt(X))+;
    if (X==) {puts(""); return ;}
    Getprime();
    DFS(,1LL,(LL)X,);
    printf("%lld\n",ans==(1LL<<)? -:ans);
    return ;
}

仔细思考一下应该是可以想到的。