[题解]（同余）POJ_3696_The Luckiest Number

还是挺难的吧......勉强看懂调了半天

首先表达式可以写成 8（10^x -1）/9，题意为求一个最小的x使L | 8（10^x -1）/9

设d=gcd（L，8）

L | 8（10^x -1）/9

<=>9L | 8（10^x -1）

<=>9L/d | 10^x -1 （因为 9L/d 和 8/d 互质了 所以 9L/d 能整除(8/d)*(10^x-1)和 8/d 无关，所以可以去掉）

<=>10^x 同余 1（mod 9L/d）

引理：

若a，n互质，则满足10^x同余1（mod n）的最小正整数x0是phi（n）的约数

反证法：

假设满足a^x 同余 1（mod n）的最小正整数x0不能整除phi（n）

设phi（n）=q*x0+r（0<r<x0)，因为a^x0 同余1（mod n），所以a^(q*x0)同余1（mod n）

根据欧拉定理a^phi(n)同余1（mod n），所以a^r同余1（mod n），与x0最小矛盾

无解的时候就是q与10不互质的时候，因为若q与10有公因子d：
1.若d=2，q=2*k，那么10^x=2^x*5^x=1%2k
   即2^x*5^x=1+2k*m，左边为偶数，右边为奇数，显然矛盾。
2.若d=5，q=5*k，那么10^x=2^x*5^x=1%5k
   即2^x*5^x=1+5k*m，左边是5的倍数，右边不是5的倍数，显然矛盾。

注意：乘的时候会爆longlong，手写乘法，要用根号的试除法求约数，不然会T

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define ll long long
using namespace std;
ll n,cnt;
ll x[];
ll gcd(ll a,ll b){
    return b?gcd(b,a%b):a;
}
ll eular(ll n){
    ll ans=n;
    for(ll i=;i*i<=n;i++){
        if(n%i==){
            ans=ans/i*(i-);
            while(n%i==)n/=i;
        }
    }
    if(n>)ans=ans/n*(n-);
    return ans;
}
ll mul(ll a,ll b,ll mod){
    ll ans=;
    while(b){
        if(b&)ans=(ans+a)%mod;
        a=(a<<)%mod;
        b>>=;
    }
    return ans;
}
ll qpow(ll a,ll b,ll mod){
    ll base=a,ans=;
    while(b){
        if(b&)ans=mul(ans,base,mod);
        base=mul(base,base,mod);
        b>>=;
    }
    return ans%mod;
}

int main(){
    int t=;
    while(){
        int fl=;cnt=;
        scanf("%lld",&n);
        if(n==)break;
        ll d=*n/gcd(n,);
        if(gcd(,d)!=){
            printf("Case %d: 0\n",++t);
        }
        else{
            ll phi=eular(d);
            for(ll i=;i*i<=phi;i++){
                if(phi%i==){
                    x[++cnt]=i;
                    if(i*i!=phi)x[++cnt]=phi/i;
                }
            }

            sort(x+,x+cnt+);
            for(int i=;i<=cnt;i++)
            if(qpow(,x[i],d)==){
                printf("Case %d: %lld\n",++t,x[i]);
                break;
            }
        }
    }
}