BIT1044 A^B mod C

题目大意：Given A,B,C, You should quickly calculate the result of A^Bmod C.(1<=A,B<2^62, 1 <= C <= 10^9)

 

考察点：快速幂

 

思路分析：这题参考了http://wenku.baidu.com/view/a1d05c8171fe910ef12df8b9.html和http://chensmiles.blog.163.com/blog/static/121463991200962851118480/

首先，有如下两个公式A^B mod C=(A^2)^(B/2) mod C,if not(odd(b))

                                   A^B mod C=(A^2)^((B div 2)*A) mod C,if odd(b)

接下来我们可以很容易想到通过不断的套公式降幂求解。但是会有一个问题，就是当A比较大的时候，A^2以及*(B div 2)*A可能会爆掉。所以我们要再处理一下如何计算A*A，

考虑将A化成二进制，如A=(a0*2^n+a1*2^n-1+......)，乘开来，再利用(a+b) mod c=((a mod c)+(b mod c)) mod c即可

 

#include<stdio.h>

unsigned long long multiply(unsigned long long a,unsigned long long b,unsigned long long c)
{
         if (a==0) return 0;
         return ((((a&1)*b)%c+(multiply(a>>1,b,c)<<1)%c)%c);
}
         
int main()
{
    unsigned long long a,b,c,ans;
    while (scanf("%llu%llu%llu",&a,&b,&c)!=EOF)
    {
          a=a%c;
          ans=1;
          while (b>0)
          {
                if (b%2==1) ans=multiply(ans,a,c);
                b=b/2;
                a=multiply(a,a,c);
          }
          printf("%llu\n",ans);
    }
    return 0;
}