2^x mod n = 1（欧拉函数）

接着上一次的题解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
#include
#include

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

结果这种方法耗时竟然比之前的暴搜长，无语了