A^B mod (大数运算）

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int MAX=;
const int mod=;
char s[MAX];
 
long long pow(long long a,long long b)
{
    long long base=a,r=;
    while(b!=)
    {
        if(b&) r=(r*base)%mod;
        base=(base*base)%mod;
        b>>=;
    }
    return r%mod;
}
 
int main()
{
    while(scanf("%s",s)!=EOF)
    {
        int len=strlen(s);
        long long num=;
        for(int i=;i<len;i++)//大数取模
            num=(num*+(int)(s[i]-''))%(mod-);
        if(num==)//说明num=mod-1
        {
            cout<<pow(,mod-)<<endl;
        }
        else
        {
            num--;
            cout<<pow(,num)<<endl;
        }
    }
    return ;
}
 
HDU4704

于任意自然数，当要求a^p%m时，就可以利用费马小定理化简，只需求(a^(p%(m-1)))%m;(p是素数)