二分求n次方根
以求x的三次方根为例,一般情况下x的三次方根一定在0到x之间,所以我们令l= 0, r = x,在l和r之间进行二分,但是当x小于1时,比如x = 0.001,很明显,x的三次方根不在0到x之间,这时我们二分的区间应该是(0,1)(这里都默认x > 0,当x为负数时方法是一样的),所以,我们应该二分的区间是(0,max(1,x))
,但是上面的方法有一点麻烦,就是我要比较x和1,那么索性我直接取l和r与x无关,即在(-10000,10000)之间进行二分(具体区间根据题意判断),这样令l=-10000,r = 10000,那么答案一定在(l, r)之间,每次令mid= (l+r)/2,然后比较midmidmid与x的大小,前者大,则令 r = mid,后者大,则令l = mid,当然每个题都会给一个要求的精度,假设精度是1e-4,那么需要使r-l的精度达到1e-6,如果要求精度是1e-6,那么r-l的精度要达到1e-8,(经验值),接下来给出求三次方根的代码。
#include<iostream>
#include<algorithm>
using namespace std;
double n;
int main()
{
cin >> n;
double l = -10000, r = 10000;
while(r - l > 1e-8)
{
double mid = (l + r) / 2;
if(mid * mid * mid >= n) r = mid;
else l = mid;
}
printf("%.6lf\n", l);
return 0;
}