uva 10006 Carmichael Numbers

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=947

打出素数表，快速幂取模。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#define ll long long
#define maxn 65000
using namespace std;

int n;
bool vis[maxn*10];


ll pow_mod(ll a,ll p,ll n)
{
    if(p==0) return 1;
    ll ans=pow_mod(a,p/2,n);
    ans=ans*ans%n;
    if(p%2==1) ans=ans*a%n;
    return ans;
}
void Get_prime()
{
    memset(vis,false,sizeof(vis));
    int m=(int)sqrt(maxn*1.0);
    for(int i=2; i<=m; i++)
    {
        if(!vis[i])
        {
            for(int j=i*i; j<=maxn; j+=i)
            {
                vis[j]=true;
            }
        }
    }
}


int main()
{
    Get_prime();
    while(scanf("%d",&n)!=EOF)
    {
        if(n==0) break;
        if(!vis[n])
        {
            printf("%d is normal.
",n);
            continue;
        }
        bool flag=true;
        for(int i=2; i<n; i++)
        {
            ll c=pow_mod((ll)i,(ll)n,(ll)n);
            if(c!=i)
            {
                flag=false;
                break;
            }
        }
        if(flag)
        {
            printf("The number %d is a Carmichael number.
",n);
        }
        else printf("%d is normal.
",n);
    }
    return 0;
}