BZOJ 3643|Phi的反函数|搜索|线性筛法

求

$$\varphi^{-1}(n)=\min\{x|\varphi(x)=n\} (n\leq 2^{31})$$

太大就输出-1的一定是搜索题

$$\begin{aligned} \varphi(n)&=n(1-\frac{1}{p_1})\cdots(1-\frac{1}{p_k})\\&=p_1^{a_1}\cdots p_k^{a_k}(1-\frac{1}{p_1})\cdots(1-\frac{1}{p_k})\\&=p_1^{a_1-1}(p_1-1)\cdots p_k^{a_k-1}(p_k-1)\end{aligned}$$

因此$(p_i-1)$是$\varphi(n)$的因数。
可以证明到$k$不会超过10，因此爆搜即可。
大于$\sqrt n$的因数要特别判断。
考虑枚举质因数，如果其-1是因数，那么就可以接受。
然后枚举质因数的幂。

#include <cstdio>
typedef long long ll;
int N, m, p[100005];
bool vis[100005]; ll ans;
void getPrimes(int n) {
    int i, j;
    for (i = 2; i <= n; ++i) {
        if (!vis[i]) p[m++] = i;
        for (j = 1; j < m && i * p[j] <= n; ++j) {
            vis[i * p[j]] = 1;
            if (i % p[j] == 0) break;
        }
    }
}
bool isPrime(int x) {
    for (int i = 1; p[i] * p[i] <= x; ++i)
        if (x % p[i] == 0) return 0;
    return 1;
}
void dfs(int k, int n, ll s) {
    if (ans != -1 && s >= ans) return;
    if (n == 1) { ans = s; return; }
    if (1ll * n * n > N && isPrime(n + 1))
        if (ans == -1 || ans > s * (n + 1))
            ans = s * (n + 1);
    for (int i = k; p[i] - 1 <= n && 1ll * (p[i] - 1) * (p[i] - 1) <= N; ++i)
        if (n % (p[i] - 1) == 0) {
            int x = n / (p[i] - 1);
            ll y = s * p[i];
            dfs(i + 1, x, y);
            while (x % p[i] == 0) {
                x /= p[i]; y *= p[i];
                dfs(i + 1, x, y);
            }
        }
}
int main() {
    scanf("%d", &N);
    getPrimes(50000);
    ans = -1; dfs(0, N, 1);
    printf("%lld\n", ans);
    return 0;
}