hdu 2462 The Luckiest number（快速乘取摸）

题意：求最小的n使得n个8组成的数mod k为 0

思路：
注意到凡是那种1111111..... 2222222..... 33333.....
之类的序列都可用这个式子来表示：
k*(10^n-1)/9
进而简化：
这个题会变成：
8*(10^n-1)/9=k*m (k是一个整数)

如果gcd(8,m)=t的话
那么原始的式子可以改为：
8/t*(10^n-1)=9m/t*k

由于8/t和9m/t完全没有公约数
所以如果8/t*(10^n-1)要整除9m/t的话，必须是
10^n-1是9m/t的整数倍

转化一下就变成了：
10^n=1(mod m)的形式

但在这样的n即是phi（m）不一定是最小的值，最小值为phi（m）的因，枚举就可以了，须注意的是直接快速模幂的话在中间会溢出，需要快速乘模幂

#include<iostream>
#include<cstdio>
#include<stack>
#include<string.h>
#include<math.h>
using namespace std;
typedef __int64 LL;
LL n;
#define N 50005
LL gcd(LL a,LL b)
{
    return b==0?a:gcd(b,a%b);
}

LL prime[N],vis[N],c;
void inint()
{
    int i,j,k;
    c=0;
    memset(vis,0,sizeof(vis));
    for(i=2; i<N; i++)
    {
        if(!vis[i]) prime[c++]=i;
        for(j=0; j<c&&(k=i*prime[j])<N; j++)
        {
            vis[k]=1;
            if(i%prime[j]==0) break;
        }
    }
}
LL euler(LL x)
{
    LL res=x,i;
    for(i=2; i*i<=x; 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;
}

LL pow(LL a,LL d)
{
    LL r=1,p=10;
    while(a)
    {
        if(a&1) r=(r*p)%d;
        p=(p*p)%d;
        a>>=1;
    }
    return r;
}
LL mmod(LL a,LL b,LL n)
{
    a=a%n;
    LL res=0;
    while(b)
    {
        if(b&1)
        {
            res=res+a;
            if(res>=n)
                res=res-n;
        }
        a=a<<1;
        if(a>=n)
            a=a-n;
        b=b>>1;
    }
    return res;
}

LL exmod(LL a,LL b,LL n)
{
    a=a%n;
    LL res=1;
    while(b>=1)
    {
        if(b&1)
            res=mmod(res,a,n);
        a=mmod(a,a,n);
        b=b>>1;
    }
    return res;
}

int judge(LL a,LL d)
{
    LL res;
    res=pow(a,d);
   // res=exmod(10,a,d);
    if(res%d==1) return 1;
    return 0;
}

stack<LL> w;
void yin(LL x,LL d)
{
    int i,j,m,f=0;
    LL a,b,res=x;
    while(!w.empty()) w.pop();
    m=sqrt(double(x)+0.5);
    for(i=1; i<=m; i++)
    {
        if(x%i==0)
        {
            a=i;
            if(judge(a,d))
            {
                f=1;
                res=a;
                break;
            }
            w.push(x/i);
        }
    }
    while(!f&&!w.empty())
    {
        a=w.top();
        w.pop();
        if(judge(a,d))
        {
            f=1;
            res=a;
            break;
        }

    }
    printf("%I64d\n",res);
}
int main()
{
    int i,j,t=1;
    LL d,x;
    inint();
    while(~scanf("%I64d",&n)&&n)
    {
        d=9*n/(gcd(8,n));
        printf("Case %d: ",t++);
        if(gcd(10,d)!=1)
            printf("0\n");
        else
        {
            x=euler(d);
            yin(x,d);
        }
    }
    return 0;
}