数论(同余方程) - The Luckiest number - POJ 3696

数论(同余方程) - The Luckiest number - POJ 3696

题意：

8是中国的幸运数字，如果一个数字的每一位都由8构成则该数字被称作是幸运数字。

现在给定一个正整数L，请问至少多少个8连在一起组成的正整数（即最小幸运数字）是L的倍数。

输入格式

输入包含多组测试用例。

每组测试用例占一行，包含一个整数L。

当输入用例L=0时，表示输入终止，该用例无需处理。

输出格式

每组测试用例输出结果占一行。

结果为“Case 1: ”+一个整数N，N代表满足条件的最小幸运数字的位数。

如果满足条件的幸运数字不存在，则N=0。

数据范围

$1\le L\le 2\times 10^9$

输入样例：

8
11
16
0

输出样例：

Case 1: 1
Case 2: 2
Case 3: 0

---

分析：

长度为x的，由同一个数字组成的数，如何表示成关于x的函数形式？

$以x个8为例：$

$\underset{x个8}{\underbrace{888...888}}=8\times \underset{x个1}{\underbrace{111...111}}=8\times \frac{999...999}{9}=\frac{8}{9}(10^x-1)。$

$本题要求888...888\equiv 0(modL)即\frac{8}{9}(10^x-1)\equiv 0(modL)的关于x的最小正整数解，$

$先求一组解(x_0,y_0)，$

$<=>8(10^x-1)\equiv 0(mod9L)，$

$<=>10^x-1\equiv 0(mod\frac{9L}{gcd(L,8)})，(因为9和8互质)$

$<=>10^x\equiv 1(mod\frac{9L}{gcd(L,8)})$

$记C=\frac{9L}{gcd(L,8)}，由于x与C不一定互质，故不满足用费马小定理使用的条件，$

$应用欧拉定理，得到可行解\varphi (C)，满足：$

$$10^{\varphi (C)}\equiv 1(modC)$$

$因为我们要求最小的正整数解，$

$可以证明：\varphi (C)的因子同样也是上式的解，$

$故我们还需枚举\varphi (C)的因子，逐个判断是否满足条件，输出最小的解。$

判断无解：

$10^x\equiv 1(modC)<=>10^x+Cy=1$

$gcd(10,C)\ne 1，则等式两边同除gcd(10,C)，发现等式右边不是整数，此时无解。$

具体步骤：

$①、计算C=\frac{9L}{gcd(L,8)}，判断是否有解。$

$②、若有解，根据定义计算C=P_1^{a_1}{P}_2^{a_2}...P_k^{a_k}的欧拉函数：$

$$\varphi (C)=C(1-\frac{1}{P_1})(1-\frac{1}{P_2})...(1-\frac{1}{P_k}),$$

$③、接着枚举\varphi (C)的所有因子i和\varphi (C)/i，对每一个质因子p，判断10^p\equiv 1(modC)是否成立，$

$输出满足条件的最小的因子p。$

注意：

$本题虽然考察同余方程，但并未使用到扩展欧几里得算法。$

$使用快速幂计算10^p对C取模时，注意乘法的过程可能会爆longlong，$

$因此，我们用快速(龟速)乘法来进行每一次乘法操作。$

代码：

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

#define ll long long

using namespace std;

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

ll quick_pow(ll a,ll b,ll mod)
{
    ll res=1;
    while(b)
    {
        if(b&1) res=quick_mul(res,a,mod);
        a=quick_mul(a,a,mod);
        b>>=1;
    }
    return res;
}

ll gcd(ll a,ll b)
{
    return b ? gcd(b,a%b) : a;
}

ll get_euler(ll C)      //定义求欧拉函数
{
    ll res=C;
    for(int i=2;i<=C/i;i++)
        if(C%i==0)
        {
            while(C%i==0) C/=i;
            res=res/i*(i-1);
        }
            
    if(C>1) res=res/C*(C-1);
    return res;
}

int main()
{
    ll L;
    int T=1;
    
    while(cin>>L,L)
    {
        ll C=L/gcd(L,8)*9;
        
        ll res=1e18;
        if(gcd(10,C)!=1) res=0;
        else
        {
            ll phi_C=get_euler(C);
            
            for(ll i=1;i<=phi_C/i;i++)
                if(phi_C%i==0)
                {
                    if(quick_pow(10,i,C)==1) res=min(res,i);
                    if(quick_pow(10,phi_C/i,C)==1) res=min(res,phi_C/i);
                }
        }
        
        printf("Case %d: %lld\n",T++,res);
    }
    
    return 0;
}