HDU—2035—人见人爱A^B—【快速幂取模】

人见人爱A^B

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

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
  
  

 

法一：暴力。。。居然可以过耶→_→

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define LL __int64
int main()
{
    LL A,B;
    while(scanf("%I64d%I64d",&A,&B)!=EOF)
    {
        if(A==0 && B==0)
            break;
        LL cj=1;
        while(B--)
            cj=(cj*A)%1000;
        printf("%I64d\n",cj);
    }
    return 0;
}

法二：二分（加速） 就是当指数是偶数的时候，底数平方，指数除以二

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define LL __int64

LL pown(LL A,LL B)
{
    LL temp=1,cj=A; //cj第一次乘的时候，值一定是A
    while(B!=1)
    {
        if(B%2==0)
        {
            cj=(cj*cj)%1000;
            B/=2;
        }
        else
        {
            temp=(temp*cj)%1000;
            B--;
        }
    }
    return cj*temp%1000;
}
int main()
{
    LL A,B;
    w
hile(scanf("%I64d%I64d",&A,&B)!=EOF)
    {
        if(A==0 && B==0)
            break;
        printf("%I64d\n",pown(A%1000,B));
    }
    return 0;
}

/*
A*A*A*A*A*A*A*A*A
=A*(AA)*(AA)*(AA)*(AA)
=A*(AA*AA)*(AA*AA)
=A*(AAAA*AAAA)
*/

法三：快速幂取模：

快速幂取模【总是会被高大上的名字吓到。。。】，一开始不知道什么是快速幂，后来网上各种找。。。原来也用到了二分的，而且很像，把上面的改改就是了。。。

附上找到的一些比较有用的资料：

百度文库的一篇，很详细：http://wenku.baidu.com/link?url=pm1rGxApZIGIL-AnyjcmerZq9Ln0dgF8JGwjlX4AaAqKi-k3Hqn5vk7sckAUu1tcvm7VmX2ZvMsNO0it5heZ4xbhovQi6RbGrNtHz_h6qO_

百度百科：快速幂 ： http://baike.baidu.com/view/4533005.htm?fr=aladdin

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define LL __int64

LL pown(LL A,LL B,LL C) //添加一个参变量，这样更加具有一般性，可以当模板了。。
{
    LL temp=1;
    A=A%C;          //依据公式：a^b mod c = (a mod c)^b mod c
                    //另外一个公式：(a*b) mod c = (a mod c)*(b mod c) mod c
    while(B>0)      //B=1的时候成到temp上，代码更简练
    {
        if(B%2==1)
            temp=(temp*A)%C;
        B/=2;       //这边可能需要好好理解一下
        A=(A*A)%C;  //最后A会多乘一次。。。但返回的是temp 嘿嘿。。。
    }
    return temp;
}
int main()
{
    LL A,B;
    while(scanf("%I64d%I64d",&A,&B)!=EOF)
    {
        if(A==0 && B==0)
            break;
        printf("%I64d\n",pown(A%1000,B,1000));
    }
    return 0;
}

法三升级：加入位运算，更快。。。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define LL __int64

LL pown(LL A,LL B,LL C) //添加一个参变量，这样更加具有一般性，可以当模板了。。
{
    LL temp=1;
    A=A%C;          //依据公式：a^b mod c = (a mod c)^b mod c
                    //另外一个公式：(a*b) mod c = (a mod c)*(b mod c) mod c
    while(B>0)      //B=1的时候成到temp上，代码更简练
    {
        if(B&1)     //if(B%2==1)，学习了。。。位运算——与运算
            temp=(temp*A)%C;
        B>>=1;       //B/=2; x<<=1就是x*2 ，x>>=1就是x/2
        A=(A*A)%C;  //最后A会多乘一次。。。但返回的是temp 嘿嘿。。。
    }
    return temp;
}
int main()
{
    LL A,B;
    while(scanf("%I64d%I64d",&A,&B)!=EOF)
    {
        if(A==0 && B==0)
            break;
        printf("%I64d\n",pown(A%1000,B,1000));
    }
    return 0;
}