【POJ3696】The Luckiest Number-欧拉定理+快速幂

测试地址：The Luckiest Number
题目大意：给出一个正整数$L(\le 2,000,000,000)$，要使正整数$\underbrace{888...8}_{K位}$能整除$L$，求最小的$K$，如果不存在这样的$K$则输出0。
做法：这题的思路很神……没看题解的时候根本不知道怎么做，调也调了好一会，我好弱啊……
这一题需要使用欧拉定理+快速幂来解决。
我们发现$\underbrace{999...9}_{K位}$可以表示为$10^K-1$，所以$\underbrace{888...8}_{K位}$可以表示为$8/9\times(10^K-1)$，那么存在一个整数$p$使得$8/9\times(10^K-1)=L\times p$，所以$10^K-1=9/8\times L\times p$。由于等式右边要是整数，那么$p$应该满足$8|(L\times p)$。因为$L$已经包含$\gcd(L,8)$这些因子，所以$p$需要包含$8/\gcd(L,8)$这些因子，所以可以将$p$写成$p=8/\gcd(L,8)\times p_1$。将其代入前面的式子，可以得到$10^K-1=9\times L/\gcd(L,8)\times p_1$，设$m=9\times L/\gcd(L,8)$，则原式可以写成$10^K-1=m\times p_1$，这样我们就把问题转化成了求同余方程$10^K\equiv 1(\mod m)$的最小的正整数解。根据欧拉定理，可以得到$10^{\varphi(m)}\equiv 1(\mod m)$，可以证明如果$K$是一个小于$\varphi(m)$的解，那么$K|\varphi(m)$。那么我们就可以按照下列步骤完成这一题：
1.求出$m$，$\varphi(m)$和$\varphi(m)$的所有质因子。这一步可以用$O(\sqrt{N})$的试除法来做。
2.令$x=\varphi(m)$，对于$\varphi(m)$的每一个质因子$p_i$，判断$10^{x/p_i}\%m$是否等于$1$，如果等于则令$x=x/p_i$，否则跳过。由于指数很大，需要使用快速幂来算出结果，而因为数字直接乘起来可能会溢出，要用快速乘代替乘法。
3.对所有质因子进行以上操作后，最后的$x$就是最小的$K$。
那么无解的情况呢？显然可知，如果$\gcd(10,m)\ne1$，那么上面的同余方程就无解。这样我们就完美解决了这一题。
犯二的地方：试除法忘记探测最后剩余是不是1了，导致可能漏掉一个质因子，以后要注意……
以下是本人代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define ll long long
using namespace std;
ll L,fac[1010]={0};

ll gcd(ll a,ll b)
{
  return (b==0)?a:gcd(b,a%b);
}

ll phi(ll x)
{
  ll p=x,s=x;
  for(ll i=2;i*i<=s;i++)
    if (!(x%i))
    {
      p=p/i*(i-1);
      while(!(x%i)) x/=i;
    }
  if (x>1) p=p/x*(x-1);
  return p;
}

void find_factor(ll x)
{
  ll s=x;
  fac[0]=0;
  for(ll i=2;i*i<=s;i++)
    if (!(x%i))
    {
      fac[++fac[0]]=i;
      while(!(x%i)) x/=i;
    }
  if (x>1) fac[++fac[0]]=x;
}

ll mult(ll a,ll b,ll mod)
{
  a%=mod,b%=mod;
  ll s=a,sum=0;
  while(b)
  {
    if (b&1)
    {
      sum+=s;
      if (sum>=mod) sum-=mod;
    }
    b>>=1;
    s<<=1;
    if (s>=mod) s-=mod;
  }
  return sum;
}

ll power(ll a,ll b,ll mod)
{
  ll s=a,sum=1;
  while(b)
  {
    if (b&1) sum=mult(sum,s,mod);
    b>>=1;s=mult(s,s,mod);
  }
  return sum;
}

int main()
{
  int t=0;
  while(scanf("%lld",&L)&&L)
  {
    t++;
    ll m=L/gcd(L,8)*9,p=phi(m),x=p;
    if (gcd(m,10)!=1) {printf("Case %d: 0\n",t);continue;}
    find_factor(p);
    for(int i=1;i<=fac[0];i++)
    {
      while(1)
      {
        x/=fac[i];
        if (power(10,x,m)!=1)
        {
          x*=fac[i];
          break;
        }
        else if (x%fac[i]) break;
      }
    }
    printf("Case %d: %lld\n",t,x);
  }

  return 0;
}