HDU 1395 2^x mod n = 1（欧拉定理+快速幂）

题目链接：https://cn.vjudge.net/contest/319236#problem/H
Sample Input

2
5

Sample Output

2^? mod 2 = 1
2^4 mod 5 = 1

翻译：
输入一个n（n>0），求书否存在一个x满足2^x mod n = 1，如果存在按格式输出最小的x，不存在按格式对应的位置输出？。
分析：
1。要使2^x mod n = 1，那么 2^x必定等于k*n+1，（k倍）。2 ^ x为偶数，易得n为奇数时，才存在x满足上式。特例：n=1时，不存在x
2。分析出来n为奇数后，可知，2和n时互质的。用欧拉定理即可结题。
欧拉定理：a，n为正整数，且两者互素，则a^phi(n) mod n=1,phi(n)为n的欧拉函数。
欧拉函数：一个正整数n，小于n且和n互质的正整数（包括1）的个数，记做phi(n)。

注：求出来的phi(n)，只是保证的是当x为phi(n)时，等式一定成立，但此时的c不一定为最小的，所以for循环遍历1到phi(n)，用快速幂找到符合情况的最小的x

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define LL long long
LL eular(LL n)/*求n的欧拉函数*/
{
    LL ans=n;
    for(int i=2; i*i<=n; i++)
    {
        if(n%i==0)
        {
            ans-=ans/i;
            while(n%i==0)
                n/=i;
        }
    }
    if(n>1)
        ans-=ans/n;
    return ans;
}
LL mod_pow(LL s,LL x,LL mod)/*求2^i%n的值*/
{
    LL res=1;
    while(x)
    {
        if(x&1)/*x为奇数*/
            res=res*s%mod;
        s=s*s%mod;
        x>>=1;
    }
    return res;
}
int main()
{
    LL n;
    while(~scanf("%lld",&n)&&n)
    {
        if(!(n&1)||n==1)/*有分析可知，n必须为奇数*/
        {
            printf("2^? mod %lld = 1\n",n);
            continue;
        }
        //奇数和2互质
        LL d=eular(n),u;//n的欧拉函数
        for(LL i=1; i<=d; i++)
        {
            LL w=mod_pow(2,i,n);//快速幂
            if(w==1)
            {
                u=i;
                break;
            }
        }
        printf("2^%lld mod %lld = 1\n",u,n);
    }
    return 0;
}

真正的勇气源于对自身的认可，以及对未来的渴望。-2019.8.13