2^x mod n = 1 HDU

欧拉定理：

如果 (a,p) 互质，那么

[a^{varphi(p)}=1 mod; p ]

但注意满足该公式的数中 (varphi(p)) 不一定是最小的。

代码：

#include <bits/stdc++.h>
using namespace std;
int getphi(int n)
{
    int res=n;
    for(int i=2;1LL*i*i<=n;i++)
    {
        if(n%i==0)
        {
            res-=res/i;
            while(n%i==0)
                n/=i;
        }
    }
    if(n>1)
        res-=res/n;
    return res;
}
int power(int a,int b,int p)
{
    int res=1;
    while(b)
    {
        if(b&1)
            res=a*res%p;
        a=a*a%p;
        b>>=1;
    }
    return res%p;
}
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        if((n&1)&&n>1)
        {
            int ans=getphi(n);
            for(int i=1;i<=ans;i++)
            {
                if(power(2,i,n)==1)
                {
                    ans=i;
                    break;
                }
            }
            printf("2^%d mod %d = 1
",ans,n);
        }
        else
            printf("2^? mod %d = 1
",n);
    }
    return 0;
}