bzoj 3643: Phi的反函数

最新推荐文章于 2019-07-04 22:58:56 发布

转载最新推荐文章于 2019-07-04 22:58:56 发布 · 97 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/WQHui/p/7536850.html

本文提供了一种求解Phi反函数的算法实现，通过深度优先搜索结合质数筛法来寻找最小的x使得φ(x)等于给定的数值。代码中详细展示了初始化质数表、检查质数、递归搜索等关键步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

3643: Phi的反函数

Time Limit: 10 Sec Memory Limit: 64 MB

Description

Input

Output

Sample Input

Sample Output

HINT

Source

Tips：

　　自己也没怎么搞懂；

Code：

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define MAXN 200008
using namespace std;
long long n,m,prime[MAXN],tot,flag[MAXN],ans;

void init(){
    for(long long i=2;i<=MAXN-8;i++){
        if(!flag[i]){
            flag[i]=1;
            prime[++tot]=i;
        }
        for(long long j=1;j<=tot&&i*prime[j]<=MAXN-8;j++){
            flag[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}

bool check(long long u){
    for(int i=1;prime[i]*prime[i]<=u;i++){
        if(u%prime[i]==0) return 0;
    }
    return 1;
}

void dfs(long long x,long long y,long long z){
    if(z>=ans) return;
    if(y==1){
        ans=min(ans,z);
        return;
    }
    if(y*y>n&&check(y+1)){
        ans=min(ans,z*(y+1));
        return;
    }
    for(int i=x;prime[i]-1<=y&&(prime[i]-1)*(prime[i]-1)<=n;i++)
    if(y%(prime[i]-1)==0){
        long long xx=y    /(prime[i]-1);
        long long yy=z*prime[i];
        dfs(i+1,xx,yy);
        while(xx%prime[i]==0){
            xx=xx/prime[i];
            yy=yy*prime[i];
            dfs(i+1,xx,yy);
        }
    }
}

int main(){
    tot=0;
    ans=(long long)1<<31;
    init();
    scanf("%lld",&n);
    dfs(1,n,1);
    if(ans==(long long)1<<31){
        printf("-1");
    }else{
        printf("%lld",ans);
    }
}

转载于:https://www.cnblogs.com/WQHui/p/7536850.html