hdu 2035 人见人爱A^B（两种解决方式：第二种快速幂取模）

人见人爱A^B

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 25752    Accepted Submission(s): 17604

Problem Description

求A^B的最后三位数表示的整数。
说明：A^B的含义是“A的B次方”

 

Input

输入数据包含多个测试实例，每个实例占一行，由两个正整数A和B组成（1<=A,B<=10000），如果A=0, B=0，则表示输入数据的结束，不做处理。

 

Output

对于每个测试实例，请输出A^B的最后三位表示的整数，每个输出占一行。

 

Sample Input

  
  

   
   2 3
12 6
6789 10000
0 0
  
  

 

Sample Output

  
  

   
   8
984
1
  
  

第一种方式：模拟求解（对于B比较小的时候可行）

#include<stdio.h>

int main()
{
    int i,j,k,t,n,m;
    while(scanf("%d %d",&n,&m)!=EOF && n || m)
    {
        t=1;k=0;
        for(i=0;i<m;i++)
        {
            t=t*n;
            if(t>999)
            {
                t=t%1000;
                k=1;
            }
        }
        
            printf("%d\n",t);
    }

    return 0;
}

第二种：快速幂取模（当B比较大，直接操作会超时时）

·二进制的用法，因为其与二叉树有着惊人的类似形式。

·算法存在基础：  a*a%Mod = (a%Mod)*(a%Mod)%Mod

快速幂理解：

1.把b写成二进制的形式：   假如说写成之后为：  1000100101011

2.现在要做的是：分别计算
     a^1
     a^10
     a^1000
     a^100000000 次方等     对Mod取余的结果，之后相乘，再对Mod取余，就是原结果
 
3.那么对于a^100000000 如何计算对Mod取余的结果？
 
                                a^100000000
    
                     a^10000000 *  a^10000000
  
         a^1000000 * a^100000   .....

a^100000 * a^100000   ............

       ......

a^1 * a^1                ...........

  
那么对于a^100000000我们只需要进行 这颗树深度h 次的计算，就可以得出结果。
 

Code:

#include <iostream>
#include <algorithm>
#include <string.h>
#include <math.h>
#include <stdio.h>
#include <queue>
using namespace std;
#define N 1000010
#define LL long long
#define Mod 1000
/*
1.取后三位： Mod 取 1000
2.快速幂： 每次对指数b判断最后一位是否为1， 之后 b /= 2;
3.a = a*a%Mod;

共三步，计算出a^b对Mod取余，最后剩余的为三位数！！

 */
int main(){
    //freopen("in.in","r",stdin);
    int a,b;
    while(scanf("%d %d",&a,&b) && a && b){
        int ans=1;
        while(b){
            if(b & 1){
                ans = ans*a%Mod;
            }
            a = a*a%Mod;
            b /= 2;
        }
        printf("%d\n",ans);
    }

    return 0;
}