我们有时可能会苦恼:为什么C++系统函数里没有立方根?我也这样苦恼过。直到有一天,我奇思妙想,用二分法求了立方根,思路如下:
0.准备工作 除了定义基本的变量,我还定义了一个双精度eps,至于它的用途,下面会讲到。另外,我定义了一个qrt,这是为了防止编译失败而定义的,可以直接把它替换成3.0这个常数值。
1.输入 输入一个数n,这很简单
2.二分 因为左指针、右指针、中间值都是double类型,所以左指针不会出现大于等于右指针的情况,只需用while(1)即可。定义左指针为1,右指针为n的平方根(一个数的立方根不可能大于其平方根),然后计算mid的立方与n的差的绝对值,如果这个值的误差小于eps就跳出循环,否则,如果其大于n,则将右指针左移,如果小于n,则将左指针右移,再次循环。
3.输出 也非常简单,直接输出mid就可以了
4.时间复杂度 因为此代码中l,r,mid都是double类型,所以很难测算,也请各位大佬多多指教。
接下来附上完整代码:
#include<bits/stdc++.h>
using namespace std;
double l,r,num,mid,eps=1e-10,qrt=3.0;//eps为精度
int main(){
cin>>num;
l=0;r=sqrt(num);
while(1){
mid=(l+r)/2;
if(fabs(num-pow(mid,qrt))<=eps) break;//fabs()函数:绝对值
else{
if(pow(mid,qrt)<num) l=mid;
else r=mid;
}
}
printf("%.9f",mid);//输出位数为(精度位数-1)
}
如有和之前博客重复,纯属巧合,此代码没有借鉴任何人的代码,另外大家不要模仿我这史一样的码风!