十九.高效率的幂运算

我们想要处理一个整数的幂

long int pow(long int X, unsigned int N)
{
   if(N == 0)
     return 1;
   if( N == 1)
     return X;
   if(N%2==0)
      return pow( X * X, N / 2);
   else
      return pow( X * X, N / 2) * X;         
}

计算X^N，可以使用N - 1次乘法自乘，用递归算法最好。
如果N是偶数则X^N = X ^N/2 * X^N/2;
如果N是奇数，则X的N次幂 = X^（N - 1）/ 2 * X^（N - 1）/ 2 * X;
例如，为了计算X^62，算法将如下进行，它只用到9次乘法：
X^3 = （X^2）*X , X^7 = (X^3) ^ 2 *X
X^15 = (X^7) ^2 X , X^31 = (X^15) ^ 2X
X^62 = (X^15) ^ 2 *X
显然，所需要的乘法次数最多是2logN,因为把问题对分最多需要两次乘法（如果N是奇数）。
(注： 因为由取幂运算得到的数一般都是相当大的，因此，我们只能在假设一台机器能够存储这样一些大整数（或有一个编译程序能够模拟它）的情况下进行分析。)

完整代码：

#include<stdio.h>
long int pow(long int X, unsigned int N)
{
   if(N == 0)
     return 1;
   if( N == 1)
     return X;
   if(N%2 ==0)
      return pow( X * X, N / 2);
   else
      return pow( X * X, N / 2) * X;         
}
int main()
{   long int X;
    unsigned int N;
	pow(X,N);
	scanf("%ld %d",&X,&N);
	printf("%ld",pow(X,N));
	return 0;
}