0.题目
实现函数double Power(double base,int exponent),求base的exponent次方,不能使用库函数,不用考虑大数问题。
1.解题思路
需要实现的函数指数为整型,不用考虑指数为浮点型,可以分为下面几类:
- 指数为正整数,不用考虑什么,直接计算结果
- 指数为0时,这个时候要看底数base,如果base为0,则无意义,因为0的0次在数学上是没有意义的;如果base不为0,则结果是1,因为任何非0的数的0次方都是1.
- 指数为负整数的时候,首先求出当指数为正整数时的结果,然后在求它的倒数
int InvalidInput = 0;//用来表示是正常返回还是异常返回(0正常,1异常)
//判断两个浮点数是否相同(不能用等号直接判断)
int Equal(double num1, double num2)
{
if ((num1 - num2)<EXP && (num2 - num1)>-EXP)
return 1;
return 0;
}
//求指数为整数的Power
double PowerWithUnsignedNum(double base, int absexponent)
{
double result = 1.0;
for (int i = 0; i < absexponent; i++)
{
result *= base;
}
return result;
}
//求base的exponent次方
double Power(double base, int exponent)
{
InvalidInput = 0;
//如果底数为0,指数小于0,则异常返回
if ((Equal(base, 0.0) == 1) && (exponent < 0))
{
InvalidInput = 1;
return 0.0;
}
unsigned int absexponent = (unsigned int)exponent;
//如果指数为负数
if (exponent < 0)
{
absexponent = (unsigned int)(-exponent);
}
//先求出当它指数为正的结果
double result = PowerWithUnsignedNum(base, absexponent);
//在求倒数
if (exponent < 0)
{
result = 1.0 / result;
}
return result;
}
如果想要更高的效率,还有一种方法处理:加入我们现在要求2的32次方,是不是可以求出2的16次方,在整体结过平方以下就得到最终结果了,而2的16次方右等于2的8次方的平方,依次类推,可以看出这是一个递归的方法。
a^n=a^(n/2)*a^(n/2) a为偶数时,
a^n=a^(n-1/2)*a^(n-1/2)*a a为奇数时
可以根据上式写出代码:
double Power_R(double base, int exponent)
{
//任何非0的数的0次方都为1
if (exponent == 0)
return 1;
//任何数的1次方是它本身
if (exponent == 1)
return base;
//求出a^(n/2)(exponent >> 1相当于除以2,但是移位运算符效率要高一些)
double result = Power_R(base, exponent >> 1);
//偶数a^(n/2)*a^(n/2)
result *= result;
//奇数a^(n-1/2)*a^(n-1/2)
if (exponent & 1 == 1)
{
result *= base;
}
return result;
}