A^B%C-- 二分，再二分_a^b%c高精度-CSDN博客

呵呵，为了解决刚才的(a/b)%m当b*m很大时引起的问题，想起了fzu1752-A^B % c这道题, 去年lvyun给自己推荐的这题。。现在才A掉……好题，不仅解决了问题，还知道了如何处理 A乘B模C的溢出问题。

===================================================================

fzu1752-A^B mod C
题目意识: Given A,B,C, You should quickly calculate the result of A^B mod C. (1<=A,B,C<2^63).

分析：这个题目的特点是C很大. 那么在快速模幂运算中 (A%C)*(A%C)就会越界 (即内部数据类型保存不了,unsigned long long也不行) 解决这个题目的好处: 如果C很大的情况有通用的解法，那么对于(a/b)%m, b*m很大也就可以解决了,从而对于(a/b)%m的公式就有用了。由于这道题目C<2^63，我们知道C+C<2^64, 所以用unsigned long long可以装的下(如何C再大一些,目前不知如何处理,高精度?) 利用上面这个特殊性，我们还是可以用一般的二分来快速求模幂运算，只要解决一个问题就可以:
(A * B) % C 如何不越界? (易知, C+C不会越界)
其实对于这个问题，我们可以再次二分来解决(A*B)%C这一问题. (A*B)%C = (A+A+...+A)%C=A%C+(2A)%C+... 呵呵，其中1,2,4取解决于B的二进制表示。两个A相加就不会越界咯 :)

代码如下(呵呵，多了两个常用东西的模版了):

#include <stdio.h>

typedef unsigned __int64 llong;

// 二分求: (A*B)%C
llong Multi(llong a,llong b,llong n)
{
        llong back=0,temp = a%n;
        while (b)
         {
               if ( (b&1) && ((back+=temp)>=n) ) back-= n;
                if ((temp<<=1) >= n) temp -= n;
               b >>= 1;
        }
         return back;
}

//二分求: A^B%C
llong Pow(llong a, llong b, llong c)
{
        llong ans = 1;
         while (b!=0)
        {
               if (b&1) ans = Multi(ans, a, c);
              a = Multi(a, a, c);
               b >>= 1;
        }
         return ans;
}

int main()
{
llong a, b, c;

        while (scanf("%I64u%I64u%I64u", &a, &b, &c)!=EOF)
                    printf("%I64u\n", Pow(a, b, c));
         return 0;
}

复杂度: (lg(c) * sigma( lg(a) + ,,,,,,,,,,,), 反正很快拉……呵呵

这道题的下面一道：Super A^B mod C

Given A,B,C, You should quickly calculate the result of A^B mod C.(1<=A,C<=1000000000,1<=B<=10^1000000).

恩,此时利用 : a^b mod c = a^(b%phi(c) + phi(c) ) mod c 其中,b>phi(c)