奇思妙想之二分法求立方根

我们有时可能会苦恼：为什么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）
}

如有和之前博客重复，纯属巧合，此代码没有借鉴任何人的代码，另外大家不要模仿我这史一样的码风！