问题描述
给一个正整数N,问是否存在N的倍数M,且M的各个位全部由数字8组成,如果存在多个取最小的M
并输出M由几个8组成。
输入格式
一行,一个整数N
输出格式
一行一个整数,表示所求的结果。
如果无解,输出0
样例输入 1
8
样例输出 1
1
样例输入 2
11
样例输出 2
2
样例输入 3
16
样例输出 3
0
样例输入 4
18
样例输出 4
9
提示
1<=N<=2000000000
注意:此题可能要使用快速乘法
我们可以知道,由k个8组成的式子我们可以表示成10^(k-1)*8/9
我们设10^(k-1)*8/9=t*n
那么移项后可以得到8*10^(k-1))= 9*t*n
那么根据题意,k要 尽量小,机智的我们要想到,先求出d=gcd(8,9*n)
(8/d)*10^(k-1)=(9*n/d)*t
令p=8/d,q=(9*n)/d则有p*10^(k-1)=q*t,并且q,p互质
我们分析:ax=by中,若a,b互质那么x一定是b的倍数
因此我们得到10^(k-1)%q=0
根据同余性质我们有10k^ ≡1 (mod q)
根据欧拉定理,当gcd(10,q)=1时,我们则有(10^φ(q)) ≡1 (mod q) 当10与q互质时才有解。
因为幂运算取模是有循环节的,即一个循环周期的长度,设为s。
如果有x^a ≡ y (mod q),那么x^(a+i*s )≡ y (mod q) i>=0且为整数
比如7^2 % 10=9 , 7^6 %10=9 , 7^10%10=9 指数循环节s=4
对于(10^k) ≡1 (mod q),设指数循环节为s
首先k=0是一个解:10^0 ≡ 1 (mod q), 则有10^(0+i*s) ≡ 1 (mod q)
也就是说10^(i*s) ≡ 1 (mod q)
而且已经确定φ(q)是一个解了,所以φ(q)一定是这个循环节s的倍数。
所以只需枚举φ(q)的因数k,找出其中最小的满足(10^k) ≡1 (mod q)的即可
#include<iostream>
#include<cstdio>
#define LL long long
using namespace std;
LL q,d;
LL gcd(LL n,LL m){
if(m==0)return n;
return gcd(m,n%m);
}
LL f(LL x){
LL i,ans=x,a=x;
for(i=2;i*i<=a;i++){
if(a%i==0){
ans=ans-ans/i;
while(a%i==0)a=a/i;
}
}
if(a>1)ans=ans-ans/a;
return ans;
}
LL ksc(LL x,LL y,LL N)
{
LL a;
for(a=0;y;x=(x<<1)%N,y=(y>>1))
if(y&1)a=(a+x)%N;
return a;
}//快速乘法函数(主要应用在二分快速幂中,防止答案过大而溢出)
LL g(LL a,LL b,LL c){
LL ans=1;
a%=c;
while(b>0){
if(b&1)ans=ksc(ans,a,c);
b=b>>1;
a=ksc(a,a,c);
}
return ans;
}//二分快速幂
bool check(LL s){
if(g(10,s,q)==1)return true;
return false;
}
int main(){
LL n,i,ans=999999999;
cin>>n;
if(n==1||n==2||n==4||n==8){
puts("1");
return 0;
}//由于算法的不完美性,因此特殊情况要特别判断
d=gcd(8,9*n);
q=9*n/d;
if(gcd(q,10)>1){
puts("0");
return 0;
}//无解的情况
LL phi=f(q);
for(i=2;i*i<=phi;i++)//不能从i=2枚举到i=phi,否则在最坏情况下会超时
if(phi%i==0){
if(check(i))
ans=min(ans,i);
if(check(phi/i))
ans=min(ans,phi/i);
}
cout<<ans;
}