對於冪運算 a^b，最原始的方法是用for迴圈一遍：

double powRaw(double base,int n)
{
	double r = 1;
	for(int i=1;i<=n;i++)
		r = r*base;
	return r;
}

考慮如下下面的事實：

對於r = a^b,若b為偶數，那麼r = a^(b/2)·a^(b/2);若b為奇數，那麼r = a^((n-1)/2)·a^((n-1)/2)·a。利用這個事實可以得到如下的遞迴解法：

double pow(double base,int n)
{
	if (n == 0)
		return 1;
	if (n == 1)
		return base;
	double t = pow(base,n >> 1);
	if (n & 0x1) 
		return t*t*base;
	else
		return t*t;
}
 

double power(double base,unsigned int n)
{
	double  r = 1;
	while(n)
	{
		if (n & 0x1)
			r *= base;
		n = n >> 1;
		base *= base;
	}
	return r;
}

http://blog.csdn.net/prstaxy/article/details/8740838