The Luckiest number （欧拉函数，问题转换，引理，快速乘）

最新推荐文章于 2022-01-20 19:11:43 发布

yezzz.

最新推荐文章于 2022-01-20 19:11:43 发布

阅读量270

点赞数 3

分类专栏：数论《算法竞赛进阶指南》好题文章标签：算法

本文链接：https://blog.csdn.net/qq_51354600/article/details/120158897

版权

好题同时被 3 个专栏收录

63 篇文章 2 订阅

订阅专栏

数论

44 篇文章 0 订阅

订阅专栏

《算法竞赛进阶指南》

12 篇文章 0 订阅

订阅专栏

本文解析了如何通过转化888连乘问题为模运算，利用欧拉函数和最大公约数技巧求解正整数LLL倍数的最小x值。重点讲解了gcd的应用以及如何找到φ(9Ld)的约数来求解。涉及算法如快速幂和Euler函数计算。

摘要由CSDN通过智能技术生成

The Luckiest number

题意：

问多少个 $8$ 连起来会是正整数 $L$ 的倍数

分析：

关键一步，巧妙转换： $x$ 个 $8$ 连在一起组成的数可以写作 $8(10^x-1)/9$
题目便是让我们求最小的x，满足： $L|\frac{8(10^x-1)}{9}$ , 设 $d = g c d (L, 8)$
$L|\frac{8(10^x-1)}{9}\Leftrightarrow9L|8(10^x-1) \Leftrightarrow\frac{9L}{d}|10^x-1 \Leftrightarrow 10^x\equiv1(mod\ \frac{9L}{d})$
- 思考：第二步到第三步的转换，为什么仅需 $g c d$ ，因为仅有约数会产生贡献
然后，问题就转换求解求解 $\varphi(\frac{9L}{d})$ 的约数
引理：若正整数 $a, n$ 互质，则满足 $a^x\equiv 1(mod\ n)$ , 的最小正整数 $x_0$ 是 $\varphi(n)$ 的约数
注意：指数比较大，快速幂里还要用快速乘，不然会炸

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

ll euler(ll n)
{
    ll ans=n;
    for(ll i=2;i*i<=n;i++)
    {
        if(n%i==0)
        {
            ans=ans/i*(i-1);
            while(n%i==0) n/=i;
        }
    }
    if(n>1) ans=ans/n*(n-1);
    return ans;
}
inline int gcd(int a,int b) { return b?gcd(b,a%b):a; }
inline ll ksc(ll a,ll b,ll p)
{
    ll res=0; a%=p; b%=p;
    while(b)
    {
        if(b&1) res=(res+a)%p;
        a=(a<<1)%p; b>>=1;
    }
	return res;
}
ll ksm(ll a,ll b,ll p)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=ksc(ans,a,p);
        a=ksc(a,a,p); b >>= 1;
    }
    return ans;
}
signed main()
{
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    ll n,cnt=0;
    while(cin>>n && n)
    {
        n=9*n/gcd(n,8);
        ll m=euler(n), ans=m;
        if(gcd(n,10)!=1) ans=0; // x与n不互质那必然为0
        else if(n==9) ans=1;
        else 
        {
            for(ll i=2;i*i<=m;i++)
            {
                if(m%i==0)
                {
                    if(ksm(10,i,n)==1) { ans=i; break; }
                    if(ksm(10,m/i,n)==1) { ans=m/i; }
                }
            }
        }
        cout<<"Case "<<++cnt<<": "<<ans<<endl;
    }

    return 0;
}