题解 [LuoguP4861]按钮

link

题意转化为：求方程 $K^x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}M)$ 的最小解。

Step 1.判无解

什么情况下会无解？方程可以转化为：$K\times K^{x-1}-aM=1$ 的形式，其中 $K,M$ 为已知。此时可以使用裴蜀定理，即 $\gcd (K,M)=1$ 时有解。

Step 2.计算

欧拉定理：若 $\gcd (a,m)=1$，则 $a^{\phi (m)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}m)$。

使用欧拉定理，可得 $K^{\phi (M)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}M)$，$\phi (M)$ 为一可行解。
而最小解肯定是 $\phi (M)$ 的一个因数。
证明：

反证。若存在一 $x$ 使得 $x\nmid \phi (M)$ 且 $K^x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}M)$，则能找到一正整数 $b$ 使得 $bx<\phi (M)<(b+1)x$。
$\because k^{\phi (M)-bx}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}M)$
$\therefore x$ 不是最小解（有更小解 $\phi (M)-bx$）。

枚举 $\phi (M)$ 的因数，用快速幂验证即可。复杂度 $O(\sqrt{n}\log n)$。

//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;
}