poj3696：同余方程，欧拉定理

感觉很不错的数学题，可惜又是看了题解才做出来的

题目大意：给定一个数n，找到8888....(x个8)这样的数中，满足能整除n的最小的x,若永远无法整除n 则输出0

做了这个题和后面的poj3358给我的感觉是这种复杂的数学题一定要哦上手去写，光想永远是想不出来的= =

做法：

基于欧拉定理：若gcd(a,m)=1 ,则满足 a^φ(m)  mod m=1， 即   a-1=k*m

88888（x个8）可以表示为 (10^x-1)/9*8，整除n

于是可以设 (10^x-1)/9*8=n*k ,移项得到 10^x-1=k*n*9/8

一看，刚好满足 a-1=k*m的形式,由于 n*9/8不一定为整数，所以我们令 m=n*9/gcd(n,8)  替代一个k=k*gcd(n,8)/8当作未知数

所以得到同余方程 10^x mod m=1

首先判断是否有解

由于 a mod m=gcd(a,m)的倍数 当gcd(10,m)>1时，显然无解，反之 则有解。

由欧拉定理只  φ(m)为此方程的一个解，但不一定是最小解

由于mod 乘法是有循环节的，由于 10^0 mod m=1成立 即对0，和φ(m)都成立，所以循环节要么是φ(m)，要么是φ(m)的约数

所以我们只需要对φ(m)进行素因子分解，判断是否满足同余方程，就可以找到最小的解

代码：

#include <iostream>
#include <stdio.h>
#include<string.h>
#include<algorithm>
#include<string>
#include<ctype.h>
using namespace std;
#define I64d lld
long long gcd(long long a,long long b)
{
    return b?gcd(b,a%b):a;
}
long long fac[];
long long nfac;
long long phi(long long n)
{
    long long res=n;
    for(long long i=;i*i<=n;i++)
    {
        if(n%i==)
        {
            res=res-res/i;
            while(n%i==)
                n/=i;
        }
    }
    if(n>)
        res=res-res/n; //可能还有大于sqrt(n)的素因子
    return res;

}
long long random(long long n)
{
    return (long long)(rand()%(n-)+);
}
long long multi(long long a,long long b,long long m)//a*b%m
{
    long long res=;
    while(b>)
    {
        if(b&)
            res=(res+a)%m;
        b>>=;
        a=(a<<)%m;
    }
    return res;
}
long long quickmod(long long a,long long b,long long m) //a^b%m
{
    long long res=;
    while(b>)
    {
        if(b&)
            res=multi(res,a,m);
        b>>=;
        a=multi(a,a,m);
    }
    return res;
}
int check(long long a,long long n,long long x,long long t)
{
    long long res=quickmod(a,x,n);
    long long last=res;
    for(int i=;i<=t;i++)
    {
        res=multi(res,res,n);
        if(res==&&last!=&&last!=n-) return ;
        last=res;
    }
    if(res!=) return ;
    return ;
}

int primetest(long long n)
{
    if(n<)return ;
    if(n==)return ;
    if((n&)==) return ;
    long long x=n-;
    long long t=;
    while((x&)==){x>>=;t++;}
    for(int i=;i<;i++)
    {
        long long a=random(n);
        if(check(a,n,x,t))
            return ;
    }
    return ;
}

long long pollardrho(long long n,long long c)
{
    long long x,y,d,i,k;
    i=;k=;
    x=random(n);
    y=x;
    while()
    {
        i++;
        x=(multi(x,x,n)+c)%n;
        long long tmp=y-x>=?y-x:x-y;
        d=gcd(tmp,n);
        if(d>&&d<n)
            return d;
        if(y==x)
            return n;
        if(i==k)
        {
            y=x;
            k+=k;
        }
    }
}
void findfac(long long n)
{
    if(n==)
        return;
    if(primetest(n))
    {
        fac[nfac++]=n;
        return;
    }
    long long p=n;
    while(p>=n)
        p=pollardrho(n,random(n-));
    findfac(p);
    findfac(n/p);
}
int main()
{
    long long n,m;
    int cas=;
    while(scanf("%I64d",&n),n)
    {
        cas++;
        m=n*/gcd(n,);
        if(gcd(m,)!=)
        {
            printf("Case %d: %d\n",cas,);
            continue;
        }
        long long p=phi(m);
        nfac=;
        findfac(p);
        for(int i=;i<nfac;i++)
        {
            p/=fac[i];
            if(quickmod(,p,m)!=)
                p*=fac[i];

        }
        printf("Case %d: %I64d\n",cas,p);
    }

    return ;
}