理科男【数论】

题目大意：

对求分数$\tfrac{A}{B}$，求其在$K$进制下是有限小数还是循环小数。如果是有限小数，求小数点后的位数；如果是循环小数，则求混循环部分和循环节的长度又分别是多少。
$Input$

3
1 8 10
17 99 10
217 990 10

$Ouput$

3 0
0 2
1 2

---

思路：

这道数论，好难啊。。。
证明我肯定是看不懂了。只能按照题解所说的打，但是什么都不懂。
在这里把 题解 的证明放一下。

首先把$\tfrac{A}{B}$约分成既约分数。设 $a[1] = A$，$r[n]$ 为原分数小数点后第$n$位的数。
显然有 $r[1] = floor(K \times \tfrac{a[1]}{B})$。
剩下来的余数 $a[2] = K \times a[1]$ $mod$ $B$。
依此类推我们有 $r[n] = floor(K \times \tfrac{a[n]}{B})$，$a[n] = K \times a[n - 1]$ $mod$ $B$。
不难发现如果 $a[p] = a[q]（p < q）$，那么小数点后第 $p$ 位到第 $q - 1$ 位这一段就可以视为一
个循环节。
暴力计算数列 $a$，找到第一个与前面重复的项，就可以找到最短循环节了。
这个重复的项前面的部分导出混循环部分。
如果最早在 $p$ 处计算到 $a[p] = 0$，那么原分数就是一个小数点后有 $p - 1$ 位的有限小数。
以上便是 50 分的解法。

---

50分解法还很容易理解，但是接下来的100分做法就真的很难理解了。。。（也许是我太菜了吧）

下面我们对 $a$ 数列的性质做一些讨论。
如果 $gcd(B, K) = 1$，对于任意的 $i$ 都有 $gcd(a[i], B) = 1$。
设 $K'$ 为 $K$ $mod$ $B$ 时的乘法逆元，即 $KK'$ $mod$ $B = 1$。由乘法逆元的性质 $K'$ 存在且唯一。
假设最早出现重复的位置是 $a[p] = a[q] (p < q)$。
如果 $p$ ≠ $1$，那么 $a[p - 1] = K' \times a[p]$ $mod$ $B = K' \times a[q]$ $mod$ $B = a[q - 1]$。
也就是出现了更早的重复，与题设矛盾。所以显然有 $p = 1$。
这时，显然原分数是一个纯循环小数，且最短循环节长度是 $q - 1$。
设 $x = q - 1$。显然 $a[q] = a[1] \times K^x$ $mod$ $B = a[1]$，于是 $K^x$ $mod$ $B = 1$。
这就转化成了求 $K$ $mod$ $B$ 的阶的问题了。
由欧拉定理 $K^{phi(B)} = 1 ($$mod$ $B)$，由阶的性质 $x | phi(B)$。
我们可以将 $phi(B)$ 分解素因数，并初始化 $x = phi(B)$。
之后考虑 $phi(B)$ 的每个素因数 $p$。如果 $K^{(x/p)} = 1 ($$mod$ $B)$，就 $x ← \tfrac{x}{p}$，并继续试除 $p$。
否则转下一个素因数。这样就可以求出 $K$ $mod$ $B$ 的阶了，这就是最短循环节的长度。
如果 $gcd(B, K) > 1$，那么 $gcd(a[2], B) > 1$。设 $gcd(a[2], B) = g$，不难发现对于任意的 $i ≥ 2$，有 $g | (a[i],B)$。
不妨设 $B' = \tfrac{B}{g}$，$a'[i] = \tfrac{a[i]}{g} (i ≥ 2)$。
若此时 $gcd(B', K) = 1$，就转化为了上面的情况。否则继续这个过程。
如果上面的转换进行了 $T$ 次，由于 $a[1]$ 到 $a[T]$ 与后面 $a$ 数列的循环无关。
卡一下范围便会知道循环节的最后一个数字与混循环部分最后一个数字一定不相等。
于是原分数的混循环长度就是 $T$ 了。
特殊地，如果在 $T$ 次转换之后得到的最后一个 $B = 1$，那么之后 $a$ 数列的值全为 0。这时
我们可以断言原分数是一个小数点后有 $T$ 位的有限小数。
以上便是 100 分的做法。

好吧这么多我知道没人会看

---

代码：

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;

long long n,a,k,sum,gcd,b,num,ans;

long long phi(long long x)
{
    long long ans=x;
    for (long long i=2;i*i<=x;i++)
     if (!(x%i))
     {
        ans=ans/i*(i-1);
        while (!(x%i)) x/=i;
     }
     if (x>1) ans=ans/x*(x-1);
    return ans;
}

long long binary(long long x,long long y,long long o)
{
    long long ans=0;
    for (;y;y>>=1)
    {
        if (y&1) ans=(ans+x)%o;
        x=(x<<1)%o;
    }
    return ans;
}

long long ksm(long long x,long long y,long long o)
{
    long long ans=1;
    for (;y;y>>=1)
    {
        if (y&1) ans=binary(ans,x,o);
        x=binary(x,x,o);
    }
    return ans;
}

int main()
{
    scanf("%lld",&n);
    while (n--)
    {
        scanf("%lld%lld%lld",&a,&b,&k);
        gcd=__gcd(a,b);
        a/=gcd;
        b/=gcd;
        sum=0;
        while (__gcd(b,k)>1)
        {
            sum++;
            gcd=__gcd(b,k);
            b/=gcd;
        }
        printf("%lld ",sum);
        if (b==1) 
        {
            printf("0\n");
            continue;
        }
        num=ans=phi(b);
        for (int i=2;i*i<=num;i++)
         if (!(num%i))
         {
            while ((!(ans%i))&&ksm(k,ans/i,b)==1) ans/=i;
            while (!(num%i)) num/=i;
         }
        if (num>1&&ksm(k,ans/num,b)==1) ans/=num;
        printf("%lld\n",ans);
    }
    return 0;
}