FZU 1759 Super A^B mod C(数论+快速幂+欧拉函数)

    思路:刚开始想直接用JAVA大数处理但是,由于再pow()中不能是大整数形的.

后来才知道有个相当于降幂公式A^(B%oural(c)+oural(c))%c=A^B%c.

#include <stdio.h>
#include <ctype.h>
#include <string.h>
#include <stdlib.h>
#include <limits.h>
#include <math.h>
#include <algorithm>
#include <bitset>
using namespace std;
typedef __int64 ll;
char b[10000000];
int oural(int x)
{
    int i,j;
    int num=x;
    for(i=2;i*i<=x;i++)
    {
        if(x%i==0)
        {
            num=(num/i)*(i-1);
            while(x%i==0)
            {
                x/=i;
            }
        }
    }
    if(x!=1)
    {
        num=(num/x)*(x-1);
    }
    return num;
}
ll quickpow(ll a,ll tmp,ll c)
{
    ll r=1,ba=a;
    while(tmp!=0)
    {
        if(tmp&1)
            r=(r*ba)%c;
        ba=(ba*ba)%c;
        tmp=tmp>>1;
    }
    return r;
}
int main()
{
    ll n,m,i,j,k,a,c;
    while(~scanf("%I64d%s%I64d",&a,b,&c))
    {
        int PHI=oural(c);
        ll tmp=0;
        for(int i=0;b[i]!='\0';i++)
        {
            tmp=(tmp*10+b[i]-'0');
            if(tmp>c)
                break;
        }
        if(tmp<=PHI)
        {
            printf("%I64d\n",quickpow(a,tmp,c));
        }
        else
        {
            tmp=0;
            for(int i=0;b[i]!='\0';i++)
            {
                tmp=(tmp*10+b[i]-'0')%PHI;
            }
            printf("%I64d\n",quickpow(a,tmp+PHI,c));
        }
    }
    return 0;
}