HDU1395+快速幂

#include<stdio.h>
int fast_pow( int a,int b,int mod ){
    int res = 1;
    while( b>=1 ){
        if( b%2==1 ){
            res = res*a%mod;
        }
        a = a*a%mod;
        b/=2;
    }
    return res;
}
int main(){
    int n;
    while( ~scanf("%d",&n) ){
        if( n==1||n%2==0 ){
            printf("2^? mod %d = 1
",n);
            continue;
        }
        int ans = 2;
        while( 1 ){
            if( fast_pow( 2,ans,n )==1 ) 
                break;
            ans++;
        }
        printf("2^%d mod %d = 1
",ans,n);
    }
    return 0;
}