题意转化为:求方程 K x ≡ 1 ( m o d M ) K^x\equiv1(\bmod M) Kx≡1(modM) 的最小解。
Step 1.判无解
什么情况下会无解?方程可以转化为: K × K x − 1 − a M = 1 K\times K^{x-1}-aM=1 K×Kx−1−aM=1 的形式,其中 K , M K,M K,M 为已知。此时可以使用裴蜀定理,即 gcd ( K , M ) = 1 \gcd(K,M)=1 gcd(K,M)=1 时有解。
Step 2.计算
欧拉定理:若 gcd ( a , m ) = 1 \gcd(a, m)=1 gcd(a,m)=1,则 a φ ( m ) ≡ 1 ( m o d m ) a^{\varphi(m)} \equiv 1(\bmod m) aφ(m)≡1(modm)。
使用欧拉定理,可得
K
φ
(
M
)
≡
1
(
m
o
d
M
)
K^{\varphi(M)}\equiv 1(\bmod M)
Kφ(M)≡1(modM),
φ
(
M
)
\varphi(M)
φ(M) 为一可行解。
而最小解肯定是
φ
(
M
)
\varphi(M)
φ(M) 的一个因数。
证明:
反证。若存在一 x x x 使得 x ∤ φ ( M ) x\not|\varphi(M) x∣φ(M) 且 K x ≡ 1 ( m o d M ) K^x\equiv1(\bmod M) Kx≡1(modM),则能找到一正整数 b b b 使得 b x < φ ( M ) < ( b + 1 ) x bx<\varphi(M)<(b+1)x bx<φ(M)<(b+1)x。
∵ k φ ( M ) − b x ≡ 1 ( m o d M ) \because k^{\varphi(M)-bx}\equiv1(\bmod M) ∵kφ(M)−bx≡1(modM)
∴ x \therefore x ∴x 不是最小解(有更小解 φ ( M ) − b x \varphi(M)-bx φ(M)−bx)。
枚举 φ ( M ) \varphi(M) φ(M) 的因数,用快速幂验证即可。复杂度 O ( n log n ) O(\sqrt n\log n) O(nlogn)。
//P4861
#include <bits/stdc++.h>
using namespace std;
//k^x === 1 mod M
//k^phi(M) === 1 mod M
int gcd(int a, int b){ return b ? gcd(b, a%b) : a; }
int phi(int k){
int ans = k;
for(int i = 2; i * i <= k; ++ i){
if(k%i==0) ans = ans / i * (i-1);
while(k%i==0) k/=i;
}
if(k) ans = ans / k * (k-1);
return ans;
}
int qp(int a, int b, int p){
int ans = 1;
while(b){
if(b&1) ans = (long long) ans * a % p;
a = (long long) a * a % p;
b >>= 1;
}
return ans;
}
int main(){
int m, k, ans = 0x7fffffff; scanf("%d%d", &m, &k);
if(gcd(m, k) != 1) return puts("Let's go Blue Jays!") & 0;
int p = phi(m);
for(int i = 1; i * i <= p; ++ i)
if(p % i == 0){
if(qp(k, i, m) == 1) ans = min(ans, i);
if(qp(k, p/i, m) == 1) ans = min(ans, p/i);
}
printf("%d\n", ans);
return 0;
}