AcWing 202. 最幸运的数字 (欧拉 + 快速幂+快速乘)

题意

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

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

输入格式

输入包含多组测试用例。

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

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

输出格式

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

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

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

数据范围

1 < L < 2*10^9

输入样例：

8
11
16
0

输出样例：

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

题解

问最少几个8组成的数字能整除L,那么把88888 转换成 (10^x-1) / 9 * 8 ,那么就是 8(10 ^ x - 1) % 9L == 0 设d = gcd(8,L), 那么(10 ^ x - 1) % (9 / d * L) == 0 ,设 9 / d * L = C ,C已知,那么就是10^x mod C = 0,由欧拉定理可知,a ^ phi(n) % n = 0 ,所以x 可以取 phi© , 引理是 Min(x) 必须是phi©的因子,那么只需要把phi©的因子中找到最小的满足解就可以了.
注意,这个题的数据有两处会超ll,所以要用龟速乘防止溢出

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


ll qpw(ll a,ll k,ll b){  //快速幂 
	ll res = 1;
	while(k){
		if(k & 1)res = qmul(res,a,b);
		a = qmul(a,a,b);
		k >>= 1;
	}
	return res; 
}

ll get_phi(ll x){  //获得x的phi 
	ll res = x;
	for(int i = 2;i <= x/i; i++){
		if(x % i == 0){
			res = res / i * (i-1);
			while(x % i == 0)x /= i; 
		}
	}
	if(x > 1)res = res / x * (x-1);
	return res; 
}

int main(){
	int T = 1;
	ll L;
	while(~scanf("%d",&L) && L){
		//gcd(8,L);
		int d = 1;
		while(d*2<=8 && L%(d*2) == 0){
			d *= 2;
		}
		ll c = L * 9 / d ;
		
		ll ans = 1e18;
		if(c % 2 == 0 || c % 5 == 0){
			ans = 0;
		} 
		else{
			ll phi = get_phi(c);
		
		//枚举因子
		for(ll i=1;i<=phi/i;i++){
			if(phi % i == 0){
				if(qpw(10,i,c) == 1)ans = min(ans,i);
			    if(qpw(10,phi/i,c) == 1)ans = min(ans,phi/i);
			}		
		} 
		}
		 printf("Case %d: %lld\n",T++,ans); 
	}
	return 0;
}