The Luckiest Number（POJ 3696，模板）

原题题面：The Luckiest Number

思路：

这题是非常好的一个模板题，考察了欧拉定理，快乘，快速幂，欧拉函数，筛质因数，
$x$ 个 $8$ 连在一起组成的正整数等价于 $\frac{8(10^x-1)}{9}$ 。其中 $x$ 就是该正整数的十进制位数。而题目就是让我们求出这样一个 $x$，满足 $L|\frac{8(10^x-1)}{9}$。
在这里，我们令 $d=gcd(L,8)$。
推导可以得出：$$L|\frac{8(10^x-1)}{9}\iff 9L|8(10^x-1)\iff \frac{9L}{d}|\frac{8}{d}(10^x-1)\iff \frac{9L}{d}|(10^x-1)\iff 10^x\equiv 1(mod\frac{9L}{d})$$

欧拉定理：
若正整数 $a,n$ 互质，则 $a^{\varphi (n)}\equiv 1(modn)$，其中 $\varphi (n)$ 为 $n$ 的欧拉函数。
由欧拉定理可以看出，$a=10$，$n=\frac{9L}{d}$，而 $x$ 集合中的最小正整数 $x_0$ 是 $\varphi (n)$ 的一个约数。
所以，我们只要依次求出欧拉函数 $\varphi (\frac{9L}{d})$，然后枚举它的所有约数，再用快速幂检查当前约数是否满足条件即可。

代码如下：

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

using namespace std;

typedef long long LL;
typedef pair<LL, int> PII;
const int N = 45000, M = 50;

LL primes[N];
int cnt;
bool st[N];

PII factor[M];
LL divider[N];
int cnt_f, cnt_d;

LL gcd(LL a, LL b)
{
    return b ? gcd(b, a % b) : a;
}
//线性筛法筛质数板子
void get_primes(LL n)
{
    for(int i = 2; i <= n; i++)
    {
        if(!st[i])  primes[cnt ++] = i;
        for(int j = 0; primes[j] <= n / i; j++)
        {
            st[i * primes[j]] = true;
            if(i % primes[j] == 0)  break;
        }
    }
}
//求欧拉函数板子
LL get_phi(LL n)
{
    LL res = n;
    for(int i = 0; i < cnt && primes[i] <= n / primes[i]; i++)
    {
        LL p = primes[i];
        if(n % p == 0)
        {
            res = res / p * (p - 1);
            while(n % p == 0)   n /= p;
        }
    }
    if(n > 1)   res = res / n * (n - 1);
    return res;
}
//筛质因数板子
void get_fac(LL n)
{
    cnt_f = 0;
    for(int i = 0; i < cnt && primes[i] <= n / primes[i]; i++)
    {
        LL p = primes[i];
        if(n % p == 0)
        {
            int s = 0;
            while(n % p == 0)
            {
                n /= p;
                s ++;
            }
            factor[++ cnt_f] = {p, s};
        }
    }
    if(n > 1)   factor[++ cnt_f] = {n, 1};
}
//深搜枚举约数板子
void dfs(int u, LL p)
{
    if(u > cnt_f)
    {
        divider[cnt_d ++] = p;
        return ;
    }
    for(int i = 0; i <= factor[u].second; i++)
    {
        dfs(u + 1, p);
        p *= factor[u].first;
    }
}
//快乘板子
LL qk_mul(LL a, LL b, LL p)
{
    LL ans = 0;
    while(b)
    {
        if(b & 1)   ans = (ans + a) % p;
        a = (a * 2) % p;
        b >>= 1;
    }
    return ans;
}
//快速幂板子
bool qmi(LL a, LL b, LL mod)
{
    LL ans = 1 % mod;
    a %= mod;
    while(b)
    {
        if(b & 1)   ans = qk_mul(ans, a, mod);
        a = qk_mul(a, a, mod);
        b >>= 1;
    }
    if(ans == 1)    return true;
    return false;
}

LL solve(LL l)
{
    LL n = l / gcd(8LL, l) * 9;
    if(gcd(10LL, n) != 1)   return 0;

    LL phi = get_phi(n);
    get_fac(phi);

    cnt_d = 0;
    dfs(1, 1);

    for(int i = 0; i < cnt_d; i++)
    {
        LL x = divider[i];
        if(qmi(10, x, n))
            return x;
    }
    return 0;
}

int main()
{
    ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);

    LL l;
    int Count = 1;
    get_primes(N);
    while(cin >> l)
    {
        if(l == 0)  break;
        cout << "Case " << Count++ << ": " << solve(l) << endl;
    }

    return 0;
}