HDU 1395

这道题目，可以使用笨一点的办法，利用蒙哥马利的快速幂取模的办法做，代码如下：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
__int64 calcute(__int64 a,__int64 b,__int64 c)//蒙哥马利快速幂取模
{
    __int64 ans=1;
    while(b)
    {
        if(b&1)
        {
            ans=(ans*a)%c;
        }
        a=(a*a)%c;
        b/=2;
    }
    return ans;
}




int main()
{
    __int64 n,m;
    while(scanf("%I64d",&n)!=EOF)
    {
        if(!(n&1)||n<=1)
        {
            printf("2^? mod %d = 1\n",n);
        }
        else
        {
            __int64 k;
            for(__int64 i=1;; i++)
            {


                m=calcute(2,i,n);
                if(m==1)
                {
                    k=i;
                    break;
                }
            }
            printf("2^%I64d mod %I64d = 1\n",k,n);
        }
    }
    return 0;
}

还有另外一种方法，暴力：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
int main()
{
    __int64 n;
    __int64 i;
    while(scanf("%I64d",&n)!=EOF)
    {
        if(!(n&1)||n<=1)
        {
            printf("2^? mod %I64d = 1\n",n);
        }
        else
        {
            __int64 sum=1;
            for(i=1; ; i++)
            {
                sum*=2;
                if(sum%n==1)
                {
                    break;
                }
                sum%=n;
            }
            printf("2^%I64d mod %I64d = 1\n",i,n);
        }
    }
    return 0;
}