NKOI 3685 8数

P3685  8数

时间限制 :  - MS    空间限制 :  65536 KB  评测说明 : 时限：1000ms

问题描述

给一个正整数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;
			
}