POJ1995(整数快速幂）

http://poj.org/problem?id=1995

题意：求(A1^B1 + A2^B2 + .....Ah^Bh)%M

直接快速幂，以前对快速幂了解不深刻，今天重新学了一遍so easy

以a^b为例：如果b是偶数那么一定可以写成 （a^2 * a^2 ....）一共是b/2个，那么其实就可以写成（a*a)^(b/2)，另a = a*a，b= b/2,此时还是求a^b,只不过a和b已经变了，但是没有问题，还是可以按照上面的方法在判断的，如果b是奇数的话就把a的一个给拿出来先与ans相成，然后现在a的个数就是偶数，就可以按照偶数方法在划分；

2^8 : 8是偶数，a = a* a = 2 * 2= 4,b = b / 2 = 4;也就变成了 4^4,4还是偶数，a = 4 * 4 = 16,b= b / 2 = 2，也就变成了16 ^ 2,2还是偶数再分 a= a * a = 16 * 16,b = b / 2 = 1也就变成了 （16 * 16）^ 1，此时1是奇数ans * a就结束了，

2^5：5是奇数,先把一个2拿出来与ans相乘,ans * a = ans * 2，然后就可以划分了，a = a * a = 2 * 2 = 4，b= b / 2 = 2，也就变成了 4 ^ 2，此时的2是偶数，又可以把a与a相乘，a= a * a = 4 * 4 ,b = b / 2 = 1，也就变成了16^1，此时的b为奇数ans * a = 2 * 16;

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
int a,b,M,sum;
int main()
{
    int t,h;
    scanf("%d", &t);
    while(t--)
    {
        scanf("%d%d",&M,&h);
        sum = 0;
        while(h--)
        {
            scanf("%d%d", &a,&b);
            int ans = 1;
            while(b)
            {
                if(b & 1)
                {
                    ans = (ans % M) * (a % M) % M;
                }
                a = (a % M ) * (a % M) % M;
                b = b >> 1;
            }
            sum = (sum + ans) % M;
        }
        printf("%d\n", sum);
    }

    return 0;
}