【题】【数论（欧拉）】NKOJ 3685 8数

NKOJ 3685 8数
时间限制 : - MS 空间限制 : 65536 KB 评测说明 : 时限：1000ms

问题描述
给一个正整数N,问是否存在N的倍数M，且M的各个位全部由数字8组成，如果存在多个取最小的M
并输出M由几个8组成。

输入格式
一行，一个整数N

输出格式
一行一个整数，表示所求的结果。如果无解，输出0

样例输入 1
8

样例输出 1
1

样例输入 2
11

样例输出 2
2

样例输入 3
16

样例输出 3
0

样例输入 4
18

样例输出 4
9

提示
1<=N<=2000000000
注意：此题可能要使用快速乘法

思路：强行配成欧拉：
x个8组成的数，可表示为（10^x-1）*8/9;
设该“8数”为n的k倍，则（10^x-1）* 8/9==k*n 显然x即为所求
移项得（10^x-1）* 8==9* k*n ，若将8直接除到右边可能右边不为整数
所以(方法：)设d=gcd(8,n) 原式为（10^x-1）* 8/d==9* k*n/d ,此时可保证右边为整数
移项得（10^x-1）* 8/d==k* 9* n/d
设p= 8/d，q= 9* n/d，则原式为（10^x-1）* p==k* q
//(结论1：)a*b==c*d，且b，d互质时，a≡0 （mod d）
因为8，9互质，所以p，q互质
所以（10^x-1）≡0 （mod q）
即 10^x≡1（mod q）
根据欧拉定理 当10与q互质时 10^(phi[q])≡1 mod(q)
所以，若gcd（10，phi[n]）！=1时无解
(结论2：)模运算取幂存在循环节，如果有x^(a) ≡ y (mod q)，那么x^(a+i * s) ≡ y (mod q) i>=0且为整数
对于 10^x≡1 (mod q)，设指数循环节为T
首先k=0是一个解:10^0 ≡ 1 (mod q), 则10^(0+i*T) ≡ 1 (mod q)
即 10^(i*s) ≡ 1 (mod q)
由于已经确定phi[q]是一个解了,所以phi[q]一定是这个循环节T的倍数。
所以只需枚举phi[q]的因数x，找出其中最小的满足 10^x ≡1 (mod q) 的x即可。

优化：
1、(注意:)枚举因数时从两端枚举，否则超时。
2、使用快速乘法。

代码：

#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long

ll yin[100000000];
ll cnt=0;

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

ll multi(ll a,ll b,ll c)
{
    ll ans=0;
    a%=c;
    while(b)
    {
        if(b&1) ans=(ans+a)%c;
        b>>=1;
        a=(a+a)%c;
    }
    return ans;
}

ll power(ll a,ll b,ll c)
{
    ll ans=1;
    a%=c;
    while(b)
    {
        if(b&1) ans=multi(ans,a,c);
        b>>=1;
        a=multi(a,a,c);
    }
    return ans;
}

ll ph(ll k)
{
    ll ans=k,a=k;
    for(int i=2;i*i<=a;i++)
    {
        if(a%i==0)
        {
            ans-=ans/i;
            while(a%i==0) a/=i;
        }
    }
    if(a>1) ans-=ans/a;
    return ans;
}

int main()
{
    ll n;scanf("%I64d",&n);
    ll d=gcd(8,n);
    ll q=9*n/d;
    if(gcd(10,q)!=1)
    {
        puts("0");
        return 0;
    }
    ll b=ph(q);
    ll m=(ll)sqrt(b+0.5),minn=b;
    for(ll i=1;i<=m;i++)
    {
        if(b%i==0)
        {
            if(power(10,i,q)==1) 
            {
                printf("%I64d",i);
                return 0;
            }
            if(power(10,b/i,q)==1)
            {
                minn=b/i;
            }
        }
    }
    printf("%I64d",minn);
}