HDOJ 1316 2^x mod n = 1

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1395

题目注意一下n小于等于1的处理……

#include

#include

using namespace std;

//蒙哥马利快速幂

int Montgomery(__int64 a,int b,int c){

    __int64 ans=1;

    while(b){

        if(b&1){

            ans=(ans*a)%c;

        }

        b>>=1;

        a=(a*a)%c;

    }

    return ans;

}

int main(){

    int n,i;

    while(scanf("%d",&n)!=EOF){

        if((n&1)==0 || n<=1){

            printf("2^? mod %d = 1\n",n);

        }

        else{

            for(i=1;;i++){

                if(Montgomery(2,i,n)==1){

                    printf("2^%d mod %d = 1\n",i,n);

                    break;

                }

            }

        }

    }

    return 0;

}