[洛谷P4861]按钮
题目
求解高次同余方程\(k^{x} \equiv 1(mod\) \(m)\),不保证\(m\)为质数
思路
我会扩展BSGS
有个欧拉定理\(k^{\varphi(m)} \equiv 1(mod\) \(m)\),可知答案一定为\(\varphi(m)\)的因子
桥豆麻袋,为什么可以使用欧拉定理?说好的\(gcd(k,m)=1\)才能用呢?
下面证明,该方程有解当且仅当\(gcd(k,m)=1\)(我比较菜所以写的有点啰嗦啦)
\(gcd(k^m,m) | k^m\) 又 \(gcd(k,m) | gcd(k^m,m)\)-->\(gcd(k,m) | k^x\)
此题中\((k^x \% m)=1\),所以\(gcd(k,m) | 1\),即\(gcd(k,m)=1\)
于是\(O(\sqrt{m})\)求出\(\varphi(m)\)后再\(O(\sqrt{\varphi(m)})\)枚举约数检验就行了
时间复杂度\((logn\sqrt{n})\)
Code
#include<bits/stdc++.h>
#define Min(x,y) ((x)<(y)?(x):(y))
#define Max(x,y) ((x)>(y)?(x):(y))
using namespace std;
typedef long long ll;
ll k,mod,phi;
template <class T>
void read(T &x)
{
char c; int sign=1;
while((c=getchar())>'9'||c<'0') if(c=='-') sign=-1; x=c-48;
while((c=getchar())>='0'&&c<='9') x=(x<<1)+(x<<3)+c-48; x*=sign;
}
ll quickpow(ll a,ll b)
{
ll ret=1;
while(b)
{
if(b&1) ret=ret*a%mod;
a=a*a%mod;
b>>=1;
}
return ret;
}
ll gcd(ll a,ll b) { return (!b) ? a : gcd(b,a%b); }
int main()
{
read(mod);read(k);
if(gcd(k,mod)!=1) { puts("Let's go Blue Jays!");return 0; }//无解
phi=mod;
ll p=mod;
for(int i=2;i*i<=p;++i)//求phi(mod)
if(p%i==0)
{
while(p%i==0) p/=i;
phi=phi*(i-1)/i;
}
if(p!=1) phi=phi*(p-1)/p;
ll ans=10000000000000;
for(int i=1;i*i<=phi;++i)
{
if(phi%i) continue;
if(quickpow(k,i)%mod==1) { ans=i; break; }
if(quickpow(k,phi/i)%mod==1) ans=phi/i;
}
cout<<ans<<endl;
return 0;
}