BJOI2014 Euler 一道简单的数论暴力

最新推荐文章于 2022-02-08 22:41:06 发布

YxuanwKeith

最新推荐文章于 2022-02-08 22:41:06 发布

阅读量1.2k

点赞数 2

分类专栏：算法-数论文章标签： phi

本文链接：https://blog.csdn.net/yxuanwkeith/article/details/51089843

版权

算法-数论专栏收录该内容

11 篇文章 1 订阅

订阅专栏

题目大意

给你一个数 $y$ 要求输出一个满足 $\varphi(x) = y$ 的最小的 $x$ ,有 $T$ 组询问。

$y <= 10^{12}$
$1 <= T <= 2$

解题思路

求 $\varphi(x)$ 时有一个众所周知的公式是 $\varphi(x) = x * \Pi\frac{p_i - 1}{p_i}$ （ $p_i$ 为 $x$ 的所有质因子）,也就是说这我们可以得到等式
$~~~~~~~~~\varphi(x) = y$
$x * \Pi\frac{p_i - 1}{p_i} = y$ （等式一）
$~~~~~~~~~~~~~~~x = y * \Pi\frac{p_i}{p_i - 1}$ （等式二）
由等式一可知，因为 $p_i$ 都是 $x$ 的质因子，所以我们可以把下面的 $p_i$ 都与 $x$ 约掉，然后我们就发现 $x$ 的质因子要不就是 $y$ 的质因子，要不就是 $y$ 的约数加一。而由等式二可知我们只需要找出 $\Pi\frac{p_i}{p_i - 1}$ 最小的符合要求的解就行了。不难发现符合要求的 x <script type="math/tex" id="MathJax-Element-557">x</script>可能的质因子只有1000多个，那么我们只需暴力判断每个质数选不选再加个全局最优解的优化就可以了。

程序

#include <cstring>
#include <cstdio>
#include <algorithm>

using namespace std;

const int MAXN = 7e3;
typedef long long LL;

LL N, Ans, Fac[MAXN], Pri[MAXN];
int tot, cnt;

void Prepare(LL N) {
    for (LL i = 1; i * i < N; i ++) {
        if (N % i != 0) continue;
        Fac[++ tot] = i;
        if (1ll * i * i == N) continue;
        Fac[++ tot] = N / i;
    }
}

bool IsPrime(LL Now) {
    for (LL i = 2; i * i <= Now; i ++) 
        if (Now % i == 0) return 0;
    return 1;
}

void GetPri() {
    for (int i = 1; i <= tot; i ++) 
        if (IsPrime(Fac[i] + 1)) Pri[++ cnt] = Fac[i] + 1;
}

bool cmp(LL a, LL b) { return a > b;}

bool Check(LL Num, LL Now) {
    if (Now == 1) return 1;
    if (Now % Pri[Num] == 0) return Check(Num, Now / Pri[Num]);
    return 0;
}

void Dfs(int Num, LL Now, LL S) {
    if (Num > cnt || S > Ans || Now == 1) return;
    if (Now % (Pri[Num] - 1) != 0) {
        Dfs(Num + 1, Now, S);
        return;
    }
    if (Now % (Pri[Num] - 1) == 0 && Check(Num, Now / (Pri[Num] - 1))) 
        Ans = min(Ans, S / (Pri[Num] - 1) * Pri[Num]);
    LL Ord = (Pri[Num] - 1);
    for (; Now % Ord == 0; Ord *= Pri[Num]) 
        Dfs(Num + 1, Now / Ord, S / (Pri[Num] - 1) * Pri[Num]);
    Dfs(Num + 1, Now, S);
}

void Solve(LL N) {
    Ans = (N == 1) ? N : N * 10;
    tot = cnt = 0;
    Prepare(N);
    GetPri();
    sort(Pri + 1, Pri + 1 + cnt, cmp);
    Dfs(1, N, N);
}

int main () {
    int Test;
    scanf("%d", &Test);
    for (; Test; Test --) {
        scanf("%lld", &N);
        Solve(N);
        printf("%lld\n", Ans);
    }
}