hdu 2837 Calculation

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

额，这两天做上了数论，搞了一个上午终于AC了，不过还是看了别人的代码。。

就是利用一个公式A^B %C  = A^( B%phi[C]  + phi[C]) %C    (B>=phi[C])

大致思路还算清晰，就是要注意中间的细节。

回头再看看吧~~

# include<stdio.h>
# include<string.h>
# include<stdlib.h>
__int64 eular(__int64 num)
{
    __int64 ret=1;
    for(int i=2;i*i<=num;i++)
    {
        if(num%i==0)
        {
            num/=i;ret*=i-1;
            while(num%i==0) {num/=i;ret*=i;}
        }
    }
    return num>1? ret*(num-1):ret;
}
__int64 mod(__int64 A,__int64 B,__int64 C) //A^B%C
{
    __int64 ans,res=1;
    int i;
    ans=1;
    for(i=1;i<=B && res<=C;i++) res=res*A; 
    while(B)
    {
        if(B%2) ans=ans*A%C;
        A=A*A%C;
        B/=2;
    }
    if(res>=C) return ans+C;
    return ans;
}

///// A^B %C = A^( B % Phi[C] + Phi[C]) % C       (B >= Phi[C])
__int64 dfs(__int64 n,__int64 m) //f[n]%m
{
    __int64 ans;
    ans=eular(m);
    if(n==0) return 1%m;
    return mod(n%10,dfs(n/10,eular(m)),m);
}

int main()
{
    int ncase;
    __int64 num,n,m;
    scanf("%d",&ncase);
    while(ncase--)
    {
        scanf("%I64d%I64d",&n,&m);
        num=dfs(n,m);
        printf("%I64d\n",num%m);
    }
    return 0;
}