接着上一次的题解https://blog.csdn.net/u010017231/article/details/84679719继续写:
首先,介绍一下思路:
易知n==1或2时,一定不存在这样的x:当n为偶数时,bn + 1(b为整数)是奇数,而2^x是偶数,故 2^x mod n = 1不可能成立
欧拉定理:若n,a为正整数,且n,a互质,则:,其中phi(n)为欧拉函数,表示小于n的正整数中与n互素的数的个数。
于是,就可以从1~phi(n)遍历,找到最小的x
附上代码:
#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
long long n;
int p[100000];
long long qpow(long long a,long long k)//快速幂
{
long long ans=1,t=a;
while(k>0)
{
ans%=n;
if(k&1)
{
ans*=(t%n);
}
t%=n;
t*=(t%n);
k>>=1;
}
return ans;
}
long long phi(long long k)
{
long long res=1;
for(long long i=2;i<=sqrt(k);i++)
{
if(k%i==0)
{
k/=i;
res*=(i-1);
while(k%i==0)
{
k/=i;
res*=i;
}
}
}
if(k>1)
res*=(k-1);
return res;
}
int main()
{
while(cin>>n)
{
if(n%2==0||n==1)
{
printf("2^? mod %lld = 1\n",n);
continue;
}
long long x=phi(n);
for(long long i=1;i<=x;i++)
if(qpow(2,i)%n==1)
{
printf("2^%lld mod %lld = 1\n",i,n);
break;
}
}
return 0;
}
结果这种方法耗时竟然比之前的暴搜长,无语了