[BZOJ3643]Phi的反函数（数学相关+搜索）

题目描述

题解

题目还是比较良心哒！给出了$\varphi$的公式。
把它拆开来看，可以知道$x=\varphi(n)={p_1}^{k_1}{p_2}^{k_2}......{p_m}^{k_m}(p_m-1)(p_m-1)......(p_m-1)$，$n={p_1}^{k_1}{p_2}^{k_2}......{p_m}^{k_m}$
我们现在知道了x，问如何求n。
看一下数据范围，质因子的种类数不会超过10。可以考虑爆搜。
先筛出所有的质数，然后用这些质数暴力凑x，每一次更新n的答案。
最后有可能会剩一个大质数，要注意特判这种情况。

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define N 100005
#define LL long long

const LL INF=1e18;
int x,m; LL ans;
int phi[N],p[N],prime[N];

void get_phi()
{
    phi[1]=1;
    for (int i=2;i<=N;++i)
    {
        if (!p[i])
        {
            prime[++prime[0]]=i;
            phi[i]=i-1;
        }
        for (int j=1;j<=prime[0]&&i*prime[j]<=N;++j)
        {
            p[i*prime[j]]=1;
            if (i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
}
bool check(LL x)
{
    for (LL i=2;i*i<=x;++i)
        if (x%i==0) return false;
    return true;
}
void dfs(int st,int now,LL n)
{
    if (n>=ans) return;
    if (now==1)
    {
        ans=min(ans,n);
        return;
    }
    if (now>m&&check(now+1))
    {
        n*=((LL)now+1);
        ans=min(ans,n);
        return;
    }
    int i;
    for (i=st+1;i<=prime[0]&&prime[i]-1<=m;++i)
    {
        if (prime[i]-1>now) break;
        if (now%(prime[i]-1)==0)
        {
            int x=now/(prime[i]-1); LL y=n*(LL)prime[i];
            dfs(i,x,y);
            while (x%prime[i]==0)
            {
                x/=prime[i]; y*=(LL)prime[i];
                dfs(i,x,y);
            }
        }
    }
} 
int main()
{
    get_phi();
    scanf("%d",&x);
    m=(int)sqrt(x);
    ans=INF;
    dfs(1,x,1);
    if (ans>2147483647) puts("-1");
    else printf("%lld\n",ans);
}