《算法竞赛进阶指南》/ poj3696 最幸运的数字 · 欧拉定理 + 同余 + 指数循环节

题解

参考这篇题解 - poj3696 快速幂的优化+欧拉函数+gcd的优化+互质

由8构成的数可以写成 $\frac{10^n-1}{9}\ast 8$，

所以题目可以改成求

首先，要知道一个性质：

当 $ax=by$ 且 $(a,b)$ 互质时，说明a的因子里没有b，那么要使等式成立，b一定在x中，即x可以乘除b

然后设 $d=gcd(8,L)$ _{我觉得gcd(8,9L)没准也可以}，则有：

此时，$\frac{8}{d}$ 与 $\frac{L}{d}$ 互质，且8里面没有9的因子，所以$\frac{8}{d}$ 与 $\frac{9L}{d}$ 也是互质

因此可以得出 $10^n-1$ 一定是 $\frac{9L}{d}$ 的倍数

即，

现在有没有觉得这个公式是不是有点像欧拉定理：

欧拉定理告诉我们 $a^{\phi (m)}$ 一定是一个满足条件的解，把公式带入一下就变成

$\phi (\frac{9L}{d})$ 是一个满足条件的解，但是不一定是题目要求的最小的解，

已知当n=0时，是一个满足条件的解，所以最小的解还可能存在 0~ $\phi (\frac{9L}{d})$ 之间，

并且，因为有模数，而余数有是有循环节的，所以一定存在一个循环r，使得 $a^r=a^{m+r}=1(modm)$，

因此，这个最小解一定是 $\phi (\frac{9L}{d})$ 的因子，

那么，遍历一遍 $\phi (\frac{9L}{d})$ 的因子就可以了

一些注意事项：

• 由于$L$可能很大，在快速幂的过程中可能会爆范围（一定的），因此必须采用快速幂+快速乘
• 拆因子的过程中，o(n)的方法会T，一定要 $\sqrt{n}$
• 本代码能过acwing，但是poj上没试过

---

---

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll L, mod;

ll phi(ll n) {
    ll res = n;
    for (ll i = 2; i * i <= n; ++i) {
        if (n % i == 0) {
            res = res / i * (i - 1);
            while (n % i == 0)n /= i;
        }
    }
    if (n > 1)res = res / n * (n - 1);
    return res;
}

ll multi(ll a, ll b) {
    ll res = 0;
    while (b) {
        if (b & 1) res = (res + a) % mod;
        a = (a * 2) % mod;//*2
        b >>= 1;
    }
    return res;
}

ll qpow(ll a, ll b) {
    ll res = 1;
    while (b) {
        if (b & 1)res = multi(res, a);
        b >>= 1;
        a = multi(a, a);
    }
    return res;
}


int main() {
    ios::sync_with_stdio(0);

    int cs = 0;

    while (cin >> L && L) {
        printf("Case %d: ", ++cs);
        ll d = __gcd(8ll, L);
        mod = L / d * 9;
        ll n = phi(mod);

        ll ans = 1e18, a;
        for (ll i = 1; i * i <= n; ++i) {
            if (n % i == 0) {
                a = qpow(10, i);
                if (a == 1) {
                    ans = i;
                    break;
                }
                a = qpow(10, n / i);
                if (a == 1) ans = min(ans, n / i);
            }
        }
        if (ans == 1e18) puts("0");
        else printf("%lld\n", ans);
    }
    return 0;
}